Lines Matching refs:sw
148 static PetscErrorCode SetupContext(DM dm, DM sw, AppCtx *user) in SetupContext() argument
180 static PetscErrorCode CheckNonNegativeWeights(DM sw, AppCtx *user) in CheckNonNegativeWeights() argument
187 PetscCall(DMSwarmGetField(sw, "w_q", NULL, NULL, (void **)&w)); in CheckNonNegativeWeights()
188 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in CheckNonNegativeWeights()
190 PetscCall(DMSwarmRestoreField(sw, "w_q", NULL, NULL, (void **)&w)); in CheckNonNegativeWeights()
225 DM sw; in MonitorEField() local
236 PetscCall(TSGetDM(ts, &sw)); in MonitorEField()
237 PetscCall(DMGetDimension(sw, &dim)); in MonitorEField()
238 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in MonitorEField()
239 PetscCall(DMSwarmGetSize(sw, &gNp)); in MonitorEField()
240 PetscCall(DMSwarmSortGetAccess(sw)); in MonitorEField()
241 PetscCall(DMSwarmGetField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&x)); in MonitorEField()
242 PetscCall(DMSwarmGetField(sw, "E_field", NULL, NULL, (void **)&E)); in MonitorEField()
243 PetscCall(DMSwarmGetField(sw, "species", NULL, NULL, (void **)&species)); in MonitorEField()
244 PetscCall(DMSwarmGetField(sw, "w_q", NULL, NULL, (void **)&weight)); in MonitorEField()
274 PetscCall(DMSwarmRestoreField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&x)); in MonitorEField()
275 PetscCall(DMSwarmRestoreField(sw, "w_q", NULL, NULL, (void **)&weight)); in MonitorEField()
276 PetscCall(DMSwarmRestoreField(sw, "E_field", NULL, NULL, (void **)&E)); in MonitorEField()
277 PetscCall(DMSwarmRestoreField(sw, "species", NULL, NULL, (void **)&species)); in MonitorEField()
284 PetscCall(DMSwarmComputeMoments(sw, "velocity", "w_q", pmoments)); in MonitorEField()
286 PetscCall(DMViewFromOptions(sw, NULL, "-sw_efield_view")); in MonitorEField()
292 DM sw; in MonitorMoments() local
297 PetscCall(TSGetDM(ts, &sw)); in MonitorMoments()
299 PetscCall(DMSwarmComputeMoments(sw, "velocity", "w_q", pmoments)); in MonitorMoments()
308 DM sw; in MonitorParticles() local
318 PetscCall(TSGetDM(ts, &sw)); in MonitorParticles()
320 PetscCall(DMGetDimension(sw, &dim)); in MonitorParticles()
321 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in MonitorParticles()
323 PetscCall(DMSwarmGetField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&x)); in MonitorParticles()
324 PetscCall(DMSwarmGetField(sw, "velocity", NULL, NULL, (void **)&v)); in MonitorParticles()
325 PetscCall(DMSwarmGetField(sw, "w_q", NULL, NULL, (void **)&w)); in MonitorParticles()
326 PetscCall(DMSwarmGetField(sw, "E_field", NULL, NULL, (void **)&E)); in MonitorParticles()
354 PetscCall(DMSwarmRestoreField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&x)); in MonitorParticles()
355 PetscCall(DMSwarmRestoreField(sw, "velocity", NULL, NULL, (void **)&v)); in MonitorParticles()
356 PetscCall(DMSwarmRestoreField(sw, "w_q", NULL, NULL, (void **)&w)); in MonitorParticles()
357 PetscCall(DMSwarmRestoreField(sw, "E_field", NULL, NULL, (void **)&E)); in MonitorParticles()
396 static PetscErrorCode CreateVelocityDM(DM sw, DM *vdm) in CreateVelocityDM() argument
405 PetscCall(DMGetDimension(sw, &dim)); in CreateVelocityDM()
406 PetscCall(DMGetApplicationContext(sw, &user)); in CreateVelocityDM()
451 static PetscErrorCode InitializeParticles_Centroid(DM sw) in InitializeParticles_Centroid() argument
453 DM_Swarm *swarm = (DM_Swarm *)sw->data; in InitializeParticles_Centroid()
465 PetscCall(PetscObjectGetComm((PetscObject)sw, &comm)); in InitializeParticles_Centroid()
468 PetscCall(DMSwarmGetNumSpecies(sw, &Ns)); in InitializeParticles_Centroid()
470 if (flg) PetscCall(DMSwarmSetNumSpecies(sw, Ns)); in InitializeParticles_Centroid()
476 PetscCall(DMGetDimension(sw, &dim)); in InitializeParticles_Centroid()
477 PetscCall(DMSwarmGetCellDM(sw, &xdm)); in InitializeParticles_Centroid()
480 PetscCall(DMSwarmGetCellDMByName(sw, "velocity", &celldm)); in InitializeParticles_Centroid()
487 PetscCall(DMSwarmSetLocalSizes(sw, Np, 0)); in InitializeParticles_Centroid()
498 PetscCall(DMSwarmGetCellDMActive(sw, &celldm)); in InitializeParticles_Centroid()
500 PetscCall(DMSwarmGetField(sw, "species", NULL, NULL, (void **)&species)); in InitializeParticles_Centroid()
501 PetscCall(DMSwarmGetField(sw, cellidname, NULL, NULL, (void **)&cellid)); in InitializeParticles_Centroid()
510 PetscCall(DMSwarmRestoreField(sw, "species", NULL, NULL, (void **)&species)); in InitializeParticles_Centroid()
511 PetscCall(DMSwarmRestoreField(sw, cellidname, NULL, NULL, (void **)&cellid)); in InitializeParticles_Centroid()
514 PetscCall(DMSwarmGetField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&x)); in InitializeParticles_Centroid()
515 PetscCall(DMSwarmGetField(sw, "velocity", NULL, NULL, (void **)&v)); in InitializeParticles_Centroid()
516 PetscCall(DMSwarmSortGetAccess(sw)); in InitializeParticles_Centroid()
525 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in InitializeParticles_Centroid()
555 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx)); in InitializeParticles_Centroid()
557 PetscCall(DMSwarmSortRestoreAccess(sw)); in InitializeParticles_Centroid()
558 PetscCall(DMSwarmRestoreField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&x)); in InitializeParticles_Centroid()
559 PetscCall(DMSwarmRestoreField(sw, "velocity", NULL, NULL, (void **)&v)); in InitializeParticles_Centroid()
577 static PetscErrorCode InitializeWeights(DM sw, PetscReal totalWeight, PetscProbFn *func, const Pets… in InitializeWeights() argument
588 …PetscInt dim, Ns, xcStart, xcEnd, vcStart, vcEnd, debug = ((DM_Swarm *)sw->data)->printWei… in InitializeWeights()
593 PetscCall(PetscObjectGetComm((PetscObject)sw, &comm)); in InitializeWeights()
595 PetscCall(DMGetDimension(sw, &dim)); in InitializeWeights()
596 PetscCall(DMSwarmGetCellDM(sw, &xdm)); in InitializeWeights()
597 PetscCall(DMSwarmGetNumSpecies(sw, &Ns)); in InitializeWeights()
599 PetscCall(DMSwarmGetCellDMByName(sw, "velocity", &celldm)); in InitializeWeights()
610 PetscCall(DMSwarmSortGetAccess(sw)); in InitializeWeights()
611 PetscCall(DMSwarmGetField(sw, "w_q", NULL, NULL, (void **)&weight)); in InitializeWeights()
621 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx)); in InitializeWeights()
676 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx)); in InitializeWeights()
678 PetscCall(DMSwarmRestoreField(sw, "w_q", NULL, NULL, (void **)&weight)); in InitializeWeights()
679 PetscCall(DMSwarmSortRestoreAccess(sw)); in InitializeWeights()
692 static PetscErrorCode InitializeParticles_PerturbedWeights(DM sw, AppCtx *user) in InitializeParticles_PerturbedWeights() argument
698 PetscCall(DMGetDimension(sw, &dim)); in InitializeParticles_PerturbedWeights()
699 PetscCall(InitializeParticles_Centroid(sw)); in InitializeParticles_PerturbedWeights()
700 …PetscCall(InitializeWeights(sw, user->totalWeight, dim == 1 ? PetscPDFCosine1D : (dim == 2 ? Petsc… in InitializeParticles_PerturbedWeights()
704 static PetscErrorCode InitializeConstants(DM sw, AppCtx *user) in InitializeConstants() argument
712 PetscCall(DMSwarmGetCellDM(sw, &dm)); in InitializeConstants()
713 PetscCall(DMGetDimension(sw, &dim)); in InitializeConstants()
714 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in InitializeConstants()
716 PetscCall(DMSwarmGetField(sw, "w_q", NULL, NULL, (void **)&weight)); in InitializeConstants()
717 PetscCall(DMSwarmGetField(sw, "species", NULL, NULL, (void **)&species)); in InitializeConstants()
722 PetscCall(DMSwarmRestoreField(sw, "w_q", NULL, NULL, (void **)&weight)); in InitializeConstants()
723 PetscCall(DMSwarmRestoreField(sw, "species", NULL, NULL, (void **)&species)); in InitializeConstants()
921 static PetscErrorCode CreateSwarm(DM dm, AppCtx *user, DM *sw) in CreateSwarm() argument
928 PetscCall(DMCreate(PetscObjectComm((PetscObject)dm), sw)); in CreateSwarm()
929 PetscCall(DMSetType(*sw, DMSWARM)); in CreateSwarm()
930 PetscCall(DMSetDimension(*sw, dim)); in CreateSwarm()
931 PetscCall(DMSwarmSetType(*sw, DMSWARM_PIC)); in CreateSwarm()
932 PetscCall(DMSetApplicationContext(*sw, user)); in CreateSwarm()
934 PetscCall(DMSwarmRegisterPetscDatatypeField(*sw, "w_q", 1, PETSC_SCALAR)); in CreateSwarm()
935 PetscCall(DMSwarmRegisterPetscDatatypeField(*sw, "velocity", dim, PETSC_REAL)); in CreateSwarm()
936 PetscCall(DMSwarmRegisterPetscDatatypeField(*sw, "species", 1, PETSC_INT)); in CreateSwarm()
937 PetscCall(DMSwarmRegisterPetscDatatypeField(*sw, "E_field", dim, PETSC_REAL)); in CreateSwarm()
942 PetscCall(DMSwarmAddCellDM(*sw, celldm)); in CreateSwarm()
948 PetscCall(CreateVelocityDM(*sw, &vdm)); in CreateSwarm()
950 PetscCall(DMSwarmAddCellDM(*sw, celldm)); in CreateSwarm()
961 PetscCall(DMSwarmAddCellDM(*sw, celldm)); in CreateSwarm()
964 PetscCall(DMSetFromOptions(*sw)); in CreateSwarm()
965 PetscCall(DMSetUp(*sw)); in CreateSwarm()
967 PetscCall(DMSwarmSetCellDMActive(*sw, "space")); in CreateSwarm()
968 user->swarm = *sw; in CreateSwarm()
970 PetscCall(InitializeParticles_PerturbedWeights(*sw, user)); in CreateSwarm()
971 PetscCall(PetscObjectSetName((PetscObject)*sw, "Particles")); in CreateSwarm()
972 PetscCall(DMViewFromOptions(*sw, NULL, "-sw_view")); in CreateSwarm()
976 static PetscErrorCode ComputeFieldAtParticles_Primal(SNES snes, DM sw, Mat M_p, PetscReal E[]) in ComputeFieldAtParticles_Primal() argument
991 PetscCall(DMGetApplicationContext(sw, (void *)&user)); in ComputeFieldAtParticles_Primal()
992 PetscCall(PetscLogEventBegin(user->ESolveEvent, snes, sw, 0, 0)); in ComputeFieldAtParticles_Primal()
993 PetscCall(DMGetDimension(sw, &dim)); in ComputeFieldAtParticles_Primal()
994 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in ComputeFieldAtParticles_Primal()
1000 PetscCall(DMSwarmCreateGlobalVectorFromField(sw, "w_q", &f)); in ComputeFieldAtParticles_Primal()
1008 PetscCall(DMSwarmDestroyGlobalVectorFromField(sw, "w_q", &f)); in ComputeFieldAtParticles_Primal()
1054 PetscCall(PetscLogEventEnd(user->ESolveEvent, snes, sw, 0, 0)); in ComputeFieldAtParticles_Primal()
1058 PetscCall(DMSwarmSortGetAccess(sw)); in ComputeFieldAtParticles_Primal()
1060 PetscCall(DMSwarmGetField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&coords)); in ComputeFieldAtParticles_Primal()
1062 PetscCall(PetscLogEventBegin(user->ETabEvent, snes, sw, 0, 0)); in ComputeFieldAtParticles_Primal()
1071 PetscCall(DMSwarmSortGetNumberOfPointsPerCell(sw, c, &Ncp)); in ComputeFieldAtParticles_Primal()
1083 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Ncp, &points)); in ComputeFieldAtParticles_Primal()
1104 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Ncp, &points)); in ComputeFieldAtParticles_Primal()
1109 PetscCall(DMSwarmRestoreField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&coords)); in ComputeFieldAtParticles_Primal()
1110 PetscCall(DMSwarmSortRestoreAccess(sw)); in ComputeFieldAtParticles_Primal()
1113 PetscCall(PetscLogEventEnd(user->ETabEvent, snes, sw, 0, 0)); in ComputeFieldAtParticles_Primal()
1117 static PetscErrorCode ComputeFieldAtParticles(SNES snes, DM sw) in ComputeFieldAtParticles() argument
1126 PetscValidHeaderSpecific(sw, DM_CLASSID, 2); in ComputeFieldAtParticles()
1127 PetscCall(DMGetDimension(sw, &dim)); in ComputeFieldAtParticles()
1128 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in ComputeFieldAtParticles()
1129 PetscCall(DMGetApplicationContext(sw, &user)); in ComputeFieldAtParticles()
1131 PetscCall(DMSwarmSetCellDMActive(sw, "moments")); in ComputeFieldAtParticles()
1133 PetscCall(DMSwarmMigrate(sw, PETSC_FALSE)); in ComputeFieldAtParticles()
1134 PetscCall(DMCreateMassMatrix(sw, user->dmPot, &M_p)); in ComputeFieldAtParticles()
1135 PetscCall(DMSwarmSetCellDMActive(sw, "space")); in ComputeFieldAtParticles()
1137 PetscCall(DMSwarmGetField(sw, "E_field", NULL, NULL, (void **)&E)); in ComputeFieldAtParticles()
1141 PetscCall(ComputeFieldAtParticles_Primal(snes, sw, M_p, E)); in ComputeFieldAtParticles()
1142 PetscCall(DMSwarmRestoreField(sw, "E_field", NULL, NULL, (void **)&E)); in ComputeFieldAtParticles()
1149 DM sw; in RHSFunction() local
1157 PetscCall(TSGetDM(ts, &sw)); in RHSFunction()
1158 PetscCall(ComputeFieldAtParticles(snes, sw)); in RHSFunction()
1160 PetscCall(DMGetDimension(sw, &dim)); in RHSFunction()
1161 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in RHSFunction()
1162 PetscCall(DMSwarmGetField(sw, "E_field", NULL, NULL, (void **)&E)); in RHSFunction()
1172 PetscCall(DMSwarmRestoreField(sw, "E_field", NULL, NULL, (void **)&E)); in RHSFunction()
1186 DM sw; in RHSJacobian() local
1191 PetscCall(TSGetDM(ts, &sw)); in RHSJacobian()
1192 PetscCall(DMGetDimension(sw, &dim)); in RHSJacobian()
1193 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in RHSJacobian()
1195 PetscCall(DMSwarmGetField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&coords)); in RHSJacobian()
1196 PetscCall(DMSwarmGetField(sw, "velocity", NULL, NULL, (void **)&vel)); in RHSJacobian()
1207 PetscCall(DMSwarmRestoreField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&coords)); in RHSJacobian()
1208 PetscCall(DMSwarmRestoreField(sw, "velocity", NULL, NULL, (void **)&vel)); in RHSJacobian()
1217 DM sw; in RHSFunctionX() local
1224 PetscCall(TSGetDM(ts, &sw)); in RHSFunctionX()
1225 PetscCall(DMGetDimension(sw, &dim)); in RHSFunctionX()
1241 DM sw; in RHSFunctionV() local
1252 PetscCall(TSGetDM(ts, &sw)); in RHSFunctionV()
1253 PetscCall(ComputeFieldAtParticles(snes, sw)); in RHSFunctionV()
1255 PetscCall(DMGetDimension(sw, &dim)); in RHSFunctionV()
1256 PetscCall(DMSwarmGetField(sw, "E_field", NULL, NULL, (void **)&E)); in RHSFunctionV()
1285 PetscCall(DMSwarmRestoreField(sw, "E_field", NULL, NULL, (void **)&E)); in RHSFunctionV()
1294 DM sw; in RHSJacobianS() local
1298 PetscCall(TSGetDM(ts, &sw)); in RHSJacobianS()
1299 PetscCall(DMGetDimension(sw, &dim)); in RHSJacobianS()
1317 DM sw; in RHSObjectiveF() local
1324 PetscCall(TSGetDM(ts, &sw)); in RHSObjectiveF()
1325 PetscCall(DMGetDimension(sw, &dim)); in RHSObjectiveF()
1333 PetscCall(DMSwarmGetField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&coords)); in RHSObjectiveF()
1334 PetscCall(DMSwarmGetField(sw, "velocity", NULL, NULL, (void **)&vel)); in RHSObjectiveF()
1335 PetscCall(DMSwarmSortGetAccess(sw)); in RHSObjectiveF()
1343 PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Ncp, &points)); in RHSObjectiveF()
1350 PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Ncp, &points)); in RHSObjectiveF()
1353 PetscCall(DMSwarmSortRestoreAccess(sw)); in RHSObjectiveF()
1354 PetscCall(DMSwarmRestoreField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&coords)); in RHSObjectiveF()
1355 PetscCall(DMSwarmRestoreField(sw, "velocity", NULL, NULL, (void **)&vel)); in RHSObjectiveF()
1362 DM sw; in RHSFunctionG() local
1371 PetscCall(TSGetDM(ts, &sw)); in RHSFunctionG()
1372 PetscCall(DMGetDimension(sw, &dim)); in RHSFunctionG()
1373 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in RHSFunctionG()
1380 PetscCall(ComputeFieldAtParticles(snes, sw)); in RHSFunctionG()
1382 PetscCall(DMSwarmGetField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&coords)); in RHSFunctionG()
1383 PetscCall(DMSwarmGetField(sw, "velocity", NULL, NULL, (void **)&vel)); in RHSFunctionG()
1384 PetscCall(DMSwarmGetField(sw, "E_field", NULL, NULL, (void **)&E)); in RHSFunctionG()
1391 PetscCall(DMSwarmRestoreField(sw, "E_field", NULL, NULL, (void **)&E)); in RHSFunctionG()
1392 PetscCall(DMSwarmRestoreField(sw, DMSwarmPICField_coor, NULL, NULL, (void **)&coords)); in RHSFunctionG()
1393 PetscCall(DMSwarmRestoreField(sw, "velocity", NULL, NULL, (void **)&vel)); in RHSFunctionG()
1401 DM sw; in CreateSolution() local
1406 PetscCall(TSGetDM(ts, &sw)); in CreateSolution()
1407 PetscCall(DMGetDimension(sw, &dim)); in CreateSolution()
1408 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in CreateSolution()
1421 DM sw; in SetProblem() local
1424 PetscCall(TSGetDM(ts, &sw)); in SetProblem()
1425 PetscCall(DMGetApplicationContext(sw, (void **)&user)); in SetProblem()
1431 PetscCall(DMGetDimension(sw, &dim)); in SetProblem()
1432 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in SetProblem()
1450 PetscCall(DMGetDimension(sw, &dim)); in SetProblem()
1451 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in SetProblem()
1479 DM sw; in DMSwarmTSRedistribute() local
1485 PetscCall(TSGetDM(ts, &sw)); in DMSwarmTSRedistribute()
1493 PetscCall(TSSetDM(ts, sw)); in DMSwarmTSRedistribute()
1509 DM sw, cdm; in line() local
1514 sw = user->swarm; in line()
1515 PetscCall(DMSwarmGetCellDM(sw, &cdm)); in line()
1518 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in line()
1541 DM sw; in InitializeSolveAndSwarm() local
1548 PetscCall(TSGetDM(ts, &sw)); in InitializeSolveAndSwarm()
1549 PetscCall(DMGetApplicationContext(sw, &user)); in InitializeSolveAndSwarm()
1550 PetscCall(DMGetDimension(sw, &dim)); in InitializeSolveAndSwarm()
1552 PetscCall(InitializeParticles_PerturbedWeights(sw, user)); in InitializeSolveAndSwarm()
1553 PetscCall(DMSwarmMigrate(sw, PETSC_TRUE)); in InitializeSolveAndSwarm()
1556 PetscCall(DMSetUp(sw)); in InitializeSolveAndSwarm()
1560 PetscCall(DMSwarmCreateGlobalVectorFromField(sw, DMSwarmPICField_coor, &gc)); in InitializeSolveAndSwarm()
1561 PetscCall(DMSwarmCreateGlobalVectorFromField(sw, "velocity", &gv)); in InitializeSolveAndSwarm()
1564 PetscCall(DMSwarmDestroyGlobalVectorFromField(sw, DMSwarmPICField_coor, &gc)); in InitializeSolveAndSwarm()
1565 PetscCall(DMSwarmDestroyGlobalVectorFromField(sw, "velocity", &gv)); in InitializeSolveAndSwarm()
1579 DM sw, cdm; in MigrateParticles() local
1584 PetscCall(TSGetDM(ts, &sw)); in MigrateParticles()
1585 PetscCall(DMGetApplicationContext(sw, &ctx)); in MigrateParticles()
1586 PetscCall(DMViewFromOptions(sw, NULL, "-migrate_view_pre")); in MigrateParticles()
1599 PetscCall(DMSwarmCreateGlobalVectorFromField(sw, DMSwarmPICField_coor, &gc)); in MigrateParticles()
1600 PetscCall(DMSwarmCreateGlobalVectorFromField(sw, "velocity", &gv)); in MigrateParticles()
1604 PetscCall(DMSwarmGetCellDM(sw, &cdm)); in MigrateParticles()
1610 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in MigrateParticles()
1635 PetscCall(DMSwarmDestroyGlobalVectorFromField(sw, "velocity", &gv)); in MigrateParticles()
1636 PetscCall(DMSwarmDestroyGlobalVectorFromField(sw, DMSwarmPICField_coor, &gc)); in MigrateParticles()
1638 PetscCall(DMSwarmMigrate(sw, PETSC_TRUE)); in MigrateParticles()
1644 PetscCall(DMSwarmVectorDefineFields(sw, 2, fieldnames)); in MigrateParticles()
1651 DM dm, sw; in main() local
1663 PetscCall(CreateSwarm(dm, &user, &sw)); in main()
1665 PetscCall(InitializeConstants(sw, &user)); in main()
1666 PetscCall(DMSetApplicationContext(sw, &user)); in main()
1670 PetscCall(TSSetDM(ts, sw)); in main()
1685 PetscCall(SetupContext(dm, sw, &user)); in main()
1691 PetscCall(CheckNonNegativeWeights(sw, &user)); in main()
1699 PetscCall(DMDestroy(&sw)); in main()