Re: [ublas] [bindings] [lapack] geev and the type of the input matrix

Subject: Re: [ublas] [bindings] [lapack] geev and the type of the input matrix
From: Thomas Klimpel (Thomas.Klimpel_at_[hidden])
Date: 2010-07-06 05:34:02

Marco Guazzone wrote:
> Just another question (sorry for abusing of your patience)
> 
> I'd like to create an "eigen" function (as a wrapper of lapack::geev)
> which accepts as input matrix both:
> * a const& matrix (since it is not modified anymore), and

What do you mean by "(since it is not modified anymore)"? lapack will definitively overwrite the original matrix during the computation and not restore it at the end.

> * a matrix_expression (in order to allow the user to call say
> eigen(trans(A),...) )

If you want that, you probably have to create a temporary matrix that lapack can modify during the computation.

Regards,
Thomas