<div dir="ltr">On Mon, Aug 4, 2008 at 6:49 AM, Techas <span dir="ltr"><<a href="mailto:techas@gmail.com">techas@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Matthew,<br>
<br>
I'm Sergio (X-FEM), we meet some weeks ago in Davis.<br>
I'm playing a little bit with petsc to evaluate how much work will<br>
take me mount my x-fem code on petsc.</blockquote><div><br>Sorry I am just replying now. I returned from Norway yesterday.<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
My first question is: can you tell me a good example of the assembly<br>
of a global finite element matrix (distributed in n procs) based on a<br>
distributed mesh?</blockquote><div><br>The answers is in two parts. If you have a structured grid, then the DA<br>construct can do everything and is very simple. You just call MatSetValuesStencil()<br>as in KSP ex2.<br><br>
If you have an unstructured grid, we have new functionality (only in the<br>development version) in a Mesh object. This will be in the upcoming release,<br>but you can see it now in src/dm/mesh/sieve/problems/Bratu.hh. The function<br>
to examine is Rhs_Unstructured() which forms the residual for this equation<br>on an unstructured mesh. Unfortunately, this is new and has little documentation.<br>Before the release, I will write some, but until then you will have to ask me<br>
questions if you want to try it out.<br><br>  Thanks,<br><br>     Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>
If you want me to write to the petsc list just tell me.<br>
<br>
thank you!<br>
Sergio.<br>
<br>
--<br>
Sergio Zlotnik, PhD<br>
Group of Dynamics of the Lithosphere -GDL-<br>
Department of Geophysics & Tectonics<br>
Institute of Earth Sciences - CSIC<br>
Sole Sabaris s/n<br>
Barcelona 08028 - SPAIN<br>
<br>
Tel: +34 93 409 54 10<br>
Fax: +34 93 411 00 12<br>
email: <a href="mailto:szlotnik@ija.csic.es">szlotnik@ija.csic.es</a><br>
<br>
Web page <a href="http://www.ija.csic.es/gt/sergioz/" target="_blank">http://www.ija.csic.es/gt/sergioz/</a><br>
</blockquote></div><br><br clear="all"><br>-- <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>

</div>