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

Fix test program

Browse Source
This commit is contained in:
sfilippone
2026-09-09 17:20:42 +02:00
parent e944da5bf7
commit 72a731ca78
4 changed files with 16 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files
samples/advanced/pdegen/amg_d_pde2d.F90
+4 -4
View File
@@ -87,7 +87,7 @@ program amg_d_pde2d
integer(psb_epk_) :: system_size
! miscellaneous
real(psb_dpk_) :: t1, t2, tprec, thier, tslv, tsmth, tpgen
real(psb_dpk_) :: t1, t2, thier, tslv, tsmth, tpgen
! sparse matrix and preconditioner
type(psb_dspmat_type) :: a
@@ -476,12 +476,12 @@ program amg_d_pde2d
end if
call psb_amx(ctxt, thier)
call psb_amx(ctxt, tprec)
call psb_amx(ctxt, tsmth)
if(iam == psb_root_) then
write(psb_out_unit,'(" ")')
write(psb_out_unit,'("Preconditioner: ",a)') trim(p_choice%descr)
write(psb_out_unit,'("Preconditioner time: ",es12.5)')thier+tprec
write(psb_out_unit,'("Preconditioner time: ",es12.5)')thier+tsmth
write(psb_out_unit,'(" ")')
end if
@@ -565,7 +565,7 @@ program amg_d_pde2d
write(psb_out_unit,'("Total preconditioner setup time : ",es12.5)') tsmth+thier
write(psb_out_unit,'("Time to solve system : ",es12.5)') tslv
write(psb_out_unit,'("Time per iteration : ",es12.5)') tslv/iter
write(psb_out_unit,'("Total time : ",es12.5)') tslv+tprec+thier
write(psb_out_unit,'("Total time : ",es12.5)') tslv+tsmth+thier
write(psb_out_unit,'("Residual 2-norm : ",es12.5)') resmx
write(psb_out_unit,'("Residual inf-norm : ",es12.5)') resmxp
write(psb_out_unit,'("Total memory occupation for X : ",i16)') vecsize
samples/advanced/pdegen/amg_d_pde3d.F90
+4 -4
View File
@@ -88,7 +88,7 @@ program amg_d_pde3d
integer(psb_epk_) :: system_size
! miscellaneous
real(psb_dpk_) :: t1, t2, tprec, thier, tslv, tsmth, tpgen
real(psb_dpk_) :: t1, t2, thier, tslv, tsmth, tpgen
! sparse matrix and preconditioner
type(psb_dspmat_type) :: a
@@ -480,12 +480,12 @@ program amg_d_pde3d
end if
call psb_amx(ctxt, thier)
call psb_amx(ctxt, tprec)
call psb_amx(ctxt, tsmth)
if(iam == psb_root_) then
write(psb_out_unit,'(" ")')
write(psb_out_unit,'("Preconditioner: ",a)') trim(p_choice%descr)
write(psb_out_unit,'("Preconditioner time: ",es12.5)')thier+tprec
write(psb_out_unit,'("Preconditioner time: ",es12.5)')thier+tsmth
write(psb_out_unit,'(" ")')
end if
@@ -569,7 +569,7 @@ program amg_d_pde3d
write(psb_out_unit,'("Total preconditioner setup time : ",es12.5)') tsmth+thier
write(psb_out_unit,'("Time to solve system : ",es12.5)') tslv
write(psb_out_unit,'("Time per iteration : ",es12.5)') tslv/iter
write(psb_out_unit,'("Total time : ",es12.5)') tslv+tprec+thier
write(psb_out_unit,'("Total time : ",es12.5)') tslv+tsmth+thier
write(psb_out_unit,'("Residual 2-norm : ",es12.5)') resmx
write(psb_out_unit,'("Residual inf-norm : ",es12.5)') resmxp
write(psb_out_unit,'("Total memory occupation for X : ",i16)') vecsize
samples/advanced/pdegen/amg_s_pde2d.F90
+4 -4
View File
@@ -87,7 +87,7 @@ program amg_s_pde2d
integer(psb_epk_) :: system_size
! miscellaneous
real(psb_dpk_) :: t1, t2, tprec, thier, tslv, tsmth, tpgen
real(psb_dpk_) :: t1, t2, thier, tslv, tsmth, tpgen
! sparse matrix and preconditioner
type(psb_sspmat_type) :: a
@@ -476,12 +476,12 @@ program amg_s_pde2d
end if
call psb_amx(ctxt, thier)
call psb_amx(ctxt, tprec)
call psb_amx(ctxt, tsmth)
if(iam == psb_root_) then
write(psb_out_unit,'(" ")')
write(psb_out_unit,'("Preconditioner: ",a)') trim(p_choice%descr)
write(psb_out_unit,'("Preconditioner time: ",es12.5)')thier+tprec
write(psb_out_unit,'("Preconditioner time: ",es12.5)')thier+tsmth
write(psb_out_unit,'(" ")')
end if
@@ -565,7 +565,7 @@ program amg_s_pde2d
write(psb_out_unit,'("Total preconditioner setup time : ",es12.5)') tsmth+thier
write(psb_out_unit,'("Time to solve system : ",es12.5)') tslv
write(psb_out_unit,'("Time per iteration : ",es12.5)') tslv/iter
write(psb_out_unit,'("Total time : ",es12.5)') tslv+tprec+thier
write(psb_out_unit,'("Total time : ",es12.5)') tslv+tsmth+thier
write(psb_out_unit,'("Residual 2-norm : ",es12.5)') resmx
write(psb_out_unit,'("Residual inf-norm : ",es12.5)') resmxp
write(psb_out_unit,'("Total memory occupation for X : ",i16)') vecsize
samples/advanced/pdegen/amg_s_pde3d.F90
+4 -4
View File
@@ -88,7 +88,7 @@ program amg_s_pde3d
integer(psb_epk_) :: system_size
! miscellaneous
real(psb_dpk_) :: t1, t2, tprec, thier, tslv, tsmth, tpgen
real(psb_dpk_) :: t1, t2, thier, tslv, tsmth, tpgen
! sparse matrix and preconditioner
type(psb_sspmat_type) :: a
@@ -480,12 +480,12 @@ program amg_s_pde3d
end if
call psb_amx(ctxt, thier)
call psb_amx(ctxt, tprec)
call psb_amx(ctxt, tsmth)
if(iam == psb_root_) then
write(psb_out_unit,'(" ")')
write(psb_out_unit,'("Preconditioner: ",a)') trim(p_choice%descr)
write(psb_out_unit,'("Preconditioner time: ",es12.5)')thier+tprec
write(psb_out_unit,'("Preconditioner time: ",es12.5)')thier+tsmth
write(psb_out_unit,'(" ")')
end if
@@ -569,7 +569,7 @@ program amg_s_pde3d
write(psb_out_unit,'("Total preconditioner setup time : ",es12.5)') tsmth+thier
write(psb_out_unit,'("Time to solve system : ",es12.5)') tslv
write(psb_out_unit,'("Time per iteration : ",es12.5)') tslv/iter
write(psb_out_unit,'("Total time : ",es12.5)') tslv+tprec+thier
write(psb_out_unit,'("Total time : ",es12.5)') tslv+tsmth+thier
write(psb_out_unit,'("Residual 2-norm : ",es12.5)') resmx
write(psb_out_unit,'("Residual inf-norm : ",es12.5)') resmxp
write(psb_out_unit,'("Total memory occupation for X : ",i16)') vecsize