Re: [Boost-users] How to do svd for a complex matrix using ublas &lapack bindings?

Subject: Re: [Boost-users] How to do svd for a complex matrix using ublas &lapack bindings?
From: Thomas Klimpel (Thomas.Klimpel_at_[hidden])
Date: 2008-12-25 07:39:52

Hi Allen,

Allen Chen wrote:
> Today I begin to learn to use boost ublas.

There is also a mailing list especially for boost ublas.

> And I have successfully done svd for a real matrix.
> To do svd for a complex matrix, I simply change codes
> from double to std::complex like the following.
> But it does not work.

The singular values "S" are still real, even for a complex matrix. So you should write "vector<double> S(2);" instead of "vector<std::complex<double>> S(2);".

Regards,
Thomas