<div class="gmail_quote">On Sun, Nov 20, 2011 at 07:00, Hui Zhang <span dir="ltr"><<a href="mailto:mike.hui.zhang@hotmail.com">mike.hui.zhang@hotmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Doest petsc has built-in *cartesian* mesh generator for finite elements?<br>
What I want is: inputs the size of a domain, outputs the geometric coordinates<br>
like (0.1, 0.2, 0.3) for each nodes, volume elements like (1,2,...,8),<br>
face elements like (1,2,...,4).<br></blockquote><div><br></div><div>You can use DMDA and then loop over the local part of the domain however you like. See DMDASetUniformCoordinates() and also finite element examples like src/snes/examples/tutorials/ex48.c or src/ksp/ksp/examples/tutorials/ex43.c or ex42.c. The function DMDAGetElements() exists, but I think it is of limited utility.</div>
</div>