| /petsc/src/ksp/ksp/tests/ |
| H A D | ex20.c | 8 PetscErrorCode FormElementStiffness(PetscReal H, PetscScalar *Ke) in FormElementStiffness() argument
10 Ke[0] = H / 6.0; in FormElementStiffness()
11 Ke[1] = -.125 * H; in FormElementStiffness()
12 Ke[2] = H / 12.0; in FormElementStiffness()
13 Ke[3] = -.125 * H; in FormElementStiffness()
14 Ke[4] = -.125 * H; in FormElementStiffness()
15 Ke[5] = H / 6.0; in FormElementStiffness()
16 Ke[6] = -.125 * H; in FormElementStiffness()
17 Ke[7] = H / 12.0; in FormElementStiffness()
18 Ke[8] = H / 12.0; in FormElementStiffness()
[all …]
|
| H A D | ex4.c | 6 PetscErrorCode FormElementStiffness(PetscReal H, PetscScalar *Ke) in FormElementStiffness() argument
8 Ke[0] = H / 6.0; in FormElementStiffness()
9 Ke[1] = -.125 * H; in FormElementStiffness()
10 Ke[2] = H / 12.0; in FormElementStiffness()
11 Ke[3] = -.125 * H; in FormElementStiffness()
12 Ke[4] = -.125 * H; in FormElementStiffness()
13 Ke[5] = H / 6.0; in FormElementStiffness()
14 Ke[6] = -.125 * H; in FormElementStiffness()
15 Ke[7] = H / 12.0; in FormElementStiffness()
16 Ke[8] = H / 12.0; in FormElementStiffness()
[all …]
|
| H A D | ex3.c | 10 PetscErrorCode FormElementStiffness(PetscReal H, PetscScalar *Ke) in FormElementStiffness() argument
13 Ke[0] = H / 6.0; in FormElementStiffness()
14 Ke[1] = -.125 * H; in FormElementStiffness()
15 Ke[2] = H / 12.0; in FormElementStiffness()
16 Ke[3] = -.125 * H; in FormElementStiffness()
17 Ke[4] = -.125 * H; in FormElementStiffness()
18 Ke[5] = H / 6.0; in FormElementStiffness()
19 Ke[6] = -.125 * H; in FormElementStiffness()
20 Ke[7] = H / 12.0; in FormElementStiffness()
21 Ke[8] = H / 12.0; in FormElementStiffness()
[all …]
|
| H A D | ex10.c | 209 PetscErrorCode Elastic20Stiff(PetscReal **Ke) in Elastic20Stiff() argument
286 for (j = 0; j < 81; j++) Ke[i][j] = 0.0; in Elastic20Stiff()
293 for (l = 0; l < 3; l++) Ke[3 * rmap[i] + k][3 * rmap[j] + l] = K[Ii + k][J + l]; in Elastic20Stiff()
302 for (j = 0; j < i; j++) Ke[i][j] = (Ke[i][j] + Ke[j][i]) / 2.0; in Elastic20Stiff()
|
| /petsc/src/ksp/pc/tests/ |
| H A D | ex6.c | 7 static PetscErrorCode FormElementStiffness(PetscReal H, PetscScalar *Ke) in FormElementStiffness() argument
10 Ke[0] = H / 6.0; in FormElementStiffness()
11 Ke[1] = -.125 * H; in FormElementStiffness()
12 Ke[2] = H / 12.0; in FormElementStiffness()
13 Ke[3] = -.125 * H; in FormElementStiffness()
14 Ke[4] = -.125 * H; in FormElementStiffness()
15 Ke[5] = H / 6.0; in FormElementStiffness()
16 Ke[6] = -.125 * H; in FormElementStiffness()
17 Ke[7] = H / 12.0; in FormElementStiffness()
18 Ke[8] = H / 12.0; in FormElementStiffness()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex19.c | 9 PetscErrorCode FormElementStiffness(PetscReal H, PetscScalar *Ke) in FormElementStiffness() argument
12 Ke[0] = H / 6.0; in FormElementStiffness()
13 Ke[1] = -.125 * H; in FormElementStiffness()
14 Ke[2] = H / 12.0; in FormElementStiffness()
15 Ke[3] = -.125 * H; in FormElementStiffness()
16 Ke[4] = -.125 * H; in FormElementStiffness()
17 Ke[5] = H / 6.0; in FormElementStiffness()
18 Ke[6] = -.125 * H; in FormElementStiffness()
19 Ke[7] = H / 12.0; in FormElementStiffness()
20 Ke[8] = H / 12.0; in FormElementStiffness()
[all …]
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex3.c | 29 PetscScalar Ke[16]; /* element matrix */ in main() local
64 PetscCall(FormElementStiffness(h * h, Ke)); in main()
71 PetscCall(MatSetValues(A, 4, idx, 4, idx, Ke, ADD_VALUES)); in main()
186 PetscErrorCode FormElementStiffness(PetscReal H, PetscScalar *Ke) in FormElementStiffness() argument
189 Ke[0] = H / 6.0; in FormElementStiffness()
190 Ke[1] = -.125 * H; in FormElementStiffness()
191 Ke[2] = H / 12.0; in FormElementStiffness()
192 Ke[3] = -.125 * H; in FormElementStiffness()
193 Ke[4] = -.125 * H; in FormElementStiffness()
194 Ke[5] = H / 6.0; in FormElementStiffness()
[all …]
|
| H A D | ex43.c | 491 static void FormStressOperatorQ1(PetscScalar Ke[], PetscScalar coords[], PetscScalar eta[]) in FormStressOperatorQ1() argument
534 Ke[i + 8 * j] = Ke[i + 8 * j] + B[k][i] * tildeD[k] * B[k][j]; in FormStressOperatorQ1()
541 static void FormGradientOperatorQ1(PetscScalar Ke[], PetscScalar coords[]) in FormGradientOperatorQ1() argument
567 Ke[IJ] = Ke[IJ] - GNx_p[di][i] * Ni_p[j] * fac; in FormGradientOperatorQ1()
591 static void FormStabilisationOperatorQ1(PetscScalar Ke[], PetscScalar coords[], PetscScalar eta[]) in FormStabilisationOperatorQ1() argument
612 …for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] = Ke[NODES_PER_EL * i + j] - fac * (Ni… in FormStabilisationOperatorQ1()
622 for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] = fac * Ke[NODES_PER_EL * i + j]; in FormStabilisationOperatorQ1()
626 static void FormScaledMassMatrixOperatorQ1(PetscScalar Ke[], PetscScalar coords[], PetscScalar eta[… in FormScaledMassMatrixOperatorQ1() argument
647 …for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] = Ke[NODES_PER_EL * i + j] - fac * Ni_… in FormScaledMassMatrixOperatorQ1()
657 for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] = fac * Ke[NODES_PER_EL * i + j]; in FormScaledMassMatrixOperatorQ1()
|
| H A D | ex70.c | 143 static void BForm_DivT(PetscScalar Ke[], PetscScalar coords[], PetscScalar eta[]) in BForm_DivT() argument
184 Ke[i + 8 * j] += B[k][i] * tildeD[k] * B[k][j]; in BForm_DivT()
191 static void BForm_Grad(PetscScalar Ke[], PetscScalar coords[]) in BForm_Grad() argument
214 Ke[IJ] -= GNx_p[di][i] * Ni_p[j] * fac; in BForm_Grad()
237 static void BForm_Stabilisation(PetscScalar Ke[], PetscScalar coords[], PetscScalar eta[]) in BForm_Stabilisation() argument
255 … for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] -= fac * (Ni_p[i] * Ni_p[j] - 0.0625); in BForm_Stabilisation()
265 for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] = fac * Ke[NODES_PER_EL * i + j]; in BForm_Stabilisation()
269 static void BForm_ScaledMassMatrix(PetscScalar Ke[], PetscScalar coords[], PetscScalar eta[]) in BForm_ScaledMassMatrix() argument
287 for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] -= fac * Ni_p[i] * Ni_p[j]; in BForm_ScaledMassMatrix()
297 for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] *= fac; in BForm_ScaledMassMatrix()
|
| H A D | ex42.c | 592 static void FormStressOperatorQ13D(PetscScalar Ke[], PetscScalar coords[], PetscScalar eta[]) in FormStressOperatorQ13D() argument
653 for (k = 0; k < 6; k++) Ke[i * nvdof + j] += B[k][i] * tildeD[k] * B[k][j]; in FormStressOperatorQ13D()
660 for (j = i; j < nvdof; j++) Ke[j * nvdof + i] = Ke[i * nvdof + j]; in FormStressOperatorQ13D()
665 static void FormGradientOperatorQ13D(PetscScalar Ke[], PetscScalar coords[]) in FormGradientOperatorQ13D() argument
691 Ke[IJ] = Ke[IJ] - GNx_p[di][i] * Ni_p[j] * fac; in FormGradientOperatorQ13D()
715 static void FormStabilisationOperatorQ13D(PetscScalar Ke[], PetscScalar coords[], PetscScalar eta[]) in FormStabilisationOperatorQ13D() argument
743 …for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] -= fac * (Ni_p[i] * Ni_p[j] - 0.015625… in FormStabilisationOperatorQ13D()
765 for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] = fac * Ke[NODES_PER_EL * i + j]; in FormStabilisationOperatorQ13D()
769 static void FormScaledMassMatrixOperatorQ13D(PetscScalar Ke[], PetscScalar coords[], PetscScalar et… in FormScaledMassMatrixOperatorQ13D() argument
798 …for (j = 0; j < NODES_PER_EL; j++) Ke[NODES_PER_EL * i + j] = Ke[NODES_PER_EL * i + j] - fac * Ni_… in FormScaledMassMatrixOperatorQ13D()
[all …]
|
| H A D | ex49.c | 417 static void FormStressOperatorQ1(PetscScalar Ke[], PetscScalar coords[], PetscScalar E[], PetscScal… in FormStressOperatorQ1() argument
473 … for (l = 0; l < 3; l++) Ke[8 * i + j] = Ke[8 * i + j] + B[k][i] * constit_D[k][l] * B[l][j]; in FormStressOperatorQ1()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex48.c | 986 PetscScalar Ke[4 * 2][4 * 2]; in THIJacobianLocal_2D() local
989 PetscCall(PetscMemzero(Ke, sizeof(Ke))); 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()
[all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex14.c | 1116 PetscScalar Ke[8 * NODE_SIZE][8 * NODE_SIZE], Kcpl[8 * NODE_SIZE][4 * PRMNODE_SIZE]; in THIJacobianLocal_Momentum() local
1121 PetscCall(PetscMemzero(Ke, sizeof(Ke))); in THIJacobianLocal_Momentum()
1143 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + dp[2] * … 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 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + dp[2] * … in THIJacobianLocal_Momentum()
1148 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d… in THIJacobianLocal_Momentum()
1149 …Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d… in THIJacobianLocal_Momentum()
1150 …Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * d… in THIJacobianLocal_Momentum()
1151 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * d… in THIJacobianLocal_Momentum()
[all …]
|
| /petsc/doc/ |
| H A D | petsc.bib | 1278 author = {DENG Lian-wen and HUANG Xiao-zhong and ZH0U Ke-sheng and YANG},
|