58 PetscErrorCode DMSwarmInterpolateFromCellToSwarm_Petsc(DM dm_swarm, const char *field, Vec U_mesh);
59 PetscErrorCode DMSwarmInterpolateFromCellToSwarm_Ceed(DM dm_swarm, const char *field, Vec U_mesh);
60 PetscErrorCode DMSwarmCheckSwarmValues(DM dm_swarm, const char *field, PetscScalar tolerance, Targe…
71   DM                 dm_mesh, dm_swarm;  in main()  local
177     PetscCall(DMCreate(comm, &dm_swarm));  in main()
178     PetscCall(DMSetType(dm_swarm, DMSWARM));  in main()
179     PetscCall(DMSetDimension(dm_swarm, dim));  in main()
180     PetscCall(DMSwarmSetType(dm_swarm, DMSWARM_PIC));  in main()
181     PetscCall(DMSwarmSetCellDM(dm_swarm, dm_mesh));  in main()
184 … PetscCall(DMSwarmRegisterPetscDatatypeField(dm_swarm, DMSwarmPICField_u, num_comp, PETSC_SCALAR));  in main()
185     PetscCall(DMSwarmFinalizeFieldRegister(dm_swarm));  in main()
186     PetscCall(DMSwarmSetLocalSizes(dm_swarm, num_points, 0));  in main()
187     PetscCall(DMSetFromOptions(dm_swarm));  in main()
190 …PetscCall(DMSwarmInitalizePointLocations(dm_swarm, point_swarm_type, num_points, num_points_per_ce…  in main()
193     PetscCall(PetscObjectSetName((PetscObject)dm_swarm, "Particle Swarm"));  in main()
194     PetscCall(DMViewFromOptions(dm_swarm, NULL, "-dm_swarm_view"));  in main()
221   PetscCall(DMSwarmCeedContextCreate(dm_swarm, ceed_resource, &swarm_ceed_context));  in main()
225     PetscCall(DMSwarmInterpolateFromCellToSwarm_Petsc(dm_swarm, DMSwarmPICField_u, U_mesh));  in main()
226     if (view_petsc_swarm) PetscCall(DMSwarmViewXDMF(dm_swarm, "swarm_petsc.xmf"));  in main()
227 …  PetscCall(DMSwarmCheckSwarmValues(dm_swarm, DMSwarmPICField_u, tolerance, target_function_proj));  in main()
232     PetscCall(DMSwarmInterpolateFromCellToSwarm_Ceed(dm_swarm, DMSwarmPICField_u, U_mesh));  in main()
233     if (view_ceed_swarm) PetscCall(DMSwarmViewXDMF(dm_swarm, "swarm_ceed.xmf"));  in main()
234 …  PetscCall(DMSwarmCheckSwarmValues(dm_swarm, DMSwarmPICField_u, tolerance, target_function_proj));  in main()
242     PetscCall(DMSwarmProjectFromSwarmToCells(dm_swarm, DMSwarmPICField_u, NULL, U_projected));  in main()
264   PetscCall(DMDestroy(&dm_swarm));  in main()
329 PetscErrorCode DMSwarmInterpolateFromCellToSwarm_Petsc(DM dm_swarm, const char *field, Vec U_mesh) {  in DMSwarmInterpolateFromCellToSwarm_Petsc()  argument
344   PetscCall(DMSwarmGetCellDM(dm_swarm, &dm_mesh));  in DMSwarmInterpolateFromCellToSwarm_Petsc()
360   PetscCall(DMSwarmSortGetAccess(dm_swarm));  in DMSwarmInterpolateFromCellToSwarm_Petsc()
362   PetscCall(DMSwarmGetField(dm_swarm, DMSwarmPICField_coor, NULL, NULL, (void **)&coords_points));  in DMSwarmInterpolateFromCellToSwarm_Petsc()
363   PetscCall(DMSwarmGetField(dm_swarm, field, &num_comp, NULL, (void **)&u_points));  in DMSwarmInterpolateFromCellToSwarm_Petsc()
375     PetscCall(DMSwarmSortGetPointsPerCell(dm_swarm, cell, &num_points_in_cell, &points_cell));  in DMSwarmInterpolateFromCellToSwarm_Petsc()
401     PetscCall(DMSwarmSortRestorePointsPerCell(dm_swarm, cell, &num_points_in_cell, &points_cell));  in DMSwarmInterpolateFromCellToSwarm_Petsc()
405 …PetscCall(DMSwarmRestoreField(dm_swarm, DMSwarmPICField_coor, NULL, NULL, (void **)&coords_points)…  in DMSwarmInterpolateFromCellToSwarm_Petsc()
406   PetscCall(DMSwarmRestoreField(dm_swarm, field, NULL, NULL, (void **)&u_points));  in DMSwarmInterpolateFromCellToSwarm_Petsc()
407   PetscCall(DMSwarmSortRestoreAccess(dm_swarm));  in DMSwarmInterpolateFromCellToSwarm_Petsc()
416 PetscErrorCode DMSwarmInterpolateFromCellToSwarm_Ceed(DM dm_swarm, const char *field, Vec U_mesh) {  in DMSwarmInterpolateFromCellToSwarm_Ceed()  argument
424   PetscCall(DMSwarmGetCellDM(dm_swarm, &dm_mesh));  in DMSwarmInterpolateFromCellToSwarm_Ceed()
434   PetscCall(DMSwarmSortGetAccess(dm_swarm));  in DMSwarmInterpolateFromCellToSwarm_Ceed()
435   PetscCall(DMSwarmPICFieldP2C(dm_swarm, field, swarm_ceed_context->u_points));  in DMSwarmInterpolateFromCellToSwarm_Ceed()
441   PetscCall(DMSwarmPICFieldC2P(dm_swarm, field, swarm_ceed_context->u_points));  in DMSwarmInterpolateFromCellToSwarm_Ceed()
442   PetscCall(DMSwarmSortRestoreAccess(dm_swarm));  in DMSwarmInterpolateFromCellToSwarm_Ceed()
451 PetscErrorCode DMSwarmCheckSwarmValues(DM dm_swarm, const char *field, PetscScalar tolerance, Targe…  in DMSwarmCheckSwarmValues()  argument
458   PetscCall(DMSwarmSortGetAccess(dm_swarm));  in DMSwarmCheckSwarmValues()
459   PetscCall(DMSwarmGetCellDM(dm_swarm, &dm_mesh));  in DMSwarmCheckSwarmValues()
462   PetscCall(DMSwarmGetField(dm_swarm, DMSwarmPICField_coor, NULL, NULL, (void **)&coords_points));  in DMSwarmCheckSwarmValues()
463   PetscCall(DMSwarmGetField(dm_swarm, field, &num_comp, NULL, (void **)&u_points));  in DMSwarmCheckSwarmValues()
470     PetscCall(DMSwarmSortGetPointsPerCell(dm_swarm, cell, &num_points_in_cell, &points));  in DMSwarmCheckSwarmValues()
480           PetscCall(PetscPrintf(PetscObjectComm((PetscObject)dm_swarm),  in DMSwarmCheckSwarmValues()
489     PetscCall(DMSwarmSortRestorePointsPerCell(dm_swarm, cell, &num_points_in_cell, &points));  in DMSwarmCheckSwarmValues()
493 …PetscCall(DMSwarmRestoreField(dm_swarm, DMSwarmPICField_coor, NULL, NULL, (void **)&coords_points)…  in DMSwarmCheckSwarmValues()
494   PetscCall(DMSwarmRestoreField(dm_swarm, field, NULL, NULL, (void **)&u_points));  in DMSwarmCheckSwarmValues()
495   PetscCall(DMSwarmSortRestoreAccess(dm_swarm));  in DMSwarmCheckSwarmValues()
496 …PetscCheck(within_tolerance, PetscObjectComm((PetscObject)dm_swarm), PETSC_ERR_USER, "Interpolatio…  in DMSwarmCheckSwarmValues()