| /petsc/src/ts/tutorials/hamiltonian/output/ |
| H A D | ex2_two_stream_rt.out | 4 c:0 [x_a,x_b] = 0.000000000000000,0.349066666666667 -> cell weight = 0.352539638258674 5 c:1 [x_a,x_b] = 0.349066666666667,0.698133333333333 -> cell weight = 0.352434113281987 6 c:2 [x_a,x_b] = 0.698133333333333,1.047200000000000 -> cell weight = 0.352226269666591 7 c:3 [x_a,x_b] = 1.047200000000000,1.396266666666667 -> cell weight = 0.351922422665025 8 c:4 [x_a,x_b] = 1.396266666666667,1.745333333333333 -> cell weight = 0.351531804557746 9 c:5 [x_a,x_b] = 1.745333333333333,2.094400000000000 -> cell weight = 0.351066284133648 10 c:6 [x_a,x_b] = 2.094400000000000,2.443466666666667 -> cell weight = 0.350540006061240 11 c:7 [x_a,x_b] = 2.443466666666667,2.792533333333333 -> cell weight = 0.349968961108035 12 c:8 [x_a,x_b] = 2.792533333333333,3.141600000000000 -> cell weight = 0.349370500266919 13 c:9 [x_a,x_b] = 3.141600000000000,3.490666666666667 -> cell weight = 0.348762807552540 [all …]
|
| H A D | ex2_two_stream_c0.out | 4 c:0 [x_a,x_b] = 0.000000000000000,0.349066666666667 -> cell weight = 0.352539638258674 5 c:1 [x_a,x_b] = 0.349066666666667,0.698133333333333 -> cell weight = 0.352434113281987 6 c:2 [x_a,x_b] = 0.698133333333333,1.047200000000000 -> cell weight = 0.352226269666591 7 c:3 [x_a,x_b] = 1.047200000000000,1.396266666666667 -> cell weight = 0.351922422665025 8 c:4 [x_a,x_b] = 1.396266666666667,1.745333333333333 -> cell weight = 0.351531804557746 9 c:5 [x_a,x_b] = 1.745333333333333,2.094400000000000 -> cell weight = 0.351066284133648 10 c:6 [x_a,x_b] = 2.094400000000000,2.443466666666667 -> cell weight = 0.350540006061240 11 c:7 [x_a,x_b] = 2.443466666666667,2.792533333333333 -> cell weight = 0.349968961108035 12 c:8 [x_a,x_b] = 2.792533333333333,3.141600000000000 -> cell weight = 0.349370500266919 13 c:9 [x_a,x_b] = 3.141600000000000,3.490666666666667 -> cell weight = 0.348762807552540 [all …]
|
| H A D | ex2_uniform_equilibrium_1d.out | 4 c:0 [x_a,x_b] = 0.000000000000000,1.256640000000000 -> cell weight = 1.256640000000000 5 c:1 [x_a,x_b] = 1.256640000000000,2.513280000000000 -> cell weight = 1.256640000000000 6 c:2 [x_a,x_b] = 2.513280000000000,3.769920000000000 -> cell weight = 1.256640000000000 7 c:3 [x_a,x_b] = 3.769920000000000,5.026560000000000 -> cell weight = 1.256640000000000 8 c:4 [x_a,x_b] = 5.026560000000000,6.283200000000000 -> cell weight = 1.256640000000000 9 c:5 [x_a,x_b] = 6.283200000000000,7.539840000000000 -> cell weight = 1.256640000000000 10 c:6 [x_a,x_b] = 7.539840000000000,8.796479999999999 -> cell weight = 1.256639999999999 11 c:7 [x_a,x_b] = 8.796479999999999,10.053120000000000 -> cell weight = 1.256640000000001 12 c:8 [x_a,x_b] = 10.053120000000000,11.309760000000001 -> cell weight = 1.256640000000001 13 c:9 [x_a,x_b] = 11.309760000000001,12.566400000000000 -> cell weight = 1.256639999999999 [all …]
|
| H A D | ex2_landau_damping_1d_dg.out | 4 c:0 [x_a,x_b] = 0.000000000000000,1.256640000000000 -> cell weight = 1.268395728819320 5 c:1 [x_a,x_b] = 1.256640000000000,2.513280000000000 -> cell weight = 1.263905419667826 6 c:2 [x_a,x_b] = 2.513280000000000,3.769920000000000 -> cell weight = 1.256639954596585 7 c:3 [x_a,x_b] = 3.769920000000000,5.026560000000000 -> cell weight = 1.249374506867984 8 c:4 [x_a,x_b] = 5.026560000000000,6.283200000000000 -> cell weight = 1.244884243120080 9 c:5 [x_a,x_b] = 6.283200000000000,7.539840000000000 -> cell weight = 1.244884299241914 10 c:6 [x_a,x_b] = 7.539840000000000,8.796479999999999 -> cell weight = 1.249374653796754 11 c:7 [x_a,x_b] = 8.796479999999999,10.053120000000000 -> cell weight = 1.256640136210245 12 c:8 [x_a,x_b] = 10.053120000000000,11.309760000000001 -> cell weight = 1.263905566595814 13 c:9 [x_a,x_b] = 11.309760000000001,12.566400000000000 -> cell weight = 1.268395784939883 [all …]
|
| H A D | ex2_uniform_primal_1d_real.out | 4 c:0 [x_a,x_b] = 0.000000000000000,1.256640000000000 -> cell weight = 1.268395728819320 5 c:1 [x_a,x_b] = 1.256640000000000,2.513280000000000 -> cell weight = 1.263905419667826 6 c:2 [x_a,x_b] = 2.513280000000000,3.769920000000000 -> cell weight = 1.256639954596585 7 c:3 [x_a,x_b] = 3.769920000000000,5.026560000000000 -> cell weight = 1.249374506867984 8 c:4 [x_a,x_b] = 5.026560000000000,6.283200000000000 -> cell weight = 1.244884243120080 9 c:5 [x_a,x_b] = 6.283200000000000,7.539840000000000 -> cell weight = 1.244884299241914 10 c:6 [x_a,x_b] = 7.539840000000000,8.796479999999999 -> cell weight = 1.249374653796754 11 c:7 [x_a,x_b] = 8.796479999999999,10.053120000000000 -> cell weight = 1.256640136210245 12 c:8 [x_a,x_b] = 10.053120000000000,11.309760000000001 -> cell weight = 1.263905566595814 13 c:9 [x_a,x_b] = 11.309760000000001,12.566400000000000 -> cell weight = 1.268395784939883 [all …]
|
| H A D | ex2_landau_damping_1d_mixed.out | 4 c:0 [x_a,x_b] = 0.000000000000000,1.256640000000000 -> cell weight = 1.268395728819320 5 c:1 [x_a,x_b] = 1.256640000000000,2.513280000000000 -> cell weight = 1.263905419667826 6 c:2 [x_a,x_b] = 2.513280000000000,3.769920000000000 -> cell weight = 1.256639954596585 7 c:3 [x_a,x_b] = 3.769920000000000,5.026560000000000 -> cell weight = 1.249374506867984 8 c:4 [x_a,x_b] = 5.026560000000000,6.283200000000000 -> cell weight = 1.244884243120080 9 c:5 [x_a,x_b] = 6.283200000000000,7.539840000000000 -> cell weight = 1.244884299241914 10 c:6 [x_a,x_b] = 7.539840000000000,8.796479999999999 -> cell weight = 1.249374653796754 11 c:7 [x_a,x_b] = 8.796479999999999,10.053120000000000 -> cell weight = 1.256640136210245 12 c:8 [x_a,x_b] = 10.053120000000000,11.309760000000001 -> cell weight = 1.263905566595814 13 c:9 [x_a,x_b] = 11.309760000000001,12.566400000000000 -> cell weight = 1.268395784939883 [all …]
|
| H A D | ex2_uniform_equilibrium_1d_real.out | 4 c:0 [x_a,x_b] = 0.000000000000000,1.256640000000000 -> cell weight = 1.256640000000000 5 c:1 [x_a,x_b] = 1.256640000000000,2.513280000000000 -> cell weight = 1.256640000000000 6 c:2 [x_a,x_b] = 2.513280000000000,3.769920000000000 -> cell weight = 1.256640000000000 7 c:3 [x_a,x_b] = 3.769920000000000,5.026560000000000 -> cell weight = 1.256640000000000 8 c:4 [x_a,x_b] = 5.026560000000000,6.283200000000000 -> cell weight = 1.256640000000000 9 c:5 [x_a,x_b] = 6.283200000000000,7.539840000000000 -> cell weight = 1.256640000000000 10 c:6 [x_a,x_b] = 7.539840000000000,8.796479999999999 -> cell weight = 1.256639999999999 11 c:7 [x_a,x_b] = 8.796479999999999,10.053120000000000 -> cell weight = 1.256640000000001 12 c:8 [x_a,x_b] = 10.053120000000000,11.309760000000001 -> cell weight = 1.256640000000001 13 c:9 [x_a,x_b] = 11.309760000000001,12.566400000000000 -> cell weight = 1.256639999999999 [all …]
|
| H A D | ex2_landau_damping_1d_c0.out | 4 c:0 [x_a,x_b] = 0.000000000000000,1.256640000000000 -> cell weight = 1.268395728819320 5 c:1 [x_a,x_b] = 1.256640000000000,2.513280000000000 -> cell weight = 1.263905419667826 6 c:2 [x_a,x_b] = 2.513280000000000,3.769920000000000 -> cell weight = 1.256639954596585 7 c:3 [x_a,x_b] = 3.769920000000000,5.026560000000000 -> cell weight = 1.249374506867984 8 c:4 [x_a,x_b] = 5.026560000000000,6.283200000000000 -> cell weight = 1.244884243120080 9 c:5 [x_a,x_b] = 6.283200000000000,7.539840000000000 -> cell weight = 1.244884299241914 10 c:6 [x_a,x_b] = 7.539840000000000,8.796479999999999 -> cell weight = 1.249374653796754 11 c:7 [x_a,x_b] = 8.796479999999999,10.053120000000000 -> cell weight = 1.256640136210245 12 c:8 [x_a,x_b] = 10.053120000000000,11.309760000000001 -> cell weight = 1.263905566595814 13 c:9 [x_a,x_b] = 11.309760000000001,12.566400000000000 -> cell weight = 1.268395784939883 [all …]
|
| H A D | ex2_uniform_primal_1d_real_pfak.out | 4 c:0 [x_a,x_b] = 0.000000000000000,1.256640000000000 -> cell weight = 1.268395728819320 5 c:1 [x_a,x_b] = 1.256640000000000,2.513280000000000 -> cell weight = 1.263905419667826 6 c:2 [x_a,x_b] = 2.513280000000000,3.769920000000000 -> cell weight = 1.256639954596585 7 c:3 [x_a,x_b] = 3.769920000000000,5.026560000000000 -> cell weight = 1.249374506867984 8 c:4 [x_a,x_b] = 5.026560000000000,6.283200000000000 -> cell weight = 1.244884243120080 9 c:5 [x_a,x_b] = 6.283200000000000,7.539840000000000 -> cell weight = 1.244884299241914 10 c:6 [x_a,x_b] = 7.539840000000000,8.796479999999999 -> cell weight = 1.249374653796754 11 c:7 [x_a,x_b] = 8.796479999999999,10.053120000000000 -> cell weight = 1.256640136210245 12 c:8 [x_a,x_b] = 10.053120000000000,11.309760000000001 -> cell weight = 1.263905566595814 13 c:9 [x_a,x_b] = 11.309760000000001,12.566400000000000 -> cell weight = 1.268395784939883 [all …]
|
| /petsc/src/dm/impls/plex/tests/output/ |
| H A D | ex1f90_0.out | 1 cell EC 0 0 2 3 4 5 2 cell pEC 0 0 2 3 4 5 3 cell EC 1 1 4 5 6 7 4 cell pEC 1 1 4 5 6 7 5 cell pEC 3 1 4 5 6 7
|
| H A D | ex97f90_0.out | 35 cell: 0 type: 6 36 cell: 1 type: 6
|
| H A D | ex97_0.out | 35 cell: 0 type: 6 36 cell: 1 type: 6
|
| /petsc/src/ts/utils/dmplexlandau/tutorials/output/ |
| H A D | ex1_re.out | 49 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=5.508970e-02, cell 0, qp 1/9 50 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=4.379594e-02, cell 0, qp 2/9 51 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=3.250658e-02, cell 0, qp 3/9 52 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=2.594060e-02, cell 4, qp 1/9 53 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=1.468071e-02, cell 4, qp 2/9 54 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=3.676391e-03, cell 4, qp 3/9 55 [0] <vec:seq> adaptToleranceFEM(): 1) Found another inner r=3.676391e-03, cell 8, qp 1/9, d=0.00… 59 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=5.508970e-02, cell 0, qp 1/9 60 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=4.379594e-02, cell 0, qp 2/9 61 [0] <vec:seq> adaptToleranceFEM(): 1) Found first inner r=3.250658e-02, cell 0, qp 3/9 [all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfvm.c | 7 …(DM dm, DM dmCell, PetscLimiter lim, PetscInt dim, PetscInt dof, PetscInt cell, PetscInt field, Pe… in DMPlexApplyLimiter_Internal() argument 20 … < fEnd) PetscCall(DMPlexApplyLimiter_Internal(dm, dmCell, lim, dim, dof, cell, field, childFace, … in DMPlexApplyLimiter_Internal() 35 ncell = cell == fcells[0] ? fcells[1] : fcells[0]; in DMPlexApplyLimiter_Internal() 66 PetscInt dim, face, cell, field, dof, cStart, cEnd, nFields; in DMPlexReconstructGradients_Internal() local 122 for (cell = (dmGrad && lim) ? cStart : cEnd; cell < cEnd; ++cell) { in DMPlexReconstructGradients_Internal() 129 PetscCall(DMPlexGetConeSize(dm, cell, &coneSize)); in DMPlexReconstructGradients_Internal() 130 PetscCall(DMPlexGetCone(dm, cell, &faces)); in DMPlexReconstructGradients_Internal() 132 PetscCall(DMPlexPointLocalFieldRead(dm, cell, field, x, &cx)); in DMPlexReconstructGradients_Internal() 134 PetscCall(DMPlexPointLocalRead(dm, cell, x, &cx)); in DMPlexReconstructGradients_Internal() 136 PetscCall(DMPlexPointLocalRead(dmCell, cell, cellgeom, &cg)); in DMPlexReconstructGradients_Internal() [all …]
|
| H A D | plexfluent.c | 285 static PetscErrorCode InsertFace(DM dm, PetscInt cell, PetscInt face, PetscInt ornt) in InsertFace() argument 291 PetscCall(DMPlexGetCone(dm, cell, &cone)); in InsertFace() 292 PetscCall(DMPlexGetConeSize(dm, cell, &coneSize)); in InsertFace() 295 …not be inserted in cone of cell %" PetscInt_FMT " with size %" PetscInt_FMT, face, cell, coneSize); in InsertFace() 296 PetscCall(DMPlexInsertCone(dm, cell, c, face)); in InsertFace() 297 PetscCall(DMPlexInsertConeOrientation(dm, cell, c, ornt)); in InsertFace() 301 static PetscErrorCode ReorderPolygon(DM dm, PetscInt cell) in ReorderPolygon() argument 307 PetscCall(DMPlexGetOrientedCone(dm, cell, &cone, &ornt)); in ReorderPolygon() 308 PetscCall(DMPlexGetConeSize(dm, cell, &coneSize)); in ReorderPolygon() 327 …G_WRONG, "Cell %" PetscInt_FMT " could not find a face match as position %" PetscInt_FMT, cell, c); in ReorderPolygon() [all …]
|
| H A D | plexgeometry.c | 421 …DMPlexLocatePoint_Simplex_1D_Internal(DM dm, const PetscScalar point[], PetscInt c, PetscInt *cell) in DMPlexLocatePoint_Simplex_1D_Internal() argument 432 if ((xi >= -eps) && (xi <= 2. + eps)) *cell = c; in DMPlexLocatePoint_Simplex_1D_Internal() 433 else *cell = DMLOCATEPOINT_POINT_NOT_FOUND; in DMPlexLocatePoint_Simplex_1D_Internal() 437 …DMPlexLocatePoint_Simplex_2D_Internal(DM dm, const PetscScalar point[], PetscInt c, PetscInt *cell) in DMPlexLocatePoint_Simplex_2D_Internal() argument 451 if ((xi[0] >= -eps) && (xi[1] >= -eps) && (xi[0] + xi[1] <= 2.0 + eps)) *cell = c; in DMPlexLocatePoint_Simplex_2D_Internal() 452 else *cell = DMLOCATEPOINT_POINT_NOT_FOUND; in DMPlexLocatePoint_Simplex_2D_Internal() 482 …exLocatePoint_Quad_2D_Linear_Internal(DM dm, const PetscScalar point[], PetscInt c, PetscInt *cell) in DMPlexLocatePoint_Quad_2D_Linear_Internal() argument 518 if (crossings % 2) *cell = c; in DMPlexLocatePoint_Quad_2D_Linear_Internal() 519 else *cell = DMLOCATEPOINT_POINT_NOT_FOUND; in DMPlexLocatePoint_Quad_2D_Linear_Internal() 524 …de DMPlexLocatePoint_Quad_2D_Internal(DM dm, const PetscScalar point[], PetscInt c, PetscInt *cell) in DMPlexLocatePoint_Quad_2D_Internal() argument [all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex52.c | 82 PetscInt cell, cStart, cEnd; in FormInitialSolution() local 92 for (cell = cStart; cell < cEnd; cell++) { in FormInitialSolution() 93 PetscCall(DMPlexComputeCellGeometryFVM(da, cell, &cellvol, centroid, normal)); in FormInitialSolution() 95 if (centroid[1] > 0.9 && centroid[1] < 0.95) u[cell] = 2.0; in FormInitialSolution() 96 } else u[cell] = 0; in FormInitialSolution() 149 PetscInt cell, cStart, cEnd, nC; in FormFunction() local 197 for (cell = cStart; cell < cEnd; cell++) { in FormFunction() 199 PetscCall(DMPlexGetConeSize(da, cell, &nF)); in FormFunction() 200 PetscCall(DMPlexGetCone(da, cell, &cellcone)); in FormFunction() 247 …flux_centre = x[cell] * ((u_plus - u_minus + 2 * user->diffusion * delta_y) / delta_x + (v_plus - … in FormFunction() [all …]
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex100.c | 79 for (PetscInt cell = cStart; cell < cEnd; cell++) { in main() local 86 … PetscCall(DMPlexVecGetClosure(dm_read, coord_section, coords, cell, &num_closure, &cell_coords)); in main() 95 …if (cell == cStart) PetscCall(PetscArraycpy(ref_cell_bounding_box_size, cell_bounding_box_size, cd… in main() 100 …tscInt_FMT " in dimension %" PetscInt_FMT ". Relative difference = %g\n", cell, cStart, d, (double… in main() 103 …PetscCall(DMPlexVecRestoreClosure(dm_read, coord_section, coords, cell, &num_closure, &cell_coords… in main()
|
| H A D | ex47.c | 420 PetscInt cell = (InitPartForRank[rank])[c], cor; in main() local 422 for (cor = 0; cor < Ncor; ++cor) cells[c * Ncor + cor] = Conn[cell][cor]; in main() 447 PetscInt cell = (InitPartForRank[rank])[c], cor; in main() local 448 PetscInt nbElemVertex = ((-1 == Conn[cell][NcorMax - 1]) ? 3 : 4); in main() 450 cells[count] = Conn[cell][cor]; in main() 473 PetscInt cell = (InitPartForRank[rank])[c], cor; in main() local 475 for (cor = 0; cor < Ncor; ++cor) cells[c * Ncor + cor] = Conn[cell][cor]; in main() 500 PetscInt cell = (InitPartForRank[rank])[c], cor; in main() local 501 PetscInt nbElemVertex = ((-1 == Conn[cell][NcorMax - 1]) ? 6 : 8); in main() 503 cells[count] = Conn[cell][cor]; in main()
|
| H A D | ex8.c | 146 static PetscErrorCode CheckFEMGeometry(DM dm, PetscInt cell, PetscInt spaceDim, PetscReal v0Ex[], P… in CheckFEMGeometry() argument 152 PetscCall(DMPlexComputeCellGeometryFEM(dm, cell, NULL, v0, J, invJ, &detJ)); in CheckFEMGeometry() 175 static PetscErrorCode CheckFVMGeometry(DM dm, PetscInt cell, PetscInt spaceDim, PetscReal centroidE… in CheckFVMGeometry() argument 182 …PetscCall(DMPlexComputeCellGeometryFVM(dm, cell, volEx ? &vol : NULL, centroidEx ? centroid : NULL… in CheckFVMGeometry() 185 …" PetscInt_FMT ", Invalid centroid[%" PetscInt_FMT "]: %g != %g diff %g", cell, d, (double)centroi… in CheckFVMGeometry() 186 …B, "Cell %" PetscInt_FMT ", Invalid normal[%" PetscInt_FMT "]: %g != %g", cell, d, (double)normal[… in CheckFVMGeometry() 188 …SC_ERR_PLIB, "Cell %" PetscInt_FMT ", Invalid volume = %g != %g diff %g", cell, (double)vol, (doub… in CheckFVMGeometry() 192 static PetscErrorCode CheckGaussLaw(DM dm, PetscInt cell) in CheckGaussLaw() argument 204 PetscCall(DMPlexGetCellType(dm, cell, &ct)); in CheckGaussLaw() 206 PetscCall(DMPlexGetConeSize(dm, cell, &coneSize)); in CheckGaussLaw() [all …]
|
| H A D | ex26.c | 152 PetscInt numCells, cell, closureSize, *closureA = NULL; in main() local 235 for (cell = 0; cell < numCells; cell++) { in main() 238 PetscCall(DMPlexGetTransitiveClosure(dm, cellID[cell], PETSC_TRUE, &closureSize, &closure)); in main() 250 … PetscCall(DMPlexRestoreTransitiveClosure(dm, cellID[cell], PETSC_TRUE, &closureSize, &closure)); in main() 420 PetscInt numCells, cell; in main() local 427 for (cell = 0; cell < numCells; cell++) { in main() 428 PetscCall(DMPlexVecGetClosure(dmS, NULL, S, cellID[cell], &clSize, &cval)); in main() 429 PetscCall(DMPlexVecGetClosure(dmS, coordSection, coord, cellID[cell], &cdimCoord, &xyz)); in main() 435 PetscCall(DMPlexVecSetClosure(dmS, NULL, S, cellID[cell], cval, INSERT_ALL_VALUES)); in main()
|
| /petsc/src/dm/impls/da/ |
| H A D | dagetelem.c | 42 PetscInt cnt = 0, cell[4], ns = 2; in DMDAGetElements_2D() local 77 cell[0] = (i - Xs) + (j - Ys) * (Xe - Xs); in DMDAGetElements_2D() 78 cell[1] = (i - Xs + 1) + (j - Ys) * (Xe - Xs); in DMDAGetElements_2D() 79 cell[2] = (i - Xs + 1) + (j - Ys + 1) * (Xe - Xs); in DMDAGetElements_2D() 80 cell[3] = (i - Xs) + (j - Ys + 1) * (Xe - Xs); in DMDAGetElements_2D() 82 for (c = 0; c < ns * nn; c++) da->e[cnt++] = cell[split[c]]; in DMDAGetElements_2D() 85 for (c = 0; c < ns * nn; c++) da->e[cnt++] = cell[c]; in DMDAGetElements_2D() 108 PetscInt cnt = 0, cell[8], ns = 6; in DMDAGetElements_3D() local 147 cell[0] = (i - Xs) + (j - Ys) * (Xe - Xs) + (k - Zs) * (Xe - Xs) * (Ye - Ys); in DMDAGetElements_3D() 148 cell[1] = (i + 1 - Xs) + (j - Ys) * (Xe - Xs) + (k - Zs) * (Xe - Xs) * (Ye - Ys); in DMDAGetElements_3D() [all …]
|
| /petsc/src/dm/impls/swarm/ |
| H A D | swarmpic_sort.c | 130 PetscErrorCode DMSwarmSortGetNumberOfPointsPerCell(DM sw, PetscInt cell, PetscInt *npoints) in DMSwarmSortGetNumberOfPointsPerCell() argument 140 …cell < ctx->ncells, PETSC_COMM_SELF, PETSC_ERR_USER, "Cell index (%" PetscInt_FMT ") is greater th… in DMSwarmSortGetNumberOfPointsPerCell() 141 …PetscCheck(cell >= 0, PETSC_COMM_SELF, PETSC_ERR_USER, "Cell index (%" PetscInt_FMT ") cannot be n… in DMSwarmSortGetNumberOfPointsPerCell() 142 *npoints = ctx->pcell_offsets[cell + 1] - ctx->pcell_offsets[cell]; in DMSwarmSortGetNumberOfPointsPerCell() 164 PetscErrorCode DMSwarmSortGetPointsPerCell(DM sw, PetscInt cell, PetscInt *npoints, PetscInt **pidl… in DMSwarmSortGetPointsPerCell() argument 174 PetscCall(DMSwarmSortGetNumberOfPointsPerCell(sw, cell, npoints)); in DMSwarmSortGetPointsPerCell() 177 pid = ctx->pcell_offsets[cell] + p; in DMSwarmSortGetPointsPerCell()
|
| /petsc/src/dm/impls/plex/tutorials/ |
| H A D | ex6.c | 121 PetscInt cStart, cEnd, cell; in ReadData2D() local 125 for (cell = cStart; cell < cEnd; ++cell) { in ReadData2D() 129 PetscCall(DMPlexVecGetClosure(dm, NULL, u, cell, &closureSize, &closure)); in ReadData2D() 130 PetscCall(PetscPrintf(PETSC_COMM_SELF, "Cell %" PetscInt_FMT "\n", cell)); in ReadData2D() 146 PetscCall(DMPlexVecRestoreClosure(dm, NULL, u, cell, &closureSize, &closure)); in ReadData2D() 154 PetscInt cStart, cEnd, cell; in ReadData3D() local 158 for (cell = cStart; cell < cEnd; ++cell) { in ReadData3D() 162 PetscCall(DMPlexVecGetClosure(dm, NULL, u, cell, &closureSize, &closure)); in ReadData3D() 163 PetscCall(PetscPrintf(PETSC_COMM_SELF, "Cell %" PetscInt_FMT "\n", cell)); in ReadData3D() 182 PetscCall(DMPlexVecRestoreClosure(dm, NULL, u, cell, &closureSize, &closure)); in ReadData3D()
|
| /petsc/doc/overview/ |
| H A D | plex_transform_table.md | 38 - Replaces each simplex cell with :math:`2^d` box cells 42 - Replaces each box cell with simplex cells 51 * - Cohesive cell extrusion
|