fdurastante/amg4psblas
1
0
Fork 0
mirror of https://github.com/sfilippone/amg4psblas.git synced 2026-10-06 06:34:57 +00:00
Code Issues Packages Projects Releases Wiki Activity

Fix sample matrix generators

Browse Source
This commit is contained in:
sfilippone
2026-05-05 15:52:03 +02:00
parent 5c055170e7
commit 3343b039e6
4 changed files with 18 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files
samples/advanced/pdegen/amg_d_genpde_mod.F90
+5 -5
View File
@@ -412,9 +412,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
!
@@ -875,8 +875,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
samples/advanced/pdegen/amg_s_genpde_mod.F90
+5 -5
View File
@@ -412,9 +412,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
!
@@ -875,8 +875,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
samples/simple/pdegen/amg_dpde_mod.f90
+4 -4
View File
@@ -144,7 +144,7 @@ contains
f_ => d_null_func_3d
end if
deltah = 1.d0/(idim+1)
deltah = 1.d0/(idim+2)
sqdeltah = deltah*deltah
deltah2 = 2.0_psb_dpk_* deltah
@@ -369,9 +369,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
!
samples/simple/pdegen/amg_spde_mod.f90
+4 -4
View File
@@ -144,7 +144,7 @@ contains
f_ => s_null_func_3d
end if
deltah = 1.d0/(idim+1)
deltah = 1.d0/(idim+2)
sqdeltah = deltah*deltah
deltah2 = 2.0_psb_spk_* deltah
@@ -369,9 +369,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
!