204 …T], PetscReal phi[PETSC_RESTRICT], PetscReal dphi[PETSC_RESTRICT][3], PetscReal *PETSC_RESTRICT jw)  in HexComputeGeometry()  argument
226   *jw = 1.0 * jdet;  in HexComputeGeometry()
791           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  in THIFunctionLocal()  local
794           HexComputeGeometry(q, hx, hy, dz, phi, dphi, &jw);  in THIFunctionLocal()
796           jw /= thi->rhog; /* scales residuals to be O(1) */  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()
827               const PetscReal jw = 0.25 * hx * hy / thi->rhog, *phi = QuadQInterp[q];  in THIFunctionLocal()  local
839                 fn[ls + l]->u += pp * jw * beta2 * u;  in THIFunctionLocal()
840                 fn[ls + l]->v += pp * jw * beta2 * v;  in THIFunctionLocal()
991         PetscReal   phi[4], dphi[4][2], jw, eta, deta, beta2, dbeta2;  in THIJacobianLocal_2D()  local
1000         jw = 0.25 * hx * hy / thi->rhog; /* rhog is only scaling */  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()
1016 … 0][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu…  in THIJacobianLocal_2D()
1017 … 0][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv…  in THIJacobianLocal_2D()
1018 … 1][ll * 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu…  in THIJacobianLocal_2D()
1019 … 1][ll * 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv…  in THIJacobianLocal_2D()
1080           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  local
1083           HexComputeGeometry(q, hx, hy, dz, phi, dphi, &jw);
1085           jw /= thi->rhog; /* residuals are scaled by this factor */
1091 …1_pd(dv[1]), dv2 = _mm_set1_pd(dv[2]), jweta = _mm_set1_pd(jw * eta), jwdeta = _mm_set1_pd(jw * de…
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] * …
1118 …][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu *…
1119 …][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv *…
1120 …][ll * 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu *…
1121 …][ll * 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv *…
1149               const PetscReal jw = 0.25 * hx * hy / thi->rhog, *phi = QuadQInterp[q];  local
1162 …               Ke[l * 2 + 0][ll * 2 + 0] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * u * u * ppl;
1163                   Ke[l * 2 + 0][ll * 2 + 1] += pp * jw * dbeta2 * u * v * ppl;
1164                   Ke[l * 2 + 1][ll * 2 + 0] += pp * jw * dbeta2 * v * u * ppl;
1165 …               Ke[l * 2 + 1][ll * 2 + 1] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * v * v * ppl;