Vector 在哪里可以找到BLAS示例代码（Fortran）？
我一直在寻找有关blas的适当文档，我发现了315页的密集材料，而ctrl-f不适用于这些材料。它提供了关于例程采用什么输入参数的所有信息，但是有很多输入参数，我真的可以使用一些示例代码。我找不到任何人。我知道必须有一些人，否则没有人能够使用这些库 具体来说，我使用MacOSX 10.5.8上通过macports安装的ATLAS，使用gcc 4.4中的gfortran（也通过macports安装）。我正在用Fortran 90编写代码。我对Fortran还是相当陌生，但我对mathematica、matlab、perl和shell脚本有相当多的经验 我希望能够初始化稠密复向量，并将其与稠密对称（但不是厄米）复矩阵相乘。矩阵的元素是通过指数的数学函数定义的——称之为f（i，j） 有人能提供一些代码或一些代码的链接吗？首先，你可以看到你正在寻找的例程是zgemv，这里是一个复杂的（'z'）矩阵（'m'）向量（'v'）乘法 如果向量只是普通数组，即它们在内存中是连续的，那么INCX和INCY参数只是1。就LDA参数而言，只需保持它等于矩阵大小即可。其他参数很简单。例如：Vector 在哪里可以找到BLAS示例代码（Fortran）？,vector,matrix,fortran,linear-algebra,blas,Vector,Matrix,Fortran,Linear Algebra,Blas,我一直在寻找有关blas的适当文档，我发现了315页的密集材料，而ctrl-f不适用于这些材料。它提供了关于例程采用什么输入参数的所有信息，但是有很多输入参数，我真的可以使用一些示例代码。我找不到任何人。我知道必须有一些人，否则没有人能够使用这些库 具体来说，我使用MacOSX 10.5.8上通过macports安装的ATLAS，使用gcc 4.4中的gfortran（也通过macports安装）。我正在用Fortran 90编写代码。我对Fortran还是相当陌生，但我对mathematica
implicit none
integer, parameter :: N=2
complex*16, parameter :: imag1 = cmplx(0.d0, 1.d0)
complex*16 :: a(N,N), x(N), y(N)
complex*16 :: alpha, beta
a(:,:)=imag1;
x(:)=1.d0
y(:)=0.d0
alpha=1.d0; beta=0.d0
call zgemv('N',N,N,alpha,a,N,x,1,beta,y,1)
print*, y
end
zsymv implicit none
integer, parameter :: N=2
complex*16, parameter :: imag1 = cmplx(0.d0, 1.d0)
complex*16 :: a(N,N), x(N), y(N)
complex*16 :: alpha, beta
a(:,:)=imag1;
x(:)=1.d0
y(:)=0.d0
alpha=1.d0; beta=0.d0
call zgemv('N',N,N,alpha,a,N,x,1,beta,y,1)
print*, y
end
编辑：上面的代码没有使用矩阵对称的事实。如果需要，请查找
zsymvzsymvztrmvzsymvzsymvztrmvzsymv