Dear people of GSL. I'm trying to adapt the example given in chapter 14 (eigensystems) of the GSL manual, on the diagonalisation of a Hilbert matrix, to the diagonalisation of a complex hermitian matrix. I can see the obvious part (changing functions "gsl_eigen_symmv" to "gsl_eigen_hermv"). Indeed, my only problem is with the first bit.
Specificaly, which would be the equivalent of array data[ ] for a matrix of complex values? I currently have my real and imaginary parts stored as two separate arrays of doubles. How can I give them the proper format to run with a code analogous to the mentioned example? What function shall I use instead of gsl_matrix_view_array()? I will really apreciate your help. Thanks in advance, Raymond _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
