Lines Matching refs:deta
374 static void THIViscosity(THI thi, PetscReal gam, PetscReal *eta, PetscReal *deta) in THIViscosity() argument
391 *deta = exponent * (*eta) / (eps + gam); in THIViscosity()
551 …ts->meter / (3e7 * units->second), gradu = u / (100 * units->meter), eta, deta, rho = 910 * units-… in THICreate() local
553 THIViscosity(thi, 0.5 * gradu * gradu, &eta, &deta); in THICreate()
559 THIViscosity(thi, 0.5 * PetscSqr(1e-3 * gradu), &eta, &deta); 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()
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()
1016 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d… in THIJacobianLocal_2D()
1017 …Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d… in THIJacobianLocal_2D()
1018 …Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * d… in THIJacobianLocal_2D()
1019 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * d… 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);
1091 …_set1_pd(dv[2]), jweta = _mm_set1_pd(jw * eta), jwdeta = _mm_set1_pd(jw * deta), dp0 = _mm_set1_pd…
1118 … 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu * (du[1…
1119 … 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv * (du[1…
1120 … 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu * (du[1…
1121 … 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv * (du[1…