I prefer to compile LAPACK and create a library. I then, create a
module with interfaces to LAPACK subroutines. For instance I have
created a function for the eigenvalues and eigenvectors of a real or a
real symmetric matrix. So when I want to calculate, say, the
eigenvectors of a real matrix A I type v=eigVectors(A), or if I want
the eigenvectors of a symmetric matrix A I type
v=eigVectors(A,sym="Y"). In this way the cumbersome argument list of
the LAPACK routines does not appear, or have to be used anywhere.

Michael