Fix matrix generation

1 month ago · d4e3d13a7d
parent 03c23d60ae
commit d4e3d13a7d
4 changed files with 6 additions and 6 deletions
--- a/samples/advanced/pdegen/amg_d_genpde_mod.F90
+++ b/samples/advanced/pdegen/amg_d_genpde_mod.F90
@ -176,7 +176,7 @@ contains
    else
      partition_ = 3
    end if
-    deltah   = done/(idim+2)
+    deltah   = done/(idim+1)
    sqdeltah = deltah*deltah
    deltah2  = 2.0_psb_dpk_* deltah
@ -643,7 +643,7 @@ contains
      f_ => d_null_func_2d
    end if
-    deltah   = done/(idim+2)
+    deltah   = done/(idim+1)
    sqdeltah = deltah*deltah
    deltah2  = 2.0_psb_dpk_* deltah
--- a/samples/advanced/pdegen/amg_s_genpde_mod.F90
+++ b/samples/advanced/pdegen/amg_s_genpde_mod.F90
@ -176,7 +176,7 @@ contains
    else
      partition_ = 3
    end if
-    deltah   = sone/(idim+2)
+    deltah   = sone/(idim+1)
    sqdeltah = deltah*deltah
    deltah2  = 2.0_psb_spk_* deltah
@ -643,7 +643,7 @@ contains
      f_ => s_null_func_2d
    end if
-    deltah   = sone/(idim+2)
+    deltah   = sone/(idim+1)
    sqdeltah = deltah*deltah
    deltah2  = 2.0_psb_spk_* deltah
--- a/samples/simple/pdegen/amg_dpde_mod.f90
+++ b/samples/simple/pdegen/amg_dpde_mod.f90
@ -144,7 +144,7 @@ contains
      f_ => d_null_func_3d
    end if
-    deltah   = 1.d0/(idim+2)
+    deltah   = 1.d0/(idim+1)
    sqdeltah = deltah*deltah
    deltah2  = 2.0_psb_dpk_* deltah
--- a/samples/simple/pdegen/amg_spde_mod.f90
+++ b/samples/simple/pdegen/amg_spde_mod.f90
@ -144,7 +144,7 @@ contains
      f_ => s_null_func_3d
    end if
-    deltah   = 1.d0/(idim+2)
+    deltah   = 1.d0/(idim+1)
    sqdeltah = deltah*deltah
    deltah2  = 2.0_psb_spk_* deltah