From d4fc54e75aceb74d8133502f6bbcfae11eab0d9f Mon Sep 17 00:00:00 2001 From: cirdans-home Date: Fri, 22 Nov 2019 09:46:53 +0100 Subject: [PATCH] Corrected error in inv_vect implementation --- base/psblas/Makefile | 3 +-- base/psblas/psb_cinv_vect.f90 | 4 ++-- base/psblas/psb_dinv_vect.f90 | 4 ++-- base/psblas/psb_sinv_vect.f90 | 4 ++-- base/psblas/psb_zinv_vect.f90 | 4 ++-- 5 files changed, 9 insertions(+), 10 deletions(-) diff --git a/base/psblas/Makefile b/base/psblas/Makefile index 781f454a..29d78acc 100644 --- a/base/psblas/Makefile +++ b/base/psblas/Makefile @@ -12,8 +12,7 @@ OBJS= psb_ddot.o psb_damax.o psb_dasum.o psb_daxpby.o\ psb_cnrm2.o psb_cnrmi.o psb_cspmm.o psb_cspsm.o \ psb_cmlt_vect.o psb_dmlt_vect.o psb_zmlt_vect.o psb_smlt_vect.o\ psb_cdiv_vect.o psb_ddiv_vect.o psb_zdiv_vect.o psb_sdiv_vect.o\ - psb_cinv_vect.o psb_dinv_vect.o psb_zinv_vect.o psb_sinv_vect.o - + psb_cinv_vect.o psb_dinv_vect.o psb_zinv_vect.o psb_sinv_vect.o\ LIBDIR=.. INCDIR=.. diff --git a/base/psblas/psb_cinv_vect.f90 b/base/psblas/psb_cinv_vect.f90 index 1df47904..f7285e0c 100644 --- a/base/psblas/psb_cinv_vect.f90 +++ b/base/psblas/psb_cinv_vect.f90 @@ -92,7 +92,7 @@ subroutine psb_cinv_vect(x,y,desc_a,info) end if if(desc_a%get_local_rows() > 0) then - call y%div(x,info) + call y%inv(x,info) end if call psb_erractionrestore(err_act) @@ -166,7 +166,7 @@ subroutine psb_cinv_vect_check(x,y,desc_a,info,flag) end if if(desc_a%get_local_rows() > 0) then - call y%div(x,info,flag) + call y%inv(x,info,flag) end if call psb_erractionrestore(err_act) diff --git a/base/psblas/psb_dinv_vect.f90 b/base/psblas/psb_dinv_vect.f90 index 0ce370b4..787f4cab 100644 --- a/base/psblas/psb_dinv_vect.f90 +++ b/base/psblas/psb_dinv_vect.f90 @@ -92,7 +92,7 @@ subroutine psb_dinv_vect(x,y,desc_a,info) end if if(desc_a%get_local_rows() > 0) then - call y%div(x,info) + call y%inv(x,info) end if call psb_erractionrestore(err_act) @@ -166,7 +166,7 @@ subroutine psb_dinv_vect_check(x,y,desc_a,info,flag) end if if(desc_a%get_local_rows() > 0) then - call y%div(x,info,flag) + call y%inv(x,info,flag) end if call psb_erractionrestore(err_act) diff --git a/base/psblas/psb_sinv_vect.f90 b/base/psblas/psb_sinv_vect.f90 index 6925d066..0be53a70 100644 --- a/base/psblas/psb_sinv_vect.f90 +++ b/base/psblas/psb_sinv_vect.f90 @@ -92,7 +92,7 @@ subroutine psb_sinv_vect(x,y,desc_a,info) end if if(desc_a%get_local_rows() > 0) then - call y%div(x,info) + call y%inv(x,info) end if call psb_erractionrestore(err_act) @@ -166,7 +166,7 @@ subroutine psb_sinv_vect_check(x,y,desc_a,info,flag) end if if(desc_a%get_local_rows() > 0) then - call y%div(x,info,flag) + call y%inv(x,info,flag) end if call psb_erractionrestore(err_act) diff --git a/base/psblas/psb_zinv_vect.f90 b/base/psblas/psb_zinv_vect.f90 index 7b36f6b0..4e61379c 100644 --- a/base/psblas/psb_zinv_vect.f90 +++ b/base/psblas/psb_zinv_vect.f90 @@ -92,7 +92,7 @@ subroutine psb_zinv_vect(x,y,desc_a,info) end if if(desc_a%get_local_rows() > 0) then - call y%div(x,info) + call y%inv(x,info) end if call psb_erractionrestore(err_act) @@ -166,7 +166,7 @@ subroutine psb_zinv_vect_check(x,y,desc_a,info,flag) end if if(desc_a%get_local_rows() > 0) then - call y%div(x,info,flag) + call y%inv(x,info,flag) end if call psb_erractionrestore(err_act)