Merge Jacobian coloring into PETSc solver

include/bout/petsc_interface.hxx

@@ -6,9 +6,9 @@
  * up a linear system.
  *
  **************************************************************************
- * Copyright 2019 C. MacMackin
+ * Copyright 2019 - 2025 BOUT++ contributors
  *
- * Contact: Ben Dudson, bd512@york.ac.uk
+ * Contact: Ben Dudson, dudson2@llnl.gov
  *
  * This file is part of BOUT++.
  *
@@ -573,7 +573,7 @@ PetscVector<T> operator*(const PetscMatrix<T>& mat, const PetscVector<T>& vec) {
 namespace bout {
 template <class T>
 constexpr auto cast_MatFDColoringFn(T func) {
-  return func;
+  return reinterpret_cast<MatFDColoringFn*>(func); // NOLINT(*-reinterpret-cast)
 }
 } // namespace bout
 #else

include/bout/petsclib.hxx

@@ -22,9 +22,9 @@
  * so it *must* be included before *any* PETSc header.
  * 
  **************************************************************************
- * Copyright 2012 B.D.Dudson, S.Farley, M.V.Umansky, X.Q.Xu
+ * Copyright 2012 - 2025 BOUT++ contributors
  *
- * Contact: Ben Dudson, bd512@york.ac.uk
+ * Contact: Ben Dudson, dudson2@llnl.gov
  * 
  * This file is part of BOUT++.
  *
@@ -60,6 +60,7 @@ class Options;
 #define PETSC_HAVE_BROKEN_RECURSIVE_MACRO
 
 #include <petsc.h> // IWYU pragma: export
+#include <petscts.h>
 #include <petscversion.h>
 
 #include "bout/boutexception.hxx"
@@ -113,6 +114,9 @@ public:
   /// was passed to the constructor.
   void setOptionsFromInputFile(SNES& snes);
 
+  /// Set options for a TS time integrator
+  void setOptionsFromInputFile(TS& ts);
+
   /*!
    * Force cleanup. This will call PetscFinalize, printing a warning
    * if any instances of PetscLib still exist