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.jo...@univ-amu.fr
<mailto:yann.jo...@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-ma...@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/>
