From: Gunter Winkler (guwi17_at_[hidden])
Date: 2005-07-28 16:35:11
Zitat von Dima Sorkin <dsorkin_at_[hidden]>:
> I have a small dense matrix (no more than 10x10),
> it is well conditioned. What can be the most portable
> way to get its inverse ? (portabe in the meaning
> of different versions of ublas).
> I'm not looking for "quick/efficient",
> just numerically stable and portable.
I would try this:
matrix<double> Ainv = identity_matrix<double>( A.size1(), A.size2() );
lu_substitute(Ainv,A); // untested