[petsc-users] Cell type for DMPlexCreateFromCellList

[Matthew Knepley]

On Fri, Feb 23, 2018 at 12:23 PM, Danyang Su <danyang.su at gmail.com> wrote:

>
> On 18-02-23 03:04 AM, Matthew Knepley wrote:
>
> On Fri, Feb 23, 2018 at 1:33 AM, Danyang Su <danyang.su at gmail.com> wrote:
>
>> Hi All,
>>
>> What cell types does DMPlexCreateFromCellList support? I test this with
>> triangle, tetrahedron and prism. Both triangle and tetrahedron work but
>> prism mesh throws error saying "Cone size 6 not supported for dimension 3".
>>
>> Could anyone tell me all the supported cell types?
>>
>
> The limitation occurs in two places:
>
>   1) Calculating edges and faces: I only know how to do this for tri, tet,
> quad, and hex. Give PETSC_FALSE to CreateFromCellList() and this error will
> go away.
>
> Passing PETSC_FALSE to CreateFromCellList() works. There is no problem in
> the distributed nodes and cells. Is there any plan to add prism(wedge) in
> the near future?
>

 Its not hard. We would just need to add one more function. However, I
would want to make sure that it conformed to what
most people want. I would like it to match the Firedrake definition for
extruded meshes.

>       You could also provide the information for interpolating prisms,
> which would depend on how they are ordered when you read in cells.
>
>   2) Even if you read them in, I have no geometric routines for prisms.
>
> That's fine for me. I just need to create dmplex based on the given
> cell/vertex list and distribute over all the processors. Then follow the
> general routines that have been used in the structured grid.
>

Okay, then that should work now.

  Thanks,

     Matt


> Thanks,
>
> Danyang
>
>
>   Thanks,
>
>      Matt
>
> Thanks,
>>
>> Danyang
>>
>>
>
>
> --
> 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
>
> https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/%7Emk51/>
>
>
>


-- 
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

https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/~mk51/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180223/bb1508b9/attachment.html>