New ICTXT arg in INIT.

MixedI8
Salvatore Filippone 6 years ago
parent 59857aee02
commit 78120dc1c0

@ -29,7 +29,7 @@
! POSSIBILITY OF SUCH DAMAGE. ! POSSIBILITY OF SUCH DAMAGE.
! !
! !
subroutine psb_cprecinit(p,ptype,info) subroutine psb_cprecinit(ictxt,p,ptype,info)
use psb_base_mod use psb_base_mod
use psb_c_prec_type, psb_protect_name => psb_cprecinit use psb_c_prec_type, psb_protect_name => psb_cprecinit
@ -37,7 +37,8 @@ subroutine psb_cprecinit(p,ptype,info)
use psb_c_diagprec, only : psb_c_diag_prec_type use psb_c_diagprec, only : psb_c_diag_prec_type
use psb_c_bjacprec, only : psb_c_bjac_prec_type use psb_c_bjacprec, only : psb_c_bjac_prec_type
implicit none implicit none
class(psb_cprec_type), intent(inout) :: p integer(psb_ipk_), intent(in) :: ictxt
class(psb_cprec_type), intent(inout) :: p
character(len=*), intent(in) :: ptype character(len=*), intent(in) :: ptype
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info
@ -49,6 +50,8 @@ subroutine psb_cprecinit(p,ptype,info)
if (info /= psb_success_) return if (info /= psb_success_) return
end if end if
p%ictxt = ictxt
select case(psb_toupper(ptype(1:len_trim(ptype)))) select case(psb_toupper(ptype(1:len_trim(ptype))))
case ('NONE','NOPREC') case ('NONE','NOPREC')

@ -29,7 +29,7 @@
! POSSIBILITY OF SUCH DAMAGE. ! POSSIBILITY OF SUCH DAMAGE.
! !
! !
subroutine psb_dprecinit(p,ptype,info) subroutine psb_dprecinit(ictxt,p,ptype,info)
use psb_base_mod use psb_base_mod
use psb_d_prec_type, psb_protect_name => psb_dprecinit use psb_d_prec_type, psb_protect_name => psb_dprecinit
@ -37,7 +37,8 @@ subroutine psb_dprecinit(p,ptype,info)
use psb_d_diagprec, only : psb_d_diag_prec_type use psb_d_diagprec, only : psb_d_diag_prec_type
use psb_d_bjacprec, only : psb_d_bjac_prec_type use psb_d_bjacprec, only : psb_d_bjac_prec_type
implicit none implicit none
class(psb_dprec_type), intent(inout) :: p integer(psb_ipk_), intent(in) :: ictxt
class(psb_dprec_type), intent(inout) :: p
character(len=*), intent(in) :: ptype character(len=*), intent(in) :: ptype
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info
@ -49,6 +50,8 @@ subroutine psb_dprecinit(p,ptype,info)
if (info /= psb_success_) return if (info /= psb_success_) return
end if end if
p%ictxt = ictxt
select case(psb_toupper(ptype(1:len_trim(ptype)))) select case(psb_toupper(ptype(1:len_trim(ptype))))
case ('NONE','NOPREC') case ('NONE','NOPREC')

@ -29,7 +29,7 @@
! POSSIBILITY OF SUCH DAMAGE. ! POSSIBILITY OF SUCH DAMAGE.
! !
! !
subroutine psb_sprecinit(p,ptype,info) subroutine psb_sprecinit(ictxt,p,ptype,info)
use psb_base_mod use psb_base_mod
use psb_s_prec_type, psb_protect_name => psb_sprecinit use psb_s_prec_type, psb_protect_name => psb_sprecinit
@ -37,7 +37,8 @@ subroutine psb_sprecinit(p,ptype,info)
use psb_s_diagprec, only : psb_s_diag_prec_type use psb_s_diagprec, only : psb_s_diag_prec_type
use psb_s_bjacprec, only : psb_s_bjac_prec_type use psb_s_bjacprec, only : psb_s_bjac_prec_type
implicit none implicit none
class(psb_sprec_type), intent(inout) :: p integer(psb_ipk_), intent(in) :: ictxt
class(psb_sprec_type), intent(inout) :: p
character(len=*), intent(in) :: ptype character(len=*), intent(in) :: ptype
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info
@ -49,6 +50,8 @@ subroutine psb_sprecinit(p,ptype,info)
if (info /= psb_success_) return if (info /= psb_success_) return
end if end if
p%ictxt = ictxt
select case(psb_toupper(ptype(1:len_trim(ptype)))) select case(psb_toupper(ptype(1:len_trim(ptype))))
case ('NONE','NOPREC') case ('NONE','NOPREC')

@ -29,7 +29,7 @@
! POSSIBILITY OF SUCH DAMAGE. ! POSSIBILITY OF SUCH DAMAGE.
! !
! !
subroutine psb_zprecinit(p,ptype,info) subroutine psb_zprecinit(ictxt,p,ptype,info)
use psb_base_mod use psb_base_mod
use psb_z_prec_type, psb_protect_name => psb_zprecinit use psb_z_prec_type, psb_protect_name => psb_zprecinit
@ -37,7 +37,8 @@ subroutine psb_zprecinit(p,ptype,info)
use psb_z_diagprec, only : psb_z_diag_prec_type use psb_z_diagprec, only : psb_z_diag_prec_type
use psb_z_bjacprec, only : psb_z_bjac_prec_type use psb_z_bjacprec, only : psb_z_bjac_prec_type
implicit none implicit none
class(psb_zprec_type), intent(inout) :: p integer(psb_ipk_), intent(in) :: ictxt
class(psb_zprec_type), intent(inout) :: p
character(len=*), intent(in) :: ptype character(len=*), intent(in) :: ptype
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info
@ -49,6 +50,8 @@ subroutine psb_zprecinit(p,ptype,info)
if (info /= psb_success_) return if (info /= psb_success_) return
end if end if
p%ictxt = ictxt
select case(psb_toupper(ptype(1:len_trim(ptype)))) select case(psb_toupper(ptype(1:len_trim(ptype))))
case ('NONE','NOPREC') case ('NONE','NOPREC')

@ -39,6 +39,7 @@ module psb_c_prec_type
use psb_c_base_prec_mod use psb_c_base_prec_mod
type psb_cprec_type type psb_cprec_type
integer(psb_ipk_) :: ictxt
class(psb_c_base_prec_type), allocatable :: prec class(psb_c_base_prec_type), allocatable :: prec
contains contains
procedure, pass(prec) :: psb_c_apply1_vect procedure, pass(prec) :: psb_c_apply1_vect
@ -60,9 +61,10 @@ module psb_c_prec_type
end interface end interface
interface psb_precinit interface psb_precinit
subroutine psb_cprecinit(prec,ptype,info) subroutine psb_cprecinit(ictxt,prec,ptype,info)
import :: psb_ipk_, psb_cprec_type import :: psb_ipk_, psb_cprec_type
implicit none implicit none
integer(psb_ipk_), intent(in) :: ictxt
class(psb_cprec_type), intent(inout) :: prec class(psb_cprec_type), intent(inout) :: prec
character(len=*), intent(in) :: ptype character(len=*), intent(in) :: ptype
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info

@ -39,6 +39,7 @@ module psb_d_prec_type
use psb_d_base_prec_mod use psb_d_base_prec_mod
type psb_dprec_type type psb_dprec_type
integer(psb_ipk_) :: ictxt
class(psb_d_base_prec_type), allocatable :: prec class(psb_d_base_prec_type), allocatable :: prec
contains contains
procedure, pass(prec) :: psb_d_apply1_vect procedure, pass(prec) :: psb_d_apply1_vect
@ -60,9 +61,10 @@ module psb_d_prec_type
end interface end interface
interface psb_precinit interface psb_precinit
subroutine psb_dprecinit(prec,ptype,info) subroutine psb_dprecinit(ictxt,prec,ptype,info)
import :: psb_ipk_, psb_dprec_type import :: psb_ipk_, psb_dprec_type
implicit none implicit none
integer(psb_ipk_), intent(in) :: ictxt
class(psb_dprec_type), intent(inout) :: prec class(psb_dprec_type), intent(inout) :: prec
character(len=*), intent(in) :: ptype character(len=*), intent(in) :: ptype
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info

@ -39,6 +39,7 @@ module psb_s_prec_type
use psb_s_base_prec_mod use psb_s_base_prec_mod
type psb_sprec_type type psb_sprec_type
integer(psb_ipk_) :: ictxt
class(psb_s_base_prec_type), allocatable :: prec class(psb_s_base_prec_type), allocatable :: prec
contains contains
procedure, pass(prec) :: psb_s_apply1_vect procedure, pass(prec) :: psb_s_apply1_vect
@ -60,9 +61,10 @@ module psb_s_prec_type
end interface end interface
interface psb_precinit interface psb_precinit
subroutine psb_sprecinit(prec,ptype,info) subroutine psb_sprecinit(ictxt,prec,ptype,info)
import :: psb_ipk_, psb_sprec_type import :: psb_ipk_, psb_sprec_type
implicit none implicit none
integer(psb_ipk_), intent(in) :: ictxt
class(psb_sprec_type), intent(inout) :: prec class(psb_sprec_type), intent(inout) :: prec
character(len=*), intent(in) :: ptype character(len=*), intent(in) :: ptype
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info

@ -39,6 +39,7 @@ module psb_z_prec_type
use psb_z_base_prec_mod use psb_z_base_prec_mod
type psb_zprec_type type psb_zprec_type
integer(psb_ipk_) :: ictxt
class(psb_z_base_prec_type), allocatable :: prec class(psb_z_base_prec_type), allocatable :: prec
contains contains
procedure, pass(prec) :: psb_z_apply1_vect procedure, pass(prec) :: psb_z_apply1_vect
@ -60,9 +61,10 @@ module psb_z_prec_type
end interface end interface
interface psb_precinit interface psb_precinit
subroutine psb_zprecinit(prec,ptype,info) subroutine psb_zprecinit(ictxt,prec,ptype,info)
import :: psb_ipk_, psb_zprec_type import :: psb_ipk_, psb_zprec_type
implicit none implicit none
integer(psb_ipk_), intent(in) :: ictxt
class(psb_zprec_type), intent(inout) :: prec class(psb_zprec_type), intent(inout) :: prec
character(len=*), intent(in) :: ptype character(len=*), intent(in) :: ptype
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info

@ -213,7 +213,7 @@ program psb_cf_sample
! !
call prec%init(ptype,info) call prec%init(ictxt,ptype,info)
! building the preconditioner ! building the preconditioner
t1 = psb_wtime() t1 = psb_wtime()

@ -213,7 +213,7 @@ program psb_df_sample
! !
call prec%init(ptype,info) call prec%init(ictxt,ptype,info)
! building the preconditioner ! building the preconditioner
t1 = psb_wtime() t1 = psb_wtime()

@ -213,7 +213,7 @@ program psb_sf_sample
! !
call prec%init(ptype,info) call prec%init(ictxt,ptype,info)
! building the preconditioner ! building the preconditioner
t1 = psb_wtime() t1 = psb_wtime()

@ -213,7 +213,7 @@ program psb_zf_sample
! !
call prec%init(ptype,info) call prec%init(ictxt,ptype,info)
! building the preconditioner ! building the preconditioner
t1 = psb_wtime() t1 = psb_wtime()

@ -609,7 +609,7 @@ program psb_d_pde2d
! prepare the preconditioner. ! prepare the preconditioner.
! !
if(iam == psb_root_) write(psb_out_unit,'("Setting preconditioner to : ",a)')ptype if(iam == psb_root_) write(psb_out_unit,'("Setting preconditioner to : ",a)')ptype
call prec%init(ptype,info) call prec%init(ictxt,ptype,info)
call psb_barrier(ictxt) call psb_barrier(ictxt)
t1 = psb_wtime() t1 = psb_wtime()

@ -650,7 +650,7 @@ program psb_d_pde3d
! prepare the preconditioner. ! prepare the preconditioner.
! !
if(iam == psb_root_) write(psb_out_unit,'("Setting preconditioner to : ",a)')ptype if(iam == psb_root_) write(psb_out_unit,'("Setting preconditioner to : ",a)')ptype
call prec%init(ptype,info) call prec%init(ictxt,ptype,info)
call psb_barrier(ictxt) call psb_barrier(ictxt)
t1 = psb_wtime() t1 = psb_wtime()

@ -609,7 +609,7 @@ program psb_s_pde2d
! prepare the preconditioner. ! prepare the preconditioner.
! !
if(iam == psb_root_) write(psb_out_unit,'("Setting preconditioner to : ",a)')ptype if(iam == psb_root_) write(psb_out_unit,'("Setting preconditioner to : ",a)')ptype
call prec%init(ptype,info) call prec%init(ictxt,ptype,info)
call psb_barrier(ictxt) call psb_barrier(ictxt)
t1 = psb_wtime() t1 = psb_wtime()

@ -650,7 +650,7 @@ program psb_s_pde3d
! prepare the preconditioner. ! prepare the preconditioner.
! !
if(iam == psb_root_) write(psb_out_unit,'("Setting preconditioner to : ",a)')ptype if(iam == psb_root_) write(psb_out_unit,'("Setting preconditioner to : ",a)')ptype
call prec%init(ptype,info) call prec%init(ictxt,ptype,info)
call psb_barrier(ictxt) call psb_barrier(ictxt)
t1 = psb_wtime() t1 = psb_wtime()

Loading…
Cancel
Save