I have noticed from other posts that better C++ support for LAPACK routines is in the works. (Thank goodness! The F2C option may produce a working program, but it’s not very nice. For example, it replaces FORTRAN goto statements with C goto statements whereas I’d prefer to have hand-written code that removes the goto statements completely).
I am curious to see how 1D and 2D dynamic arrays are implemented in the proposed C++ programs. Is the STL <vector> library being used as is, or are classes defined? If so, could the code for the class definitions be viewed?