<div class="gmail_quote">On Wed, Jan 4, 2012 at 13:18, TAY wee-beng <span dir="ltr"><<a href="mailto:zonexo@gmail.com">zonexo@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":6cy">So the 1st step should be checking the load balancing. If it's more or less balanced, will slicing it in 3 directions further improve the speed?<br></div></blockquote><div><br></div><div>You want some combination of balancing and small surface area. Slicing in 3 directions usually improves this. Depending on the anisotropy in the physics, it could be better or worse for solver convergence rates.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":6cy">
<br>
Another thing is that I hope to do some form of adaptive mesh refinement.<br></div></blockquote><div><br></div><div>On a Cartesian mesh? With what sort of discretization. There are lots of packages for this, each one targeting some class of discretizations and problems.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":6cy">
<br>
I'm a bit confused. Are partitioning software like ParMETIS, Zoltan or Isorropia also used for adaptive mesh refinement?<br>
<br>
Or which open source software can do that with PETSc and in Fortran? I searched and got libMesh, for use with PETSc and paramesh, which is in Fortran.</div></blockquote></div><br>