fdurastante/amg4psblas
1
0
Fork 0
mirror of https://github.com/sfilippone/amg4psblas.git synced 2026-10-06 14:44:55 +00:00
Code Issues Packages Projects Releases Wiki Activity

Merged MUMPS fixes

Browse Source
This commit is contained in:
Salvatore Filippone
2026-04-19 12:49:11 +02:00
parent f99b563f83
commit d45fffe482
4 changed files with 28 additions and 48 deletions
Download Patch File Download Diff File Expand all files Collapse all files
amgprec/impl/solver/amg_c_mumps_solver_bld.F90
+7 -12
View File
@@ -80,24 +80,19 @@ subroutine c_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
call psb_info(ctxt, iam, np)
if (sv%ipar(1) == amg_local_solver_ ) then
call psb_init(ctxt1,np=1,basectxt=ctxt,ids=(/iam/))
icomm = psb_get_mpi_comm(ctxt1)
allocate(sv%local_ctxt,stat=info)
sv%local_ctxt = ctxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm,sv%local_ctxt
call psb_info(ctxt1, me, np)
npr = np
else if (sv%ipar(1) == amg_global_solver_ ) then
icomm = psb_get_mpi_comm(ctxt)
!write(*,*)iam,'mumps_bld: global +++++>',icomm,ctxt
call psb_info(ctxt, iam, np)
me = iam
npr = np
call psb_init(ctxt1,basectxt=ctxt)
else
info = psb_err_internal_error_
call psb_errpush(info,name,&
& a_err='Invalid local/global solver in MUMPS')
goto 9999
end if
!allocate(sv%local_ctxt,stat=info)
icomm = ctxt1%get_mpic()
sv%local_ctxt = ctxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm!,sv%local_ctxt%ctxt
call psb_info(ctxt1, me, npr)
npc = 1
if (debug_level >= psb_debug_outer_) &
& write(debug_unit,*) me,' ',trim(name),' start'
@@ -115,8 +110,8 @@ subroutine c_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
end if
end if
call mpi_comm_dup(icomm,sv%id%comm,info)
sv%id%comm = icomm
sv%id%job = -1
sv%id%par = 1
if (sv%ipar(3) == 2) then
amgprec/impl/solver/amg_d_mumps_solver_bld.F90
+7 -12
View File
@@ -80,24 +80,19 @@ subroutine d_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
call psb_info(ctxt, iam, np)
if (sv%ipar(1) == amg_local_solver_ ) then
call psb_init(ctxt1,np=1,basectxt=ctxt,ids=(/iam/))
icomm = psb_get_mpi_comm(ctxt1)
allocate(sv%local_ctxt,stat=info)
sv%local_ctxt = ctxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm,sv%local_ctxt
call psb_info(ctxt1, me, np)
npr = np
else if (sv%ipar(1) == amg_global_solver_ ) then
icomm = psb_get_mpi_comm(ctxt)
!write(*,*)iam,'mumps_bld: global +++++>',icomm,ctxt
call psb_info(ctxt, iam, np)
me = iam
npr = np
call psb_init(ctxt1,basectxt=ctxt)
else
info = psb_err_internal_error_
call psb_errpush(info,name,&
& a_err='Invalid local/global solver in MUMPS')
goto 9999
end if
!allocate(sv%local_ctxt,stat=info)
icomm = ctxt1%get_mpic()
sv%local_ctxt = ctxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm!,sv%local_ctxt%ctxt
call psb_info(ctxt1, me, npr)
npc = 1
if (debug_level >= psb_debug_outer_) &
& write(debug_unit,*) me,' ',trim(name),' start'
@@ -115,8 +110,8 @@ subroutine d_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
end if
end if
call mpi_comm_dup(icomm,sv%id%comm,info)
sv%id%comm = icomm
sv%id%job = -1
sv%id%par = 1
if (sv%ipar(3) == 2) then
amgprec/impl/solver/amg_s_mumps_solver_bld.F90
+7 -12
View File
@@ -80,24 +80,19 @@ subroutine s_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
call psb_info(ctxt, iam, np)
if (sv%ipar(1) == amg_local_solver_ ) then
call psb_init(ctxt1,np=1,basectxt=ctxt,ids=(/iam/))
icomm = psb_get_mpi_comm(ctxt1)
allocate(sv%local_ctxt,stat=info)
sv%local_ctxt = ctxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm,sv%local_ctxt
call psb_info(ctxt1, me, np)
npr = np
else if (sv%ipar(1) == amg_global_solver_ ) then
icomm = psb_get_mpi_comm(ctxt)
!write(*,*)iam,'mumps_bld: global +++++>',icomm,ctxt
call psb_info(ctxt, iam, np)
me = iam
npr = np
call psb_init(ctxt1,basectxt=ctxt)
else
info = psb_err_internal_error_
call psb_errpush(info,name,&
& a_err='Invalid local/global solver in MUMPS')
goto 9999
end if
!allocate(sv%local_ctxt,stat=info)
icomm = ctxt1%get_mpic()
sv%local_ctxt = ctxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm!,sv%local_ctxt%ctxt
call psb_info(ctxt1, me, npr)
npc = 1
if (debug_level >= psb_debug_outer_) &
& write(debug_unit,*) me,' ',trim(name),' start'
@@ -115,8 +110,8 @@ subroutine s_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
end if
end if
call mpi_comm_dup(icomm,sv%id%comm,info)
sv%id%comm = icomm
sv%id%job = -1
sv%id%par = 1
if (sv%ipar(3) == 2) then
amgprec/impl/solver/amg_z_mumps_solver_bld.F90
+7 -12
View File
@@ -80,24 +80,19 @@ subroutine z_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
call psb_info(ctxt, iam, np)
if (sv%ipar(1) == amg_local_solver_ ) then
call psb_init(ctxt1,np=1,basectxt=ctxt,ids=(/iam/))
icomm = psb_get_mpi_comm(ctxt1)
allocate(sv%local_ctxt,stat=info)
sv%local_ctxt = ctxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm,sv%local_ctxt
call psb_info(ctxt1, me, np)
npr = np
else if (sv%ipar(1) == amg_global_solver_ ) then
icomm = psb_get_mpi_comm(ctxt)
!write(*,*)iam,'mumps_bld: global +++++>',icomm,ctxt
call psb_info(ctxt, iam, np)
me = iam
npr = np
call psb_init(ctxt1,basectxt=ctxt)
else
info = psb_err_internal_error_
call psb_errpush(info,name,&
& a_err='Invalid local/global solver in MUMPS')
goto 9999
end if
!allocate(sv%local_ctxt,stat=info)
icomm = ctxt1%get_mpic()
sv%local_ctxt = ctxt1
!write(*,*)iam,'mumps_bld: local +++++>',icomm!,sv%local_ctxt%ctxt
call psb_info(ctxt1, me, npr)
npc = 1
if (debug_level >= psb_debug_outer_) &
& write(debug_unit,*) me,' ',trim(name),' start'
@@ -115,8 +110,8 @@ subroutine z_mumps_solver_bld(a,desc_a,sv,info,b,amold,vmold,imold)
end if
end if
call mpi_comm_dup(icomm,sv%id%comm,info)
sv%id%comm = icomm
sv%id%job = -1
sv%id%par = 1
if (sv%ipar(3) == 2) then