RF Eigenproblem

firedrake/__init__.py

@@ -83,6 +83,7 @@ def init_petsc():
 from firedrake.utility_meshes import *
 from firedrake.variational_solver import *
 from firedrake.eigensolver import *
+from firedrake.rf_eigensolver import *
 from firedrake.ensemble import *
 from firedrake.randomfunctiongen import *
 from firedrake.external_operators import *

firedrake/adjoint/__init__.py

@@ -27,6 +27,11 @@
     continue_disk_checkpointing, stop_disk_checkpointing, \
     checkpointable_mesh  # noqa F401
 from firedrake.adjoint_utils import get_solve_blocks  # noqa F401
+from pyadjoint.optimization.tao_solver import (
+    ReducedFunctionalMat, ReducedFunctionalMatCtx, HessianAction,
+    TLMAction, AdjointAction, RFAction, check_rf_action,
+    flatten_parameters, PETScVecInterface, get_valid_comm
+)
 
 from pyadjoint.verification import taylor_test, taylor_to_dict  # noqa F401
 from pyadjoint.drivers import compute_gradient, compute_derivative, compute_hessian  # noqa F401