524   const PetscReal   *coords, *vel;  in RHSFunction()  local
534   PetscCall(DMSwarmGetField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSFunction()
545     const PetscReal vy0   = vel[p * dim + 1];  in RHSFunction()
554   PetscCall(DMSwarmRestoreField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSFunction()
570   const PetscReal *coords, *vel;  in RHSJacobian()  local
579   PetscCall(DMSwarmGetField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSJacobian()
583     const PetscReal vy0     = vel[p * dim + 1];  in RHSJacobian()
593   PetscCall(DMSwarmRestoreField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSJacobian()
626   const PetscReal   *coords, *vel;  in RHSFunctionV()  local
635   PetscCall(DMSwarmGetField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSFunctionV()
647     const PetscReal vy0   = vel[p * dim + 1];  in RHSFunctionV()
655   PetscCall(DMSwarmRestoreField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSFunctionV()
690   PetscReal         *vel, *coords, m_p = 1., q_p = -1.;  in RHSObjectiveF()  local
706   PetscCall(DMSwarmGetField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSObjectiveF()
720       const PetscReal vy0   = vel[p * dim + 1];  in RHSObjectiveF()
732   PetscCall(DMSwarmRestoreField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSObjectiveF()
745   const PetscReal   *coords, *vel;  in RHSFunctionG()  local
755   PetscCall(DMSwarmGetField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSFunctionG()
768     const PetscReal vy0   = vel[p * dim + 1];  in RHSFunctionG()
776   PetscCall(DMSwarmRestoreField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in RHSFunctionG()
966   const PetscReal   *coords, *vel;  in ComputeError()  local
981   PetscCall(DMSwarmGetField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in ComputeError()
987     const PetscReal    v0    = DMPlex_NormD_Internal(dim, &vel[p * dim]);  in ComputeError()
1008   PetscCall(DMSwarmRestoreField(sw, "initVelocity", NULL, NULL, (void **)&vel));  in ComputeError()