[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