Add a function for evaluating K_rtheta

1 files changed, 6 insertions, 0 deletions

diff --git a/td_constraints.c b/td_constraints.c
index fe6bc18..663b83a 100644
--- a/td_constraints.c
+++ b/td_constraints.c
@@ -192,6 +192,12 @@ constraint_eval_mom_t_confflat(TDConstraintEvalContext *ctx,
     }
 }
 
+double tdi_constraint_eval_k_rtheta(TDConstraintEvalContext *ctx, double r, double theta)
+{
+    const TDFamilyDef *fd = td_families[ctx->family];
+    return fd->eval_krt(ctx, r, theta);
+}
+
 int tdi_constraint_eval_init(TDConstraintEvalContext *ctx)
 {
     TDConstraintEvalPriv *priv = ctx->priv;