Re: [Help-gsl] Invert a complex matrix ?

[Brian Gough]

Mathieu Lenoble writes:
 > I have a problem to find how invert a complex matrix! I read the
 > How To of the GSL and i followed some programs found on the web but
 > it doesn't still working !!  So, to invert my complex matrix i use
 > the gsl_linalg_complex_LU_decomp(...)  following by
 > gsl_linalg_complex_LU_invert(....) functions.
 > 
 > and here is the results that i obtained :
 > 
 > esinv(0,0)=5.20715e-141+1.06685e-194*j esinv(0,1)=1.17094e-321+0*j
 > esinv(1,0)=2.10994e305+2.06841e149*j esinv(1,1)=1+0*j
 > 

Your program looks ok. It works for me:

bjg|debian> ./a.out

esinv(0,0)=0.197106+-0.30952*j 
esinv(0,1)=0.295506+-0.185834*j
esinv(1,0)=0.317136+0.274791*j 
esinv(1,1)=-0.0654989+-0.268088*j

If you still have problems please send a detailed bug report to
address@hidden with all the information about your system and version
of GSL, etc.

-- 
Brian Gough