| /petsc/src/dm/impls/swarm/ |
| H A D | swarmpic.c | 1168 PetscInt *pidx, Npc, q; in DMSwarmInitializeCoordinates() local
1170 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in DMSwarmInitializeCoordinates()
1172 for (q = 0; q < Npc; ++q) { in DMSwarmInitializeCoordinates()
1182 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx)); in DMSwarmInitializeCoordinates()
1355 PetscInt *pindices, Npc; in DMProjectFieldLocal_Swarm() local
1357 PetscCall(DMSwarmSortGetPointsPerCell(dm, cell, &Npc, &pindices)); in DMProjectFieldLocal_Swarm()
1358 maxC = PetscMax(maxC, Npc); in DMProjectFieldLocal_Swarm()
1359 PetscCall(DMSwarmSortRestorePointsPerCell(dm, cell, &Npc, &pindices)); in DMProjectFieldLocal_Swarm()
1365 PetscInt *pindices, Npc; in DMProjectFieldLocal_Swarm() local
1369 PetscCall(DMSwarmSortGetPointsPerCell(dm, cell, &Npc, &pindices)); in DMProjectFieldLocal_Swarm()
[all …]
|
| H A D | swarm.c | 1002 PetscInt Npc; in DMSwarmComputeGradientMatrix_Private() local
1004 PetscCall(DMSwarmSortGetPointsPerCell(sw, cell, &Npc, &pind)); in DMSwarmComputeGradientMatrix_Private()
1005 maxNpc = PetscMax(maxNpc, Npc); in DMSwarmComputeGradientMatrix_Private()
1006 PetscCall(DMSwarmSortRestorePointsPerCell(sw, cell, &Npc, &pind)); in DMSwarmComputeGradientMatrix_Private()
1018 PetscInt numFIndices, Npc; in DMSwarmComputeGradientMatrix_Private() local
1023 PetscCall(DMSwarmSortGetPointsPerCell(sw, cell, &Npc, &pind)); in DMSwarmComputeGradientMatrix_Private()
1024 for (PetscInt j = 0; j < Npc; ++j) { in DMSwarmComputeGradientMatrix_Private()
1034 PetscCall(PetscFECreateTabulation(fe, 1, Npc, xi, 1, &Tcoarse)); in DMSwarmComputeGradientMatrix_Private()
1036 PetscCall(PetscArrayzero(elemMat, Npc * cdim * totDim)); in DMSwarmComputeGradientMatrix_Private()
1038 for (PetscInt j = 0; j < Npc; ++j) { in DMSwarmComputeGradientMatrix_Private()
[all …]
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex11.c | 12 PetscInt Npc; // The number of partices per cell member
65 options->Npc = 1; in ProcessOptions()
74 …oundedInt("-npc", "Number of particles per cell", __FILE__, options->Npc, &options->Npc, NULL, 0)); in ProcessOptions()
186 PetscInt Npc = user->Npc, dim, Np; in CreateSwarm() local
200 PetscCall(DMSwarmInsertPointsUsingCellDM(*sw, DMSWARMPIC_LAYOUT_GAUSS, Npc)); in CreateSwarm()
|
| /petsc/src/ts/tutorials/hamiltonian/ |
| H A D | ex2.c | 646 PetscInt dim, Npc, cStart, cEnd, cell = ctx->velocity_monitor; in MonitorVelocity() local
687 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in MonitorVelocity()
690 PetscCall(VecSetSizes(cellv, Npc * dim, Npc)); in MonitorVelocity()
693 PetscCall(VecSetSizes(cellw, Npc, Npc)); in MonitorVelocity()
697 for (PetscInt q = 0; q < Npc; ++q) { in MonitorVelocity()
710 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx)); in MonitorVelocity()
762 PetscInt *pidx, Npc, q; in MonitorPositions_2D() local
763 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in MonitorPositions_2D()
764 for (q = 0; q < Npc; ++q) { in MonitorPositions_2D()
779 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx)); in MonitorPositions_2D()
[all …]
|
| H A D | ex4.c | 1162 PetscInt dim, Npc, cStart, cEnd, cell = ctx->velocity_monitor; in MonitorVelocity() local
1203 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in MonitorVelocity()
1206 PetscCall(VecSetSizes(cellv, Npc * dim, Npc)); in MonitorVelocity()
1209 PetscCall(VecSetSizes(cellw, Npc, Npc)); in MonitorVelocity()
1213 for (PetscInt q = 0; q < Npc; ++q) { in MonitorVelocity()
1226 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx)); in MonitorVelocity()
1278 PetscInt *pidx, Npc, q; in MonitorPositions_2D() local
1279 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in MonitorPositions_2D()
1280 for (q = 0; q < Npc; ++q) { in MonitorPositions_2D()
1295 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx)); in MonitorPositions_2D()
[all …]
|
| H A D | ex3.c | 459 PetscInt dim, xcStart, xcEnd, vcStart, vcEnd, Ns, Np, Npc, debug; in InitializeParticles_Centroid() local
485 Npc = (vcEnd - vcStart) * Ns; in InitializeParticles_Centroid()
486 Np = (xcEnd - xcStart) * Npc; in InitializeParticles_Centroid()
504 for (PetscInt q = 0; q < Npc / Ns; ++q, ++p) { in InitializeParticles_Centroid()
522 PetscInt *pidx, Npc; in InitializeParticles_Centroid() local
525 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in InitializeParticles_Centroid()
528 for (PetscInt q = 0; q < Npc / Ns; ++q) { in InitializeParticles_Centroid()
555 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx)); in InitializeParticles_Centroid()
614 PetscInt *pidx, Npc; in InitializeWeights() local
621 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in InitializeWeights()
[all …]
|
| /petsc/src/ksp/ksp/utils/dm/ |
| H A D | dmproject.c | 855 PetscInt dim, Ns, Npc, Np, cStart, cEnd, debug; in InitializeParticles_Regular() local
876 Npc = Ns * PetscPowInt(n, dim); in InitializeParticles_Regular()
877 Np = (cEnd - cStart) * Npc; in InitializeParticles_Regular()
884 PetscCall(PetscPrintf(comm, "Regular layout using %" PetscInt_FMT " particles per cell\n", Npc)); in InitializeParticles_Regular()
897 for (PetscInt q = 0; q < Npc / Ns; ++q, ++p) { in InitializeParticles_Regular()
929 …PetscCheck(Npc == npc, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Invalid number of points per cell %"… in InitializeParticles_Regular()
942 for (PetscInt q = 0; q < Npc / Ns; ++q) { in InitializeParticles_Regular()
971 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx)); in InitializeParticles_Regular()
1075 PetscInt *pidx, Npc; in DMSwarmRemap_Colella_Internal() local
1078 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in DMSwarmRemap_Colella_Internal()
[all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex77.c | 59 PetscInt Npc; /* The initial number of particles per cell */ member
297 options->Npc = 1; in ProcessOptions()
308 …("-Npc", "The initial number of particles per cell", "ex77.c", options->Npc, &options->Npc, NULL)); in ProcessOptions()
498 PetscCall(DMSwarmSetPointCoordinatesRandom(dm, user->Npc)); in SetInitialParticleConditions()
607 PetscCall(DMSwarmSetLocalSizes(sdm, (cEnd - cStart) * user->Npc, 0)); in SetupDiscretization()
611 for (p = 0; p < user->Npc; ++p) { in SetupDiscretization()
612 const PetscInt n = c * user->Npc + p; in SetupDiscretization()
618 PetscCall(DMSwarmSetPointCoordinatesRandom(sdm, user->Npc)); in SetupDiscretization()
|