[Help-gsl] Inversion of a symmetric matrix

[Filip Westarp]

Hi there,

I'm using gsl to compute the inverse of a symmetric matrix - first calling gsl_linalg_LU_decomp() to get the LU decomposition and then gsl_linalg_LU_invert() to get the matrix inverse. 
However, the resulting matrix isn't symmetric, as I expect it to be. 
Am I doing something wrong - or have I misunderstood something?

Thanks for your help,
Filip