Dear Mohammad,

As a user of PETSc for the last 8 years, since my days as an
undergraduate, and now as a professional staff scientist at a
supercomputing center, I can say with some confidence that there are
no codes like PETSc in C++ or any other language in terms of quality
of implementation, documentation, and support.  Can you tell us a
little more about your current implementation?  It is true that PETSc
does not support multiple types in the same build, but you do get your
choice of floating-point values and real or complex types.  Also,
PETSc has several C++ components within it, and one of the supported
ways of building it is in 'C++' mode, see -c-language in the configure
options.

If you insist on departing us (we'll miss you), I suggest you look at
Sandia's Trilinos package: http://trilinos.sandia.gov/

Good Luck,
Aron

On Sat, Oct 30, 2010 at 7:54 AM, Mohammad Mirzadeh
<m.mirzadeh@engineering.ucsb.edu> wrote:
> Dear all,
> I was wondering if anyone can refer me to a package similar to PETSc but
> that is written in c++? Right now I have a large code written in c++ for
> doing CFD simulations that I need to transform from serial to parallel.
> Initially I was thinking of PETSc and tried using it but found that PETSc is
> written in C and thus does not allow to have arrays (in parallel) of
> arbitrary type. I have a big data structure and it is much easier for me to
> retain the current structure and form of the code. As a result I was
> wondering if you guys know of any similar package in C++ ?(in the sense that
> it can provide with efficient linear solvers in parallel while hiding most
> of MPI from the user)
> I could think of HYPRE but then again I am not sure it is written in C++.
> Thanks,
> Mohammad