diff --git a/mlprec/mld_d_sludist_solver.F90 b/mlprec/mld_d_sludist_solver.F90 index 264c287a..db97e483 100644 --- a/mlprec/mld_d_sludist_solver.F90 +++ b/mlprec/mld_d_sludist_solver.F90 @@ -96,7 +96,7 @@ module mld_d_sludist_solver & bind(c,name='mld_dsludist_solve') result(info) use iso_c_binding integer(c_int) :: info - integer(c_int), value :: itrans,n,ldb + integer(c_int), value :: itrans,n,nrhs,ldb real(c_double) :: b(ldb,*) type(c_ptr), value :: lufactors end function mld_dsludist_solve @@ -218,10 +218,10 @@ contains info = psb_success_ - write(0,*) 'SLUDIST INTERFACE IS CURRENTLY BROKEN. TO BE FIXED' - info=psb_err_internal_error_ - call psb_errpush(info,name) - goto 9999 +!!$ write(0,*) 'SLUDIST INTERFACE IS CURRENTLY BROKEN. TO BE FIXED' +!!$ info=psb_err_internal_error_ +!!$ call psb_errpush(info,name) +!!$ goto 9999 call x%v%sync() call y%v%sync() call sv%apply(alpha,x%v%v,beta,y%v%v,desc_data,trans,work,info) @@ -275,10 +275,10 @@ contains if (debug_level >= psb_debug_outer_) & & write(debug_unit,*) me,' ',trim(name),' start' - write(0,*) 'SLUDIST INTERFACE IS CURRENTLY BROKEN. TO BE FIXED' - info=psb_err_internal_error_ - call psb_errpush(info,name) - goto 9999 +!!$ write(0,*) 'SLUDIST INTERFACE IS CURRENTLY BROKEN. TO BE FIXED' +!!$ info=psb_err_internal_error_ +!!$ call psb_errpush(info,name) +!!$ goto 9999 if (psb_toupper(upd) == 'F') then @@ -298,6 +298,7 @@ contains call psb_loc_to_glob(acsr%ja(1:nztota),desc_a,info,iact='I') acsr%ja(:) = acsr%ja(:) - 1 acsr%irp(:) = acsr%irp(:) - 1 + write(0,*) 'ACSR ',maxval(acsr%ja),minval(acsr%ja),nrow_a,nztota ifrst = ifrst - 1 info = mld_dsludist_fact(nglob,nrow_a,nztota,ifrst,& & acsr%val,acsr%irp,acsr%ja,sv%lufactors,&