--- trunk/doc/probe/paper.tex 2010/08/11 16:21:39 274 +++ trunk/doc/probe/paper.tex 2010/08/11 16:45:29 275 @@ -235,11 +235,11 @@ } double vx = 0.0 - for (int i = 1-s; i < s; i++) { + for (int i = 1-s; i <= s; i++) { double vy = 0.0 - for (int j = 1-s; j < s; j++) { + for (int j = 1-s; j <= s; j++) { double vz = 0.0; - for (int k = 1-s; k < s; k++) { + for (int k = 1-s; k <= s; k++) { vz += V[nx+i][ny+j][nz+k] * hz[k+s-1]; } vy += vz * hy[j+s-1]; @@ -280,8 +280,8 @@ hz = a + t*(b + t*(c + t*d)); double vy[4], vz[4]; - for (int i = 1-s; i < s; i++) { - for (int j = 1-s; j < s; j++) { + for (int i = 1-s; i <= s; i++) { + for (int j = 1-s; j <= s; j++) { double4 v = (double4) ( V[nx+i][ny+j][nz-1], V[nx+i][ny+j][nz],