Corrected error in inv_vect implementation

merge-paraggr-newops
cirdans-home 5 years ago
parent c1f3b2d9d4
commit d4fc54e75a

@ -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=..

@ -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)

@ -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)

@ -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)

@ -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)

Loading…
Cancel
Save