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!