From 70d2537dfd98ff67ed7678f422a6888f6eeb3191 Mon Sep 17 00:00:00 2001
From: Phillip Alday <me@phillipalday.com>
Date: Fri, 20 May 2022 10:41:52 -0500
Subject: [PATCH] add comment

---
 src/linearmixedmodel.jl | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/src/linearmixedmodel.jl b/src/linearmixedmodel.jl
index 4c4c7ac6c..92ed27fda 100644
--- a/src/linearmixedmodel.jl
+++ b/src/linearmixedmodel.jl
@@ -444,6 +444,11 @@ function fit!(
         val = try
             objective(updateL!(setθ!(m, x)))
         catch ex
+            # This can happen when the optimizer drifts into an area where
+            # there isn't enough shrinkage. Why finitial? Generally, it will
+            # be the (near) worst case scenario value, so the optimizer won't
+            # view it as an optimum. Using Inf messes up the quadratic
+            # approximation in BOBYQA.
             ex isa PosDefException || rethrow()
             iter == 0 && rethrow()
             m.optsum.finitial