On Fri, Oct 14, 2011 at 10:21 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<br>
  Why have all this stuff?<br>
<br>
  <a href="http://petsc.cs.iit.edu/petsc/petsc-dev/rev/49c718f99df2" target="_blank">http://petsc.cs.iit.edu/petsc/petsc-dev/rev/49c718f99df2</a><br>
<br>
   Why not build into the PETSc C code and structures any (small) additional stuff needed to do all this? Why create another entirely new duplicate structure in python to do it?</blockquote><div><br></div><div>Alright reasons:</div>
<div><br></div><div>1) C sucks for this kind of processing. I use the right language for the right job. If not, we could do everything with C++ templates.</div><div><br></div><div>2) It meshes with the build system, so I can build several versions of the code. You do this with make, which is quaint.</div>
<div><br></div><div>Good reasons:</div><div><br></div><div>3) This is for runs across many versions of the code. I can't even switch MPI_COMM_WORLD sizes in your model.</div><div><br></div><div>    Matt</div><div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><font color="#888888"><br>
    Barry<br></font></blockquote></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>