mld2p4-2:

tests/newslv/mld_d_tlu_solver.f90
 tests/newslv/mld_d_tlu_solver_impl.f90

Fix NEWSLV sample dir.
stopcriterion
Salvatore Filippone 12 years ago
parent c0507d8eb2
commit 71dd991b9b

@ -107,7 +107,7 @@ module mld_d_tlu_solver
import :: psb_desc_type, mld_d_tlu_solver_type, psb_d_vect_type, psb_dpk_, & import :: psb_desc_type, mld_d_tlu_solver_type, psb_d_vect_type, psb_dpk_, &
& psb_dspmat_type, psb_d_base_sparse_mat, psb_d_base_vect_type & psb_dspmat_type, psb_d_base_sparse_mat, psb_d_base_vect_type
type(psb_desc_type), intent(in) :: desc_data type(psb_desc_type), intent(in) :: desc_data
class(mld_d_tlu_solver_type), intent(in) :: sv class(mld_d_tlu_solver_type), intent(inout) :: sv
real(psb_dpk_),intent(inout) :: x(:) real(psb_dpk_),intent(inout) :: x(:)
real(psb_dpk_),intent(inout) :: y(:) real(psb_dpk_),intent(inout) :: y(:)
real(psb_dpk_),intent(in) :: alpha,beta real(psb_dpk_),intent(in) :: alpha,beta
@ -253,15 +253,14 @@ contains
info = psb_success_ info = psb_success_
call psb_erractionsave(err_act) call psb_erractionsave(err_act)
ival = sv%stringval(val)
call mld_stringval(val,ival,info) if (ival >=0) then
if (info == psb_success_) call sv%set(what,ival,info) call sv%set(what,ival,info)
if (info /= psb_success_) then
info = psb_err_from_subroutine_
call psb_errpush(info, name)
goto 9999
end if end if
if (info /= psb_success_) goto 9999
call psb_erractionrestore(err_act) call psb_erractionrestore(err_act)
return return

@ -189,7 +189,7 @@ subroutine mld_d_tlu_solver_apply(alpha,sv,x,beta,y,desc_data,trans,work,info)
use mld_d_tlu_solver, mld_protect_name => mld_d_tlu_solver_apply use mld_d_tlu_solver, mld_protect_name => mld_d_tlu_solver_apply
implicit none implicit none
type(psb_desc_type), intent(in) :: desc_data type(psb_desc_type), intent(in) :: desc_data
class(mld_d_tlu_solver_type), intent(in) :: sv class(mld_d_tlu_solver_type), intent(inout) :: sv
real(psb_dpk_),intent(inout) :: x(:) real(psb_dpk_),intent(inout) :: x(:)
real(psb_dpk_),intent(inout) :: y(:) real(psb_dpk_),intent(inout) :: y(:)
real(psb_dpk_),intent(in) :: alpha,beta real(psb_dpk_),intent(in) :: alpha,beta

Loading…
Cancel
Save