bugfixes from PETSc's developer Junchao Zhang

See merge request !14

Update code for improved error handling, caching, etc

Closes #34

See merge request !13

stefano zampini authored 1 year ago and Mark OLESEN committed 1 year ago

COMP: noexcept for casting petscWrappedVector to Vec

stefano zampini authored 1 year ago and Mark OLESEN committed 1 year ago

Two macros:
  AssertPETSc is meant to be called within OpenFOAM call that calls PETSc
  PetscCall instead is meant to be called within callbacks used by PETSc

- previously had UPstream::waitRequests() instead of
  UPstream::waitRequests(startOfRequests), which means that it would
  implicitly wait for *all* requests to finished (not just the local
  scope).

- also remove mpi-specific locations

- use setter method for context initialized

- fix constness on updateMat

COMP: use std::max_element instead of Foam::max(List..)

- PetscInt may be an unusually sized type and the Foam::max(List...)
  convenience field function may fail.
  Eg, if PetscInt is something like 'long long' and there is
  no pTraits<long long>::min, max defined.

COMP: add additional library link dependencies

See merge request !9

Support for fast GPU matrix assembly

See merge request !8

Fix caching

See merge request !7

with the code before, the periodic iters for the preconditioner were
updated only if the matrix was changed. Now, periodic updates are decoupled
Added PC setup stage to better distinguish low level PETSc events

- for OpenFOAM-v2006 shifted the location of flange.ans into tutorial
  resources specifically to make tutorial tests easier for this add-on
  solver.

  Only attempt to copy flange.ans for older versions to avoid spurious
  error message for new versions.

- unused fall-through value, but should default to "Always"