mirror of
https://github.com/sfilippone/amg4psblas.git
synced 2026-10-06 22:55:12 +00:00
Modified IPK/LPK/EPK/MPK.
This commit is contained in:
Salvatore Filippone
@@ -81,7 +81,7 @@ program mld_cexample_1lev
|
||||
|
||||
! other variables
|
||||
integer :: i,info,j,m_problem
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: ierr, ircode
|
||||
real(psb_spk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -96,7 +96,7 @@ program mld_cexample_ml
|
||||
! other variables
|
||||
integer :: choice
|
||||
integer :: i,info,j,m_problem
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: ierr, ircode
|
||||
real(psb_spk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -81,7 +81,7 @@ program mld_dexample_1lev
|
||||
|
||||
! other variables
|
||||
integer :: i,info,j,m_problem
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: ierr, ircode
|
||||
real(psb_dpk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -96,7 +96,7 @@ program mld_dexample_ml
|
||||
! other variables
|
||||
integer :: choice
|
||||
integer :: i,info,j,m_problem
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: ierr, ircode
|
||||
real(psb_dpk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -81,7 +81,7 @@ program mld_sexample_1lev
|
||||
|
||||
! other variables
|
||||
integer :: i,info,j,m_problem
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: ierr, ircode
|
||||
real(psb_spk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -96,7 +96,7 @@ program mld_sexample_ml
|
||||
! other variables
|
||||
integer :: choice
|
||||
integer :: i,info,j,m_problem
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: ierr, ircode
|
||||
real(psb_spk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -81,7 +81,7 @@ program mld_zexample_1lev
|
||||
|
||||
! other variables
|
||||
integer :: i,info,j,m_problem
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: ierr, ircode
|
||||
real(psb_dpk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -96,7 +96,7 @@ program mld_zexample_ml
|
||||
! other variables
|
||||
integer :: choice
|
||||
integer :: i,info,j,m_problem
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: ierr, ircode
|
||||
real(psb_dpk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -88,7 +88,7 @@ program mld_dexample_1lev
|
||||
|
||||
! other variables
|
||||
integer :: i,info,j
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: idim, nlev, ierr, ircode
|
||||
real(psb_dpk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -108,7 +108,7 @@ program mld_dexample_ml
|
||||
! other variables
|
||||
integer :: choice
|
||||
integer :: i,info,j
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: idim, ierr, ircode
|
||||
real(psb_dpk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -88,7 +88,7 @@ program mld_sexample_1lev
|
||||
|
||||
! other variables
|
||||
integer :: i,info,j
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: idim, nlev, ierr, ircode
|
||||
real(psb_spk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -108,7 +108,7 @@ program mld_sexample_ml
|
||||
! other variables
|
||||
integer :: choice
|
||||
integer :: i,info,j
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
integer :: idim, ierr, ircode
|
||||
real(psb_spk_) :: resmx, resmxp
|
||||
real(psb_dpk_) :: t1, t2, tprec
|
||||
|
||||
@@ -346,7 +346,7 @@ contains
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me, ncol
|
||||
integer(psb_mpk_) :: ictxt, np, me, ncol
|
||||
integer(psb_ipk_) :: err_act,ntaggr,nzl
|
||||
integer(psb_ipk_), allocatable :: ilaggr(:), nlaggr(:)
|
||||
type(psb_cspmat_type) :: ac, am2, am3, am4
|
||||
|
||||
@@ -92,7 +92,7 @@ subroutine mld_c_lev_aggrmap_bld(p,a,desc_a,ilaggr,nlaggr,op_prol,info)
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me
|
||||
integer(psb_mpk_) :: ictxt, np, me
|
||||
integer(psb_ipk_) :: err_act
|
||||
integer(psb_ipk_) :: nzl, ntaggr
|
||||
integer(psb_ipk_) :: debug_level, debug_unit
|
||||
|
||||
@@ -102,7 +102,7 @@ subroutine mld_c_lev_aggrmat_asb(p,a,desc_a,ilaggr,nlaggr,op_prol,info)
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me
|
||||
integer(psb_mpk_) :: ictxt, np, me
|
||||
integer(psb_ipk_) :: err_act
|
||||
type(psb_cspmat_type) :: ac, op_restr
|
||||
type(psb_c_coo_sparse_mat) :: acoo, bcoo
|
||||
|
||||
@@ -346,7 +346,7 @@ contains
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me, ncol
|
||||
integer(psb_mpk_) :: ictxt, np, me, ncol
|
||||
integer(psb_ipk_) :: err_act,ntaggr,nzl
|
||||
integer(psb_ipk_), allocatable :: ilaggr(:), nlaggr(:)
|
||||
type(psb_dspmat_type) :: ac, am2, am3, am4
|
||||
|
||||
@@ -92,7 +92,7 @@ subroutine mld_d_lev_aggrmap_bld(p,a,desc_a,ilaggr,nlaggr,op_prol,info)
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me
|
||||
integer(psb_mpk_) :: ictxt, np, me
|
||||
integer(psb_ipk_) :: err_act
|
||||
integer(psb_ipk_) :: nzl, ntaggr
|
||||
integer(psb_ipk_) :: debug_level, debug_unit
|
||||
|
||||
@@ -102,7 +102,7 @@ subroutine mld_d_lev_aggrmat_asb(p,a,desc_a,ilaggr,nlaggr,op_prol,info)
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me
|
||||
integer(psb_mpk_) :: ictxt, np, me
|
||||
integer(psb_ipk_) :: err_act
|
||||
type(psb_dspmat_type) :: ac, op_restr
|
||||
type(psb_d_coo_sparse_mat) :: acoo, bcoo
|
||||
|
||||
@@ -346,7 +346,7 @@ contains
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me, ncol
|
||||
integer(psb_mpk_) :: ictxt, np, me, ncol
|
||||
integer(psb_ipk_) :: err_act,ntaggr,nzl
|
||||
integer(psb_ipk_), allocatable :: ilaggr(:), nlaggr(:)
|
||||
type(psb_sspmat_type) :: ac, am2, am3, am4
|
||||
|
||||
@@ -92,7 +92,7 @@ subroutine mld_s_lev_aggrmap_bld(p,a,desc_a,ilaggr,nlaggr,op_prol,info)
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me
|
||||
integer(psb_mpk_) :: ictxt, np, me
|
||||
integer(psb_ipk_) :: err_act
|
||||
integer(psb_ipk_) :: nzl, ntaggr
|
||||
integer(psb_ipk_) :: debug_level, debug_unit
|
||||
|
||||
@@ -102,7 +102,7 @@ subroutine mld_s_lev_aggrmat_asb(p,a,desc_a,ilaggr,nlaggr,op_prol,info)
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me
|
||||
integer(psb_mpk_) :: ictxt, np, me
|
||||
integer(psb_ipk_) :: err_act
|
||||
type(psb_sspmat_type) :: ac, op_restr
|
||||
type(psb_s_coo_sparse_mat) :: acoo, bcoo
|
||||
|
||||
@@ -346,7 +346,7 @@ contains
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me, ncol
|
||||
integer(psb_mpk_) :: ictxt, np, me, ncol
|
||||
integer(psb_ipk_) :: err_act,ntaggr,nzl
|
||||
integer(psb_ipk_), allocatable :: ilaggr(:), nlaggr(:)
|
||||
type(psb_zspmat_type) :: ac, am2, am3, am4
|
||||
|
||||
@@ -92,7 +92,7 @@ subroutine mld_z_lev_aggrmap_bld(p,a,desc_a,ilaggr,nlaggr,op_prol,info)
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me
|
||||
integer(psb_mpk_) :: ictxt, np, me
|
||||
integer(psb_ipk_) :: err_act
|
||||
integer(psb_ipk_) :: nzl, ntaggr
|
||||
integer(psb_ipk_) :: debug_level, debug_unit
|
||||
|
||||
@@ -102,7 +102,7 @@ subroutine mld_z_lev_aggrmat_asb(p,a,desc_a,ilaggr,nlaggr,op_prol,info)
|
||||
|
||||
! Local variables
|
||||
character(len=20) :: name
|
||||
integer(psb_mpik_) :: ictxt, np, me
|
||||
integer(psb_mpk_) :: ictxt, np, me
|
||||
integer(psb_ipk_) :: err_act
|
||||
type(psb_zspmat_type) :: ac, op_restr
|
||||
type(psb_z_coo_sparse_mat) :: acoo, bcoo
|
||||
|
||||
@@ -65,10 +65,10 @@ module mld_base_prec_type
|
||||
use psb_const_mod
|
||||
use psb_base_mod, only :&
|
||||
& psb_desc_type, psb_i_vect_type, psb_i_base_vect_type,&
|
||||
& psb_ipk_, psb_dpk_, psb_spk_, psb_long_int_k_, &
|
||||
& psb_ipk_, psb_dpk_, psb_spk_, psb_epk_, &
|
||||
& psb_cdfree, psb_halo_, psb_none_, psb_sum_, psb_avg_, &
|
||||
& psb_nohalo_, psb_square_root_, psb_toupper, psb_root_,&
|
||||
& psb_sizeof_int, psb_sizeof_long_int, psb_sizeof_sp, &
|
||||
& psb_sizeof_ip, psb_sizeof_lp, psb_sizeof_sp, &
|
||||
& psb_sizeof_dp, psb_sizeof,&
|
||||
& psb_cd_get_context, psb_info, psb_min, psb_sum, psb_bcast,&
|
||||
& psb_sizeof, psb_free, psb_cdfree, &
|
||||
|
||||
@@ -111,7 +111,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_check(sm,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -121,7 +121,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_restr_v(sm,x,trans,work,info,data)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
@@ -136,7 +136,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_restr_a(sm,x,trans,work,info,data)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
@@ -151,7 +151,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_prol_v(sm,x,trans,work,info,data)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
@@ -166,7 +166,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_prol_a(sm,x,trans,work,info,data)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
@@ -183,7 +183,7 @@ module mld_c_as_smoother
|
||||
subroutine mld_c_as_smoother_apply_vect(alpha,sm,x,beta,y,desc_data,&
|
||||
& trans,sweeps,work,wv,info,init,initu)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
type(psb_desc_type), intent(in) :: desc_data
|
||||
@@ -205,7 +205,7 @@ module mld_c_as_smoother
|
||||
subroutine mld_c_as_smoother_apply(alpha,sm,x,beta,y,desc_data,&
|
||||
& trans,sweeps,work,info,init,initu)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_,&
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_,&
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
type(psb_desc_type), intent(in) :: desc_data
|
||||
@@ -225,7 +225,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_c_base_sparse_mat, psb_ipk_,&
|
||||
& psb_i_base_vect_type
|
||||
implicit none
|
||||
@@ -242,7 +242,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_cnv(sm,info,amold,vmold,imold)
|
||||
import :: psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, &
|
||||
& psb_c_base_sparse_mat, psb_ipk_, psb_i_base_vect_type
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
@@ -256,7 +256,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_seti(sm,what,val,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -268,7 +268,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_setc(sm,what,val,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -280,7 +280,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_setr(sm,what,val,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -292,7 +292,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_cseti(sm,what,val,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -304,7 +304,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_csetc(sm,what,val,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -316,7 +316,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_csetr(sm,what,val,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -328,7 +328,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_free(sm,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -338,7 +338,7 @@ module mld_c_as_smoother
|
||||
interface
|
||||
subroutine mld_c_as_smoother_dmp(sm,ictxt,level,info,prefix,head,smoother,solver)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_long_int_k_, psb_desc_type, &
|
||||
& psb_spk_, mld_c_as_smoother_type, psb_epk_, psb_desc_type, &
|
||||
& psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(in) :: sm
|
||||
@@ -366,10 +366,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_as_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
if (allocated(sm%sv)) val = val + sm%sv%sizeof()
|
||||
val = val + sm%nd%sizeof()
|
||||
|
||||
@@ -379,7 +379,7 @@ contains
|
||||
function c_as_smoother_get_nzeros(sm) result(val)
|
||||
implicit none
|
||||
class(mld_c_as_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(sm%sv)) &
|
||||
|
||||
@@ -60,7 +60,7 @@
|
||||
module mld_c_base_smoother_mod
|
||||
|
||||
use mld_c_base_solver_mod
|
||||
use psb_base_mod, only : psb_desc_type, psb_cspmat_type, psb_long_int_k_,&
|
||||
use psb_base_mod, only : psb_desc_type, psb_cspmat_type, psb_epk_,&
|
||||
& psb_c_vect_type, psb_c_base_vect_type, psb_c_base_sparse_mat, &
|
||||
& psb_spk_, psb_i_base_vect_type, psb_erractionsave, psb_error_handler
|
||||
|
||||
@@ -355,7 +355,7 @@ contains
|
||||
function c_base_smoother_get_nzeros(sm) result(val)
|
||||
implicit none
|
||||
class(mld_c_base_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(sm%sv)) &
|
||||
@@ -366,7 +366,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_base_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -368,7 +368,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_base_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
|
||||
@@ -378,7 +378,7 @@ contains
|
||||
function c_base_solver_get_nzeros(sv) result(val)
|
||||
implicit none
|
||||
class(mld_c_base_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
end function c_base_solver_get_nzeros
|
||||
|
||||
@@ -237,7 +237,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_diag_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -250,7 +250,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_diag_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -543,7 +543,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_gs_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -558,10 +558,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_gs_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
val = val + sv%l%sizeof()
|
||||
val = val + sv%u%sizeof()
|
||||
|
||||
|
||||
@@ -517,7 +517,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_ilu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -533,10 +533,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_ilu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 2*psb_sizeof_int + (2*psb_sizeof_sp)
|
||||
val = 2*psb_sizeof_ip + (2*psb_sizeof_sp)
|
||||
val = val + sv%dv%sizeof()
|
||||
val = val + sv%l%sizeof()
|
||||
val = val + sv%u%sizeof()
|
||||
|
||||
@@ -156,7 +156,7 @@ module mld_c_jac_smoother
|
||||
interface
|
||||
subroutine mld_c_jac_smoother_dmp(sm,ictxt,level,info,prefix,head,smoother,solver)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_spk_, mld_c_jac_smoother_type, psb_long_int_k_, psb_desc_type, &
|
||||
& psb_spk_, mld_c_jac_smoother_type, psb_epk_, psb_desc_type, &
|
||||
& psb_ipk_
|
||||
implicit none
|
||||
class(mld_c_jac_smoother_type), intent(in) :: sm
|
||||
@@ -230,10 +230,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_jac_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
if (allocated(sm%sv)) val = val + sm%sv%sizeof()
|
||||
val = val + sm%nd%sizeof()
|
||||
|
||||
@@ -245,7 +245,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_jac_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -468,14 +468,14 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_mumps_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
#if defined(HAVE_MUMPS_)
|
||||
val = (sv%id%INFOG(22)+sv%id%INFOG(32))*1d+6
|
||||
#else
|
||||
val = 0
|
||||
#endif
|
||||
! val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
! val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
! val = val + sv%symbsize
|
||||
! val = val + sv%numsize
|
||||
return
|
||||
|
||||
+16
-16
@@ -57,7 +57,7 @@ module mld_c_onelev_mod
|
||||
use mld_c_base_smoother_mod
|
||||
use psb_base_mod, only : psb_cspmat_type, psb_c_vect_type, &
|
||||
& psb_c_base_vect_type, psb_clinmap_type, psb_spk_, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type, psb_i_base_vect_type, &
|
||||
& psb_ipk_, psb_epk_, psb_desc_type, psb_i_base_vect_type, &
|
||||
& psb_erractionsave, psb_error_handler
|
||||
!
|
||||
!
|
||||
@@ -194,7 +194,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_build(lv,info,amold,vmold,imold)
|
||||
import :: psb_c_base_sparse_mat, psb_c_base_vect_type, &
|
||||
& psb_i_base_vect_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
class(mld_c_onelev_type), target, intent(inout) :: lv
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -208,7 +208,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_descr(lv,il,nl,ilmin,info,iout)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_c_onelev_type), intent(in) :: lv
|
||||
@@ -235,7 +235,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_free(lv,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
|
||||
class(mld_c_onelev_type), intent(inout) :: lv
|
||||
@@ -247,7 +247,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_check(lv,info)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_c_onelev_type), intent(inout) :: lv
|
||||
@@ -259,7 +259,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_seti(lv,what,val,info,pos)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -274,7 +274,7 @@ module mld_c_onelev_mod
|
||||
interface
|
||||
subroutine mld_c_base_onelev_setsm(lv,val,info,pos)
|
||||
import :: psb_spk_, mld_c_onelev_type, mld_c_base_smoother_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -288,7 +288,7 @@ module mld_c_onelev_mod
|
||||
interface
|
||||
subroutine mld_c_base_onelev_setsv(lv,val,info,pos)
|
||||
import :: psb_spk_, mld_c_onelev_type, mld_c_base_solver_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -303,7 +303,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_setc(lv,what,val,info,pos)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_c_onelev_type), intent(inout) :: lv
|
||||
@@ -318,7 +318,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_setr(lv,what,val,info,pos)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
class(mld_c_onelev_type), intent(inout) :: lv
|
||||
@@ -334,7 +334,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_cseti(lv,what,val,info,pos)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -350,7 +350,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_csetc(lv,what,val,info,pos)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_c_onelev_type), intent(inout) :: lv
|
||||
@@ -365,7 +365,7 @@ module mld_c_onelev_mod
|
||||
subroutine mld_c_base_onelev_csetr(lv,what,val,info,pos)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
class(mld_c_onelev_type), intent(inout) :: lv
|
||||
@@ -381,7 +381,7 @@ module mld_c_onelev_mod
|
||||
& solver,global_num)
|
||||
import :: psb_cspmat_type, psb_c_vect_type, psb_c_base_vect_type, &
|
||||
& psb_clinmap_type, psb_spk_, mld_c_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
class(mld_c_onelev_type), intent(in) :: lv
|
||||
integer(psb_ipk_), intent(in) :: level
|
||||
@@ -400,7 +400,7 @@ contains
|
||||
function c_base_onelev_get_nzeros(lv) result(val)
|
||||
implicit none
|
||||
class(mld_c_onelev_type), intent(in) :: lv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(lv%sm)) &
|
||||
@@ -412,7 +412,7 @@ contains
|
||||
function c_base_onelev_sizeof(lv) result(val)
|
||||
implicit none
|
||||
class(mld_c_onelev_type), intent(in) :: lv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -425,7 +425,7 @@ contains
|
||||
function mld_c_get_nzeros(prec) result(val)
|
||||
implicit none
|
||||
class(mld_cprec_type), intent(in) :: prec
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(prec%precv)) then
|
||||
@@ -438,10 +438,10 @@ contains
|
||||
function mld_cprec_sizeof(prec) result(val)
|
||||
implicit none
|
||||
class(mld_cprec_type), intent(in) :: prec
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
val = val + psb_sizeof_int
|
||||
val = val + psb_sizeof_ip
|
||||
if (allocated(prec%precv)) then
|
||||
do i=1, size(prec%precv)
|
||||
val = val + prec%precv(i)%sizeof()
|
||||
|
||||
@@ -402,10 +402,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_c_slu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
|
||||
val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
val = val + sv%symbsize
|
||||
val = val + sv%numsize
|
||||
return
|
||||
|
||||
@@ -111,7 +111,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_check(sm,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -121,7 +121,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_restr_v(sm,x,trans,work,info,data)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
@@ -136,7 +136,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_restr_a(sm,x,trans,work,info,data)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
@@ -151,7 +151,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_prol_v(sm,x,trans,work,info,data)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
@@ -166,7 +166,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_prol_a(sm,x,trans,work,info,data)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
@@ -183,7 +183,7 @@ module mld_d_as_smoother
|
||||
subroutine mld_d_as_smoother_apply_vect(alpha,sm,x,beta,y,desc_data,&
|
||||
& trans,sweeps,work,wv,info,init,initu)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
type(psb_desc_type), intent(in) :: desc_data
|
||||
@@ -205,7 +205,7 @@ module mld_d_as_smoother
|
||||
subroutine mld_d_as_smoother_apply(alpha,sm,x,beta,y,desc_data,&
|
||||
& trans,sweeps,work,info,init,initu)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_,&
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_,&
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
type(psb_desc_type), intent(in) :: desc_data
|
||||
@@ -225,7 +225,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_d_base_sparse_mat, psb_ipk_,&
|
||||
& psb_i_base_vect_type
|
||||
implicit none
|
||||
@@ -242,7 +242,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_cnv(sm,info,amold,vmold,imold)
|
||||
import :: psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, &
|
||||
& psb_d_base_sparse_mat, psb_ipk_, psb_i_base_vect_type
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
@@ -256,7 +256,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_seti(sm,what,val,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -268,7 +268,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_setc(sm,what,val,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -280,7 +280,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_setr(sm,what,val,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -292,7 +292,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_cseti(sm,what,val,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -304,7 +304,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_csetc(sm,what,val,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -316,7 +316,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_csetr(sm,what,val,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -328,7 +328,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_free(sm,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -338,7 +338,7 @@ module mld_d_as_smoother
|
||||
interface
|
||||
subroutine mld_d_as_smoother_dmp(sm,ictxt,level,info,prefix,head,smoother,solver)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_long_int_k_, psb_desc_type, &
|
||||
& psb_dpk_, mld_d_as_smoother_type, psb_epk_, psb_desc_type, &
|
||||
& psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(in) :: sm
|
||||
@@ -366,10 +366,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_as_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
if (allocated(sm%sv)) val = val + sm%sv%sizeof()
|
||||
val = val + sm%nd%sizeof()
|
||||
|
||||
@@ -379,7 +379,7 @@ contains
|
||||
function d_as_smoother_get_nzeros(sm) result(val)
|
||||
implicit none
|
||||
class(mld_d_as_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(sm%sv)) &
|
||||
|
||||
@@ -60,7 +60,7 @@
|
||||
module mld_d_base_smoother_mod
|
||||
|
||||
use mld_d_base_solver_mod
|
||||
use psb_base_mod, only : psb_desc_type, psb_dspmat_type, psb_long_int_k_,&
|
||||
use psb_base_mod, only : psb_desc_type, psb_dspmat_type, psb_epk_,&
|
||||
& psb_d_vect_type, psb_d_base_vect_type, psb_d_base_sparse_mat, &
|
||||
& psb_dpk_, psb_i_base_vect_type, psb_erractionsave, psb_error_handler
|
||||
|
||||
@@ -355,7 +355,7 @@ contains
|
||||
function d_base_smoother_get_nzeros(sm) result(val)
|
||||
implicit none
|
||||
class(mld_d_base_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(sm%sv)) &
|
||||
@@ -366,7 +366,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_base_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -368,7 +368,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_base_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
|
||||
@@ -378,7 +378,7 @@ contains
|
||||
function d_base_solver_get_nzeros(sv) result(val)
|
||||
implicit none
|
||||
class(mld_d_base_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
end function d_base_solver_get_nzeros
|
||||
|
||||
@@ -237,7 +237,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_diag_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -250,7 +250,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_diag_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -543,7 +543,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_gs_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -558,10 +558,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_gs_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
val = val + sv%l%sizeof()
|
||||
val = val + sv%u%sizeof()
|
||||
|
||||
|
||||
@@ -517,7 +517,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_ilu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -533,10 +533,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_ilu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
val = val + sv%dv%sizeof()
|
||||
val = val + sv%l%sizeof()
|
||||
val = val + sv%u%sizeof()
|
||||
|
||||
@@ -156,7 +156,7 @@ module mld_d_jac_smoother
|
||||
interface
|
||||
subroutine mld_d_jac_smoother_dmp(sm,ictxt,level,info,prefix,head,smoother,solver)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dpk_, mld_d_jac_smoother_type, psb_long_int_k_, psb_desc_type, &
|
||||
& psb_dpk_, mld_d_jac_smoother_type, psb_epk_, psb_desc_type, &
|
||||
& psb_ipk_
|
||||
implicit none
|
||||
class(mld_d_jac_smoother_type), intent(in) :: sm
|
||||
@@ -230,10 +230,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_jac_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
if (allocated(sm%sv)) val = val + sm%sv%sizeof()
|
||||
val = val + sm%nd%sizeof()
|
||||
|
||||
@@ -245,7 +245,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_jac_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -468,14 +468,14 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_mumps_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
#if defined(HAVE_MUMPS_)
|
||||
val = (sv%id%INFOG(22)+sv%id%INFOG(32))*1d+6
|
||||
#else
|
||||
val = 0
|
||||
#endif
|
||||
! val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
! val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
! val = val + sv%symbsize
|
||||
! val = val + sv%numsize
|
||||
return
|
||||
|
||||
+16
-16
@@ -57,7 +57,7 @@ module mld_d_onelev_mod
|
||||
use mld_d_base_smoother_mod
|
||||
use psb_base_mod, only : psb_dspmat_type, psb_d_vect_type, &
|
||||
& psb_d_base_vect_type, psb_dlinmap_type, psb_dpk_, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type, psb_i_base_vect_type, &
|
||||
& psb_ipk_, psb_epk_, psb_desc_type, psb_i_base_vect_type, &
|
||||
& psb_erractionsave, psb_error_handler
|
||||
!
|
||||
!
|
||||
@@ -194,7 +194,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_build(lv,info,amold,vmold,imold)
|
||||
import :: psb_d_base_sparse_mat, psb_d_base_vect_type, &
|
||||
& psb_i_base_vect_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
class(mld_d_onelev_type), target, intent(inout) :: lv
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -208,7 +208,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_descr(lv,il,nl,ilmin,info,iout)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_d_onelev_type), intent(in) :: lv
|
||||
@@ -235,7 +235,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_free(lv,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
|
||||
class(mld_d_onelev_type), intent(inout) :: lv
|
||||
@@ -247,7 +247,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_check(lv,info)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_d_onelev_type), intent(inout) :: lv
|
||||
@@ -259,7 +259,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_seti(lv,what,val,info,pos)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -274,7 +274,7 @@ module mld_d_onelev_mod
|
||||
interface
|
||||
subroutine mld_d_base_onelev_setsm(lv,val,info,pos)
|
||||
import :: psb_dpk_, mld_d_onelev_type, mld_d_base_smoother_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -288,7 +288,7 @@ module mld_d_onelev_mod
|
||||
interface
|
||||
subroutine mld_d_base_onelev_setsv(lv,val,info,pos)
|
||||
import :: psb_dpk_, mld_d_onelev_type, mld_d_base_solver_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -303,7 +303,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_setc(lv,what,val,info,pos)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_d_onelev_type), intent(inout) :: lv
|
||||
@@ -318,7 +318,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_setr(lv,what,val,info,pos)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
class(mld_d_onelev_type), intent(inout) :: lv
|
||||
@@ -334,7 +334,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_cseti(lv,what,val,info,pos)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -350,7 +350,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_csetc(lv,what,val,info,pos)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_d_onelev_type), intent(inout) :: lv
|
||||
@@ -365,7 +365,7 @@ module mld_d_onelev_mod
|
||||
subroutine mld_d_base_onelev_csetr(lv,what,val,info,pos)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
class(mld_d_onelev_type), intent(inout) :: lv
|
||||
@@ -381,7 +381,7 @@ module mld_d_onelev_mod
|
||||
& solver,global_num)
|
||||
import :: psb_dspmat_type, psb_d_vect_type, psb_d_base_vect_type, &
|
||||
& psb_dlinmap_type, psb_dpk_, mld_d_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
class(mld_d_onelev_type), intent(in) :: lv
|
||||
integer(psb_ipk_), intent(in) :: level
|
||||
@@ -400,7 +400,7 @@ contains
|
||||
function d_base_onelev_get_nzeros(lv) result(val)
|
||||
implicit none
|
||||
class(mld_d_onelev_type), intent(in) :: lv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(lv%sm)) &
|
||||
@@ -412,7 +412,7 @@ contains
|
||||
function d_base_onelev_sizeof(lv) result(val)
|
||||
implicit none
|
||||
class(mld_d_onelev_type), intent(in) :: lv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -425,7 +425,7 @@ contains
|
||||
function mld_d_get_nzeros(prec) result(val)
|
||||
implicit none
|
||||
class(mld_dprec_type), intent(in) :: prec
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(prec%precv)) then
|
||||
@@ -438,10 +438,10 @@ contains
|
||||
function mld_dprec_sizeof(prec) result(val)
|
||||
implicit none
|
||||
class(mld_dprec_type), intent(in) :: prec
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
val = val + psb_sizeof_int
|
||||
val = val + psb_sizeof_ip
|
||||
if (allocated(prec%precv)) then
|
||||
do i=1, size(prec%precv)
|
||||
val = val + prec%precv(i)%sizeof()
|
||||
|
||||
@@ -402,10 +402,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_slu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
|
||||
val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
val = val + sv%symbsize
|
||||
val = val + sv%numsize
|
||||
return
|
||||
|
||||
@@ -411,10 +411,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_sludist_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
|
||||
val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
val = val + sv%symbsize
|
||||
val = val + sv%numsize
|
||||
return
|
||||
|
||||
@@ -406,10 +406,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_d_umf_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
|
||||
val = 2*psb_sizeof_long_int
|
||||
val = 2*psb_sizeof_lp
|
||||
val = val + sv%symbsize
|
||||
val = val + sv%numsize
|
||||
return
|
||||
|
||||
@@ -111,7 +111,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_check(sm,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -121,7 +121,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_restr_v(sm,x,trans,work,info,data)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
@@ -136,7 +136,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_restr_a(sm,x,trans,work,info,data)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
@@ -151,7 +151,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_prol_v(sm,x,trans,work,info,data)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
@@ -166,7 +166,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_prol_a(sm,x,trans,work,info,data)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
@@ -183,7 +183,7 @@ module mld_s_as_smoother
|
||||
subroutine mld_s_as_smoother_apply_vect(alpha,sm,x,beta,y,desc_data,&
|
||||
& trans,sweeps,work,wv,info,init,initu)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
type(psb_desc_type), intent(in) :: desc_data
|
||||
@@ -205,7 +205,7 @@ module mld_s_as_smoother
|
||||
subroutine mld_s_as_smoother_apply(alpha,sm,x,beta,y,desc_data,&
|
||||
& trans,sweeps,work,info,init,initu)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_,&
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_,&
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
type(psb_desc_type), intent(in) :: desc_data
|
||||
@@ -225,7 +225,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_s_base_sparse_mat, psb_ipk_,&
|
||||
& psb_i_base_vect_type
|
||||
implicit none
|
||||
@@ -242,7 +242,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_cnv(sm,info,amold,vmold,imold)
|
||||
import :: psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, &
|
||||
& psb_s_base_sparse_mat, psb_ipk_, psb_i_base_vect_type
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
@@ -256,7 +256,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_seti(sm,what,val,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -268,7 +268,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_setc(sm,what,val,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -280,7 +280,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_setr(sm,what,val,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -292,7 +292,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_cseti(sm,what,val,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -304,7 +304,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_csetc(sm,what,val,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -316,7 +316,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_csetr(sm,what,val,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -328,7 +328,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_free(sm,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -338,7 +338,7 @@ module mld_s_as_smoother
|
||||
interface
|
||||
subroutine mld_s_as_smoother_dmp(sm,ictxt,level,info,prefix,head,smoother,solver)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_long_int_k_, psb_desc_type, &
|
||||
& psb_spk_, mld_s_as_smoother_type, psb_epk_, psb_desc_type, &
|
||||
& psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(in) :: sm
|
||||
@@ -366,10 +366,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_as_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
if (allocated(sm%sv)) val = val + sm%sv%sizeof()
|
||||
val = val + sm%nd%sizeof()
|
||||
|
||||
@@ -379,7 +379,7 @@ contains
|
||||
function s_as_smoother_get_nzeros(sm) result(val)
|
||||
implicit none
|
||||
class(mld_s_as_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(sm%sv)) &
|
||||
|
||||
@@ -60,7 +60,7 @@
|
||||
module mld_s_base_smoother_mod
|
||||
|
||||
use mld_s_base_solver_mod
|
||||
use psb_base_mod, only : psb_desc_type, psb_sspmat_type, psb_long_int_k_,&
|
||||
use psb_base_mod, only : psb_desc_type, psb_sspmat_type, psb_epk_,&
|
||||
& psb_s_vect_type, psb_s_base_vect_type, psb_s_base_sparse_mat, &
|
||||
& psb_spk_, psb_i_base_vect_type, psb_erractionsave, psb_error_handler
|
||||
|
||||
@@ -355,7 +355,7 @@ contains
|
||||
function s_base_smoother_get_nzeros(sm) result(val)
|
||||
implicit none
|
||||
class(mld_s_base_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(sm%sv)) &
|
||||
@@ -366,7 +366,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_base_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -368,7 +368,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_base_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
|
||||
@@ -378,7 +378,7 @@ contains
|
||||
function s_base_solver_get_nzeros(sv) result(val)
|
||||
implicit none
|
||||
class(mld_s_base_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
end function s_base_solver_get_nzeros
|
||||
|
||||
@@ -237,7 +237,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_diag_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -250,7 +250,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_diag_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -543,7 +543,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_gs_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -558,10 +558,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_gs_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
val = val + sv%l%sizeof()
|
||||
val = val + sv%u%sizeof()
|
||||
|
||||
|
||||
@@ -517,7 +517,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_ilu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -533,10 +533,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_ilu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 2*psb_sizeof_int + psb_sizeof_sp
|
||||
val = 2*psb_sizeof_ip + psb_sizeof_sp
|
||||
val = val + sv%dv%sizeof()
|
||||
val = val + sv%l%sizeof()
|
||||
val = val + sv%u%sizeof()
|
||||
|
||||
@@ -156,7 +156,7 @@ module mld_s_jac_smoother
|
||||
interface
|
||||
subroutine mld_s_jac_smoother_dmp(sm,ictxt,level,info,prefix,head,smoother,solver)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_spk_, mld_s_jac_smoother_type, psb_long_int_k_, psb_desc_type, &
|
||||
& psb_spk_, mld_s_jac_smoother_type, psb_epk_, psb_desc_type, &
|
||||
& psb_ipk_
|
||||
implicit none
|
||||
class(mld_s_jac_smoother_type), intent(in) :: sm
|
||||
@@ -230,10 +230,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_jac_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
if (allocated(sm%sv)) val = val + sm%sv%sizeof()
|
||||
val = val + sm%nd%sizeof()
|
||||
|
||||
@@ -245,7 +245,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_jac_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -468,14 +468,14 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_mumps_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
#if defined(HAVE_MUMPS_)
|
||||
val = (sv%id%INFOG(22)+sv%id%INFOG(32))*1d+6
|
||||
#else
|
||||
val = 0
|
||||
#endif
|
||||
! val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
! val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
! val = val + sv%symbsize
|
||||
! val = val + sv%numsize
|
||||
return
|
||||
|
||||
+16
-16
@@ -57,7 +57,7 @@ module mld_s_onelev_mod
|
||||
use mld_s_base_smoother_mod
|
||||
use psb_base_mod, only : psb_sspmat_type, psb_s_vect_type, &
|
||||
& psb_s_base_vect_type, psb_slinmap_type, psb_spk_, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type, psb_i_base_vect_type, &
|
||||
& psb_ipk_, psb_epk_, psb_desc_type, psb_i_base_vect_type, &
|
||||
& psb_erractionsave, psb_error_handler
|
||||
!
|
||||
!
|
||||
@@ -194,7 +194,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_build(lv,info,amold,vmold,imold)
|
||||
import :: psb_s_base_sparse_mat, psb_s_base_vect_type, &
|
||||
& psb_i_base_vect_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
class(mld_s_onelev_type), target, intent(inout) :: lv
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -208,7 +208,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_descr(lv,il,nl,ilmin,info,iout)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_s_onelev_type), intent(in) :: lv
|
||||
@@ -235,7 +235,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_free(lv,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
|
||||
class(mld_s_onelev_type), intent(inout) :: lv
|
||||
@@ -247,7 +247,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_check(lv,info)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_s_onelev_type), intent(inout) :: lv
|
||||
@@ -259,7 +259,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_seti(lv,what,val,info,pos)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -274,7 +274,7 @@ module mld_s_onelev_mod
|
||||
interface
|
||||
subroutine mld_s_base_onelev_setsm(lv,val,info,pos)
|
||||
import :: psb_spk_, mld_s_onelev_type, mld_s_base_smoother_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -288,7 +288,7 @@ module mld_s_onelev_mod
|
||||
interface
|
||||
subroutine mld_s_base_onelev_setsv(lv,val,info,pos)
|
||||
import :: psb_spk_, mld_s_onelev_type, mld_s_base_solver_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -303,7 +303,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_setc(lv,what,val,info,pos)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_s_onelev_type), intent(inout) :: lv
|
||||
@@ -318,7 +318,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_setr(lv,what,val,info,pos)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
class(mld_s_onelev_type), intent(inout) :: lv
|
||||
@@ -334,7 +334,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_cseti(lv,what,val,info,pos)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -350,7 +350,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_csetc(lv,what,val,info,pos)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_s_onelev_type), intent(inout) :: lv
|
||||
@@ -365,7 +365,7 @@ module mld_s_onelev_mod
|
||||
subroutine mld_s_base_onelev_csetr(lv,what,val,info,pos)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
class(mld_s_onelev_type), intent(inout) :: lv
|
||||
@@ -381,7 +381,7 @@ module mld_s_onelev_mod
|
||||
& solver,global_num)
|
||||
import :: psb_sspmat_type, psb_s_vect_type, psb_s_base_vect_type, &
|
||||
& psb_slinmap_type, psb_spk_, mld_s_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
class(mld_s_onelev_type), intent(in) :: lv
|
||||
integer(psb_ipk_), intent(in) :: level
|
||||
@@ -400,7 +400,7 @@ contains
|
||||
function s_base_onelev_get_nzeros(lv) result(val)
|
||||
implicit none
|
||||
class(mld_s_onelev_type), intent(in) :: lv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(lv%sm)) &
|
||||
@@ -412,7 +412,7 @@ contains
|
||||
function s_base_onelev_sizeof(lv) result(val)
|
||||
implicit none
|
||||
class(mld_s_onelev_type), intent(in) :: lv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -425,7 +425,7 @@ contains
|
||||
function mld_s_get_nzeros(prec) result(val)
|
||||
implicit none
|
||||
class(mld_sprec_type), intent(in) :: prec
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(prec%precv)) then
|
||||
@@ -438,10 +438,10 @@ contains
|
||||
function mld_sprec_sizeof(prec) result(val)
|
||||
implicit none
|
||||
class(mld_sprec_type), intent(in) :: prec
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
val = val + psb_sizeof_int
|
||||
val = val + psb_sizeof_ip
|
||||
if (allocated(prec%precv)) then
|
||||
do i=1, size(prec%precv)
|
||||
val = val + prec%precv(i)%sizeof()
|
||||
|
||||
@@ -402,10 +402,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_s_slu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
|
||||
val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
val = val + sv%symbsize
|
||||
val = val + sv%numsize
|
||||
return
|
||||
|
||||
@@ -111,7 +111,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_check(sm,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -121,7 +121,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_restr_v(sm,x,trans,work,info,data)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
@@ -136,7 +136,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_restr_a(sm,x,trans,work,info,data)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
@@ -151,7 +151,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_prol_v(sm,x,trans,work,info,data)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
@@ -166,7 +166,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_prol_a(sm,x,trans,work,info,data)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
@@ -183,7 +183,7 @@ module mld_z_as_smoother
|
||||
subroutine mld_z_as_smoother_apply_vect(alpha,sm,x,beta,y,desc_data,&
|
||||
& trans,sweeps,work,wv,info,init,initu)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
type(psb_desc_type), intent(in) :: desc_data
|
||||
@@ -205,7 +205,7 @@ module mld_z_as_smoother
|
||||
subroutine mld_z_as_smoother_apply(alpha,sm,x,beta,y,desc_data,&
|
||||
& trans,sweeps,work,info,init,initu)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_,&
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_,&
|
||||
& psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
type(psb_desc_type), intent(in) :: desc_data
|
||||
@@ -225,7 +225,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, &
|
||||
& psb_desc_type, psb_z_base_sparse_mat, psb_ipk_,&
|
||||
& psb_i_base_vect_type
|
||||
implicit none
|
||||
@@ -242,7 +242,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_cnv(sm,info,amold,vmold,imold)
|
||||
import :: psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, &
|
||||
& psb_z_base_sparse_mat, psb_ipk_, psb_i_base_vect_type
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
@@ -256,7 +256,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_seti(sm,what,val,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -268,7 +268,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_setc(sm,what,val,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -280,7 +280,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_setr(sm,what,val,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(in) :: what
|
||||
@@ -292,7 +292,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_cseti(sm,what,val,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -304,7 +304,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_csetc(sm,what,val,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -316,7 +316,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_csetr(sm,what,val,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
character(len=*), intent(in) :: what
|
||||
@@ -328,7 +328,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_free(sm,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, psb_desc_type, psb_ipk_
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, psb_desc_type, psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(inout) :: sm
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -338,7 +338,7 @@ module mld_z_as_smoother
|
||||
interface
|
||||
subroutine mld_z_as_smoother_dmp(sm,ictxt,level,info,prefix,head,smoother,solver)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_long_int_k_, psb_desc_type, &
|
||||
& psb_dpk_, mld_z_as_smoother_type, psb_epk_, psb_desc_type, &
|
||||
& psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(in) :: sm
|
||||
@@ -366,10 +366,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_as_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
if (allocated(sm%sv)) val = val + sm%sv%sizeof()
|
||||
val = val + sm%nd%sizeof()
|
||||
|
||||
@@ -379,7 +379,7 @@ contains
|
||||
function z_as_smoother_get_nzeros(sm) result(val)
|
||||
implicit none
|
||||
class(mld_z_as_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(sm%sv)) &
|
||||
|
||||
@@ -60,7 +60,7 @@
|
||||
module mld_z_base_smoother_mod
|
||||
|
||||
use mld_z_base_solver_mod
|
||||
use psb_base_mod, only : psb_desc_type, psb_zspmat_type, psb_long_int_k_,&
|
||||
use psb_base_mod, only : psb_desc_type, psb_zspmat_type, psb_epk_,&
|
||||
& psb_z_vect_type, psb_z_base_vect_type, psb_z_base_sparse_mat, &
|
||||
& psb_dpk_, psb_i_base_vect_type, psb_erractionsave, psb_error_handler
|
||||
|
||||
@@ -355,7 +355,7 @@ contains
|
||||
function z_base_smoother_get_nzeros(sm) result(val)
|
||||
implicit none
|
||||
class(mld_z_base_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(sm%sv)) &
|
||||
@@ -366,7 +366,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_base_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -368,7 +368,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_base_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
|
||||
@@ -378,7 +378,7 @@ contains
|
||||
function z_base_solver_get_nzeros(sv) result(val)
|
||||
implicit none
|
||||
class(mld_z_base_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
end function z_base_solver_get_nzeros
|
||||
|
||||
@@ -237,7 +237,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_diag_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -250,7 +250,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_diag_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -543,7 +543,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_gs_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -558,10 +558,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_gs_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
val = val + sv%l%sizeof()
|
||||
val = val + sv%u%sizeof()
|
||||
|
||||
|
||||
@@ -517,7 +517,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_ilu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
@@ -533,10 +533,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_ilu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 2*psb_sizeof_int + (2*psb_sizeof_dp)
|
||||
val = 2*psb_sizeof_ip + (2*psb_sizeof_dp)
|
||||
val = val + sv%dv%sizeof()
|
||||
val = val + sv%l%sizeof()
|
||||
val = val + sv%u%sizeof()
|
||||
|
||||
@@ -156,7 +156,7 @@ module mld_z_jac_smoother
|
||||
interface
|
||||
subroutine mld_z_jac_smoother_dmp(sm,ictxt,level,info,prefix,head,smoother,solver)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_dpk_, mld_z_jac_smoother_type, psb_long_int_k_, psb_desc_type, &
|
||||
& psb_dpk_, mld_z_jac_smoother_type, psb_epk_, psb_desc_type, &
|
||||
& psb_ipk_
|
||||
implicit none
|
||||
class(mld_z_jac_smoother_type), intent(in) :: sm
|
||||
@@ -230,10 +230,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_jac_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = psb_sizeof_int
|
||||
val = psb_sizeof_ip
|
||||
if (allocated(sm%sv)) val = val + sm%sv%sizeof()
|
||||
val = val + sm%nd%sizeof()
|
||||
|
||||
@@ -245,7 +245,7 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_jac_smoother_type), intent(in) :: sm
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -468,14 +468,14 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_mumps_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
#if defined(HAVE_MUMPS_)
|
||||
val = (sv%id%INFOG(22)+sv%id%INFOG(32))*1d+6
|
||||
#else
|
||||
val = 0
|
||||
#endif
|
||||
! val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
! val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
! val = val + sv%symbsize
|
||||
! val = val + sv%numsize
|
||||
return
|
||||
|
||||
+16
-16
@@ -57,7 +57,7 @@ module mld_z_onelev_mod
|
||||
use mld_z_base_smoother_mod
|
||||
use psb_base_mod, only : psb_zspmat_type, psb_z_vect_type, &
|
||||
& psb_z_base_vect_type, psb_zlinmap_type, psb_dpk_, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type, psb_i_base_vect_type, &
|
||||
& psb_ipk_, psb_epk_, psb_desc_type, psb_i_base_vect_type, &
|
||||
& psb_erractionsave, psb_error_handler
|
||||
!
|
||||
!
|
||||
@@ -194,7 +194,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_build(lv,info,amold,vmold,imold)
|
||||
import :: psb_z_base_sparse_mat, psb_z_base_vect_type, &
|
||||
& psb_i_base_vect_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
class(mld_z_onelev_type), target, intent(inout) :: lv
|
||||
integer(psb_ipk_), intent(out) :: info
|
||||
@@ -208,7 +208,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_descr(lv,il,nl,ilmin,info,iout)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_z_onelev_type), intent(in) :: lv
|
||||
@@ -235,7 +235,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_free(lv,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
|
||||
class(mld_z_onelev_type), intent(inout) :: lv
|
||||
@@ -247,7 +247,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_check(lv,info)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_z_onelev_type), intent(inout) :: lv
|
||||
@@ -259,7 +259,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_seti(lv,what,val,info,pos)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -274,7 +274,7 @@ module mld_z_onelev_mod
|
||||
interface
|
||||
subroutine mld_z_base_onelev_setsm(lv,val,info,pos)
|
||||
import :: psb_dpk_, mld_z_onelev_type, mld_z_base_smoother_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -288,7 +288,7 @@ module mld_z_onelev_mod
|
||||
interface
|
||||
subroutine mld_z_base_onelev_setsv(lv,val,info,pos)
|
||||
import :: psb_dpk_, mld_z_onelev_type, mld_z_base_solver_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -303,7 +303,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_setc(lv,what,val,info,pos)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_z_onelev_type), intent(inout) :: lv
|
||||
@@ -318,7 +318,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_setr(lv,what,val,info,pos)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
class(mld_z_onelev_type), intent(inout) :: lv
|
||||
@@ -334,7 +334,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_cseti(lv,what,val,info,pos)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
! Arguments
|
||||
@@ -350,7 +350,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_csetc(lv,what,val,info,pos)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
! Arguments
|
||||
class(mld_z_onelev_type), intent(inout) :: lv
|
||||
@@ -365,7 +365,7 @@ module mld_z_onelev_mod
|
||||
subroutine mld_z_base_onelev_csetr(lv,what,val,info,pos)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
Implicit None
|
||||
|
||||
class(mld_z_onelev_type), intent(inout) :: lv
|
||||
@@ -381,7 +381,7 @@ module mld_z_onelev_mod
|
||||
& solver,global_num)
|
||||
import :: psb_zspmat_type, psb_z_vect_type, psb_z_base_vect_type, &
|
||||
& psb_zlinmap_type, psb_dpk_, mld_z_onelev_type, &
|
||||
& psb_ipk_, psb_long_int_k_, psb_desc_type
|
||||
& psb_ipk_, psb_epk_, psb_desc_type
|
||||
implicit none
|
||||
class(mld_z_onelev_type), intent(in) :: lv
|
||||
integer(psb_ipk_), intent(in) :: level
|
||||
@@ -400,7 +400,7 @@ contains
|
||||
function z_base_onelev_get_nzeros(lv) result(val)
|
||||
implicit none
|
||||
class(mld_z_onelev_type), intent(in) :: lv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(lv%sm)) &
|
||||
@@ -412,7 +412,7 @@ contains
|
||||
function z_base_onelev_sizeof(lv) result(val)
|
||||
implicit none
|
||||
class(mld_z_onelev_type), intent(in) :: lv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
|
||||
val = 0
|
||||
|
||||
@@ -425,7 +425,7 @@ contains
|
||||
function mld_z_get_nzeros(prec) result(val)
|
||||
implicit none
|
||||
class(mld_zprec_type), intent(in) :: prec
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
if (allocated(prec%precv)) then
|
||||
@@ -438,10 +438,10 @@ contains
|
||||
function mld_zprec_sizeof(prec) result(val)
|
||||
implicit none
|
||||
class(mld_zprec_type), intent(in) :: prec
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer(psb_ipk_) :: i
|
||||
val = 0
|
||||
val = val + psb_sizeof_int
|
||||
val = val + psb_sizeof_ip
|
||||
if (allocated(prec%precv)) then
|
||||
do i=1, size(prec%precv)
|
||||
val = val + prec%precv(i)%sizeof()
|
||||
|
||||
@@ -402,10 +402,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_slu_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
|
||||
val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
val = val + sv%symbsize
|
||||
val = val + sv%numsize
|
||||
return
|
||||
|
||||
@@ -411,10 +411,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_sludist_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
|
||||
val = 2*psb_sizeof_int + psb_sizeof_dp
|
||||
val = 2*psb_sizeof_ip + psb_sizeof_dp
|
||||
val = val + sv%symbsize
|
||||
val = val + sv%numsize
|
||||
return
|
||||
|
||||
@@ -406,10 +406,10 @@ contains
|
||||
implicit none
|
||||
! Arguments
|
||||
class(mld_z_umf_solver_type), intent(in) :: sv
|
||||
integer(psb_long_int_k_) :: val
|
||||
integer(psb_epk_) :: val
|
||||
integer :: i
|
||||
|
||||
val = 2*psb_sizeof_long_int
|
||||
val = 2*psb_sizeof_lp
|
||||
val = val + sv%symbsize
|
||||
val = val + sv%numsize
|
||||
return
|
||||
|
||||
@@ -138,7 +138,7 @@ program mld_cf_sample
|
||||
|
||||
! solver paramters
|
||||
integer(psb_ipk_) :: iter, ircode, nlv
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
real(psb_spk_) :: err
|
||||
|
||||
character(len=5) :: afmt
|
||||
|
||||
@@ -138,7 +138,7 @@ program mld_df_sample
|
||||
|
||||
! solver paramters
|
||||
integer(psb_ipk_) :: iter, ircode, nlv
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
real(psb_dpk_) :: err
|
||||
|
||||
character(len=5) :: afmt
|
||||
|
||||
@@ -138,7 +138,7 @@ program mld_sf_sample
|
||||
|
||||
! solver paramters
|
||||
integer(psb_ipk_) :: iter, ircode, nlv
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
real(psb_spk_) :: err
|
||||
|
||||
character(len=5) :: afmt
|
||||
|
||||
@@ -138,7 +138,7 @@ program mld_zf_sample
|
||||
|
||||
! solver paramters
|
||||
integer(psb_ipk_) :: iter, ircode, nlv
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
real(psb_dpk_) :: err
|
||||
|
||||
character(len=5) :: afmt
|
||||
|
||||
@@ -153,7 +153,7 @@ program mld_d_pde3d
|
||||
|
||||
! solver parameters
|
||||
integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, nlv
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
real(psb_dpk_) :: err, eps
|
||||
|
||||
! other variables
|
||||
|
||||
@@ -445,7 +445,7 @@ program mld_d_pde2d
|
||||
|
||||
! solver parameters
|
||||
integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, nlv
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
real(psb_dpk_) :: err, resmx, resmxp
|
||||
|
||||
! Krylov solver data
|
||||
|
||||
@@ -498,7 +498,7 @@ program mld_d_pde3d
|
||||
|
||||
! solver parameters
|
||||
integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, nlv
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
real(psb_dpk_) :: err, resmx, resmxp
|
||||
|
||||
! Krylov solver data
|
||||
|
||||
@@ -445,7 +445,7 @@ program mld_s_pde2d
|
||||
|
||||
! solver parameters
|
||||
integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, nlv
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
real(psb_spk_) :: err, resmx, resmxp
|
||||
|
||||
! Krylov solver data
|
||||
|
||||
@@ -498,7 +498,7 @@ program mld_s_pde3d
|
||||
|
||||
! solver parameters
|
||||
integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst, nlv
|
||||
integer(psb_long_int_k_) :: amatsize, precsize, descsize
|
||||
integer(psb_epk_) :: amatsize, precsize, descsize
|
||||
real(psb_spk_) :: err, resmx, resmxp
|
||||
|
||||
! Krylov solver data
|
||||
|
||||
Reference in New Issue
Block a user