这篇博客简要介绍使用MKL中的blas计算向量的元素之和
代码如下:
program MKL_asum
use blas95
implicit none
integer, parameter :: n = 10
real*8 :: x(n) = 0, res
call random_seed()
call random_number(x)
write(*,*) 'The vector x is...'
write(*,*) x
res = asum(x)
write(*,*) 'The sum of the vector x is...'
write(*,*) res
end program MKL_asum
结果如下:
The vector x is...
0.398639747353446 0.798604414277417 8.118048864432666E-002
4.452416290741119E-002 8.636616046592745E-002 0.226938333031608
0.153932851778479 0.107662005420434 0.849353654400938
0.339383395317750
The sum of the vector x is...
3.08658521359774