Lines Matching refs:l
712 PetscInt l, ll; in PointwiseNonlinearity() local
719 for (l = 0; l < 8; l++) { in PointwiseNonlinearity()
720 *u += phi[l] * n[l].u; in PointwiseNonlinearity()
721 *v += phi[l] * n[l].v; in PointwiseNonlinearity()
723 du[ll] += dphi[l][ll] * n[l].u; in PointwiseNonlinearity()
724 dv[ll] += dphi[l][ll] * n[l].v; in PointwiseNonlinearity()
733 PetscInt l, ll; in PointwiseNonlinearity2D() local
740 for (l = 0; l < 4; l++) { in PointwiseNonlinearity2D()
741 *u += phi[l] * n[l].u; in PointwiseNonlinearity2D()
742 *v += phi[l] * n[l].v; in PointwiseNonlinearity2D()
744 du[ll] += dphi[l][ll] * n[l].u; in PointwiseNonlinearity2D()
745 dv[ll] += dphi[l][ll] * n[l].v; in PointwiseNonlinearity2D()
754 PetscInt xs, ys, xm, ym, zm, i, j, k, q, l; in THIFunctionLocal() local
786 for (l = 0; l < 4; l++) n[l].u = n[l].v = 0; in THIFunctionLocal()
799 for (l = ls; l < 8; l++) { /* test functions */ in THIFunctionLocal()
801 const PetscReal pp = phi[l], *dp = dphi[l]; 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()
830 for (l = 0; l < 4; l++) { in THIFunctionLocal()
831 u += phi[l] * n[l].u; in THIFunctionLocal()
832 v += phi[l] * n[l].v; in THIFunctionLocal()
833 rbeta2 += phi[l] * pn[l].beta2; in THIFunctionLocal()
837 for (l = 0; l < 4; l++) { in THIFunctionLocal()
838 const PetscReal pp = phi[l]; in THIFunctionLocal()
839 fn[ls + l]->u += pp * jw * beta2 * u; in THIFunctionLocal()
840 fn[ls + l]->v += pp * jw * beta2 * v; in THIFunctionLocal()
967 PetscInt xs, ys, xm, ym, i, j, q, l, ll; in THIJacobianLocal_2D() local
993 for (l = 0; l < 4; l++) { in THIJacobianLocal_2D()
994 phi[l] = QuadQInterp[q][l]; in THIJacobianLocal_2D()
995 dphi[l][0] = QuadQDeriv[q][l][0] * 2. / hx; in THIJacobianLocal_2D()
996 dphi[l][1] = QuadQDeriv[q][l][1] * 2. / hy; in THIJacobianLocal_2D()
997 h += phi[l] * pn[l].h; in THIJacobianLocal_2D()
998 rbeta2 += phi[l] * pn[l].beta2; in THIJacobianLocal_2D()
1003 for (l = 0; l < 4; l++) { in THIJacobianLocal_2D()
1004 const PetscReal pp = phi[l], *dp = dphi[l]; 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 …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()
1045 PetscInt xs, ys, xm, ym, zm, i, j, k, q, l, ll; local
1076 for (l = 0; l < 4; l++) n[l].u = n[l].v = 0;
1087 for (l = ls; l < 8; l++) { /* test functions */
1088 const PetscReal *PETSC_RESTRICT dp = dphi[l];
1102 for (ll = l; ll < 8; ll++) {
1105 …if (amode == THIASSEMBLY_TRIDIAGONAL && (l - ll) % 4) continue; /* these entries would not be inse…
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 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d…
1119 …Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d…
1120 …Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * d…
1121 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * d…
1128 …__m128d keu = _mm_loadu_pd(&Ke[l * 2 + 0][ll * 2 + 0]), kev = _mm_loadu_pd(&Ke[l * 2 + 1][ll * 2 +…
1134 …_mm_storeu_pd(&Ke[l * 2 + 0][ll * 2 + 0], _mm_add_pd(keu, _mm_add_pd(_mm_mul_pd(t1, t0), _mm_add_p…
1135 …_mm_storeu_pd(&Ke[l * 2 + 1][ll * 2 + 0], _mm_add_pd(kev, _mm_add_pd(_mm_mul_pd(t2, t0), _mm_add_p…
1152 for (l = 0; l < 4; l++) {
1153 u += phi[l] * n[l].u;
1154 v += phi[l] * n[l].v;
1155 rbeta2 += phi[l] * pn[l].beta2;
1158 for (l = 0; l < 4; l++) {
1159 const PetscReal pp = phi[l];
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;
1183 … for (l = 0; l < 4; l++) { /* Copy out each of the blocks, discarding horizontal coupling */
1184 const PetscInt l4 = l + 4;
1186 {rc[l].k, rc[l].j, rc[l].i, 0},
1191 …{Ke[2 * l + 0][2 * l + 0], Ke[2 * l + 0][2 * l + 1], Ke[2 * l + 0][2 * l4 + 0], Ke[2 * l + 0][2…
1192 …{Ke[2 * l + 1][2 * l + 0], Ke[2 * l + 1][2 * l + 1], Ke[2 * l + 1][2 * l4 + 0], Ke[2 * l + 1][2…
1193 …{Ke[2 * l4 + 0][2 * l + 0], Ke[2 * l4 + 0][2 * l + 1], Ke[2 * l4 + 0][2 * l4 + 0], Ke[2 * l4 + 0][…
1194 …{Ke[2 * l4 + 1][2 * l + 0], Ke[2 * l4 + 1][2 * l + 1], Ke[2 * l4 + 1][2 * l4 + 0], Ke[2 * l4 + 1][…
1199 …{Ke[2 * l + 0][2 * l + 0], Ke[2 * l + 0][2 * l + 1], Ke[2 * l + 0][2 * l4 + 0], Ke[2 * l + 0][2…
1200 …{Ke[2 * l + 1][2 * l + 0], Ke[2 * l + 1][2 * l + 1], Ke[2 * l + 1][2 * l4 + 0], Ke[2 * l + 1][2…
1201 …{Ke[2 * l + 0][2 * l4 + 0], Ke[2 * l + 1][2 * l4 + 0], Ke[2 * l4 + 0][2 * l4 + 0], Ke[2 * l4 + 0][…
1202 …{Ke[2 * l + 0][2 * l4 + 1], Ke[2 * l + 1][2 * l4 + 1], Ke[2 * l4 + 1][2 * l4 + 0], Ke[2 * l4 + 1][…
1209 for (l = 0; l < 8; l++) {
1210 for (ll = l + 1; ll < 8; ll++) {
1211 Ke[ll * 2 + 0][l * 2 + 0] = Ke[l * 2 + 0][ll * 2 + 0];
1212 Ke[ll * 2 + 1][l * 2 + 0] = Ke[l * 2 + 0][ll * 2 + 1];
1213 Ke[ll * 2 + 0][l * 2 + 1] = Ke[l * 2 + 1][ll * 2 + 0];
1214 Ke[ll * 2 + 1][l * 2 + 1] = Ke[l * 2 + 1][ll * 2 + 1];