Searched refs:dpl (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/snes/tutorials/ |
| H A D | ex48.c | 1006 const PetscReal ppl = phi[ll], *dpl = dphi[ll]; in THIJacobianLocal_2D() local 1008 dgdu = 2. * du[0] * dpl[0] + dv[1] * dpl[0] + 0.5 * (du[1] + dv[0]) * dpl[1]; in THIJacobianLocal_2D() 1009 dgdv = 2. * dv[1] * dpl[1] + du[0] * dpl[1] + 0.5 * (du[1] + dv[0]) * dpl[0]; 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() 1104 const PetscReal *PETSC_RESTRICT dpl = dphi[ll]; local 1110 …dgdu = 2. * du[0] * dpl[0] + dv[1] * dpl[0] + 0.5 * (du[1] + dv[0]) * dpl[1] + 0.5 * du[2] * dpl[2… 1111 …dgdv = 2. * dv[1] * dpl[1] + du[0] * dpl[1] + 0.5 * (du[1] + dv[0]) * dpl[0] + 0.5 * dv[2] * dpl[2… [all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex14.c | 1138 const PetscReal *restrict dpl = dphi[ll]; in THIJacobianLocal_Momentum() local 1140 …dgdu = 2. * du[0] * dpl[0] + dv[1] * dpl[0] + 0.5 * (du[1] + dv[0]) * dpl[1] + 0.5 * du[2] * dpl[2… in THIJacobianLocal_Momentum() 1141 …dgdv = 2. * dv[1] * dpl[1] + du[0] * dpl[1] + 0.5 * (du[1] + dv[0]) * dpl[0] + 0.5 * dv[2] * dpl[2… in THIJacobianLocal_Momentum() 1143 … 2 + 0][ll * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + dp[2] * jw * e… in THIJacobianLocal_Momentum() 1144 … Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * eta * 2. * dpl[1] + dp[1] * jw * eta * dpl[0]; in THIJacobianLocal_Momentum() 1145 … Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * eta * 2. * dpl[0] + dp[0] * jw * eta * dpl[1]; in THIJacobianLocal_Momentum() 1146 … 2 + 1][ll * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + dp[2] * jw * e… in THIJacobianLocal_Momentum() 1157 …const PetscReal dpl[] = {QuadQDeriv[q % 4][ll][0] / hx, QuadQDeriv[q % 4][ll][1] / hy}; /* surface… in THIJacobianLocal_Momentum() local 1158 … Kcpl[FieldIndex(Node, l, u)][FieldIndex(PrmNode, ll, h)] += pp * jw * thi->rhog * dpl[0]; in THIJacobianLocal_Momentum() 1159 … Kcpl[FieldIndex(Node, l, u)][FieldIndex(PrmNode, ll, b)] += pp * jw * thi->rhog * dpl[0]; in THIJacobianLocal_Momentum() [all …]
|