249   for (int tid = 0; tid < numthreads; tid++) {  in go()  local
253     PetscCall(DMCreate(PETSC_COMM_SELF, &dm_t[tid]));  in go()
254     PetscCall(DMSetType(dm_t[tid], DMPLEX));  in go()
255     PetscCall(DMSetFromOptions(dm_t[tid]));  in go()
259     PetscCall(DMSetField(dm_t[tid], field, NULL, (PetscObject)fe));  in go()
260     PetscCall(DMCreateDS(dm_t[tid]));  in go()
263     PetscCall(DMCreateGlobalVector(dm_t[tid], &rho_t[tid]));  in go()
264     PetscCall(DMCreateGlobalVector(dm_t[tid], &rhs_t[tid]));  in go()
266     PetscCall(DMGetBoundingBox(dm_t[tid], lo, hi));  in go()
267     if (tid == target) {  in go()
268       PetscCall(DMViewFromOptions(dm_t[tid], NULL, "-dm_view"));  in go()
273 …PetscCall(PetscInfo(dm_t[tid], " lo = %g hi = %g n = %" PetscInt_FMT " h = %g hp = %g\n", (double)…  in go()
282   for (int tid = 0; tid < numthreads; tid++) { // change size of particle list a little  in go()  local
283     Np_t[tid] = Np2[0] * Np2[1];  in go()
284     PetscCall(PetscMalloc3(Np_t[tid], &xx_t[tid], Np_t[tid], &yy_t[tid], Np_t[tid], &wp_t[tid]));  in go()
285     if (tid == target) moments_0[0] = moments_0[1] = moments_0[2] = 0;  in go()
288         xx_t[tid][pp] = lo[0] + hp[0] / 2. + pi * hp[0];  in go()
289         yy_t[tid][pp] = lo[1] + hp[1] / 2. + pj * hp[1];  in go()
291           PetscReal x[] = {xx_t[tid][pp], yy_t[tid][pp]};  in go()
292           PetscCall(maxwellian(2, x, 1.0, vol / (PetscReal)Np_t[tid], &wp_t[tid][pp]));  in go()
294 …if (tid == target) { //energy_0 += wp_t[tid][pp]*(PetscSqr(xx_t[tid][pp])+PetscSqr(yy_t[tid][pp]));  in go()
295           moments_0[0] += wp_t[tid][pp];  in go()
296           moments_0[1] += wp_t[tid][pp] * xx_t[tid][pp]; // x-momentum  in go()
297           moments_0[2] += wp_t[tid][pp] * (PetscSqr(xx_t[tid][pp]) + PetscSqr(yy_t[tid][pp]));  in go()
308 …for (int tid = 0; tid < numthreads; tid++) PetscCallAbort(PETSC_COMM_SELF, createSwarm(dm_t[tid], …  in go()  local
310 …tid = 0; tid < numthreads; tid++) PetscCallAbort(PETSC_COMM_SELF, particlesToGrid(dm_t[tid], sw_t[…  in go()  local
314 …for (int tid = 0; tid < numthreads; tid++) PetscCallAbort(PETSC_COMM_SELF, VecCopy(rho_t[tid], rhs…  in go()  local
316 …tid = 0; tid < numthreads; tid++) PetscCallAbort(PETSC_COMM_SELF, gridToParticles(dm_t[tid], sw_t[…  in go()  local
318   for (int tid = 0; tid < numthreads; tid++) {  in go()  local
319     PetscCall(MatDestroy(&M_p_t[tid]));  in go()
320     PetscCall(DMDestroy(&sw_t[tid]));  in go()
326   for (int tid = 0; tid < numthreads; tid++) {  in go()  local
327     PetscCall(VecDestroy(&rho_t[tid]));  in go()
328     PetscCall(VecDestroy(&rhs_t[tid]));  in go()
329     PetscCall(DMDestroy(&dm_t[tid]));  in go()
330     PetscCall(PetscFree3(xx_t[tid], yy_t[tid], wp_t[tid]));  in go()