<div dir="ltr"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div class="im"><div>Thanks for the hint. I was wondering whether I should start with a working prototype based on numpy, then do parallelization with PETSc on demand - or start with PETSc wherever matrices and vectors are needed.</div>
</div></div></blockquote><div><br></div><div>You need to write parallel code in a Single Program Multiple Data fashion if you want it to scale well to hundreds or thousands of processors.  This is the best model supported by mpi4py and petsc4py.</div>
<div><br></div><div>A</div></div></div>