247   PRange    eta;  member
424 static void THIViscosity(THI thi, PetscReal gam, PetscReal *eta, PetscReal *deta)  in THIViscosity()  argument
440   *eta     = Bd2 * PetscPowReal(eps + gam, exponent);  in THIViscosity()
441   *deta    = exponent * (*eta) / (eps + gam);  in THIViscosity()
626   PRangeClear(&thi->eta);  in THICreate()
630 … * units->meter / (3e7 * units->second), gradu = u / (100 * units->meter), eta, deta, rho = 910 * …  in THICreate()  local
632     THIViscosity(thi, 0.5 * gradu * gradu, &eta, &deta);  in THICreate()
637 …%8.2g, ratio %8.2g\n", (double)u, (double)gradu, (double)eta, (double)(2 * eta * gradu, 2 * eta * …  in THICreate()
638       THIViscosity(thi, 0.5 * PetscSqr(1e-3 * gradu), &eta, &deta);  in THICreate()
639 …, (double)(1e-3 * u), (double)(1e-3 * gradu), (double)eta, (double)(2 * eta * 1e-3 * gradu, 2 * et…  in THICreate()
729 …estrict v, PetscScalar du[restrict 3], PetscScalar dv[restrict 3], PetscReal *eta, PetscReal *deta)  in PointwiseNonlinearity()  argument
747   THIViscosity(thi, PetscRealPart(gam), eta, deta);  in PointwiseNonlinearity()
789           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  in THIFunctionLocal_3D()  local
797           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIFunctionLocal_3D()
799           if (q == 0) etabase = eta;  in THIFunctionLocal_3D()
800           RangeUpdate(&etamin, &etamax, eta);  in THIFunctionLocal_3D()
804 …fn[l]->u += dp[0] * jw * eta * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp…  in THIFunctionLocal_3D()
805 …fn[l]->v += dp[1] * jw * eta * (2. * du[0] + 4. * dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp…  in THIFunctionLocal_3D()
853   PetscCall(PRangeMinMax(&thi->eta, etamin, etamax));  in THIFunctionLocal_3D()
980 …", m, (double)nrm, (double)PetscRealPart(val0), (double)PetscRealPart(val2), (double)thi->eta.cmin,  in THIMatrixStatistics()
981 …                         (double)thi->eta.cmax, (double)thi->beta2.cmin, (double)thi->beta2.cmax));  in THIMatrixStatistics()
1077 …onverged range [%g, %g]\n", (double)thi->eta.min, (double)thi->eta.max, (double)thi->eta.cmin, (do…  in THISolveStatistics()
1128           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  in THIJacobianLocal_Momentum()  local
1132           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIJacobianLocal_Momentum()
1134           if (q == 0) etabase = eta;  in THIJacobianLocal_Momentum()
1143 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + dp[2] * …  in THIJacobianLocal_Momentum()
1144 …           Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * eta * 2. * dpl[1] + dp[1] * jw * eta * dpl[0];  in THIJacobianLocal_Momentum()
1145 …           Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * eta * 2. * dpl[0] + dp[0] * jw * eta * dpl[1];  in THIJacobianLocal_Momentum()
1146 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + dp[2] * …  in THIJacobianLocal_Momentum()