260   PRange    eta;  member
374 static void THIViscosity(THI thi, PetscReal gam, PetscReal *eta, PetscReal *deta)  in THIViscosity()  argument
390   *eta     = Bd2 * PetscPowReal(eps + gam, exponent);  in THIViscosity()
391   *deta    = exponent * (*eta) / (eps + gam);  in THIViscosity()
547   PRangeClear(&thi->eta);  in THICreate()
551 … * units->meter / (3e7 * units->second), gradu = u / (100 * units->meter), eta, deta, rho = 910 * …  in THICreate()  local
553     THIViscosity(thi, 0.5 * gradu * gradu, &eta, &deta);  in THICreate()
558 …, ratio %8.2g\n", (double)u, (double)gradu, (double)eta, (double)(2 * eta * gradu), (double)(2 * e…  in THICreate()
559       THIViscosity(thi, 0.5 * PetscSqr(1e-3 * gradu), &eta, &deta);  in THICreate()
560 …uble)(1e-3 * u), (double)(1e-3 * gradu), (double)eta, (double)(2 * eta * 1e-3 * gradu), (double)(2…  in THICreate()
710 …v, PetscScalar du[PETSC_RESTRICT], PetscScalar dv[PETSC_RESTRICT], PetscReal *eta, PetscReal *deta)  in PointwiseNonlinearity()  argument
728   THIViscosity(thi, PetscRealPart(gam), eta, deta);  in PointwiseNonlinearity()
731 …etscScalar *u, PetscScalar *v, PetscScalar du[], PetscScalar dv[], PetscReal *eta, PetscReal *deta)  in PointwiseNonlinearity2D()  argument
749   THIViscosity(thi, PetscRealPart(gam), eta, deta);  in PointwiseNonlinearity2D()
791           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  in THIFunctionLocal()  local
795           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIFunctionLocal()
797           if (q == 0) etabase = eta;  in THIFunctionLocal()
798           RangeUpdate(&etamin, &etamax, eta);  in THIFunctionLocal()
802 …fn[l]->u += dp[0] * jw * eta * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp…  in THIFunctionLocal()
803 …fn[l]->v += dp[1] * jw * eta * (2. * du[0] + 4. * dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp…  in THIFunctionLocal()
851   PetscCall(PRangeMinMax(&thi->eta, etamin, etamax));  in THIFunctionLocal()
871 …  (double)thi->eta.cmin, (double)thi->eta.cmax, (double)thi->beta2.cmin, (double)thi->beta2.cmax));  in THIMatrixStatistics()
959 …onverged range %g to %g\n", (double)thi->eta.min, (double)thi->eta.max, (double)thi->eta.cmin, (do…  in THISolveStatistics()
991         PetscReal   phi[4], dphi[4][2], jw, eta, deta, beta2, dbeta2;  in THIJacobianLocal_2D()  local
1001         PointwiseNonlinearity2D(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIJacobianLocal_2D()
1011 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + pp * jw …  in THIJacobianLocal_2D()
1012             Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * eta * 2. * dpl[1] + dp[1] * jw * eta * dpl[0];  in THIJacobianLocal_2D()
1013             Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * eta * 2. * dpl[0] + dp[0] * jw * eta * dpl[1];  in THIJacobianLocal_2D()
1014 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + pp * jw …  in THIJacobianLocal_2D()
1080           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  local
1084           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);
1086           if (q == 0) etabase = eta;
1091 …1 = _mm_set1_pd(dv[1]), dv2 = _mm_set1_pd(dv[2]), jweta = _mm_set1_pd(jw * eta), jwdeta = _mm_set1…
1113 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + dp[2] * …
1114 …           Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * eta * 2. * dpl[1] + dp[1] * jw * eta * dpl[0];
1115 …           Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * eta * 2. * dpl[0] + dp[0] * jw * eta * dpl[1];
1116 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + dp[2] * …