Lines Matching refs:ll
712 PetscInt l, ll; in PointwiseNonlinearity() local
722 for (ll = 0; ll < 3; ll++) { 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
743 for (ll = 0; ll < 2; ll++) { in PointwiseNonlinearity2D()
744 du[ll] += dphi[l][ll] * n[l].u; in PointwiseNonlinearity2D()
745 dv[ll] += dphi[l][ll] * n[l].v; in PointwiseNonlinearity2D()
967 PetscInt xs, ys, xm, ym, i, j, q, l, ll; in THIJacobianLocal_2D() local
1005 for (ll = 0; ll < 4; ll++) { in THIJacobianLocal_2D()
1006 const PetscReal ppl = phi[ll], *dpl = dphi[ll]; 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
1100 for (ll = ls; ll < 8; ll++) { /* trial functions */
1102 for (ll = l; ll < 8; ll++) {
1104 const PetscReal *PETSC_RESTRICT dpl = dphi[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…
1160 for (ll = 0; ll < 4; ll++) {
1161 const PetscReal ppl = phi[ll];
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;
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];