Lines Matching refs:jacobian
110 … *phypts, PetscReal *jxw, PetscReal *phi, PetscReal *dphidx, PetscReal *jacobian, PetscReal *ijaco… in Compute_Lagrange_Basis_1D_Internal() argument
115 PetscAssertPointer(jacobian, 9); in Compute_Lagrange_Basis_1D_Internal()
133 jacobian[0] = ijacobian[0] = volume[0] = 0.0; in Compute_Lagrange_Basis_1D_Internal()
136 jacobian[0] += dNi_dxi[i] * vertices[0]; in Compute_Lagrange_Basis_1D_Internal()
141 *volume = jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
142 ijacobian[0] = 1.0 / jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
162 jacobian[0] = ijacobian[0] = volume[0] = 0.0; in Compute_Lagrange_Basis_1D_Internal()
165 jacobian[0] += dNi_dxi[i] * vertices[0]; in Compute_Lagrange_Basis_1D_Internal()
170 *volume = jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
171 ijacobian[0] = 1.0 / jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
220 … *jxw, PetscReal *phi, PetscReal *dphidx, PetscReal *dphidy, PetscReal *jacobian, PetscReal *ijaco… in Compute_Lagrange_Basis_2D_Internal() argument
225 PetscAssertPointer(jacobian, 10); in Compute_Lagrange_Basis_2D_Internal()
249 PetscCall(PetscArrayzero(jacobian, 4)); in Compute_Lagrange_Basis_2D_Internal()
253 jacobian[0] += dNi_dxi[i] * vertices[0]; in Compute_Lagrange_Basis_2D_Internal()
254 jacobian[2] += dNi_dxi[i] * vertices[1]; in Compute_Lagrange_Basis_2D_Internal()
255 jacobian[1] += dNi_deta[i] * vertices[0]; in Compute_Lagrange_Basis_2D_Internal()
256 jacobian[3] += dNi_deta[i] * vertices[1]; in Compute_Lagrange_Basis_2D_Internal()
265 PetscCall(DMatrix_Invert_2x2_Internal(jacobian, ijacobian, volume)); in Compute_Lagrange_Basis_2D_Internal()
283 PetscCall(PetscArrayzero(jacobian, 4)); in Compute_Lagrange_Basis_2D_Internal()
287 jacobian[0] = (coords[0 * 3 + 0] - x2); in Compute_Lagrange_Basis_2D_Internal()
288 jacobian[1] = (coords[1 * 3 + 0] - x2); in Compute_Lagrange_Basis_2D_Internal()
289 jacobian[2] = (coords[0 * 3 + 1] - y2); in Compute_Lagrange_Basis_2D_Internal()
290 jacobian[3] = (coords[1 * 3 + 1] - y2); in Compute_Lagrange_Basis_2D_Internal()
293 PetscCall(DMatrix_Invert_2x2_Internal(jacobian, ijacobian, volume)); in Compute_Lagrange_Basis_2D_Internal()
307 const PetscReal phipts_x = coords[2 * 3 + 0] + jacobian[0] * r + jacobian[1] * s; in Compute_Lagrange_Basis_2D_Internal()
308 const PetscReal phipts_y = coords[2 * 3 + 1] + jacobian[2] * r + jacobian[3] * s; in Compute_Lagrange_Basis_2D_Internal()
376 …hi, PetscReal *dphidx, PetscReal *dphidy, PetscReal *dphidz, PetscReal *jacobian, PetscReal *ijaco… in Compute_Lagrange_Basis_3D_Internal() argument
381 PetscAssertPointer(jacobian, 11); in Compute_Lagrange_Basis_3D_Internal()
416 PetscCall(PetscArrayzero(jacobian, 9)); in Compute_Lagrange_Basis_3D_Internal()
420 jacobian[0] += dNi_dxi[i] * vertex[0]; in Compute_Lagrange_Basis_3D_Internal()
421 jacobian[3] += dNi_dxi[i] * vertex[1]; in Compute_Lagrange_Basis_3D_Internal()
422 jacobian[6] += dNi_dxi[i] * vertex[2]; in Compute_Lagrange_Basis_3D_Internal()
423 jacobian[1] += dNi_deta[i] * vertex[0]; in Compute_Lagrange_Basis_3D_Internal()
424 jacobian[4] += dNi_deta[i] * vertex[1]; in Compute_Lagrange_Basis_3D_Internal()
425 jacobian[7] += dNi_deta[i] * vertex[2]; in Compute_Lagrange_Basis_3D_Internal()
426 jacobian[2] += dNi_dzeta[i] * vertex[0]; in Compute_Lagrange_Basis_3D_Internal()
427 jacobian[5] += dNi_dzeta[i] * vertex[1]; in Compute_Lagrange_Basis_3D_Internal()
428 jacobian[8] += dNi_dzeta[i] * vertex[2]; in Compute_Lagrange_Basis_3D_Internal()
437 PetscCall(DMatrix_Invert_3x3_Internal(jacobian, ijacobian, volume)); in Compute_Lagrange_Basis_3D_Internal()
455 PetscCall(PetscArrayzero(jacobian, 9)); in Compute_Lagrange_Basis_3D_Internal()
459 jacobian[0] = coords[1 * 3 + 0] - x0; in Compute_Lagrange_Basis_3D_Internal()
460 jacobian[1] = coords[2 * 3 + 0] - x0; in Compute_Lagrange_Basis_3D_Internal()
461 jacobian[2] = coords[3 * 3 + 0] - x0; in Compute_Lagrange_Basis_3D_Internal()
462 jacobian[3] = coords[1 * 3 + 1] - y0; in Compute_Lagrange_Basis_3D_Internal()
463 jacobian[4] = coords[2 * 3 + 1] - y0; in Compute_Lagrange_Basis_3D_Internal()
464 jacobian[5] = coords[3 * 3 + 1] - y0; in Compute_Lagrange_Basis_3D_Internal()
465 jacobian[6] = coords[1 * 3 + 2] - z0; in Compute_Lagrange_Basis_3D_Internal()
466 jacobian[7] = coords[2 * 3 + 2] - z0; in Compute_Lagrange_Basis_3D_Internal()
467 jacobian[8] = coords[3 * 3 + 2] - z0; in Compute_Lagrange_Basis_3D_Internal()
470 PetscCall(DMatrix_Invert_3x3_Internal(jacobian, ijacobian, volume)); in Compute_Lagrange_Basis_3D_Internal()
560 PetscReal jacobian[9], ijacobian[9], volume; in DMMoabFEMComputeBasis() local
575 …weight_product, fe_basis, compute_der ? fe_basis_derivatives[0] : NULL, jacobian, ijacobian, &volu… in DMMoabFEMComputeBasis()
578 …is_derivatives[0] : NULL, compute_der ? fe_basis_derivatives[1] : NULL, jacobian, ijacobian, &volu… in DMMoabFEMComputeBasis()
581 …is_derivatives[1] : NULL, compute_der ? fe_basis_derivatives[2] : NULL, jacobian, ijacobian, &volu… in DMMoabFEMComputeBasis()
678 …tscReal *quadrature, PetscReal *phypts, PetscReal *phibasis, PetscReal *jacobian, PetscReal *ijaco… in FEMComputeBasis_JandF() argument
683 …ernal(nverts, coordinates, 1, quadrature, phypts, NULL, phibasis, NULL, jacobian, ijacobian, volum… in FEMComputeBasis_JandF()
686 …nverts, coordinates, 1, quadrature, phypts, NULL, phibasis, NULL, NULL, jacobian, ijacobian, volum… in FEMComputeBasis_JandF()
689 …, coordinates, 1, quadrature, phypts, NULL, phibasis, NULL, NULL, NULL, jacobian, ijacobian, volum… in FEMComputeBasis_JandF()
725 PetscReal phibasis[8], jacobian[9], ijacobian[9], volume; in DMMoabPToRMapping() local
736 PetscCall(PetscArrayzero(jacobian, dim * dim)); in DMMoabPToRMapping()
744 …PetscCall(FEMComputeBasis_JandF(dim, nverts, coordinates, natparam, phypts, phibasis, jacobian, ij… in DMMoabPToRMapping()
780 …PetscCall(FEMComputeBasis_JandF(dim, nverts, coordinates, natparam, phypts, phibasis, jacobian, ij… in DMMoabPToRMapping()