Update example_blas_gemv.f90

memory aliasing error
Euler-37 · Apr 13, 2024 · 4e60fb2 · 4e60fb2
1 parent c26253a
commit 4e60fb2
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/example/linalg/example_blas_gemv.f90 b/example/linalg/example_blas_gemv.f90
@@ -2,13 +2,13 @@ program example_gemv
   use stdlib_linalg, only: eye
   use stdlib_linalg_blas, only: sp,gemv
   implicit none(type,external)
-  real(sp) :: A(2, 2), B(2)
+  real(sp) :: A(2, 2), B(2), C(2)
   B = [1.0,2.0]
   A = eye(2)
 
   ! Use legacy BLAS interface 
-  call gemv('No transpose',m=size(A,1),n=size(A,2),alpha=1.0,a=A,lda=size(A,1),x=B,incx=1,beta=0.0,y=B,incy=1)
+  call gemv('No transpose',m=size(A,1),n=size(A,2),alpha=1.0,a=A,lda=size(A,1),x=B,incx=1,beta=0.0,y=C,incy=1)
 
-  print *, B ! returns 1.0 2.0
+  print *, C ! returns 1.0 2.0
 
 end program example_gemv