Lines Matching refs:sw

43 PetscErrorCode createSwarm(const DM dm, DM *sw)  in createSwarm()  argument
48 PetscCall(DMCreate(PETSC_COMM_SELF, sw)); in createSwarm()
49 PetscCall(DMSetType(*sw, DMSWARM)); in createSwarm()
50 PetscCall(DMSetDimension(*sw, dim)); in createSwarm()
51 PetscCall(DMSwarmSetType(*sw, DMSWARM_PIC)); in createSwarm()
52 PetscCall(DMSwarmSetCellDM(*sw, dm)); in createSwarm()
53 PetscCall(DMSwarmRegisterPetscDatatypeField(*sw, "w_q", Nc, PETSC_SCALAR)); in createSwarm()
54 PetscCall(DMSwarmFinalizeFieldRegister(*sw)); in createSwarm()
55 PetscCall(DMSetFromOptions(*sw)); in createSwarm()
59 PetscErrorCode gridToParticles(const DM dm, DM sw, PetscReal *moments, Vec rhs, Mat M_p) in gridToParticles() argument
117 PetscCall(DMSwarmCreateMassMatrixSquare(sw, dm, &PM_p)); in gridToParticles()
123 PetscCall(DMSwarmCreateGlobalVectorFromField(sw, "w_q", &ff)); // this grabs access !!!!! in gridToParticles()
138 PetscCall(DMSetOutputSequenceNumber(sw, timestep, time)); in gridToParticles()
140 PetscCall(DMSwarmDestroyGlobalVectorFromField(sw, "w_q", &ff)); in gridToParticles()
145 PetscCall(DMSwarmGetLocalSize(sw, &Np)); in gridToParticles()
146 PetscCall(DMSwarmGetField(sw, "w_q", &bs, &dtype, (void **)&wq)); in gridToParticles()
147 PetscCall(DMSwarmGetField(sw, "DMSwarmPIC_coor", &bs, &dtype, (void **)&coords)); in gridToParticles()
154 PetscCall(DMSwarmRestoreField(sw, "DMSwarmPIC_coor", &bs, &dtype, (void **)&coords)); in gridToParticles()
155 PetscCall(DMSwarmRestoreField(sw, "w_q", &bs, &dtype, (void **)&wq)); in gridToParticles()
161 PetscErrorCode particlesToGrid(const DM dm, DM sw, const PetscInt Np, const PetscInt a_tid, const P… in particlesToGrid() argument
171 PetscCall(DMSwarmSetLocalSizes(sw, Np, zero)); in particlesToGrid()
172 PetscCall(DMSwarmGetField(sw, "w_q", &bs, &dtype, (void **)&wq)); in particlesToGrid()
173 PetscCall(DMSwarmGetField(sw, "DMSwarmPIC_coor", &bs, &dtype, (void **)&coords)); in particlesToGrid()
179 PetscCall(DMSwarmRestoreField(sw, "DMSwarmPIC_coor", &bs, &dtype, (void **)&coords)); in particlesToGrid()
180 PetscCall(DMSwarmRestoreField(sw, "w_q", &bs, &dtype, (void **)&wq)); in particlesToGrid()
181 PetscCall(DMSwarmMigrate(sw, removePoints)); in particlesToGrid()
182 PetscCall(PetscObjectSetName((PetscObject)sw, "Particle Grid")); in particlesToGrid()
183 if (a_tid == target) PetscCall(DMViewFromOptions(sw, NULL, "-swarm_view")); in particlesToGrid()
187 PetscCall(DMCreateMassMatrix(sw, dm, &M_p)); in particlesToGrid()
190 PetscCall(DMSwarmCreateGlobalVectorFromField(sw, "w_q", &ff)); // this grabs access !!!!! in particlesToGrid()
193 PetscCall(DMSwarmDestroyGlobalVectorFromField(sw, "w_q", &ff)); in particlesToGrid()