fdurastante
/
amg4psblas
mirror of https://github.com/sfilippone/amg4psblas.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Emergency fix for silly interface mistake with MUMPS.

Browse Source
stopcriterion
Salvatore Filippone 6 years ago
parent 0b57535762
commit b0d1b09021
4 changed files with 10 additions and 10 deletions
Show Stats Download Patch File Download Diff File

6
mlprec/impl/solver/mld_c_mumps_solver_bld.F90
Unescape Escape View File

@ -104,7 +104,7 @@
allocate(sv%id,stat=info) allocate(sv%id,stat=info)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_alloc_dealloc_ info=psb_err_alloc_dealloc_
call psb_errpush(info,name,a_err='mld_dmumps_default') call psb_errpush(info,name,a_err='mld_cmumps_default')
goto 9999 goto 9999
end if end if
end if end if
@ -114,7 +114,7 @@
sv%id%job = -1 sv%id%job = -1
sv%id%par = 1 sv%id%par = 1
call cmumps(sv%id) call cmumps(sv%id)
!WARNING: CALLING dMUMPS WITH JOB=-1 DESTROY THE SETTING OF DEFAULT:TO FIX !WARNING: CALLING cmumps WITH JOB=-1 DESTROY THE SETTING OF DEFAULT:TO FIX
if (allocated(sv%icntl)) then if (allocated(sv%icntl)) then
do i=1,mld_mumps_icntl_size do i=1,mld_mumps_icntl_size
if (allocated(sv%icntl(i)%item)) then if (allocated(sv%icntl(i)%item)) then
@ -169,7 +169,7 @@
end if end if
!call psb_barrier(ictxt) !call psb_barrier(ictxt)
write(*,*)iam, ' calling mumps N,nz,nz_loc',sv%id%n,sv%id%nnz,sv%id%nnz_loc write(*,*)iam, ' calling mumps N,nz,nz_loc',sv%id%n,sv%id%nnz,sv%id%nnz_loc
call dmumps(sv%id) call cmumps(sv%id)
!call psb_barrier(ictxt) !call psb_barrier(ictxt)
info = sv%id%infog(1) info = sv%id%infog(1)
if (info /= psb_success_) then if (info /= psb_success_) then

2
mlprec/impl/solver/mld_d_mumps_solver_bld.F90
Unescape Escape View File

@ -114,7 +114,7 @@
sv%id%job = -1 sv%id%job = -1
sv%id%par = 1 sv%id%par = 1
call dmumps(sv%id) call dmumps(sv%id)
!WARNING: CALLING dMUMPS WITH JOB=-1 DESTROY THE SETTING OF DEFAULT:TO FIX !WARNING: CALLING dmumps WITH JOB=-1 DESTROY THE SETTING OF DEFAULT:TO FIX
if (allocated(sv%icntl)) then if (allocated(sv%icntl)) then
do i=1,mld_mumps_icntl_size do i=1,mld_mumps_icntl_size
if (allocated(sv%icntl(i)%item)) then if (allocated(sv%icntl(i)%item)) then

6
mlprec/impl/solver/mld_s_mumps_solver_bld.F90
Unescape Escape View File

@ -104,7 +104,7 @@
allocate(sv%id,stat=info) allocate(sv%id,stat=info)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_alloc_dealloc_ info=psb_err_alloc_dealloc_
call psb_errpush(info,name,a_err='mld_dmumps_default') call psb_errpush(info,name,a_err='mld_smumps_default')
goto 9999 goto 9999
end if end if
end if end if
@ -114,7 +114,7 @@
sv%id%job = -1 sv%id%job = -1
sv%id%par = 1 sv%id%par = 1
call smumps(sv%id) call smumps(sv%id)
!WARNING: CALLING dMUMPS WITH JOB=-1 DESTROY THE SETTING OF DEFAULT:TO FIX !WARNING: CALLING smumps WITH JOB=-1 DESTROY THE SETTING OF DEFAULT:TO FIX
if (allocated(sv%icntl)) then if (allocated(sv%icntl)) then
do i=1,mld_mumps_icntl_size do i=1,mld_mumps_icntl_size
if (allocated(sv%icntl(i)%item)) then if (allocated(sv%icntl(i)%item)) then
@ -169,7 +169,7 @@
end if end if
!call psb_barrier(ictxt) !call psb_barrier(ictxt)
write(*,*)iam, ' calling mumps N,nz,nz_loc',sv%id%n,sv%id%nnz,sv%id%nnz_loc write(*,*)iam, ' calling mumps N,nz,nz_loc',sv%id%n,sv%id%nnz,sv%id%nnz_loc
call dmumps(sv%id) call smumps(sv%id)
!call psb_barrier(ictxt) !call psb_barrier(ictxt)
info = sv%id%infog(1) info = sv%id%infog(1)
if (info /= psb_success_) then if (info /= psb_success_) then

6
mlprec/impl/solver/mld_z_mumps_solver_bld.F90
Unescape Escape View File

@ -104,7 +104,7 @@
allocate(sv%id,stat=info) allocate(sv%id,stat=info)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_alloc_dealloc_ info=psb_err_alloc_dealloc_
call psb_errpush(info,name,a_err='mld_dmumps_default') call psb_errpush(info,name,a_err='mld_zmumps_default')
goto 9999 goto 9999
end if end if
end if end if
@ -114,7 +114,7 @@
sv%id%job = -1 sv%id%job = -1
sv%id%par = 1 sv%id%par = 1
call zmumps(sv%id) call zmumps(sv%id)
!WARNING: CALLING dMUMPS WITH JOB=-1 DESTROY THE SETTING OF DEFAULT:TO FIX !WARNING: CALLING zmumps WITH JOB=-1 DESTROY THE SETTING OF DEFAULT:TO FIX
if (allocated(sv%icntl)) then if (allocated(sv%icntl)) then
do i=1,mld_mumps_icntl_size do i=1,mld_mumps_icntl_size
if (allocated(sv%icntl(i)%item)) then if (allocated(sv%icntl(i)%item)) then
@ -169,7 +169,7 @@
end if end if
!call psb_barrier(ictxt) !call psb_barrier(ictxt)
write(*,*)iam, ' calling mumps N,nz,nz_loc',sv%id%n,sv%id%nnz,sv%id%nnz_loc write(*,*)iam, ' calling mumps N,nz,nz_loc',sv%id%n,sv%id%nnz,sv%id%nnz_loc
call dmumps(sv%id) call zmumps(sv%id)
!call psb_barrier(ictxt) !call psb_barrier(ictxt)
info = sv%id%infog(1) info = sv%id%infog(1)
if (info /= psb_success_) then if (info /= psb_success_) then

Write Preview
Loading…
Cancel
Save