[Lapackpp-devel] Appreciation and suggestion

[Aosheng R. <ro...@gm...>]

Many thanks go to Dr.Stimming for adding the inverse of complex-value matrix
to lapackpp.  As Dr. Stimming said that  many  times  we could  ignore  the
inverse  of  a matrix  and   re-arrange algorithms for the same  purpose.
But for some  cases, we are lazy to do so and rely on a direct inverse of
the matrix.

LAPACKPP is definitely one of the most useful packages  for engineering
applications.

Right now,  I have a random thought that may increase the significant
efforts of Dr. Stimming. Also I am not sure if such a thought makes sense to
the other people.  My thought is Just for discussion.

Thought:  I wonder If the COMPLEX type in LAPACKPP  is  identical with the
complex<double> in GNU C++.   Right now,  the Complex type in LAPACKPP is
the structure - type, instead of class - type, as in C++.  Within the
LAPACKPP class, Class-Type COMPLEX in C++ is converted to the Structure-Type
COMPLEX, compatible with  Fortran  LAPACK. The direct benifit from such a
definition of Class-Type COMPLEX is that one can do operations on the
product and division using  complex<double> type in GNU C++,  once results
are obtained from the LAPACKPP. With the present definition of COMPLEX, we
have to do term-by-term calculation based  on the decomposition of real and
imaginary parts.

I appreciate the wonderful job of Dr. Stimming who brings us  an easy to
implement  interface  of  LAPACK  with  C++.

Regards,
-------------
Rong