From 0bc654a449de9308beaad98ee00861338654a447 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Fri, 9 Feb 2018 14:55:29 +0100 Subject: Attempt to enforce regularity at origin. probably useless, committing half-done in case it's ever needed again --- nlsolve.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'nlsolve.c') diff --git a/nlsolve.c b/nlsolve.c index 5d691e0..bb7e083 100644 --- a/nlsolve.c +++ b/nlsolve.c @@ -97,7 +97,7 @@ struct NLSolvePriv { }; int tdi_nlsolve_solve(NLSolveContext *ctx, NLEqCallback eq_eval, - NLEqJacobianCallback eq_jac_eval, void *opaque, double *coeffs) + NLEqJacobianCallback eq_jac_eval, void *opaque, double *coeffs, double L) { NLSolvePriv *s = ctx->priv; int64_t start, totaltime_start; @@ -189,7 +189,7 @@ int tdi_nlsolve_solve(NLSolveContext *ctx, NLEqCallback eq_eval, } // solve for delta - ret = tdi_pssolve_solve(s->ps_ctx, s->eq_coeffs, s->rhs, s->delta); + ret = tdi_pssolve_solve(s->ps_ctx, s->eq_coeffs, s->rhs, s->delta, L); if (ret < 0) return ret; -- cgit v1.2.3