[petsc-dev] coordinates stored twice in snes ex12?

Matthew Knepley knepley at gmail.com
Tue Nov 19 15:14:41 CST 2013


On Tue, Nov 19, 2013 at 1:14 PM, Geoffrey Irving <irving at naml.us> wrote:

> It looks like when ex12 clones the DM, it reuses the mesh data
> structure but copies the coordinate vector.  Is there a reason it
> can't reuse that storage?  Given that the PetscFEM routines accept the
> same position array for dependent and auxiliary fields, does it even
> make sense to have distinct coordinate arrays?


Yes, this is lazy. It could be replaced by

DMGetCoordinateSection(dm, &s)
DMGetCoordinatesLocal(dm, &c)
DMSetCoordinateSection(dmAux, s)
DMSetCoordinatesLocal(dmAux, c)

    Matt


>
> Geoffrey
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20131119/b7930e19/attachment.html>


More information about the petsc-dev mailing list