[Help-gsl] Re: eigenvectors of non symmetric matrix?

[David Doria]

Ah, it seems I do in fact have 1.8 because I was using a version ported to
windows because I spend half the day now developing some windows stuff so
its just easier to stay in it.  Do you think I can just steal the
gsl_eigen.h files from the new release ? or will that break everything?

On Thu, Mar 6, 2008 at 5:08 PM, David Doria <address@hidden> wrote:

> I am taking an outer product:
>
> a b^T
> where a and b are column vectors.  Then I want the eigen values and
> vectors of the resulting matrix (called mat3).
>
> I tried to use:
>     gsl_eigen_symmv_workspace * EigenWorkspace = gsl_eigen_symmv_alloc
> (2);
>     gsl_eigen_symmv (mat3, EigenValues, EigenVectors, EigenWorkspace);
>
> but it gave the wrong results.  I guess this is because it was expecting a
> symmetric matrix? Is the only other choice to use:
>     gsl_eigen_hermv_workspace * EigenWorkspace = gsl_eigen_hermv_alloc
> (2);
>     gsl_eigen_hermv (mat3, EigenValues, EigenVectors, EigenWorkspace);
>
> but for that, I'd have to first make mat3 a complex matrix (or so says the
> error haha)?
>
> Please let me know.
>
> --
> Thanks,
>
> David




-- 
Thanks,

David