Fix matrix generation

Makefile

@@ -1,5 +1,5 @@
 include Make.inc
-
+ 
 all: dirs mods objs libd
 	@echo "====================================="
 	@echo "PSBLAS libraries Compilation Successful."

test/pdegen/psb_d_pde2d.F90

@@ -467,8 +467,8 @@ contains
           ! compute gridpoint coordinates
           call idx2ijk(ix,iy,glob_row,idim,idim)
           ! x, y coordinates
-          x = (ix-1)*deltah
-          y = (iy-1)*deltah
+          x = (ix)*deltah
+          y = (iy)*deltah
           
           zt(k) = f_(x,y)
           ! internal point: build discretization

test/pdegen/psb_d_pde3d.F90

@@ -496,9 +496,9 @@ contains
           ! compute gridpoint coordinates
           call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim)
           ! x, y, z coordinates
-          x = (ix-1)*deltah
-          y = (iy-1)*deltah
-          z = (iz-1)*deltah
+          x = (ix)*deltah
+          y = (iy)*deltah
+          z = (iz)*deltah
           zt(k) = f_(x,y,z)
           ! internal point: build discretization
           !

test/pdegen/psb_s_pde2d.F90

@@ -467,8 +467,8 @@ contains
           ! compute gridpoint coordinates
           call idx2ijk(ix,iy,glob_row,idim,idim)
           ! x, y coordinates
-          x = (ix-1)*deltah
-          y = (iy-1)*deltah
+          x = (ix)*deltah
+          y = (iy)*deltah
           
           zt(k) = f_(x,y)
           ! internal point: build discretization

test/pdegen/psb_s_pde3d.F90

@@ -496,9 +496,9 @@ contains
           ! compute gridpoint coordinates
           call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim)
           ! x, y, z coordinates
-          x = (ix-1)*deltah
-          y = (iy-1)*deltah
-          z = (iz-1)*deltah
+          x = (ix)*deltah
+          y = (iy)*deltah
+          z = (iz)*deltah
           zt(k) = f_(x,y,z)
           ! internal point: build discretization
           !