Older versions of PETSc have generally less error checking and debug
versions have more error checking than optimized version. Once you update the
code for newer versions it should work clearly for both debug and optimized
versions of PETSc.
Barry
> On Jun 24, 2023, at 12:08 PM, Paweł
Look for places where you call VecSetValues(), make sure that after the call
before you use the vector for some other use you call
VecAssemblyBegin/VecAssemblyEnd
> On Jun 23, 2023, at 6:57 PM, Paweł Stebliński wrote:
>
> I'm using petsc 3.19.2 --with-debuging=yes (generated executable a
Search through all your code looking for calls to
VecGetArray()/VecRestoreArray(). For all uses where you only need to read from
the array, replace
the calls with VecGetArrayRead()/VecRestoreArrayRead(). And add the const
modifier to the array declaration.
> On Jun 24, 2023, at 5:32
On Fri, Jun 23, 2023 at 10:12 PM neil liu wrote:
> Thanks, Matt.
> It seems DMPlexComputeCellGeometryFEM works well with the quadrature
> points to deliver Jacobian and inverse one.
> Will it be a good choice ?
>
Yes, it is intended to compute these quantities for you. Let me know if it
does
