[Trilinos-Users] Trilinos interface to LAPACK

[Jeremy Conlin]

Heidi,
Thank you for your response.  I received another private response also
suggesting I use the SerialDenseMatrix class.  I'm going to investigate this
further.

Thanks,
Jeremy


On Jan 15, 2008 4:15 PM, hkthorn <hkthorn at sandia.gov> wrote:

>  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/64dc1cf4/attachment.html