[petsc-users] ex42 question

[Jed Brown]

Sang pham van <pvsang002 at gmail.com> writes:

> Yes, It does converge toward 0-flow when I refine the mesh !
>
> In the code, coordinates are set by
>   ierr =
> DMDASetUniformCoordinates(da_prop,0.0+0.5*dx,1.0-0.5*dx,0.0+0.5*dy,1.0-0.5*dy,0.,0);CHKERRQ(ierr);
> Does it set coordinate for elements center here?
> When I want to try non-uniform grid, should I just modify coordinates
> attached with the DM, or I need more implementation?

You can modify coordinates.  It looks like the code does the correct
coordinate transformations.  (Dave's real code runs on non-affine meshes
for science reasons.)
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140416/6af2bf12/attachment.pgp>