Regarding<br>&quot;...Maybe PETSc is an overkill for you, but if the problem is
that simple, why not use Fortran and regular LAPACK/PLAPACK?,,,&quot;<br>I haven&#39;t written fortran since IV and 77 in the early 80&#39;s. The matrix multiplication part is simple but the ambient application is squirrelly and I&quot;m more facile in C or C++. So I just want basic linear algebra library I can link my C to, but maybe all y&#39;all use Matlab for simple things.<br>
Probably I should just do it with, say, Armadillo, and learn something without spilling too much blood from fingertips.<br>Peter<br>
<br><br><div class="gmail_quote">On Mon, Nov 22, 2010 at 12:38 PM, Gus Correa <span dir="ltr">&lt;<a href="mailto:gus@ldeo.columbia.edu">gus@ldeo.columbia.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div class="im">Peter St. John wrote:<br>
&gt; Gus,<br>
&gt; In my case, I&quot;m not concerned with PDEs but the &quot;linear system solvers&quot;<br>
&gt; themselves,which PETSci gets from other packages?<br>
</div>Hi Peter<br>
<br>
Yes, you can compile PETSC with support of a variety of linear algebra<br>
packages.  See their installation instructions and documentation.<br>
I think there is support for PLAPACK, please check their<br>
install and docs pages:<br>
<br>
<a href="http://www.mcs.anl.gov/petsc/petsc-2/documentation/installation.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-2/documentation/installation.html</a><br>
<a href="http://www.mcs.anl.gov/petsc/petsc-2/documentation/index.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-2/documentation/index.html</a><br>
<div class="im"><br>
&gt; So for example,<br>
&gt; looking over their table at:<br>
&gt;  <a href="http://www.mcs.anl.gov/petsc/petsc-as/documentation/linearsolvertable.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/documentation/linearsolvertable.html</a><br>
&gt;  I&#39;d want the &quot;dense&quot; case and therefore PLAPACK, so now I know to<br>
&gt; google PLAPACK.<br>
&gt; I just want to solve a system of simultaneous linear equations within<br>
&gt; the vernacular of C/C++.<br>
<br>
</div>AFAIK, PETSc has C, Fortran, and C++ interfaces.<br>
If I remember right, you can choose between C++ or C when<br>
you compile it.<br>
You can choose between real (for most linear algebra )<br>
or complex (a must for FFTs, and FFTW support) scalars.<br>
You can also choose which<br>
linear algebra packages you will use (which you can compile separately,<br>
or let PETSc handle - it handles many of them),<br>
which BLAS/LAPACK to use (say, Goto BLAS,<br>
Intel MKL, AMD ACML, or let PETSc install its own version).<br>
<br>
I am not a PETSc user, but some people here used PETSc to solve inverse<br>
and  forward problems, problems with sparse and dense matrices, etc.<br>
As it often happens when science (not computer science)<br>
is the main goal, once the &quot;it works&quot; stage is reached,<br>
the code starts to be used to produce papers,<br>
more and more features are added,<br>
and it is never again improved for efficiency,<br>
never goes beyond the prototype phase in this regard.<br>
<br>
In terms of efficiency, the results didn&#39;t shine,<br>
in terms of science/results output they were very good.<br>
However, it may well be that the way the problem<br>
(not PETSc) was programmed led to the inefficient code.<br>
<br>
PETSc experts in the list:<br>
For the benefit of Peter, Ivan, and the list,<br>
would you kindly jump in and say something about PETSc<br>
efficiency and other features?<br>
Txs.<br>
<br>
Maybe PETSc is an overkill for you, but if the problem is<br>
that simple, why not use Fortran and regular LAPACK/PLAPACK?<br>
<br>
My two cents.<br>
Gus Correa<br>
<div class="im"><br>
&gt; Thanks,<br>
&gt; Peter<br>
&gt;<br>
&gt; On Mon, Nov 22, 2010 at 11:42 AM, Gus Correa &lt;<a href="mailto:gus@ldeo.columbia.edu">gus@ldeo.columbia.edu</a><br>
</div><div class="im">&gt; &lt;mailto:<a href="mailto:gus@ldeo.columbia.edu">gus@ldeo.columbia.edu</a>&gt;&gt; wrote:<br>
&gt;<br>
&gt;     Would PETSc be of use for you?<br>
&gt;<br>
&gt;      <a href="http://www.mcs.anl.gov/petsc/petsc-as/" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/</a><br>
&gt;<br>
&gt;     Gus Correa<br>
&gt;<br>
&gt;     Ivan Marin wrote:<br>
&gt;      &gt; I&#39;m interested also in the bindings for C++ and LAPACK/ScaLAPACK. I<br>
&gt;      &gt; did some tests in the past for simple stuff, and worked, but very<br>
&gt;      &gt; simple linear system solvers.<br>
&gt;      &gt;<br>
&gt;      &gt; Ivan Marin<br>
&gt;      &gt;<br>
&gt;      &gt; Civil Engineering Dept<br>
&gt;      &gt; University of Minnesota<br>
&gt;      &gt; 500 Pillsbury Dr. S.E.<br>
&gt;      &gt; Minneapolis, MN 55455<br>
&gt;      &gt;<br>
&gt;      &gt; Laboratório de Hidráulica Computacional - LHC<br>
&gt;      &gt; Departamento de Hidráulica e Saneamento - SHS<br>
&gt;      &gt; Escola de Engenharia de São Carlos - EESC<br>
&gt;      &gt; Universidade de São Paulo - USP<br>
&gt;      &gt;<br>
&gt;      &gt; <a href="http://albatroz.shs.eesc.usp.br" target="_blank">http://albatroz.shs.eesc.usp.br</a><br>
&gt;      &gt; +55 16 3373 8270<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; 2010/11/17 Peter St. John &lt;<a href="mailto:peter.st.john@gmail.com">peter.st.john@gmail.com</a><br>
</div>&gt;     &lt;mailto:<a href="mailto:peter.st.john@gmail.com">peter.st.john@gmail.com</a>&gt;&gt;:<br>
<div class="im">&gt;      &gt;&gt; I&#39;m going to dabble a bit with writing some C/C++ with a linear<br>
&gt;     algebra<br>
&gt;      &gt;&gt; library. My fortran is rusty (although I did some work<br>
&gt;     translating fortran<br>
&gt;      &gt;&gt; to C in this millennium :-).<br>
&gt;      &gt;&gt; The wiki suggests: LAPACK++, clapack, and Armadillo. Would this<br>
&gt;     choice have<br>
&gt;      &gt;&gt; any real bearing on compiling against MPI later?<br>
&gt;      &gt;&gt; Just if anyone is happy with any combination (LAPACK++ and<br>
&gt;     OpenMPI, say) I&#39;d<br>
&gt;      &gt;&gt; be glad to hear it.<br>
&gt;      &gt;&gt;<br>
&gt;      &gt;&gt; I just installed (k)ubuntu, which seems fine for my purposes so<br>
&gt;     far. I do<br>
&gt;      &gt;&gt; miss fvwm, and the default &quot;single click to invoke an<br>
&gt;     application&quot; is<br>
&gt;      &gt;&gt; totally nuts but I fixed it.<br>
&gt;      &gt;&gt; Thanks,<br>
&gt;      &gt;&gt; Peter<br>
&gt;      &gt;&gt;<br>
&gt;      &gt;&gt; _______________________________________________<br>
&gt;      &gt;&gt; Beowulf mailing list, <a href="mailto:Beowulf@beowulf.org">Beowulf@beowulf.org</a><br>
</div>&gt;     &lt;mailto:<a href="mailto:Beowulf@beowulf.org">Beowulf@beowulf.org</a>&gt; sponsored by Penguin Computing<br>
<div class="im">&gt;      &gt;&gt; To change your subscription (digest mode or unsubscribe) visit<br>
&gt;      &gt;&gt; <a href="http://www.beowulf.org/mailman/listinfo/beowulf" target="_blank">http://www.beowulf.org/mailman/listinfo/beowulf</a><br>
&gt;      &gt;&gt;<br>
&gt;      &gt;&gt;<br>
&gt;      &gt; _______________________________________________<br>
&gt;      &gt; Beowulf mailing list, <a href="mailto:Beowulf@beowulf.org">Beowulf@beowulf.org</a><br>
</div>&gt;     &lt;mailto:<a href="mailto:Beowulf@beowulf.org">Beowulf@beowulf.org</a>&gt; sponsored by Penguin Computing<br>
<div class="im">&gt;      &gt; To change your subscription (digest mode or unsubscribe) visit<br>
&gt;     <a href="http://www.beowulf.org/mailman/listinfo/beowulf" target="_blank">http://www.beowulf.org/mailman/listinfo/beowulf</a><br>
&gt;<br>
&gt;     _______________________________________________<br>
&gt;     Beowulf mailing list, <a href="mailto:Beowulf@beowulf.org">Beowulf@beowulf.org</a><br>
</div>&gt;     &lt;mailto:<a href="mailto:Beowulf@beowulf.org">Beowulf@beowulf.org</a>&gt; sponsored by Penguin Computing<br>
<div><div></div><div class="h5">&gt;     To change your subscription (digest mode or unsubscribe) visit<br>
&gt;     <a href="http://www.beowulf.org/mailman/listinfo/beowulf" target="_blank">http://www.beowulf.org/mailman/listinfo/beowulf</a><br>
&gt;<br>
&gt;<br>
<br>
_______________________________________________<br>
Beowulf mailing list, <a href="mailto:Beowulf@beowulf.org">Beowulf@beowulf.org</a> sponsored by Penguin Computing<br>
To change your subscription (digest mode or unsubscribe) visit <a href="http://www.beowulf.org/mailman/listinfo/beowulf" target="_blank">http://www.beowulf.org/mailman/listinfo/beowulf</a><br>
</div></div></blockquote></div><br>