[petsc-dev] DMPlexComputeGeometryFEM

[Matthew Knepley]

On Fri, Nov 25, 2022 at 7:14 AM Yann Jobic <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

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/ 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/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20221126/ad460dd7/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: testComputeGeomFEM.c
Type: application/octet-stream
Size: 2565 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20221126/ad460dd7/attachment.obj>