Implement evaluating derivatives of k_rtheta.

1 files changed, 12 insertions, 0 deletions

diff --git a/td_constraints.c b/td_constraints.c
index dd62150..41c1c5b 100644
--- a/td_constraints.c
+++ b/td_constraints.c
@@ -218,6 +218,18 @@ double tdi_constraint_eval_k_rtheta(TDConstraintEvalContext *ctx, double r, doub
     return fd->eval_krt(ctx, r, theta);
 }
 
+double tdi_constraint_eval_dk_rtheta_r(TDConstraintEvalContext *ctx, double r, double theta)
+{
+    const TDFamilyDef *fd = td_families[ctx->family];
+    return fd->eval_krt_dr(ctx, r, theta);
+}
+
+double tdi_constraint_eval_dk_rtheta_t(TDConstraintEvalContext *ctx, double r, double theta)
+{
+    const TDFamilyDef *fd = td_families[ctx->family];
+    return fd->eval_krt_dt(ctx, r, theta);
+}
+
 int tdi_constraint_eval_init(TDConstraintEvalContext *ctx)
 {
     TDConstraintEvalPriv *priv = ctx->priv;