Fix fh tol (#74)

* add tolerance parameter in fh

Author: nathanemac

src/fh_model.jl

@@ -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