fdurastante
/
psblas3
mirror of https://github.com/sfilippone/psblas3.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Fix sample matrix generators

Browse Source
merge-maint-par
sfilippone 2 months ago
parent 6d4ed9575a
commit 6c9337c728
4 changed files with 14 additions and 14 deletions
Show Stats Download Patch File Download Diff File

6
test/pdegen/psb_d_pde2d.F90
Unescape Escape View File

@ -230,7 +230,7 @@ contains
f_ => d_null_func_2d f_ => d_null_func_2d
end if end if
deltah = done/(idim+1) deltah = done/(idim+2)
sqdeltah = deltah*deltah sqdeltah = deltah*deltah
deltah2 = (2*done)* deltah deltah2 = (2*done)* deltah
@ -467,8 +467,8 @@ contains
! compute gridpoint coordinates ! compute gridpoint coordinates
call idx2ijk(ix,iy,glob_row,idim,idim) call idx2ijk(ix,iy,glob_row,idim,idim)
! x, y coordinates ! x, y coordinates
x = (ix-1)*deltah x = (ix)*deltah
y = (iy-1)*deltah y = (iy)*deltah
zt(k) = f_(x,y) zt(k) = f_(x,y)
! internal point: build discretization ! internal point: build discretization

8
test/pdegen/psb_d_pde3d.F90
Unescape Escape View File

@ -246,7 +246,7 @@ contains
f_ => d_null_func_3d f_ => d_null_func_3d
end if end if
deltah = done/(idim+1) deltah = done/(idim+2)
sqdeltah = deltah*deltah sqdeltah = deltah*deltah
deltah2 = (2*done)* deltah deltah2 = (2*done)* deltah
@ -496,9 +496,9 @@ contains
! compute gridpoint coordinates ! compute gridpoint coordinates
call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim) call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim)
! x, y, z coordinates ! x, y, z coordinates
x = (ix-1)*deltah x = (ix)*deltah
y = (iy-1)*deltah y = (iy)*deltah
z = (iz-1)*deltah z = (iz)*deltah
zt(k) = f_(x,y,z) zt(k) = f_(x,y,z)
! internal point: build discretization ! internal point: build discretization
! !

6
test/pdegen/psb_s_pde2d.F90
Unescape Escape View File

@ -230,7 +230,7 @@ contains
f_ => s_null_func_2d f_ => s_null_func_2d
end if end if
deltah = sone/(idim+1) deltah = sone/(idim+2)
sqdeltah = deltah*deltah sqdeltah = deltah*deltah
deltah2 = (2*sone)* deltah deltah2 = (2*sone)* deltah
@ -467,8 +467,8 @@ contains
! compute gridpoint coordinates ! compute gridpoint coordinates
call idx2ijk(ix,iy,glob_row,idim,idim) call idx2ijk(ix,iy,glob_row,idim,idim)
! x, y coordinates ! x, y coordinates
x = (ix-1)*deltah x = (ix)*deltah
y = (iy-1)*deltah y = (iy)*deltah
zt(k) = f_(x,y) zt(k) = f_(x,y)
! internal point: build discretization ! internal point: build discretization

8
test/pdegen/psb_s_pde3d.F90
Unescape Escape View File

@ -246,7 +246,7 @@ contains
f_ => s_null_func_3d f_ => s_null_func_3d
end if end if
deltah = sone/(idim+1) deltah = sone/(idim+2)
sqdeltah = deltah*deltah sqdeltah = deltah*deltah
deltah2 = (2*sone)* deltah deltah2 = (2*sone)* deltah
@ -496,9 +496,9 @@ contains
! compute gridpoint coordinates ! compute gridpoint coordinates
call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim) call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim)
! x, y, z coordinates ! x, y, z coordinates
x = (ix-1)*deltah x = (ix)*deltah
y = (iy-1)*deltah y = (iy)*deltah
z = (iz-1)*deltah z = (iz)*deltah
zt(k) = f_(x,y,z) zt(k) = f_(x,y,z)
! internal point: build discretization ! internal point: build discretization
! !

Write Preview
Loading…
Cancel
Save