diff --git a/mlprec/impl/mld_c_bld_mlhier_array.f90 b/mlprec/impl/mld_c_bld_mlhier_array.f90 index 5cdabd24..da0d2279 100644 --- a/mlprec/impl/mld_c_bld_mlhier_array.f90 +++ b/mlprec/impl/mld_c_bld_mlhier_array.f90 @@ -139,7 +139,7 @@ subroutine mld_c_bld_mlhier_array(nplevs,a,desc_a,precv,info) ! precv(1)%base_a => a precv(1)%base_desc => desc_a - + newsz = 0 array_build_loop: do i=2, iszv ! ! Check on the iprcparm contents: they should be the same diff --git a/mlprec/impl/mld_caggrmat_smth_asb.f90 b/mlprec/impl/mld_caggrmat_smth_asb.f90 index e40665ab..a54f5484 100644 --- a/mlprec/impl/mld_caggrmat_smth_asb.f90 +++ b/mlprec/impl/mld_caggrmat_smth_asb.f90 @@ -404,7 +404,7 @@ subroutine mld_caggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest & 'starting spspmm 3' call psb_spspmm(op_restr,am3,ac,info) 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 call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3') goto 9999 diff --git a/mlprec/impl/mld_d_bld_mlhier_array.f90 b/mlprec/impl/mld_d_bld_mlhier_array.f90 index 73002567..3618de1b 100644 --- a/mlprec/impl/mld_d_bld_mlhier_array.f90 +++ b/mlprec/impl/mld_d_bld_mlhier_array.f90 @@ -139,7 +139,7 @@ subroutine mld_d_bld_mlhier_array(nplevs,a,desc_a,precv,info) ! precv(1)%base_a => a precv(1)%base_desc => desc_a - + newsz = 0 array_build_loop: do i=2, iszv ! ! Check on the iprcparm contents: they should be the same diff --git a/mlprec/impl/mld_daggrmat_smth_asb.f90 b/mlprec/impl/mld_daggrmat_smth_asb.f90 index e86f7244..01ae8661 100644 --- a/mlprec/impl/mld_daggrmat_smth_asb.f90 +++ b/mlprec/impl/mld_daggrmat_smth_asb.f90 @@ -404,7 +404,7 @@ subroutine mld_daggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest & 'starting spspmm 3' call psb_spspmm(op_restr,am3,ac,info) 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 call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3') goto 9999 diff --git a/mlprec/impl/mld_s_bld_mlhier_array.f90 b/mlprec/impl/mld_s_bld_mlhier_array.f90 index 9a516028..3f249d34 100644 --- a/mlprec/impl/mld_s_bld_mlhier_array.f90 +++ b/mlprec/impl/mld_s_bld_mlhier_array.f90 @@ -139,7 +139,7 @@ subroutine mld_s_bld_mlhier_array(nplevs,a,desc_a,precv,info) ! precv(1)%base_a => a precv(1)%base_desc => desc_a - + newsz = 0 array_build_loop: do i=2, iszv ! ! Check on the iprcparm contents: they should be the same diff --git a/mlprec/impl/mld_saggrmat_smth_asb.f90 b/mlprec/impl/mld_saggrmat_smth_asb.f90 index 2719e166..595a1580 100644 --- a/mlprec/impl/mld_saggrmat_smth_asb.f90 +++ b/mlprec/impl/mld_saggrmat_smth_asb.f90 @@ -404,7 +404,7 @@ subroutine mld_saggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest & 'starting spspmm 3' call psb_spspmm(op_restr,am3,ac,info) 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 call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3') goto 9999 diff --git a/mlprec/impl/mld_z_bld_mlhier_array.f90 b/mlprec/impl/mld_z_bld_mlhier_array.f90 index b8a20df4..206cdf30 100644 --- a/mlprec/impl/mld_z_bld_mlhier_array.f90 +++ b/mlprec/impl/mld_z_bld_mlhier_array.f90 @@ -139,7 +139,7 @@ subroutine mld_z_bld_mlhier_array(nplevs,a,desc_a,precv,info) ! precv(1)%base_a => a precv(1)%base_desc => desc_a - + newsz = 0 array_build_loop: do i=2, iszv ! ! Check on the iprcparm contents: they should be the same diff --git a/mlprec/impl/mld_zaggrmat_smth_asb.f90 b/mlprec/impl/mld_zaggrmat_smth_asb.f90 index 30544b9c..9e1065c3 100644 --- a/mlprec/impl/mld_zaggrmat_smth_asb.f90 +++ b/mlprec/impl/mld_zaggrmat_smth_asb.f90 @@ -404,7 +404,7 @@ subroutine mld_zaggrmat_smth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_rest & 'starting spspmm 3' call psb_spspmm(op_restr,am3,ac,info) 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 call psb_errpush(psb_err_internal_error_,name,a_err='Build ac = op_restr x am3') goto 9999