mlprec/impl/mld_c_bld_mlhier_array.f90
 mlprec/impl/mld_caggrmat_smth_asb.f90
 mlprec/impl/mld_d_bld_mlhier_array.f90
 mlprec/impl/mld_daggrmat_smth_asb.f90
 mlprec/impl/mld_s_bld_mlhier_array.f90
 mlprec/impl/mld_saggrmat_smth_asb.f90
 mlprec/impl/mld_z_bld_mlhier_array.f90
 mlprec/impl/mld_zaggrmat_smth_asb.f90

Initialize newsz in mlhier_array.
Why was ac in 'coo'?
stopcriterion
Salvatore Filippone 9 years ago
parent 0749aec572
commit a4b0a64647

@ -139,7 +139,7 @@ subroutine mld_c_bld_mlhier_array(nplevs,a,desc_a,precv,info)
! !
precv(1)%base_a => a precv(1)%base_a => a
precv(1)%base_desc => desc_a precv(1)%base_desc => desc_a
newsz = 0
array_build_loop: do i=2, iszv array_build_loop: do i=2, iszv
! !
! Check on the iprcparm contents: they should be the same ! Check on the iprcparm contents: they should be the same

@ -404,7 +404,7 @@ subroutine mld_caggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest
& 'starting spspmm 3' & 'starting spspmm 3'
call psb_spspmm(op_restr,am3,ac,info) call psb_spspmm(op_restr,am3,ac,info)
if (info == psb_success_) call am3%free() if (info == psb_success_) call am3%free()
if (info == psb_success_) call ac%cscnv(info,type='coo',dupl=psb_dupl_add_) if (info == psb_success_) call ac%cscnv(info,type='csr',dupl=psb_dupl_add_)
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3') call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3')
goto 9999 goto 9999

@ -139,7 +139,7 @@ subroutine mld_d_bld_mlhier_array(nplevs,a,desc_a,precv,info)
! !
precv(1)%base_a => a precv(1)%base_a => a
precv(1)%base_desc => desc_a precv(1)%base_desc => desc_a
newsz = 0
array_build_loop: do i=2, iszv array_build_loop: do i=2, iszv
! !
! Check on the iprcparm contents: they should be the same ! Check on the iprcparm contents: they should be the same

@ -404,7 +404,7 @@ subroutine mld_daggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest
& 'starting spspmm 3' & 'starting spspmm 3'
call psb_spspmm(op_restr,am3,ac,info) call psb_spspmm(op_restr,am3,ac,info)
if (info == psb_success_) call am3%free() if (info == psb_success_) call am3%free()
if (info == psb_success_) call ac%cscnv(info,type='coo',dupl=psb_dupl_add_) if (info == psb_success_) call ac%cscnv(info,type='csr',dupl=psb_dupl_add_)
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3') call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3')
goto 9999 goto 9999

@ -139,7 +139,7 @@ subroutine mld_s_bld_mlhier_array(nplevs,a,desc_a,precv,info)
! !
precv(1)%base_a => a precv(1)%base_a => a
precv(1)%base_desc => desc_a precv(1)%base_desc => desc_a
newsz = 0
array_build_loop: do i=2, iszv array_build_loop: do i=2, iszv
! !
! Check on the iprcparm contents: they should be the same ! Check on the iprcparm contents: they should be the same

@ -404,7 +404,7 @@ subroutine mld_saggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest
& 'starting spspmm 3' & 'starting spspmm 3'
call psb_spspmm(op_restr,am3,ac,info) call psb_spspmm(op_restr,am3,ac,info)
if (info == psb_success_) call am3%free() if (info == psb_success_) call am3%free()
if (info == psb_success_) call ac%cscnv(info,type='coo',dupl=psb_dupl_add_) if (info == psb_success_) call ac%cscnv(info,type='csr',dupl=psb_dupl_add_)
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3') call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3')
goto 9999 goto 9999

@ -139,7 +139,7 @@ subroutine mld_z_bld_mlhier_array(nplevs,a,desc_a,precv,info)
! !
precv(1)%base_a => a precv(1)%base_a => a
precv(1)%base_desc => desc_a precv(1)%base_desc => desc_a
newsz = 0
array_build_loop: do i=2, iszv array_build_loop: do i=2, iszv
! !
! Check on the iprcparm contents: they should be the same ! Check on the iprcparm contents: they should be the same

@ -404,7 +404,7 @@ subroutine mld_zaggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest
& 'starting spspmm 3' & 'starting spspmm 3'
call psb_spspmm(op_restr,am3,ac,info) call psb_spspmm(op_restr,am3,ac,info)
if (info == psb_success_) call am3%free() if (info == psb_success_) call am3%free()
if (info == psb_success_) call ac%cscnv(info,type='coo',dupl=psb_dupl_add_) if (info == psb_success_) call ac%cscnv(info,type='csr',dupl=psb_dupl_add_)
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3') call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3')
goto 9999 goto 9999

Loading…
Cancel
Save