Thanks for the quick answer. We already have a wrapper class for BLAS/Lapack and it's currently being used for small problems as well as solving the - possibly huge - equation system. We want to use PETSc to solve this problem in parallel and would have hoped to be able to simultaneously drop the wrapper. I guess a mixed approach would be perfect then?
On Fri, Nov 11, 2011 at 3:05 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote: > On Fri, Nov 11, 2011 at 03:35, markus.sons at gmail.com < > markus.sons at gmail.com> wrote: > >> So, what do you recommend? Using PETSc for the large-scale computations >> and some simple Vec and Mat class for small, local stuff? > > > Don't use PETSc Mat/Vec for very small local problems like the 3x3 or 4x4 > matrices. If you are already addicted to templates and overloading, you > might check out a library like Eigen which is competitive with BLAS/Lapack > for some problems on some architectures. > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111111/90fd099b/attachment.htm>
