78   PetscScalar    Hx, Hy, Hz, HxHydHz, HyHzdHx, HxHzdHy;  in ComputeRHS()  local
87   Hz      = 1.0 / (PetscReal)(mz - 1);  in ComputeRHS()
88   HxHydHz = Hx * Hy / Hz;  in ComputeRHS()
89   HxHzdHy = Hx * Hz / Hy;  in ComputeRHS()
90   HyHzdHx = Hy * Hz / Hx;  in ComputeRHS()
100           barray[k][j][i] = Hx * Hy * Hz;  in ComputeRHS()
121   PetscScalar v[7], Hx, Hy, Hz, HxHydHz, HyHzdHx, HxHzdHy;  in ComputeMatrix()  local
130   Hz      = 1.0 / (PetscReal)(mz - 1);  in ComputeMatrix()
131   HxHydHz = Hx * Hy / Hz;  in ComputeMatrix()
132   HxHzdHy = Hx * Hz / Hy;  in ComputeMatrix()
133   HyHzdHx = Hy * Hz / Hx;  in ComputeMatrix()