Boost logo

Ublas :

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?


>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: 505-284-4291
Fax: 505-844-0092
E-Mail: sdpautz_at_[hidden]