Home
last modified time | relevance | path

Searched refs:NSD (Results 1 – 4 of 4) sorted by relevance

/petsc/src/ksp/ksp/tutorials/
H A Dex42.c18 #define NSD 3 /* number of spatial dimensions */ macro
26 PetscScalar gp_coords[NSD * GAUSS_POINTS];
189 PetscScalar cx = coords[NSD * n + 0]; in ShapeFunctionQ13D_Evaluate_dx()
190 PetscScalar cy = coords[NSD * n + 1]; in ShapeFunctionQ13D_Evaluate_dx()
191 PetscScalar cz = coords[NSD * n + 2]; in ShapeFunctionQ13D_Evaluate_dx()
228 static void ConstructGaussQuadrature3D(PetscInt *ngp, PetscScalar gp_xi[][NSD], PetscScalar gp_weig… in ConstructGaussQuadrature3D() argument
566 II = u_eqn[NSD * n].i; in DMDASetValuesLocalStencil3D_ADD_VALUES()
567 J = u_eqn[NSD * n].j; in DMDASetValuesLocalStencil3D_ADD_VALUES()
568 K = u_eqn[NSD * n].k; in DMDASetValuesLocalStencil3D_ADD_VALUES()
570 fields_F[K][J][II].u_dof = fields_F[K][J][II].u_dof + Fe_u[NSD * n]; in DMDASetValuesLocalStencil3D_ADD_VALUES()
[all …]
H A Dex70.c38 #define NSD 2 /* number of spatial dimensions */ macro
124 s_u[NSD * i + 0] = 3 * element[i]; in DMDAGetElementEqnums_up()
125 s_u[NSD * i + 1] = 3 * element[i] + 1; in DMDAGetElementEqnums_up()
145 PetscScalar gp_xi[GAUSS_POINTS][NSD], gp_weight[GAUSS_POINTS]; in BForm_DivT()
146 PetscScalar GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in BForm_DivT()
193 PetscScalar gp_xi[GAUSS_POINTS][NSD], gp_weight[GAUSS_POINTS]; in BForm_Grad()
194 PetscScalar Ni_p[NODES_PER_EL], GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in BForm_Grad()
209 for (di = 0; di < NSD; di++) { /* u dofs */ in BForm_Grad()
239 PetscScalar gp_xi[GAUSS_POINTS][NSD], gp_weight[GAUSS_POINTS]; in BForm_Stabilisation()
240 PetscScalar Ni_p[NODES_PER_EL], GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in BForm_Stabilisation()
[all …]
H A Dex43.c63 #define NSD 2 /* number of spatial dimensions */ macro
497 PetscScalar GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in FormStressOperatorQ1()
548 PetscScalar GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in FormGradientOperatorQ1()
562 for (di = 0; di < NSD; di++) { /* u dofs */ in FormGradientOperatorQ1()
598 PetscScalar GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in FormStabilisationOperatorQ1()
633 PetscScalar GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in FormScaledMassMatrixOperatorQ1()
668 PetscScalar GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in FormMomentumRhsQ1()
682 Fe[NSD * i] += fac * Ni_p[i] * fx[p]; in FormMomentumRhsQ1()
683 Fe[NSD * i + 1] += fac * Ni_p[i] * fy[p]; in FormMomentumRhsQ1()
692 el_coords[NSD * 0] = _coords[ej][ei].x; in GetElementCoords()
[all …]
H A Dex49.c77 #define NSD 2 /* number of spatial dimensions */ macro
423 PetscScalar GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in FormStressOperatorQ1()
488 PetscScalar GNi_p[NSD][NODES_PER_EL], GNx_p[NSD][NODES_PER_EL]; in FormMomentumRhsQ1()
502 Fe[NSD * i] += fac * Ni_p[i] * fx[p]; in FormMomentumRhsQ1()
503 Fe[NSD * i + 1] += fac * Ni_p[i] * fy[p]; in FormMomentumRhsQ1()
555 el_coords[NSD * 0 + 0] = _coords[ej][ei].x; in GetElementCoords()
556 el_coords[NSD * 0 + 1] = _coords[ej][ei].y; in GetElementCoords()
557 el_coords[NSD * 1 + 0] = _coords[ej + 1][ei].x; in GetElementCoords()
558 el_coords[NSD * 1 + 1] = _coords[ej + 1][ei].y; in GetElementCoords()
559 el_coords[NSD * 2 + 0] = _coords[ej + 1][ei + 1].x; in GetElementCoords()
[all …]