Starting with netlib.org/blas you see that the routine you're looking for is zgemv, here netlib.org/blas/zgemv. F it's a complex ('z') matrix ('m') vector ('v') multiply If your vectors are just normal arrays, i.e. They are contiguous in memory, then INCX and INCY arguments are just 1.As far as LDA parameter is concerned, just keep it equal to the matrix size.
Other parameters are straightforward. For example: 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 In general, every time I need a BLAS or LAPACK routine, I look up the parameters on netlib EDIT: the code above doesn't use the fact that your matrix is symmetric. If you want that, then look up the ztrmv routine.
Starting with netlib.org/blas/, you see that the routine you're looking for is zgemv, here netlib.org/blas/zgemv. F --- it's a complex ('z') matrix ('m') vector ('v') multiply. If your vectors are just normal arrays, i.e.
They are contiguous in memory, then INCX and INCY arguments are just 1. As far as LDA parameter is concerned, just keep it equal to the matrix size. Other parameters are straightforward.
For example: 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 In general, every time I need a BLAS or LAPACK routine, I look up the parameters on netlib. EDIT: the code above doesn't use the fact that your matrix is symmetric. If you want that, then look up the ztrmv routine.
I have been searching for decent documentation on blas, and I have found some 315 pages of dense material that ctrl-f does not work on. It provides all the information regarding what input arguments the routines take, but there are a LOT of input arguments and I could really use some example code. I am unable to locate any.
I know there has to be some or no one would be able to use these libraries!
I cant really gove you an answer,but what I can give you is a way to a solution, that is you have to find the anglde that you relate to or peaks your interest. A good paper is one that people get drawn into because it reaches them ln some way.As for me WW11 to me, I think of the holocaust and the effect it had on the survivors, their families and those who stood by and did nothing until it was too late.