代码如下:
program MKL_axpy
use blas95
implicit none
integer, parameter :: n = 10
real*8 :: x(n), y(n)
real*8 :: a = 2.d0
call random_seed()
call random_number(x)
write(*,*) 'x is...'
write(*,*) x
call random_number(y)
write(*,*) 'y is...'
write(*,*) y
!// compute y:= a*x + y
call axpy(x, y, a)
write(*,*) 'a*x + y is...'
write(*,*) y
end program MKL_axpy
结果如下:
x is...
0.711440635133745 0.239022631345746 0.531877574142810
0.408197494082518 0.205041923294106 0.193456430660466
0.551417066655331 8.433182591954472E-002 0.104865774006392
0.691408261549519
y is...
0.421106025015009 0.550881903537065 2.059467218376093E-002
0.765671392475285 0.243057506000570 0.894994875450880
0.851445092527582 0.620371592106103 0.818479128913359
0.326524560691130
a*x + y is...
1.84398729528250 1.02892716622856 1.08434982046938
1.58206638064032 0.653141352588782 1.28190773677181
1.95427922583824 0.789035243945193 1.02821067692614
1.70934108379017