| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex35.cxx | 289 PetscInt i, q, nconn, nc, npoints; in ComputeRHS() local
317 PetscCall(DMMoabGetElementConnectivity(dm, ehandle, &nconn, &connect)); in ComputeRHS()
318 …nconn == 3 || nconn == 4, PETSC_COMM_WORLD, PETSC_ERR_ARG_WRONG, "Only TRI3/QUAD4 element bases ar… in ComputeRHS()
320 PetscCall(PetscArrayzero(localv, nconn)); in ComputeRHS()
323 PetscCall(DMMoabGetVertexCoordinates(dm, nconn, connect, vpos)); in ComputeRHS()
327 PetscCall(DMMoabGetFieldDofsLocal(dm, nconn, connect, 0, dof_indices)); in ComputeRHS()
329 PetscCall(DMMoabGetFieldDofs(dm, nconn, connect, 0, dof_indices)); 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()
349 PetscCall(DMMoabCheckBoundaryVertices(dm, nconn, connect, dbdry)); in ComputeRHS()
[all …]
|
| H A D | ex36.cxx | 260 PetscInt i, q, nconn; in ComputeRHS_MOAB() local
288 PetscCall(DMMoabGetElementConnectivity(dm, ehandle, &nconn, &connect)); in ComputeRHS_MOAB()
289 …nconn == 4 || nconn == 8, PETSC_COMM_WORLD, PETSC_ERR_ARG_WRONG, "Only HEX8/TET4 element bases are… in ComputeRHS_MOAB()
292 PetscCall(DMMoabGetVertexCoordinates(dm, nconn, connect, vpos)); in ComputeRHS_MOAB()
295 PetscCall(DMMoabGetFieldDofsLocal(dm, nconn, connect, 0, dof_indices)); in ComputeRHS_MOAB()
299 PetscCall(DMMoabFEMComputeBasis(user->dim, nconn, vpos, quadratureObj, phypts, jxw, phi, NULL)); in ComputeRHS_MOAB()
301 PetscCall(PetscArrayzero(localv, nconn)); in ComputeRHS_MOAB()
305 const PetscInt offset = q * nconn; in ComputeRHS_MOAB()
307 for (i = 0; i < nconn; ++i) localv[i] += jxw[q] * phi[offset + i] * ff; in ComputeRHS_MOAB()
316 PetscCall(DMMoabCheckBoundaryVertices(dm, nconn, connect, dbdry)); in ComputeRHS_MOAB()
[all …]
|
| /petsc/src/dm/impls/moab/ |
| H A D | dmmoab.cxx | 663 PetscErrorCode DMMoabGetVertexCoordinates(DM dm, PetscInt nconn, const moab::EntityHandle *conn, Pe… in DMMoabGetVertexCoordinates() argument
676 …merr = dmmoab->hierarchy->get_coordinates(const_cast<moab::EntityHandle *>(conn), nconn, dmmoab->h… in DMMoabGetVertexCoordinates()
679 merr = dmmoab->mbiface->get_coords(conn, nconn, vpos); in DMMoabGetVertexCoordinates()
702 PetscErrorCode DMMoabGetVertexConnectivity(DM dm, moab::EntityHandle vhandle, PetscInt *nconn, moab… in DMMoabGetVertexConnectivity() argument
723 if (nconn) *nconn = connect.size(); in DMMoabGetVertexConnectivity()
742 PetscErrorCode DMMoabRestoreVertexConnectivity(DM dm, moab::EntityHandle ehandle, PetscInt *nconn, … in DMMoabRestoreVertexConnectivity() argument
749 if (nconn) *nconn = 0; in DMMoabRestoreVertexConnectivity()
770 PetscErrorCode DMMoabGetElementConnectivity(DM dm, moab::EntityHandle ehandle, PetscInt *nconn, con… in DMMoabGetElementConnectivity() argument
787 if (nconn) *nconn = nnodes; in DMMoabGetElementConnectivity()
853 PetscErrorCode DMMoabCheckBoundaryVertices(DM dm, PetscInt nconn, const moab::EntityHandle *cnt, Pe… in DMMoabCheckBoundaryVertices() argument
[all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexvtu.c | 7 PetscInt nconn; /* number of entries in cell->vertex connectivity array */ member
59 PetscCall(PetscMalloc3(piece->nconn, &conn, piece->ncells, &offsets, piece->ncells, &types)); in DMPlexGetVTKConnectivity()
115 …PetscCheck(countconn == piece->nconn, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Inconsistent connectivity … in DMPlexGetVTKConnectivity()
206 piece.nconn = 0; in DMPlexVTKWriteAll_VTU()
224 piece.nconn++; in DMPlexVTKWriteAll_VTU()
231 piece.nconn += dof / dimEmbed; in DMPlexVTKWriteAll_VTU()
254 boffset += gpiece[r].nconn * sizeof(PetscVTKInt) + (gpiece[r].nconn ? sizeof(PetscInt64) : 0); in DMPlexVTKWriteAll_VTU()
439 maxsize = PetscMax(maxsize, (PetscInt)(gpiece[r].nconn * sizeof(PetscVTKInt))); in DMPlexVTKWriteAll_VTU()
519 … PetscCall(TransferWrite(comm, viewer, fp, r, 0, connectivity, buffer, piece.nconn, MPI_INT, tag)); in DMPlexVTKWriteAll_VTU()
759 …PetscCall(TransferWrite(comm, viewer, fp, r, 0, NULL, buffer, gpiece[r].nconn, MPI_INT, tag)); … in DMPlexVTKWriteAll_VTU()
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 545 PetscInt nconn = iit[v + 1] - iit[v]; in PCBDDCNedelecSupport() local
547 if (!PetscBTLookup(btee, jjt[k])) nconn--; in PCBDDCNedelecSupport()
548 if (nconn > order) ends++; in PCBDDCNedelecSupport()
549 else if (nconn == order) ints++; in PCBDDCNedelecSupport()
|