get_mpicomm is now a function

stopcriterion
Salvatore Filippone 7 years ago
parent 8588c0e5c6
commit f4543fc8b5

@ -78,14 +78,14 @@ subroutine c_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
call psb_info(ictxt, iam, np)
if (sv%ipar(1) == mld_local_solver_ ) then
call psb_init(ictxt1,np=1,basectxt=ictxt,ids=(/iam/))
call psb_get_mpicomm(ictxt1, icomm)
icomm = psb_get_mpicomm(ictxt1)
allocate(sv%local_ictxt,stat=info)
sv%local_ictxt = ictxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm,sv%local_ictxt
call psb_info(ictxt1, me, np)
npr = np
else if (sv%ipar(1) == mld_global_solver_ ) then
call psb_get_mpicomm(ictxt,icomm)
icomm = psb_get_mpicomm(ictxt)
!write(*,*)iam,'mumps_bld: global +++++>',icomm,ictxt
call psb_info(ictxt, iam, np)
me = iam

@ -78,14 +78,14 @@ subroutine d_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
call psb_info(ictxt, iam, np)
if (sv%ipar(1) == mld_local_solver_ ) then
call psb_init(ictxt1,np=1,basectxt=ictxt,ids=(/iam/))
call psb_get_mpicomm(ictxt1, icomm)
icomm = psb_get_mpicomm(ictxt1)
allocate(sv%local_ictxt,stat=info)
sv%local_ictxt = ictxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm,sv%local_ictxt
call psb_info(ictxt1, me, np)
npr = np
else if (sv%ipar(1) == mld_global_solver_ ) then
call psb_get_mpicomm(ictxt,icomm)
icomm = psb_get_mpicomm(ictxt)
!write(*,*)iam,'mumps_bld: global +++++>',icomm,ictxt
call psb_info(ictxt, iam, np)
me = iam

@ -78,14 +78,14 @@ subroutine s_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
call psb_info(ictxt, iam, np)
if (sv%ipar(1) == mld_local_solver_ ) then
call psb_init(ictxt1,np=1,basectxt=ictxt,ids=(/iam/))
call psb_get_mpicomm(ictxt1, icomm)
icomm = psb_get_mpicomm(ictxt1)
allocate(sv%local_ictxt,stat=info)
sv%local_ictxt = ictxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm,sv%local_ictxt
call psb_info(ictxt1, me, np)
npr = np
else if (sv%ipar(1) == mld_global_solver_ ) then
call psb_get_mpicomm(ictxt,icomm)
icomm = psb_get_mpicomm(ictxt)
!write(*,*)iam,'mumps_bld: global +++++>',icomm,ictxt
call psb_info(ictxt, iam, np)
me = iam

@ -78,14 +78,14 @@ subroutine z_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
call psb_info(ictxt, iam, np)
if (sv%ipar(1) == mld_local_solver_ ) then
call psb_init(ictxt1,np=1,basectxt=ictxt,ids=(/iam/))
call psb_get_mpicomm(ictxt1, icomm)
icomm = psb_get_mpicomm(ictxt1)
allocate(sv%local_ictxt,stat=info)
sv%local_ictxt = ictxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm,sv%local_ictxt
call psb_info(ictxt1, me, np)
npr = np
else if (sv%ipar(1) == mld_global_solver_ ) then
call psb_get_mpicomm(ictxt,icomm)
icomm = psb_get_mpicomm(ictxt)
!write(*,*)iam,'mumps_bld: global +++++>',icomm,ictxt
call psb_info(ictxt, iam, np)
me = iam

Loading…
Cancel
Save