67   PetscReal     *centroid, *coords, *velocity, *xi0, *v0, *J, *invJ, detJ, *vals;  in SetInitialCoordinates()  local
89   PetscCall(DMSwarmGetField(sw, "velocity", NULL, NULL, (void **)&velocity));  in SetInitialCoordinates()
120         velocity[p * dim + d] = v_val;  in SetInitialCoordinates()
125   PetscCall(DMSwarmRestoreField(sw, "velocity", NULL, NULL, (void **)&velocity));  in SetInitialCoordinates()
138   PetscReal   *velocity;  in SetInitialConditions()  local
149   PetscCall(DMSwarmGetField(dmSw, "velocity", NULL, NULL, (void **)&velocity));  in SetInitialConditions()
154       for (d = 0; d < dim; d++) initialConditions[n * dim + d] = velocity[n * dim + d];  in SetInitialConditions()
158   PetscCall(DMSwarmRestoreField(dmSw, "velocity", NULL, NULL, (void **)&velocity));  in SetInitialConditions()
287 …eGradS(PetscInt dim, PetscInt Np, const PetscReal vp[], const PetscReal velocity[], PetscReal inte…  in ComputeGradS()  argument
304         for (q = 0; q < Np; ++q) sum += Gaussian(dim, &velocity[q * dim], epsilon, vc_l);  in ComputeGradS()
339   PetscReal         *velocity;  in RHSFunctionParticles()  local
348   PetscCall(VecGetArray(sol, &velocity));  in RHSFunctionParticles()
356     PetscCall(ComputeGradS(dim, Np, &velocity[p * dim], velocity, gradS_p, user));  in RHSFunctionParticles()
361       PetscCall(ComputeGradS(dim, Np, &velocity[q * dim], velocity, gradS_q, user));  in RHSFunctionParticles()
379   PetscCall(VecRestoreArray(sol, &velocity));  in RHSFunctionParticles()
393   PetscScalar       *velocity;  in UpdateSwarm()  local
399   PetscCall(DMSwarmGetField(sw, "velocity", NULL, NULL, (void **)&velocity));  in UpdateSwarm()
403   for (idx = 0; idx < n; ++idx) velocity[idx] = u[idx];  in UpdateSwarm()
405   PetscCall(DMSwarmRestoreField(sw, "velocity", NULL, NULL, (void **)&velocity));  in UpdateSwarm()