Will there be a way to set coordinates if I am using a homemade mesh generator?
What about having coordinates for the edges/faces/cells as well (to be interpreted as we want)? On Wed, Sep 26, 2012 at 3:34 PM, Matthew Knepley <knepley at gmail.com> wrote: > On Wed, Sep 26, 2012 at 5:24 PM, Chris Eldred <chris.eldred at gmail.com> > wrote: >> >> On the subject: >> >> How do coordinates work for user-defined DMComplex objects? What >> determines the Coordinate section (and which points does it contain?)? >> How is the coordinate vector set? > > > I am just rewriting this. Now it will be > > DMGetCoordinateDM(DM, DM *) > DMGetCoordinates(DM, Vec *) > DMGetCoordinatesLocal(DM, Vec *) > > and there is a convenience > > DMComplexGetCoordinateSection(DM, PetscSection *) > > By default, I read in the vertex coordinates from the mesh generator. > However, > the idea is that you can interpret the coordinates however you want. We just > give you a DM and a vector. You can construct the section however you want, > although without vertex coordinates some default things will not work, like > VTK > output. > > Matt > >> >> >> -- >> Chris Eldred >> DOE Computational Science Graduate Fellow >> Graduate Student, Atmospheric Science, Colorado State University >> B.S. Applied Computational Physics, Carnegie Mellon University, 2009 >> chris.eldred at gmail.com > > > > > -- > 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 -- Chris Eldred DOE Computational Science Graduate Fellow Graduate Student, Atmospheric Science, Colorado State University B.S. Applied Computational Physics, Carnegie Mellon University, 2009 chris.eldred at gmail.com
