<div class="gmail_quote">On Fri, Nov 11, 2011 at 03:35, <a href="mailto:markus.sons@gmail.com">markus.sons@gmail.com</a> <span dir="ltr">&lt;<a href="mailto:markus.sons@gmail.com">markus.sons@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">So, what do you recommend? Using PETSc for the large-scale computations and some simple Vec and Mat class for small, local stuff?</blockquote>
</div><br><div>Don&#39;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.</div>