Re: Who uses clapack?

From: Tino (tino52_at_yahoo.com)
Date: 12/13/04


Date: 13 Dec 2004 06:23:28 -0800


Bill Shortall wrote:
> "Victor Eijkhout" <see.sig@for.addy> wrote in message
> news:1gokv1s.16zyrkl1ytczkbN%see.sig@for.addy...
> > Tino <tino52@yahoo.com> wrote:
> >
> > > On a slightly related note, is there any effort (or even thought)
being
> > > put into the development of a "standardized" C++ object-oriented
LAPACK
> > > library? Something to wrap around the BLAS and LAPACK, perhaps?
> >
> > The thought has occurred to many people. I mean, passing LDA into a
C++
> > routine is the programming equivalent of turning a crank handle on
the
> > space shuttle. Now, what to do about it....
> >
> > I've started looking at Sidl/Babel to get a more native looking
> > interface in other languages than F77. Not sure to what extent that
will
> > satisfy you.
> >
> > V.
> > --
> > email: lastname at cs utk edu
> > homepage: www cs utk edu tilde lastname
>
> Hi Victor,
>
> For the last few years I have been working on a
> set of C++ classes to do linear algebra. I was trying to
> acheive most of the functionality of Lapack/Blas but
> in a very user friendly fashion. The classes center
> arround a general purpose vector class along with a dense rectangular
Matrix
> class and then some more compilcated classes like a vector of vectors
and
> vector of matrices.
> All of them use operator overloading so you can write
> A = B + C * D where A,B,C, D are matrices and all of the
operators
> work for both real and complex matrices. Akll of the basic Linear
algebra
> functions are supported
> ie solution of equations, decompositions, SVD eigenvectors etc. etc.
> The entire library is less than 5000 lines of code
> and compiles into a static library of ~~ 3 megabytes
> You don't have to compile it but you will need either a
> Miccrosoft VC6 or Linux GCC compiler to use it.
> I have demo version available and if anyone in the group wants to
play
> with it and has a compatible compiler send me a line. Once you've
used it
> you will never go back to LAPACK
> pecos@cminet.net
> regards...Bill

I have also developed a C++ library which acts as nothing more than a
wrapper for the BLAS/LAPACK or Intel MKL libraries. In my case, I
claim only to support a subset of the BLAS/LAPACK functionality (though
it is 100% of the functionality that I use). I also support all of the
different matrix formats (dense, banded, symmetric), though only in
double precision (no complex matrices). The reason I was interested in
a more "official" effort is, that I believe that LAPACK's popularity
and utility is in its generality and I would welcome a sort of
standardization of the C++ interface as there seem to be dozens if not
hundereds of them floating around.

Ryan



Relevant Pages

  • Re: It Hurts When I Do This
    ... Being an old libraries guy (I used to do compiler run-time libraries before working on compilers themselves), I'd have suggested moving this code to a library routine - just pass it the descriptors of the LHS and RHS. ... But adding a library routine isn't without cost - it creates support headaches when users try to link objects with older libraries. ... I was a compiler developer for 15 years. ...
    (comp.lang.fortran)
  • Re: [FAQ] mini-FAQ V1.68 - essential reading for those new to the newsgroup
    ... Free Pascal is a 32/64 Object Pascal Compiler ... existing libraries, and a large set of non-visual classes in the ... The compiler is licensed under the GNU General Public License. ... incomplete support for Linux. ...
    (comp.lang.pascal.borland)
  • Re: STM32 ARM toolset advice?
    ... compilers work better than a compiler vendor. ... abortive attempts to fix them themselves. ... I don't think the open-source libraries have worse quality than the ... Commercial support can be ...
    (comp.arch.embedded)
  • Re: Usage of complex numbers?
    ... > into libraries. ... is the array, which C and C++ don't even support in the base language except ... The C++ committee placed the support ... OpenMP, where you get help from the compiler, than to make a bunch of calls ...
    (comp.lang.fortran)
  • Re: Usage of complex numbers?
    ... > into libraries. ... is the array, which C and C++ don't even support in the base language except ... The C++ committee placed the support ... OpenMP, where you get help from the compiler, than to make a bunch of calls ...
    (comp.programming)