323   CharacteristicPointDA2D Qi;  in CharacteristicSolve()  local
378   for (Qi.j = js; Qi.j < je; Qi.j++) {  in CharacteristicSolve()
379     for (Qi.i = is; Qi.i < ie; Qi.i++) {  in CharacteristicSolve()
380       interpIndices[0] = Qi.i;  in CharacteristicSolve()
381       interpIndices[1] = Qi.j;  in CharacteristicSolve()
384       Qi.x = Qi.i - velocityValues[0] * dt / 2.0;  in CharacteristicSolve()
385       Qi.y = Qi.j - velocityValues[1] * dt / 2.0;  in CharacteristicSolve()
388       Qi.proc = DMDAGetNeighborRelative(da, Qi.x, Qi.y);  in CharacteristicSolve()
391       PetscCall(DMDAMapCoordsToPeriodicDomain(da, &Qi.x, &Qi.y));  in CharacteristicSolve()
392       PetscCall(CharacteristicAddPoint(c, &Qi));  in CharacteristicSolve()
405     Qi = c->queue[n];  in CharacteristicSolve()
406     if (c->neighbors[Qi.proc] == rank) {  in CharacteristicSolve()
407       interpIndices[0] = Qi.x;  in CharacteristicSolve()
408       interpIndices[1] = Qi.y;  in CharacteristicSolve()
416       Qi.x = 0.5 * (velocityValues[0] + velocityValuesOld[0]);  in CharacteristicSolve()
417       Qi.y = 0.5 * (velocityValues[1] + velocityValuesOld[1]);  in CharacteristicSolve()
419     c->queue[n] = Qi;  in CharacteristicSolve()
431     Qi               = c->queueRemote[n];  in CharacteristicSolve()
432     interpIndices[0] = Qi.x;  in CharacteristicSolve()
433     interpIndices[1] = Qi.y;  in CharacteristicSolve()
441     Qi.x              = 0.5 * (velocityValues[0] + velocityValuesOld[0]);  in CharacteristicSolve()
442     Qi.y              = 0.5 * (velocityValues[1] + velocityValuesOld[1]);  in CharacteristicSolve()
443     c->queueRemote[n] = Qi;  in CharacteristicSolve()
465     Qi   = c->queue[n];  in CharacteristicSolve()
466     Qi.x = Qi.i - Qi.x * dt;  in CharacteristicSolve()
467     Qi.y = Qi.j - Qi.y * dt;  in CharacteristicSolve()
470     Qi.proc = DMDAGetNeighborRelative(da, Qi.x, Qi.y);  in CharacteristicSolve()
473     PetscCall(DMDAMapCoordsToPeriodicDomain(da, &Qi.x, &Qi.y));  in CharacteristicSolve()
475     c->queue[n] = Qi;  in CharacteristicSolve()
542     Qi = c->queue[n];  in CharacteristicSolve()
543 …or (comp = 0; comp < c->numFieldComp; comp++) solArray[Qi.j][Qi.i * dof + c->fieldComp[comp]] = Qi…  in CharacteristicSolve()