fix tolerance choice in fh

Author: nathanemac

src/fh_model.jl

@@ -1,6 +1,6 @@
 export fh_model, fh_nls_model
 
-function FH_smooth_term(; abstol = 1e-14, reltol = 1e-14)
+function FH_smooth_term()
   # FH model = van-der-Pol oscillator when I = b = c = 0
   # x' = μ(x - x^3/3 - y)
   # y' = x/μ -> here μ = 12.5
@@ -29,7 +29,7 @@ function FH_smooth_term(; abstol = 1e-14, reltol = 1e-14)
   data = noise .+ b
 
   # solve FH with parameters p
-  function simulate(p)
+  function simulate(p, abstol = 1e-14, reltol = 1e-14)
     temp_prob = DifferentialEquations.remake(prob_FH, p = p)
     sol = DifferentialEquations.solve(
       temp_prob,
@@ -46,15 +46,15 @@ function FH_smooth_term(; abstol = 1e-14, reltol = 1e-14)
   end
 
   # define residual vector
-  function residual(p)
-    F = simulate(p)
+  function residual(p, args...)
+    F = simulate(p, args...)
     F .-= data
     return F
   end
 
   # misfit = ‖residual‖² / 2
-  function misfit(p)
-    F = residual(p)
+  function misfit(p, args...)
+    F = residual(p, args...)
     return dot(F, F) / 2
   end