diff --git a/samples/advanced/pdegen/amg_d_pde2d.F90 b/samples/advanced/pdegen/amg_d_pde2d.F90 index 0451cd59..877b5cbf 100644 --- a/samples/advanced/pdegen/amg_d_pde2d.F90 +++ b/samples/advanced/pdegen/amg_d_pde2d.F90 @@ -102,7 +102,7 @@ program amg_d_pde2d ! solver parameters integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, nlv - integer(psb_epk_) :: amatsize, precsize, descsize + integer(psb_epk_) :: amatsize, precsize, descsize, vecsize real(psb_dpk_) :: err, resmx, resmxp ! Solver data @@ -537,10 +537,12 @@ program amg_d_pde2d resmx = psb_genrm2(r,desc_a,info) resmxp = psb_geamax(r,desc_a,info) + vecsize = x%sizeof() amatsize = a%sizeof() descsize = desc_a%sizeof() precsize = prec%sizeof() system_size = desc_a%get_global_rows() + call psb_sum(ctxt,vecsize) call psb_sum(ctxt,amatsize) call psb_sum(ctxt,descsize) call psb_sum(ctxt,precsize) @@ -564,9 +566,12 @@ program amg_d_pde2d write(psb_out_unit,'("Total time : ",es12.5)') tslv+tprec+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 : ",i12)') vecsize write(psb_out_unit,'("Total memory occupation for A : ",i12)') amatsize write(psb_out_unit,'("Total memory occupation for DESC_A : ",i12)') descsize write(psb_out_unit,'("Total memory occupation for PREC : ",i12)') precsize + write(psb_out_unit,'("Total memory occupation : ",i12)') & + & amatsize + descsize+precsize+2*vecsize write(psb_out_unit,'("Storage format for A : ",a )') a%get_fmt() write(psb_out_unit,'("Storage format for DESC_A : ",a )') desc_a%get_fmt() diff --git a/samples/advanced/pdegen/amg_d_pde3d.F90 b/samples/advanced/pdegen/amg_d_pde3d.F90 index 33860515..a59d2ed3 100644 --- a/samples/advanced/pdegen/amg_d_pde3d.F90 +++ b/samples/advanced/pdegen/amg_d_pde3d.F90 @@ -103,7 +103,7 @@ program amg_d_pde3d ! solver parameters integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, nlv - integer(psb_epk_) :: amatsize, precsize, descsize + integer(psb_epk_) :: amatsize, precsize, descsize, vecsize real(psb_dpk_) :: err, resmx, resmxp ! Solver data @@ -541,10 +541,12 @@ program amg_d_pde3d resmx = psb_genrm2(r,desc_a,info) resmxp = psb_geamax(r,desc_a,info) + vecsize = x%sizeof() amatsize = a%sizeof() descsize = desc_a%sizeof() precsize = prec%sizeof() system_size = desc_a%get_global_rows() + call psb_sum(ctxt,vecsize) call psb_sum(ctxt,amatsize) call psb_sum(ctxt,descsize) call psb_sum(ctxt,precsize) @@ -568,9 +570,12 @@ program amg_d_pde3d write(psb_out_unit,'("Total time : ",es12.5)') tslv+tprec+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 : ",i12)') vecsize write(psb_out_unit,'("Total memory occupation for A : ",i12)') amatsize write(psb_out_unit,'("Total memory occupation for DESC_A : ",i12)') descsize write(psb_out_unit,'("Total memory occupation for PREC : ",i12)') precsize + write(psb_out_unit,'("Total memory occupation : ",i12)') & + & amatsize + descsize+precsize+2*vecsize write(psb_out_unit,'("Storage format for A : ",a )') a%get_fmt() write(psb_out_unit,'("Storage format for DESC_A : ",a )') desc_a%get_fmt() diff --git a/samples/advanced/pdegen/amg_s_pde2d.F90 b/samples/advanced/pdegen/amg_s_pde2d.F90 index 9ce86cc5..2241637d 100644 --- a/samples/advanced/pdegen/amg_s_pde2d.F90 +++ b/samples/advanced/pdegen/amg_s_pde2d.F90 @@ -102,7 +102,7 @@ program amg_s_pde2d ! solver parameters integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, nlv - integer(psb_epk_) :: amatsize, precsize, descsize + integer(psb_epk_) :: amatsize, precsize, descsize, vecsize real(psb_spk_) :: err, resmx, resmxp ! Solver data @@ -537,10 +537,12 @@ program amg_s_pde2d resmx = psb_genrm2(r,desc_a,info) resmxp = psb_geamax(r,desc_a,info) + vecsize = x%sizeof() amatsize = a%sizeof() descsize = desc_a%sizeof() precsize = prec%sizeof() system_size = desc_a%get_global_rows() + call psb_sum(ctxt,vecsize) call psb_sum(ctxt,amatsize) call psb_sum(ctxt,descsize) call psb_sum(ctxt,precsize) @@ -564,9 +566,12 @@ program amg_s_pde2d write(psb_out_unit,'("Total time : ",es12.5)') tslv+tprec+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 : ",i12)') vecsize write(psb_out_unit,'("Total memory occupation for A : ",i12)') amatsize write(psb_out_unit,'("Total memory occupation for DESC_A : ",i12)') descsize write(psb_out_unit,'("Total memory occupation for PREC : ",i12)') precsize + write(psb_out_unit,'("Total memory occupation : ",i12)') & + & amatsize + descsize+precsize+2*vecsize write(psb_out_unit,'("Storage format for A : ",a )') a%get_fmt() write(psb_out_unit,'("Storage format for DESC_A : ",a )') desc_a%get_fmt() diff --git a/samples/advanced/pdegen/amg_s_pde3d.F90 b/samples/advanced/pdegen/amg_s_pde3d.F90 index 4c140cca..234e91e7 100644 --- a/samples/advanced/pdegen/amg_s_pde3d.F90 +++ b/samples/advanced/pdegen/amg_s_pde3d.F90 @@ -103,7 +103,7 @@ program amg_s_pde3d ! solver parameters integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, nlv - integer(psb_epk_) :: amatsize, precsize, descsize + integer(psb_epk_) :: amatsize, precsize, descsize, vecsize real(psb_spk_) :: err, resmx, resmxp ! Solver data @@ -541,10 +541,12 @@ program amg_s_pde3d resmx = psb_genrm2(r,desc_a,info) resmxp = psb_geamax(r,desc_a,info) + vecsize = x%sizeof() amatsize = a%sizeof() descsize = desc_a%sizeof() precsize = prec%sizeof() system_size = desc_a%get_global_rows() + call psb_sum(ctxt,vecsize) call psb_sum(ctxt,amatsize) call psb_sum(ctxt,descsize) call psb_sum(ctxt,precsize) @@ -568,9 +570,12 @@ program amg_s_pde3d write(psb_out_unit,'("Total time : ",es12.5)') tslv+tprec+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 : ",i12)') vecsize write(psb_out_unit,'("Total memory occupation for A : ",i12)') amatsize write(psb_out_unit,'("Total memory occupation for DESC_A : ",i12)') descsize write(psb_out_unit,'("Total memory occupation for PREC : ",i12)') precsize + write(psb_out_unit,'("Total memory occupation : ",i12)') & + & amatsize + descsize+precsize+2*vecsize write(psb_out_unit,'("Storage format for A : ",a )') a%get_fmt() write(psb_out_unit,'("Storage format for DESC_A : ",a )') desc_a%get_fmt()