Searched refs:Ke1D (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex51.c | 31 PetscScalar *Ke1D, *Ke2D, *Me1D, *Me2D; in main() local
75 PetscCall(PetscMalloc1(num1Dnodes * num1Dnodes, &Ke1D)); in main()
87 PetscCall(Form1DElementStiffness(h, p, gllNode, gllWgts, Ke1D)); in main()
89 PetscCall(Form2DElementStiffness(p, Ke1D, Me1D, Ke2D)); in main()
106 PetscCall(PetscFree(Ke1D)); in main()
244 … Form1DElementStiffness(PetscReal H, PetscInt P, PetscReal *gqn, PetscReal *gqw, PetscScalar *Ke1D) in Form1DElementStiffness() argument
253 Ke1D[indx] = 0.0; in Form1DElementStiffness()
254 …for (k = 0; k < (P + 1); ++k) Ke1D[indx] += (1. / H) * gqw[k] * derivPolyBasisFunc(P, i, gqn, gqn[… in Form1DElementStiffness()
287 static PetscErrorCode Form2DElementStiffness(PetscInt P, PetscScalar *Ke1D, PetscScalar *Me1D, Pets… in Form2DElementStiffness() argument
300 Ke2D[indx3] = Ke1D[indx1] * Me1D[indx2] + Me1D[indx1] * Ke1D[indx2]; in Form2DElementStiffness()
|