[petsc-dev] DMPlexComputeGeometryFEM

[Yann Jobic]

Hi Matt,
Thanks a lot, the new way to do it is clear and efficent.
Yann

Le 11/26/2022 à 10:57 PM, Matthew Knepley a écrit :
> On Fri, Nov 25, 2022 at 7:14 AM Yann Jobic <yann.jobic at univ-amu.fr 
> <mailto:yann.jobic at univ-amu.fr>> wrote:
> 
>     Hi all,
> 
>     I made a very simple example which uses the DMPlexComputeGeometryFEM
>     function.
> 
> 
> You are right. We changed PetscFEGeom, but that broke this function. I 
> have updated your code
> with the new way to get this information. I also removed the function in 
> this MR:
> 
> https://gitlab.com/petsc/petsc/-/merge_requests/5879 
> <https://gitlab.com/petsc/petsc/-/merge_requests/5879>
> 
> I am also rewriting this as we speak in order to make it compatible with 
> libCEED, so it will change
> a little in the next release unfortunately.
> 
>    Thanks,
> 
>       Matt
> 
>     However, i've got an error :
>     ./testComputeGeomFEM -heat_petscspace_degree 0
>         Quadrature of order 1 on 1 points (dim 2)
>     [0]PETSC ERROR: --------------------- Error Message
>     --------------------------------------------------------------
>     [0]PETSC ERROR: Null argument, when expecting valid pointer
>     [0]PETSC ERROR: Null Pointer: Parameter # 7
>     [0]PETSC ERROR: See https://petsc.org/release/faq/
>     <https://petsc.org/release/faq/> for trouble shooting.
>     [0]PETSC ERROR: Petsc Release Version 3.18.1, unknown
>     [0]PETSC ERROR:
>     /home/jobic/projet/fe-utils/petsc/fetools/cmake-build-debug/testComputeGeomFEM
>     on a  named luke by jobic Fri Nov 25 14:08:34 2022
>     [0]PETSC ERROR: Configure options
>     --prefix=/local/lib/petsc/3.18/p0/gcc/nompi --with-mpi=0
>     --with-debugging=1 --with-blacs=1 --download-zlib,--download-p4est
>     --download-triangle=1 --with-single-library=0 --with-large-file-io=1
>     --with-shared-libraries=0 -CFLAGS=" -g -O0" -CXXFLAGS=" -g -O0"
>     -FFLAGS=" -g -O0" PETSC_ARCH=nompi_gcc
>     [0]PETSC ERROR: #1 DMPlexComputeCellGeometryFEM() at
>     /home/devel/src_linux/petsc-3.18.0/src/dm/impls/plex/plexgeometry.c:2158
>     [0]PETSC ERROR: #2 DMPlexComputeGeometryFEM() at
>     /home/devel/src_linux/petsc-3.18.0/src/dm/impls/plex/plexgeometry.c:2524
>     [0]PETSC ERROR: #3 main() at
>     /home/jobic/projet/fe-utils/petsc/fetools/src/testComputeGeomFEM.c:31
>     [0]PETSC ERROR: PETSc Option Table entries:
>     [0]PETSC ERROR: -heat_petscspace_degree 0
>     [0]PETSC ERROR: ----------------End of Error Message -------send entire
>     error message to petsc-maint at mcs.anl.gov----------
> 
>     Do you know where it comes from ?
> 
>     Thanks,
> 
>     Yann
> 
> 
> 
> -- 
> 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/>