Lines Matching refs:uarray
161 PetscScalar u, uxx, uyy, **uarray, **f, **udot; in FormIFunction() local
184 PetscCall(DMDAVecGetArrayRead(da, localU, &uarray)); in FormIFunction()
197 f[j][i] = uarray[j][i]; /* F = U */ in FormIFunction()
200 f[j][i] = uarray[j][i] - uarray[j + 1][i + 1]; in FormIFunction()
202 f[j][i] = uarray[j][i] - uarray[j + 1][i - 1]; in FormIFunction()
204 f[j][i] = uarray[j][i] - uarray[j - 1][i + 1]; in FormIFunction()
206 f[j][i] = uarray[j][i] - uarray[j - 1][i - 1]; in FormIFunction()
208 f[j][i] = uarray[j][i] - uarray[j][i + 1]; in FormIFunction()
210 f[j][i] = uarray[j][i] - uarray[j][i - 1]; in FormIFunction()
212 f[j][i] = uarray[j][i] - uarray[j + 1][i]; in FormIFunction()
214 f[j][i] = uarray[j][i] - uarray[j - 1][i]; in FormIFunction()
218 u = uarray[j][i]; in FormIFunction()
220 uxx = (-2.0 * u + uarray[j][i - 1] + uarray[j][i + 1]); in FormIFunction()
221 uyy = (-2.0 * u + uarray[j - 1][i] + uarray[j + 1][i]); in FormIFunction()
224 …uxx = 2.0 * uxx / 3.0 + (0.5 * (uarray[j - 1][i - 1] + uarray[j - 1][i + 1] + uarray[j + 1][i - 1]… in FormIFunction()
225 …uyy = 2.0 * uyy / 3.0 + (0.5 * (uarray[j - 1][i - 1] + uarray[j - 1][i + 1] + uarray[j + 1][i - 1]… in FormIFunction()
233 PetscCall(DMDAVecRestoreArrayRead(da, localU, &uarray)); in FormIFunction()