Bugfix to stop crash when initialized at minimizer (fixes #24)

Author: lindonroberts

dfols/__init__.py

@@ -39,7 +39,7 @@
 from __future__ import absolute_import, division, print_function, unicode_literals
 
 # DFO-LS version
-__version__ = '1.5.2'
+__version__ = '1.5.3'
 
 # Main solver & exit flags
 from .solver import *

dfols/solver.py

@@ -148,7 +148,9 @@ def solve_main(objfun, x0, argsf, xl, xu, projections, npt, rhobeg, rhoend, maxf
                 exit_info = ExitInformation(EXIT_SUCCESS, "Objective is sufficiently small")
 
         if exit_info is not None:
-            return x0, r0_avg, sumsq(r0_avg), None, num_samples_run, nf, nx, nruns_so_far+1, exit_info, diagnostic_info
+            xmin_eval_num = 0
+            jacmin_eval_nums = np.array([0], dtype=int)
+            return x0, r0_avg, sumsq(r0_avg), None, num_samples_run, nf, nx, nruns_so_far+1, exit_info, diagnostic_info, xmin_eval_num, jacmin_eval_nums
 
     else:  # have old r0 information (e.g. from previous restart), use this instead
 

dfols/tests/test_solver.py

@@ -210,3 +210,22 @@ def runTest(self):
         self.assertTrue(array_compare(soln.resid, rosenbrock(soln.x), thresh=1e-10), "Wrong resid")
         self.assertTrue(array_compare(soln.jacobian, rosenbrock_jacobian(soln.x), thresh=1e-2), "Wrong Jacobian")
         self.assertTrue(abs(soln.obj - fmin) < 1e-4, "Wrong fmin")
+
+class TestDiagnosticInto(unittest.TestCase):
+    # Return diagnostic info - minimal working example for issue #23
+    def runTest(self):
+        ## Minimal working example for DFO-LS bug #23
+        objfun = lambda x: x
+        x0 = np.arange(3)
+        xmin = np.zeros(x0.shape)
+        soln = dfols.solve(objfun=objfun, x0=x0, user_params={"logging.save_diagnostic_info": True})
+        self.assertTrue(array_compare(soln.x, xmin, thresh=1e-2), "Wrong xmin")
+
+class TestStartAtMin(unittest.TestCase):
+    # Start solver at solution - minimal working example for issue #24
+    def runTest(self):
+        objfun = lambda x: x
+        x0 = np.zeros(3)
+        xmin = x0.copy()
+        soln = dfols.solve(objfun=objfun, x0=x0, user_params={"logging.save_diagnostic_info": True})
+        self.assertTrue(array_compare(soln.x, xmin, thresh=1e-8), "Wrong xmin")

docs/build/doctrees/environment.pickle

@@ -1,4 +1,4 @@
 # Sphinx build info version 1
 # This file hashes the configuration used when building these files. When it is not found, a full rebuild will be done.
-config: 865700494aceade7d85c2a3d7b404b00
+config: e27aac59f14e76a4142e0b9f9d374216
 tags: 645f666f9bcd5a90fca523b33c5a78b7

docs/build/doctrees/history.doctree

@@ -75,3 +75,7 @@ Version 1.5.1 (10 Oct 2024)
 Version 1.5.2 (28 Oct 2024)
 ---------------------------
 * Bugfix for saving diagnostic info (bug introduced in v1.5.1)
+
+Version 1.5.3 (30 Oct 2024)
+---------------------------
+* Bugfix when starting solver at problem minimizer (bug introduced in v1.5.1)

docs/build/doctrees/index.doctree

@@ -1,5 +1,5 @@
 const DOCUMENTATION_OPTIONS = {
-    VERSION: '1.5.2',
+    VERSION: '1.5.3',
     LANGUAGE: 'en',
     COLLAPSE_INDEX: false,
     BUILDER: 'html',

docs/build/html/.buildinfo

@@ -4,7 +4,7 @@
   <meta charset="utf-8" /><meta name="viewport" content="width=device-width, initial-scale=1" />
 
   <meta name="viewport" content="width=device-width, initial-scale=1.0" />
-  <title>Advanced Usage &mdash; DFO-LS v1.5.2 documentation</title>
+  <title>Advanced Usage &mdash; DFO-LS v1.5.3 documentation</title>
       <link rel="stylesheet" type="text/css" href="_static/pygments.css?v=fa44fd50" />
       <link rel="stylesheet" type="text/css" href="_static/css/theme.css?v=19f00094" />
 
@@ -15,7 +15,7 @@
 
         <script src="_static/jquery.js?v=5d32c60e"></script>
         <script src="_static/_sphinx_javascript_frameworks_compat.js?v=2cd50e6c"></script>
-        <script src="_static/documentation_options.js?v=73275c37"></script>
+        <script src="_static/documentation_options.js?v=d75fd82e"></script>
         <script src="_static/doctools.js?v=888ff710"></script>
         <script src="_static/sphinx_highlight.js?v=dc90522c"></script>
         <script async="async" src="https://cdn.jsdelivr.net/npm/mathjax@3/es5/tex-mml-chtml.js"></script>
@@ -38,7 +38,7 @@
             DFO-LS
           </a>
               <div class="version">
-                1.5.2
+                1.5.3
               </div>
 <div role="search">
   <form id="rtd-search-form" class="wy-form" action="search.html" method="get">