[petsc-users] specifying vertex coordinates using DMPlexCreateFromCellListParallel

[Adrian Croucher]

hi

> That would great if it allows elements with mixed faces,

Yes, as I've already mentioned to Matt, I need this too for 6-node wedge elements (probably the same thing as your 'prisms').


- Adrian
>
> Thank you
> -Hassan
>
> From: Matthew Knepley [mailto:knepley at gmail.com]
> Sent: Thursday, June 29, 2017 10:36 AM
> To: Hassan Raiesi <Hassan.Raiesi at aero.bombardier.com>
> Cc: petsc-users at mcs.anl.gov
> Subject: Re: [petsc-users] specifying vertex coordinates using DMPlexCreateFromCellListParallel
>
> On Thu, Jun 29, 2017 at 9:32 AM, Hassan Raiesi <Hassan.Raiesi at aero.bombardier.com<mailto:Hassan.Raiesi at aero.bombardier.com>> wrote:
> Ok,
>
> The problem is that the faceSizes varies for pyramids and prisms, they have faces with variable number of nodes (i.e for pyramids, four triangular faces and one quad face, this I don?t see how to specify?
> What if I report all faces as quad faces, for triangular faces, there would be a repeated vertex, would it make any problem?
>
> No, you are right, that will not work. The whole interpolation code would have to be rewritten to allow faces with
> variable size. Its not hard, but it is time-consuming. Since I am moving this summer, this will not be done in a timely
> manner unfortunately.