424 static void THIViscosity(THI thi, PetscReal gam, PetscReal *eta, PetscReal *deta)  in THIViscosity()  argument
441   *deta    = exponent * (*eta) / (eps + gam);  in THIViscosity()
630 …ts->meter / (3e7 * units->second), gradu = u / (100 * units->meter), eta, deta, rho = 910 * units-…  in THICreate()  local
632     THIViscosity(thi, 0.5 * gradu * gradu, &eta, &deta);  in THICreate()
638       THIViscosity(thi, 0.5 * PetscSqr(1e-3 * gradu), &eta, &deta);  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()
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()
1148 … 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu * (du[1…  in THIJacobianLocal_Momentum()
1149 … 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv * (du[1…  in THIJacobianLocal_Momentum()
1150 … 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu * (du[1…  in THIJacobianLocal_Momentum()
1151 … 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv * (du[1…  in THIJacobianLocal_Momentum()