← Back to team overview

dolfin team mailing list archive

Re: LAPACK wrapper

 



On 22/05/10 10:48, Anders Logg wrote:
On Sat, May 22, 2010 at 09:28:55AM +0100, Garth N. Wells wrote:
Is anyone using the LAPACK linear algebra wrappers? Since they are
not subclasses of GenericFoo, I suggest that we remove them and
users use Armadillo (http://arma.sourceforge.net/) instead when
wanting to access LAPACK/ATLAS functions. Armadillo is very easy to
use. Any objections to removing the LAPACK wrappers?

We have a Blueprint for wrapping Armadillo, but this would only be
worthwhile if we want to assemble into an Armadillo-based object.
Since Armadillo will most commonly be used for small matrices when
performance is critical, and expression templates are one of its
very attractive features, I don't see much point in wrapping it.

Garth

I agree it would be better to wrap Armadillo (or use it directly), but
the LAPACK wrappers are used heavily in the Extrapolation class so
they can't be removed until that has been ported to use Armadillo.


Turned out that only a few lines needed to be changed.

Non-const integers are passed by reference as few times in Extrapolation. Is it the intention that they change and are then used?

Garth

--
Anders



Follow ups

References