295   DMSwarmCellDM b;  in DMSwarmCellDMCreate()  local
305 …PetscCall(PetscHeaderCreate(b, DMSWARMCELLDM_CLASSID, "DMSwarmCellDM", "Background DM for a Swarm"…  in DMSwarmCellDMCreate()
307   PetscCall(PetscObjectSetName((PetscObject)b, name));  in DMSwarmCellDMCreate()
309   b->dm  = dm;  in DMSwarmCellDMCreate()
310   b->Nf  = Nf;  in DMSwarmCellDMCreate()
311   b->Nfc = Nfc;  in DMSwarmCellDMCreate()
312   PetscCall(PetscMalloc1(b->Nf, &b->dmFields));  in DMSwarmCellDMCreate()
313   for (PetscInt f = 0; f < b->Nf; ++f) PetscCall(PetscStrallocpy(dmFields[f], &b->dmFields[f]));  in DMSwarmCellDMCreate()
314   PetscCall(PetscMalloc1(b->Nfc, &b->coordFields));  in DMSwarmCellDMCreate()
315 …for (PetscInt f = 0; f < b->Nfc; ++f) PetscCall(PetscStrallocpy(coordFields[f], &b->coordFields[f]…  in DMSwarmCellDMCreate()
317   PetscCall(PetscStrallocpy(cellid, &b->cellid));  in DMSwarmCellDMCreate()
318   *celldm = b;  in DMSwarmCellDMCreate()
348   PetscInt           i, j, k, bs, b, n_estimate, n_curr, n_new_est, p, n_found, Nfc;  in DMSwarmSetPointsUniformCoordinates()  local
372   for (b = 0; b < bs; b++) {  in DMSwarmSetPointsUniformCoordinates()
373     if (npoints[b] > 1) {  in DMSwarmSetPointsUniformCoordinates()
374       dx[b] = (max[b] - min[b]) / ((PetscReal)(npoints[b] - 1));  in DMSwarmSetPointsUniformCoordinates()
376       dx[b] = 0.0;  in DMSwarmSetPointsUniformCoordinates()
378     _npoints[b] = npoints[b];  in DMSwarmSetPointsUniformCoordinates()
393         for (b = 0; b < bs; b++) xp[b] = min[b] + ijk[b] * dx[b];  in DMSwarmSetPointsUniformCoordinates()
394         for (b = 0; b < bs; b++) {  in DMSwarmSetPointsUniformCoordinates()
395           if (xp[b] < lmin[b]) point_inside = PETSC_FALSE;  in DMSwarmSetPointsUniformCoordinates()
396           if (xp[b] > lmax[b]) point_inside = PETSC_FALSE;  in DMSwarmSetPointsUniformCoordinates()
421         for (b = 0; b < bs; b++) xp[b] = min[b] + ijk[b] * dx[b];  in DMSwarmSetPointsUniformCoordinates()
422         for (b = 0; b < bs; b++) {  in DMSwarmSetPointsUniformCoordinates()
423           if (xp[b] < lmin[b]) point_inside = PETSC_FALSE;  in DMSwarmSetPointsUniformCoordinates()
424           if (xp[b] > lmax[b]) point_inside = PETSC_FALSE;  in DMSwarmSetPointsUniformCoordinates()
427           for (b = 0; b < bs; b++) _pos[bs * n_estimate + b] = xp[b];  in DMSwarmSetPointsUniformCoordinates()
463 …for (b = 0; b < bs; b++) swarm_coor[bs * (n_curr + n_found) + b] = PetscRealPart(_coor[bs * p + b]…  in DMSwarmSetPointsUniformCoordinates()
502   PetscInt           i, N, bs, b, n_estimate, n_curr, n_new_est, p, n_found;  in DMSwarmSetPointCoordinates()  local
535     for (b = 0; b < bs; b++) {  in DMSwarmSetPointCoordinates()
536       gmin[b] = PetscMin(gmin[b], PetscRealPart(_coor[bs * i + b]));  in DMSwarmSetPointCoordinates()
537       gmax[b] = PetscMax(gmax[b], PetscRealPart(_coor[bs * i + b]));  in DMSwarmSetPointCoordinates()
566     for (b = 0; b < bs; b++) {  in DMSwarmSetPointCoordinates()
567       if (my_coor[bs * i + b] < gmin[b]) point_inside = PETSC_FALSE;  in DMSwarmSetPointCoordinates()
568       if (my_coor[bs * i + b] > gmax[b]) point_inside = PETSC_FALSE;  in DMSwarmSetPointCoordinates()
584     for (b = 0; b < bs; b++) {  in DMSwarmSetPointCoordinates()
585       if (my_coor[bs * i + b] < gmin[b]) point_inside = PETSC_FALSE;  in DMSwarmSetPointCoordinates()
586       if (my_coor[bs * i + b] > gmax[b]) point_inside = PETSC_FALSE;  in DMSwarmSetPointCoordinates()
589       for (b = 0; b < bs; b++) _pos[bs * n_estimate + b] = my_coor[bs * i + b];  in DMSwarmSetPointCoordinates()
624 …for (b = 0; b < bs; b++) swarm_coor[bs * (n_curr + n_found) + b] = PetscRealPart(_coor[bs * p + b]…  in DMSwarmSetPointCoordinates()
1381 …     for (PetscInt b = 0; b < cbs[i]; ++b, ++off) xr[off] = coordVals[i][pindices[p] * cbs[i] + b];  in DMProjectFieldLocal_Swarm()  local
1387 …for (PetscInt b = 0; b < bs[i]; ++b, ++off) val[p * totbs + off] = vals[i][pindices[p] * bs[i] + b…  in DMProjectFieldLocal_Swarm()  local