Re: [petsc-users] DMPlex and CGNS

2023-01-16 Thread Jed Brown
Matthew Knepley writes: > On Mon, Jan 16, 2023 at 6:15 PM Jed Brown wrote: > >> How soon do you need this? I understand the grumbling about CGNS, but it's >> easy to build, uses HDF5 parallel IO in a friendly way, supports high order >> elements, and is generally pretty expressive. I wrote a

Re: [petsc-users] DMPlex and CGNS

2023-01-16 Thread Matthew Knepley
On Mon, Jan 16, 2023 at 6:15 PM Jed Brown wrote: > How soon do you need this? I understand the grumbling about CGNS, but it's > easy to build, uses HDF5 parallel IO in a friendly way, supports high order > elements, and is generally pretty expressive. I wrote a parallel writer > (with some

Re: [petsc-users] DMPlex and CGNS

2023-01-16 Thread Jed Brown
How soon do you need this? I understand the grumbling about CGNS, but it's easy to build, uses HDF5 parallel IO in a friendly way, supports high order elements, and is generally pretty expressive. I wrote a parallel writer (with some limitations that I'll remove) and plan to replace the current

[petsc-users] DMPlex and CGNS

2023-01-16 Thread Ferrand, Jesus A.
Dear PETSc team: I would like to use DMPlex to partition a mesh stored as a CGNS file. I configured my installation with --download_cgns = 1, got me a .cgns file and called DMPlexCreateCGNSFromFile() on it. Doing so got me this error: [0]PETSC ERROR: - Error Message

Re: [petsc-users] Nonconforming object sizes using TAO (TAOBNTR)

2023-01-16 Thread Alexis Marboeuf
Hi Matt, After investigation, it fails because, at some point, the boolean needH is set to PETSC_FALSE when initializing the BNK method with TAOBNKInitialize (line 103 of $PETSC_DIR/src/tao/bound/impls/bnk/bntr.c). The Hessian and the precondtitioner are thus not updated throughout the TAO