| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex36.cxx | 265 PetscReal *phi, *phypts, *jxw; in ComputeRHS_MOAB() local
277 PetscCall(PetscMalloc3(user->VPERE * npoints, &phi, npoints * 3, &phypts, npoints, &jxw)); in ComputeRHS_MOAB()
299 PetscCall(DMMoabFEMComputeBasis(user->dim, nconn, vpos, quadratureObj, phypts, jxw, phi, NULL)); in ComputeRHS_MOAB()
307 for (i = 0; i < nconn; ++i) localv[i] += jxw[q] * phi[offset + i] * ff; in ComputeRHS_MOAB()
343 PetscCall(PetscFree3(phi, phypts, jxw)); in ComputeRHS_MOAB()
361 PetscReal *phi, *dphi[3], *phypts, *jxw; in ComputeMatrix_MOAB() local
376 …>VPERE * npoints, &dphi[1], user->VPERE * npoints, &dphi[2], npoints * 3, &phypts, npoints, &jxw)); in ComputeMatrix_MOAB()
394 PetscCall(DMMoabFEMComputeBasis(user->dim, nconn, vpos, quadratureObj, phypts, jxw, phi, dphi)); in ComputeMatrix_MOAB()
409 …array[i * nconn + j] += jxw[q] * (rho * (dphi[0][offset + i] * dphi[0][offset + j] + dphi[1][offse… in ComputeMatrix_MOAB()
448 PetscCall(PetscFree6(phi, dphi[0], dphi[1], dphi[2], phypts, jxw)); in ComputeMatrix_MOAB()
|
| H A D | ex35.cxx | 294 PetscReal *phi, *phypts, *jxw; in ComputeRHS() local
306 PetscCall(PetscMalloc3(user->VPERE * npoints, &phi, npoints * 3, &phypts, npoints, &jxw)); in ComputeRHS()
334 PetscCall(DMMoabFEMComputeBasis(2, nconn, vpos, quadratureObj, phypts, jxw, phi, NULL)); in ComputeRHS()
340 for (i = 0; i < nconn; ++i) localv[i] += jxw[q] * phi[q * nconn + i] * ff; in ComputeRHS()
379 PetscCall(PetscFree3(phi, phypts, jxw)); in ComputeRHS()
397 PetscReal *phi, *dphi[2], *phypts, *jxw; in ComputeMatrix() local
412 …>VPERE * npoints, &dphi[0], user->VPERE * npoints, &dphi[1], npoints * 3, &phypts, npoints, &jxw)); in ComputeMatrix()
434 PetscCall(DMMoabFEMComputeBasis(2, nconn, vpos, quadratureObj, phypts, jxw, phi, dphi)); in ComputeMatrix()
446 …for (j = 0; j < nconn; ++j) array[i * nconn + j] += jxw[q] * rho * (dphi[0][q * nconn + i] * dphi[… in ComputeMatrix()
488 PetscCall(PetscFree5(phi, dphi[0], dphi[1], phypts, jxw)); in ComputeMatrix()
|
| /petsc/src/dm/impls/moab/ |
| H A D | dmmbfem.cxx | 110 … const PetscInt npts, const PetscReal *quad, PetscReal *phypts, PetscReal *jxw, PetscReal *phi, Pe… in Compute_Lagrange_Basis_1D_Internal() argument
143 jxw[j] *= *volume; in Compute_Lagrange_Basis_1D_Internal()
172 if (jxw) jxw[j] *= *volume; in Compute_Lagrange_Basis_1D_Internal()
220 … const PetscInt npts, const PetscReal *quad, PetscReal *phypts, PetscReal *jxw, PetscReal *phi, Pe… in Compute_Lagrange_Basis_2D_Internal() argument
268 if (jxw) jxw[j] *= *volume; in Compute_Lagrange_Basis_2D_Internal()
304 if (jxw) jxw[j] *= 0.5; in Compute_Lagrange_Basis_2D_Internal()
376 … const PetscInt npts, const PetscReal *quad, PetscReal *phypts, PetscReal *jxw, PetscReal *phi, Pe… in Compute_Lagrange_Basis_3D_Internal() argument
440 if (jxw) jxw[j] *= (*volume); in Compute_Lagrange_Basis_3D_Internal()
494 if (jxw) jxw[j] *= *volume; in Compute_Lagrange_Basis_3D_Internal()
|