[Trilinos-Users] Trilinos interface to LAPACK

[hkthorn]

Hi Jeremy,

Trilinos provides wrappers to LAPACK functionality in the Epetra and Teuchos
packages.  The Amesos-LAPACK interface is solely for solving linear systems
using DGETRS/DGETRF.  The LAPACK wrappers in Epetra and Teuchos are not all
inclusive, there is not a wrapper for every single LAPACK routine.  The
Epetra interface to LAPACK gives access to some single and double precision
floating point routines.  The Teuchos interface to LAPACK is templated and
gives access to some LAPACK routines for all four native types (float,
double, complex<float>, complex<double>).

At this time the wrappers in Teuchos have more eigenvalue/eigenvector LAPACK
routines and are still being actively enhanced with more LAPACK routines as
requested.  Yes, the interfaces still use pointers instead of a std::vector.
However, in using these wrappers you can create std::vectors for the storage
and just pass in the pointer to the first element of the vector to the
LAPACK wrapper.  
In my experience, a perfect compliment to the LAPACK wrappers in Teuchos is
the SerialDenseMatrix class in Teuchos.  This templated class has been
designed to provide serial dense matrix storage that is compatible with
LAPACK routines.

Let me know if you have any further questions.

Thanks,
Heidi

-- 
Heidi K. Thornquist
Electrical & Microsystems Modeling
Sandia National Laboratories
P.O. Box 5800, MS 0316
Albuquerque, NM  87185-0316
Office:  505-284-8426
Fax:     505-284-5451


On 1/14/08 7:36 PM, "Jeremy Conlin" <jeremit0 at gmail.com> wrote:

> I am trying to use the Trilinos interface to LAPACK and I'm not sure what is
> the best way to do this.  I am primarily using it to calculate eigenvalues and
> eigenvectors, nothing fancy.  First of all there are several packages that
> each have their own interface to LAPACK including: Epetra, Amesos, and
> Teuchos.  When I first learned of Trilinos (I'm new to the scene) I was hoping
> I could pass either a C++ std::vector or some Trilinos objects to the
> functions instead of merely arrays.
> 
> So I guess my question is, what is the best way to interface with LAPACK?  Can
> I pass Trilinos objects (or std::vector) instead of arrays?  Which package
> offers the best interface?  Why are there so many?
> 
> Thanks,
> Jeremy Conlin
> 
> 
> _______________________________________________
> Trilinos-Users mailing list
> Trilinos-Users at software.sandia.gov
> http://software.sandia.gov/mailman/listinfo/trilinos-users


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20080115/8636a377/attachment.html