Global memory breaking BLAS/LAPACK routines?

[w0ss4g3]

Hi,

I've got a code in fortran where I've used the variable N as shared memory - its the order of a polynomial approximation and appears in pretty much every routine I write.

Would this assignment have an effect on the routines which use N as a variable in the LAPACK/BLAS libraries?

For example DGEMV:
SUBROUTINE DGEMV(TRANS,M,N,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)

One call is using: CALL DGEMV(trans, mat_rows, mat_cols, alpha, tempmat, mat_rows, tempvec_in, 1, beta, tempvec_out, 1)

Is my value, mat_cols being passed to the subroutine, or is the globally valued N being used instead?

Thanks for any help - I'm sure its a dumb question, but I'm at a loose end as to why my solver isn't agreeing with a known solution!

[Julien Langou]

Having N as a global variable in your program should not be a problem. N is input only in the BLAS interface.
The problem might be in your code. If you strongly suspect the BLAS DGEMM call, then link with the reference BLAS library which is provided within LAPACK (which you can get from http://www.netlib.org/lapack/lapack.tgz ).
(It is slower than optimized BLAS library but safer. So you can assess problems like this.)
Then if the problem is still there and you still strongly suspect the BLAS DGEMM call, simply remove the BLAS DGEMM call from your code and have three hand-made loops for the MM and start debugging from there.
I do not think the problem comes from the BLAS.
Julien.