333   PetscScalar            *velocityValues, *velocityValuesOld;  in CharacteristicSolve()  local
357   PetscCall(PetscMalloc1(c->numVelocityComp, &velocityValues));  in CharacteristicSolve()
382 …l(velocityArray, interpIndices, c->numVelocityComp, c->velocityComp, velocityValues, c->velocityCt…  in CharacteristicSolve()
383 …erp(c->velocity, interpIndices, c->numVelocityComp, c->velocityComp, velocityValues, c->velocityCt…  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()
410 …l(velocityArray, interpIndices, c->numVelocityComp, c->velocityComp, velocityValues, c->velocityCt…  in CharacteristicSolve()
413 …erp(c->velocity, interpIndices, c->numVelocityComp, c->velocityComp, velocityValues, c->velocityCt…  in CharacteristicSolve()
416       Qi.x = 0.5 * (velocityValues[0] + velocityValuesOld[0]);  in CharacteristicSolve()
417       Qi.y = 0.5 * (velocityValues[1] + velocityValuesOld[1]);  in CharacteristicSolve()
435 …l(velocityArray, interpIndices, c->numVelocityComp, c->velocityComp, velocityValues, c->velocityCt…  in CharacteristicSolve()
438 …erp(c->velocity, interpIndices, c->numVelocityComp, c->velocityComp, velocityValues, c->velocityCt…  in CharacteristicSolve()
441     Qi.x              = 0.5 * (velocityValues[0] + velocityValuesOld[0]);  in CharacteristicSolve()
442     Qi.y              = 0.5 * (velocityValues[1] + velocityValuesOld[1]);  in CharacteristicSolve()
551   PetscCall(PetscFree(velocityValues));  in CharacteristicSolve()