Re: [ublas] matrix expressions (was find2 required by Matrix)

From: Shawn D. Pautz (sdpautz_at_[hidden])
Date: 2005-10-19 11:57:29

I think I understand now, thanks.  But now that "Matrix" no longer 
appears in the template machinery, should I write

E LU = A;

rather than

Matrix LU = A;

Perhaps an equivalent question: Is E always a concrete matrix class?

Shawn

>
>uBLAS has a much better way to write this kind of function. You should use...
> template <class E>
> typename matrix_expression<E>::value_type determinant (const matrix_expression<E>& A)
> {
>
>This will match any Matrix Expression. In this case you will need to create a temporary
>  Matrix LU = A;
>
>  
>
-- 
Shawn D. Pautz
Radiation Transport Department, MS 1179
Sandia National Laboratories
P.O. Box 5800
Albuquerque, NM 87185
Phone: [phone]
Fax: [phone]
E-Mail: sdpautz_at_[hidden]