Solve the discretized system exactly on the coarsest level.

1 files changed, 7 insertions, 0 deletions

diff --git a/ell_grid_solve.h b/ell_grid_solve.h
index 70782c8..a12e9c3 100644
--- a/ell_grid_solve.h
+++ b/ell_grid_solve.h
@@ -52,6 +52,13 @@ enum EGSType {
      * mg2di_egs_solve() does a single relaxation step
      */
     EGS_SOLVER_RELAXATION,
+    /**
+     * Solve the equation exactly by contructing a linear system and solving it with LAPACK.
+     *
+     * solver_data is NULL
+     * mg2di_egs_solve() solves the discretized system exactly (up to roundoff error)
+     */
+    EGS_SOLVER_EXACT,
 };
 
 typedef struct EGSInternal EGSInternal;