- 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

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

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

- remnants from older code prior to wrapping PETSc vectors directly.

- leave gAverage() and gSum(). Still potentially useful.

- might have PREFIX/ARCH/include PREFIX/include

- change very misleading enum from `cachingTypes` to `updateTypes`
  and fix enumeration order.

- default constructed cache was "None" (meant to be no caching)
  but should have been "Always" (ie, always update).

- consistent with other modules

- support mixed-precision (SPDP) mode

stefano zampini authored 4 years ago and Mark OLESEN committed 4 years ago

- use KSP_NORM_NONE if using foam convergence criterion

- support monitoring OpenFOAM residual reduction irrespective of the
  convergence testing

- remove occurences of PETSC_NULL and nullptr in favor of NULL

stefano zampini authored 4 years ago and Mark OLESEN committed 4 years ago

- call SetUseInitialGuessNonzero before KSPSetFromOptions, so that
  users can still customize it.

- We no longer need a patched petsc version that compute the L1 residual
  for us.

- Use a monitor to report the initial residual from PETSc when using
  PETSc convergence testing.

- Cleanup source code from unused functions/variables

- Use PetscFunctionBeginUser/PetscFunctionReturn and PETSc error
  handling if called from PETSc

stefano zampini authored 4 years ago and Mark OLESEN committed 4 years ago

- the general MatSetUp call is a dummy call if we already preallocated
  with the above preallocators this is intended to allow for a
  different matrix type, only for the purpose of testing

stefano zampini authored 4 years ago and Mark OLESEN committed 4 years ago

- need VecResetArray otherwise the CPU array placed will never be updated

STYLE: make PetscWrappedVector constructor explicit, use default constructor

- compile into FOAM_USER_LIBBIN for consistency with other
  modules

- normally not defined, but can be used to define an alternative
  installation location. More refinements in the future, pending
  changes in OpenFOAM.