First round of fixes for precdescr

tspmm
Salvatore Filippone 3 years ago
parent 2f45f8631b
commit 485a94765b

@ -735,20 +735,40 @@ contains
& matrix_names(pm%coarse_mat) & matrix_names(pm%coarse_mat)
select case(pm%coarse_solve) select case(pm%coarse_solve)
case (amg_bjac_,amg_as_) case (amg_bjac_,amg_as_)
write(iout,*) ' Number of sweeps : ',&
& pm%sweeps_pre
write(iout,*) ' Coarse solver: ',& write(iout,*) ' Coarse solver: ',&
& 'Block Jacobi' & 'Block Jacobi'
case (amg_l1_bjac_)
write(iout,*) ' Number of sweeps : ',& write(iout,*) ' Number of sweeps : ',&
& pm%sweeps_pre & pm%sweeps_pre
case (amg_l1_bjac_)
write(iout,*) ' Coarse solver: ',& write(iout,*) ' Coarse solver: ',&
& 'L1-Block Jacobi' & 'L1-Block Jacobi'
case (amg_jac_)
write(iout,*) ' Number of sweeps : ',& write(iout,*) ' Number of sweeps : ',&
& pm%sweeps_pre & pm%sweeps_pre
case (amg_jac_)
write(iout,*) ' Coarse solver: ',& write(iout,*) ' Coarse solver: ',&
& 'Point Jacobi' & 'Point Jacobi'
write(iout,*) ' Number of sweeps : ',&
& pm%sweeps_pre
case (amg_l1_jac_)
write(iout,*) ' Coarse solver: ',&
& 'L1-Jacobi'
write(iout,*) ' Number of sweeps : ',&
& pm%sweeps_pre
case (amg_l1_fbgs_)
write(iout,*) ' Coarse solver: ',&
& 'L1 Forward-Backward Gauss-Seidel (Hybrid)'
write(iout,*) ' Number of sweeps : ',&
& pm%sweeps_pre
case (amg_l1_gs_)
write(iout,*) ' Coarse solver: ',&
& 'L1 Gauss-Seidel (Hybrid)'
write(iout,*) ' Number of sweeps : ',&
& pm%sweeps_pre
case (amg_fbgs_)
write(iout,*) ' Coarse solver: ',&
& 'Forward-Backward Gauss-Seidel (Hybrid)'
write(iout,*) ' Number of sweeps : ',&
& pm%sweeps_pre
case default case default
write(iout,*) ' Coarse solver: ',& write(iout,*) ' Coarse solver: ',&
& amg_fact_names(pm%coarse_solve) & amg_fact_names(pm%coarse_solve)

@ -466,15 +466,9 @@ contains
end if end if
write(iout_,*) ' method: ',sv%method write(iout_,*) ' method: ',sv%method
write(iout_,*) ' kprec: ',sv%kprec write(iout_,*) ' kprec: ',sv%kprec
if (sv%i_sub_solve > 0) then call sv%prec%descr(iout_,info)
write(iout_,*) ' sub_solve: ',amg_fact_names(sv%i_sub_solve)
else
write(iout_,*) ' sub_solve: ',sv%sub_solve
end if
write(iout_,*) ' itmax: ',sv%itmax write(iout_,*) ' itmax: ',sv%itmax
write(iout_,*) ' eps: ',sv%eps write(iout_,*) ' eps: ',sv%eps
write(iout_,*) ' fillin: ',sv%fillin
call psb_erractionrestore(err_act) call psb_erractionrestore(err_act)
return return

@ -466,15 +466,9 @@ contains
end if end if
write(iout_,*) ' method: ',sv%method write(iout_,*) ' method: ',sv%method
write(iout_,*) ' kprec: ',sv%kprec write(iout_,*) ' kprec: ',sv%kprec
if (sv%i_sub_solve > 0) then call sv%prec%descr(iout_,info)
write(iout_,*) ' sub_solve: ',amg_fact_names(sv%i_sub_solve)
else
write(iout_,*) ' sub_solve: ',sv%sub_solve
end if
write(iout_,*) ' itmax: ',sv%itmax write(iout_,*) ' itmax: ',sv%itmax
write(iout_,*) ' eps: ',sv%eps write(iout_,*) ' eps: ',sv%eps
write(iout_,*) ' fillin: ',sv%fillin
call psb_erractionrestore(err_act) call psb_erractionrestore(err_act)
return return

@ -466,15 +466,9 @@ contains
end if end if
write(iout_,*) ' method: ',sv%method write(iout_,*) ' method: ',sv%method
write(iout_,*) ' kprec: ',sv%kprec write(iout_,*) ' kprec: ',sv%kprec
if (sv%i_sub_solve > 0) then call sv%prec%descr(iout_,info)
write(iout_,*) ' sub_solve: ',amg_fact_names(sv%i_sub_solve)
else
write(iout_,*) ' sub_solve: ',sv%sub_solve
end if
write(iout_,*) ' itmax: ',sv%itmax write(iout_,*) ' itmax: ',sv%itmax
write(iout_,*) ' eps: ',sv%eps write(iout_,*) ' eps: ',sv%eps
write(iout_,*) ' fillin: ',sv%fillin
call psb_erractionrestore(err_act) call psb_erractionrestore(err_act)
return return

@ -466,15 +466,9 @@ contains
end if end if
write(iout_,*) ' method: ',sv%method write(iout_,*) ' method: ',sv%method
write(iout_,*) ' kprec: ',sv%kprec write(iout_,*) ' kprec: ',sv%kprec
if (sv%i_sub_solve > 0) then call sv%prec%descr(iout_,info)
write(iout_,*) ' sub_solve: ',amg_fact_names(sv%i_sub_solve)
else
write(iout_,*) ' sub_solve: ',sv%sub_solve
end if
write(iout_,*) ' itmax: ',sv%itmax write(iout_,*) ' itmax: ',sv%itmax
write(iout_,*) ' eps: ',sv%eps write(iout_,*) ' eps: ',sv%eps
write(iout_,*) ' fillin: ',sv%fillin
call psb_erractionrestore(err_act) call psb_erractionrestore(err_act)
return return

Loading…
Cancel
Save