Hi Matthew, That is ok. May I ask you one more question?
Sorry. I just start to learn how to use DMPLex. If I have an msh file including the triangular elements (t3). After I use DMPlex to read this mesh, how can I know the element type (actually which is t3) by using the functions inside DMPlex? Best, Jian ________________________________ From: Matthew Knepley <[email protected]> Sent: 21 August 2019 14:24:56 To: Jian Zhang - 3ME Cc: Jed Brown; [email protected] Subject: Re: [petsc-users] Getting the connectivity from DMPlex On Wed, Aug 21, 2019 at 8:21 AM Jian Zhang - 3ME <[email protected]<mailto:[email protected]>> wrote: Hi Matthew, Thanks for your answer. I just looked at the pestc4py 3.11.0 (I think it is the newest one). Actually, I did not find DMPlexGetClosure function in the DMPlex.pyx, also not in DM.pyx. Sorry, I messed up the name in my mail. Its getTransitiveClosure() https://bitbucket.org/petsc/petsc4py/src/9598d8fb2bb0baa79f99ebf76991ddecd1fd173e/src/PETSc/DMPlex.pyx#lines-332 Thanks, Matt Best, Jian ________________________________ From: Matthew Knepley <[email protected]<mailto:[email protected]>> Sent: 21 August 2019 14:03:19 To: Jian Zhang - 3ME Cc: Jed Brown; [email protected]<mailto:[email protected]> Subject: Re: [petsc-users] Getting the connectivity from DMPlex On Wed, Aug 21, 2019 at 7:48 AM Jian Zhang - 3ME via petsc-users <[email protected]<mailto:[email protected]>> wrote: Hi Jed, Thank you very much. I tried to use DMPlexGetCone, but the output is the edge ids, not the vertice ids. For the function DMPlexGetClosureIndices, I can not find it in petsc4py. Do you know any alternative ways to solve this? DMPlexGetClosure() should be in petsc4py. It also returns orientations, which you can just filter out. Thanks, Matt Best, Jian ________________________________ From: Jed Brown <[email protected]<mailto:[email protected]>> Sent: 21 August 2019 13:42:41 To: Jian Zhang - 3ME; [email protected]<mailto:[email protected]> Subject: Re: [petsc-users] Getting the connectivity from DMPlex Jian Zhang - 3ME via petsc-users <[email protected]<mailto:[email protected]>> writes: > Hi guys, > > I am trying to get the element connectivity from DMPlex. The input is the > element id, and the output should be the vertice ids. Which function should I > use to achieve this? Thanks in advance. See DMPlexGetCone or DMPlexGetClosureIndices. -- 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.cse.buffalo.edu/~knepley/> -- 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.cse.buffalo.edu/~knepley/>
