... the letters stand for. I'm sure I left a couple out, which is something my compiler will not likely miss. Given that I've herded the dependencies into this directory, how do I best build this module? Netlib has a nice interface for fetching the driver routine you need plus all dependencies - e.g. "dsyev plus dependencies" from http://www.netlib.org/lapack/double/

...> be hoping to take a eal matrix and use lapack to give me the Hessenberg form, which, I can then use for bigger and better things. Hi, Lapack has both driver routines and expert routines. If you just want to get started with for example eigenvalue systems then drivers like DGEEV or DSYEV are a good choice. They will do the transformations to triangular form and other stuff for ...

... to write the caller in f90 or f95 but certainly in free form. sgehrd.f looks like fixed form f77. My first question is how do I deal with calling an f77 routine with main compiled in contemporary gfortran? As others explained, that's no problem. Though LAPACK95 is nice in that it provides explicit interfaces, that among other things helps protect the programmer from his own ...

...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

... library? I've had zero luck building large ANYTHING on Windows from source, so I add pieces until it works. I hadn't noticed that there's a half dozen calls to other routines, not all of which are in the same folder. EXTERNAL SAXPY, SGEHD2, SGEMM, SLAHR2, SLARFB, STRMM, $ XERBLA We'll see if it cascades much farther. I suppose I could put ...

....  I would like to write the caller in f90 or f95 but certainly in free form.  sgehrd.f looks like fixed form f77. My first question is how do I deal with calling an f77 routine with main compiled in contemporary gfortran? gfortran (like g95) will compile both types of source from the same command line: C:\Users\epc\temp>type main.f90 k = 1 print *, k call sub(k) print *, k...

...ILO, IHI, A, LDA, TAU, WORK, LWORK, INFO ) * * -- LAPACK routine (version 3.1) -- * Univ. of Tennessee, Univ. of California Berkeley...LWORK = -1, then a workspace query is assumed; the routine * only calculates the optimal size of the WORK array, returns ... question is how do I deal with calling an f77 routine with main compiled in contemporary gfortran? -- A good politician ...

... am, m...@iae.nl (Marcel Hendrix) wrote: "dkel...@hotmail.com" <dkel...@hotmail.com> writes Re: Sorting routine On Sep 8, 10:14 am, m...@iae.nl (Marcel Hendrix) wrote: "dkel...@hotmail.com" <dkel...@hotmail.com> writes Re: Sorting routine [..] I don't understand it. Not only isn't there a difference but the results are almost ...

"dkelvey@hotmail.com" <dkelvey@hotmail.com> writes Re: Sorting routine On Sep 8, 10:14 am, m...@iae.nl (Marcel Hendrix) wrote: "dkel...@hotmail.com" <dkel...@hotmail.com> writes Re: Sorting routine On Sep 7, 8:25= am, m...@iae.nl (Marcel Hendrix) wrote: ---snip--- 8 =: nbits PRIVATE ( 16 is 261 ms, 4 ...

"dkelvey@hotmail.com" <dkelvey@hotmail.com> writes Re: Sorting routine On Sep 8, 10:14 am, m...@iae.nl (Marcel Hendrix) wrote: "dkel...@hotmail.com" <dkel...@hotmail.com> writes Re: Sorting routine [..] I don't understand it. Not only isn't there a difference but the results are almost identical. I've looked at your code and I...