Add public API for evaluating the residual at an arbitrary location.

1 files changed, 9 insertions, 0 deletions

diff --git a/teukolsky_data.h b/teukolsky_data.h
index 786554e..861613d 100644
--- a/teukolsky_data.h
+++ b/teukolsky_data.h
@@ -207,4 +207,13 @@ int td_eval_lapse(const TDContext *td,
                   const unsigned int diff_order[2],
                   double *out);
 
+/**
+ * Evaluate the residual of the three constraint equations at the specified
+ * coordinates.
+ */
+int td_eval_residual(const TDContext *td,
+                     size_t nb_coords_r,     const double *r,
+                     size_t nb_coords_theta, const double *theta,
+                     double *out[3]);
+
 #endif /* TEUKOLSKY_DATA_H */