inverse_wlc now handles F=HUGE_VAL.

diff --git a/src/sawsim.nw b/src/sawsim.nw
index cc506972b42ba0871e313e0daa70e3348d138c96..f1e7f2c18256fffc84721d5f47fa491064fc2d1b 100644 (file)
--- a/src/sawsim.nw
+++ b/src/sawsim.nw
@@ -3945,6 +3945,8 @@ static double inverse_wlc(double F, double T, double p, double L)
   int num_roots;
   assert(F >= 0);
   assert(T > 0); assert(p > 0); assert(L > 0);
+  if (F == HUGE_VAL)
+    return L;
   num_roots = gsl_poly_solve_cubic(-(FT+2.25),2*FT+1.5,-FT, &xL0, &xL1, &xL2);
   assert(num_roots > 0); assert(xL0 >= -DOUBLE_PRECISION); assert(xL0 < 1);
   if (xL0 < 0) xL0 = 0.0;