Corrected bug in out-of-place axpby

5 years ago · ea6ad3a85f
parent 177bc6ecb3
commit ea6ad3a85f
4 changed files with 4 additions and 4 deletions
--- a/base/psblas/psb_caxpby.f90
+++ b/base/psblas/psb_caxpby.f90
@ -263,7 +263,7 @@ subroutine psb_caxpby_vect_out(alpha, x, beta, y,&

  if(desc_a%get_local_rows() > 0) then
    call z%axpby(desc_a%get_local_rows(),&
-         & alpha,x,beta,z,info)
+         & alpha,x,beta,y,info)
  end if

  call psb_erractionrestore(err_act)
--- a/base/psblas/psb_daxpby.f90
+++ b/base/psblas/psb_daxpby.f90
@ -263,7 +263,7 @@ subroutine psb_daxpby_vect_out(alpha, x, beta, y,&

  if(desc_a%get_local_rows() > 0) then
    call z%axpby(desc_a%get_local_rows(),&
-         & alpha,x,beta,z,info)
+         & alpha,x,beta,y,info)
  end if

  call psb_erractionrestore(err_act)
--- a/base/psblas/psb_saxpby.f90
+++ b/base/psblas/psb_saxpby.f90
@ -263,7 +263,7 @@ subroutine psb_saxpby_vect_out(alpha, x, beta, y,&

  if(desc_a%get_local_rows() > 0) then
    call z%axpby(desc_a%get_local_rows(),&
-         & alpha,x,beta,z,info)
+         & alpha,x,beta,y,info)
  end if

  call psb_erractionrestore(err_act)
--- a/base/psblas/psb_zaxpby.f90
+++ b/base/psblas/psb_zaxpby.f90
@ -263,7 +263,7 @@ subroutine psb_zaxpby_vect_out(alpha, x, beta, y,&

  if(desc_a%get_local_rows() > 0) then
    call z%axpby(desc_a%get_local_rows(),&
-         & alpha,x,beta,z,info)
+         & alpha,x,beta,y,info)
  end if

  call psb_erractionrestore(err_act)