[petsc-dev] DMComplex and Coordinates

[Chris Eldred]

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