fdurastante/psblas3
1
0
Fork 0
mirror of https://github.com/sfilippone/psblas3.git synced 2026-10-15 18:39:10 +00:00
Code Issues Packages Projects Releases Wiki Activity

cosmetic fixes

Browse Source
This commit is contained in:
Salvatore Filippone
2020-11-02 21:27:34 +01:00
parent 288e3cd341
commit 4a274e43ff
1 changed files with 89 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files
base/serial/psi_d_serial_impl.f90
+89 -5
View File
@@ -42,7 +42,7 @@ subroutine psi_daxpby(m,n,alpha, x, beta, y, info)
integer(psb_ipk_) :: err_act
integer(psb_ipk_) :: lx, ly
integer(psb_ipk_) :: ierr(5)
character(len=20) :: name, ch_err
character(len=20) :: name, ch_err
name='psb_geaxpby'
info=psb_success_
@@ -78,8 +78,9 @@ subroutine psi_daxpby(m,n,alpha, x, beta, y, info)
goto 9999
end if
if ((m>0).and.(n>0)) call daxpby(m,n,alpha,x,lx,beta,y,ly,info)
if (m>0) call daxpby(m,ione,alpha,x,lx,beta,y,ly,info)
call psb_erractionrestore(err_act)
return
@@ -99,9 +100,9 @@ subroutine psi_daxpbyv(m,alpha, x, beta, y, info)
real(psb_dpk_), intent (in) :: alpha, beta
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_) :: err_act
integer(psb_ipk_) :: lx, ly
integer(psb_ipk_) :: lx, ly, i
integer(psb_ipk_) :: ierr(5)
character(len=20) :: name, ch_err
character(len=20) :: name, ch_err
name='psb_geaxpby'
info=psb_success_
@@ -131,7 +132,90 @@ subroutine psi_daxpbyv(m,alpha, x, beta, y, info)
goto 9999
end if
if (m>0) call daxpby(m,ione,alpha,x,lx,beta,y,ly,info)
if (alpha.eq.@XZERO@) then
if (beta.eq.@XZERO@) then
do i=1,m
y(i) = @XZERO@
enddo
else if (beta.eq.@XONE@) then
!
! Do nothing!
!
else if (beta.eq.-@XONE@) then
do i=1,m
y(i) = - y(i)
enddo
else
do i=1,m
y(i) = beta*y(i)
enddo
endif
else if (alpha.eq.@XONE@) then
if (beta.eq.@XZERO@) then
do i=1,m
y(i) = x(i)
enddo
else if (beta.eq.@XONE@) then
do i=1,m
y(i) = x(i) + y(i)
enddo
else if (beta.eq.-@XONE@) then
do i=1,m
y(i) = x(i) - y(i)
enddo
else
do i=1,m
y(i) = x(i) + beta*y(i)
enddo
endif
else if (alpha.eq.-@XONE@) then
if (beta.eq.@XZERO@) then
do i=1,m
y(i) = -x(i)
enddo
else if (beta.eq.@XONE@) then
do i=1,m
y(i) = -x(i) + y(i)
enddo
else if (beta.eq.-@XONE@) then
do i=1,m
y(i) = -x(i) - y(i)
enddo
else
do i=1,m
y(i) = -x(i) + beta*y(i)
enddo
endif
else
if (beta.eq.@XZERO@) then
do i=1,m
y(i) = alpha*x(i)
enddo
else if (beta.eq.@XONE@) then
do i=1,m
y(i) = alpha*x(i) + y(i)
enddo
else if (beta.eq.-@XONE@) then
do i=1,m
y(i) = alpha*x(i) - y(i)
enddo
else
do i=1,m
y(i) = alpha*x(i) + beta*y(i)
enddo
endif
endif
call psb_erractionrestore(err_act)
return