[Help-gsl] SVD and complex matrix !

[Mathieu Lenoble]

Hi,

I need to find the singular value decomposition of a complex matrix but I
dont find the appropriate fonction in the GSL.
So i was wondering if there is a fonction for that ?

Thanks,

Mathieu.


----- Original Message ----- 
From: "Brian Gough" <address@hidden>
To: "Mathieu Lenoble" <address@hidden>
Cc: <address@hidden>
Sent: Tuesday, March 30, 2004 10:11 PM
Subject: Re: [Help-gsl] Invert a complex matrix ?


> 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
>