229   PetscInt numthreads = PetscNumOMPThreads;  in go()  local
231   PetscInt numthreads = 1;  in go()  local
236 …PetscCheck(numthreads <= MAX_NUM_THRDS, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Too many thre…  in go()
237 …PetscCheck(numthreads > 0, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "No threads %" PetscInt_FMT…  in go()
239   if (target >= numthreads) target = numthreads - 1;  in go()
249   for (int tid = 0; tid < numthreads; tid++) {  in go()
282   for (int tid = 0; tid < numthreads; tid++) { // change size of particle list a little  in go()
308 …for (int tid = 0; tid < numthreads; tid++) PetscCallAbort(PETSC_COMM_SELF, createSwarm(dm_t[tid], …  in go()
310 …for (int tid = 0; tid < numthreads; tid++) PetscCallAbort(PETSC_COMM_SELF, particlesToGrid(dm_t[ti…  in go()
314 …for (int tid = 0; tid < numthreads; tid++) PetscCallAbort(PETSC_COMM_SELF, VecCopy(rho_t[tid], rhs…  in go()
316 …for (int tid = 0; tid < numthreads; tid++) PetscCallAbort(PETSC_COMM_SELF, gridToParticles(dm_t[ti…  in go()
318   for (int tid = 0; tid < numthreads; tid++) {  in go()
324 …)moments_1[2], (double)((moments_1[2] - moments_0[2]) / moments_0[2]), Np[0] * Np[1], numthreads));  in go()
326   for (int tid = 0; tid < numthreads; tid++) {  in go()