Re: [Help-gsl] Inversion of a symmetric matrix

[Brian Gough]

Filip Westarp writes:
 > 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?

The procedure sounds reasonable (check the reference manual for actual
usage of the functions ) .... if you still get incorrect (unsymmetric)
results please submit a bug report to address@hidden with an example
program that shows the problem.

Thanks,

-- 
Brian Gough

Network Theory Ltd
15 Royal Park
Bristol BS8 3AL
United Kingdom

Tel: +44 (0)117 3179309
Fax: +44 (0)117 9048108
Web: http://www.network-theory.co.uk/