199   PetscScalar *uarray, *f, *uoldarray, h, uave, c;  in IFunction_LaxFriedrichs()  local
219   PetscCall(DMDAVecGetArrayRead(da, localUold, &uoldarray));  in IFunction_LaxFriedrichs()
226     uave = 0.5 * (uoldarray[i - 1] + uoldarray[i + 1]);  in IFunction_LaxFriedrichs()
227     f[i] = uarray[i] - uave + c * 0.5 * (uoldarray[i + 1] - uoldarray[i - 1]);  in IFunction_LaxFriedrichs()
232   PetscCall(DMDAVecRestoreArrayRead(da, localUold, &uoldarray));  in IFunction_LaxFriedrichs()
247   PetscScalar *uarray, *f, *uoldarray, h, RFlux, LFlux, lambda;  in IFunction_LaxWendroff()  local
267   PetscCall(DMDAVecGetArrayRead(da, localUold, &uoldarray));  in IFunction_LaxWendroff()
275 …RFlux = 0.5 * a * (uoldarray[i + 1] + uoldarray[i]) - a * a * 0.5 * lambda * (uoldarray[i + 1] - u…  in IFunction_LaxWendroff()
276 …LFlux = 0.5 * a * (uoldarray[i - 1] + uoldarray[i]) - a * a * 0.5 * lambda * (uoldarray[i] - uolda…  in IFunction_LaxWendroff()
277     f[i]  = uarray[i] - uoldarray[i] + lambda * (RFlux - LFlux);  in IFunction_LaxWendroff()
282   PetscCall(DMDAVecRestoreArrayRead(da, localUold, &uoldarray));  in IFunction_LaxWendroff()