C/C++ Minpack

This is a C version of the minpack minimization package. It has been derived from the fortran code using f2c and some limited manual editing. Note that you need to link against libf2c to use this version of minpack. Extern "C" linkage permits the package routines to be called from C++. Check ftp://netlib.bell-labs.com/netlib/f2c for the latest f2c version. For general minpack info and test programs, see the accompanying readme.txt and http://www.netlib.org/minpack/.

CMake is the standard build system:

cmake -B build && cmake --build build

and ctest --test-dir build runs the tests. See the project home page for build options (precision variants, BLAS/LAPACK, shared libraries). A plain Makefile is also provided for backward compatibility: type make to compile and make install to install in /usr/local, or modify the Makefile to suit your needs.

Manolis Lourakis -- lourakis at ics forth gr, July 2002 Institute of Computer Science, Foundation for Research and Technology - Hellas Heraklion, Crete, Greece

Repackaging by Frederic Devernay -- frederic dot devernay at m4x dot org

The project home page is at http://devernay.github.io/cminpack

Testing

Two test routes are available:

• CMake / CTest (recommended): after building, run ctest --test-dir build. This runs the standard example tests (compared against examples/ref/*.ref with the dependency-free C tool cmpfiles), the self-checking regression tests, the intensive driver programs as smoke tests (run to completion, no NaN), and — when Python 3 is available — the pure-C-vs-f2c cross-check. Set -DCMINPACK_CROSSCHECK=OFF to disable that cross-check.
• Makefile: make check runs the standard tests for the double, long double and float builds, then make crosscheck. The double standard tests are strict (a failure fails the build); the long double and float tests are informational.

The cross-check (examples/crosscheck.py) verifies the one strict, deterministic invariant — the pure-C and f2c implementations are byte-identical at double precision. Differences against the original FORTRAN MINPACK, and single/extended-precision iteration-path differences, are reported for information only (see the next section).

Python 3 is optional. It is not needed to build cminpack, nor for the standard or smoke tests (which use cmpfiles); only the cross-check uses it (via examples/compare.py and examples/driver_check.py). When Python 3 is missing, both build systems skip only that cross-check — printing a disclaimer that coverage is incomplete — and run everything else.

Numerical differences from FORTRAN MINPACK

On the difficult test problems (the Moré/Garbow/Hillstrom functions exercised by the intensive driver programs examples/*drv*), cminpack and the original FORTRAN MINPACK can report different iteration, function- and Jacobian-evaluation counts, and last-digit differences in the results. This is expected and harmless: the problems still converge to equally valid solutions. Test messages reporting such differences are normal, not a sign of a broken build.

The difficult test problems themselves are from Moré, Garbow & Hillstrom, Testing Unconstrained Optimization Software, ACM TOMS 7(1):17-41 (1981) (doi:10.1145/355934.355936; free open-access technical report ANL-AMD-TM-324: https://www.osti.gov/servlets/purl/6650344), and the companion Algorithm 566, ACM TOMS 7(1):136-140 (1981) (doi:10.1145/355934.355943). They are defined in examples/ssqfcn.c (least squares) and examples/vecfcn.c (systems of equations); the driver programs run each one from progressively harder starting points (x0, 10*x0, 100*x0), so the hardest problems (e.g. problem 10, the Meyer function) are deliberately pushed to non-convergence.

The differences do not come from the algorithm:

• dpmpar returns identical machine constants on both sides, so the convergence tolerance is the same.
• cminpack's own pure-C and f2c builds are bit-identical to each other.
• The dominant cause is floating-point contraction: gcc and gfortran fuse a*b + c into a fused multiply-add (FMA) at different places, so intermediate values differ by one unit in the last place (ULP). On ill-conditioned problems a single ULP early in the iteration can flip a trust-region accept/reject decision, and the two runs then follow different paths to (equally valid) results. Compiling both sides with -ffp-contract=off removes most of the divergence.

(The enorm scaling constants rdwarf/rgiant, which cminpack tunes to the IEEE range rather than MINPACK's 1980 values, are a separate and minor effect: they only change the norm of out-of-range vectors in the last bit. See src/enorm.c for the full explanation.)

Because of this, examples/crosscheck.py and make check treat pure-C == f2c as the strict, deterministic invariant — a mismatch there is a real bug and fails the build — while differences against the original FORTRAN, and single/extended-precision driver differences, are reported for information only and never fail the build.

History

• version 1.3.13 (16/07/2026):

  • Fix a division by zero in covar1 when the Jacobian rank equals the number of residuals (m == rank, e.g. square full-rank problems), which produced Inf/NaN throughout the covariance matrix
  • Make the banded finite-difference branch of fdjac1 call the user function with iflag=2, like the dense branch, fdjac2 and the FORTRAN version
  • Fix the USE_BLAS Newton correction in lmpar to use all n components (it was truncated to the original Jacobian rank, giving a wrong step for rank-deficient problems)
  • Fix the jpvt memset size in the USE_LAPACK qrfac, and make the work-array size checks in hybrd, hybrj, hybrd1, hybrj1 and lmdif1 overflow-safe (also in the f2c versions)
  • Guard a harmless transient infinity in dogleg for rank-deficient Jacobians (also in the f2c version)
  • Add a USE_LAPACK CMake option to build the LAPACK-based QR factorization (qrfac), which was previously only reachable through the Makefile
  • Exercise the BLAS and LAPACK builds in CI and run the full test suite (including tlmdifc) for them; add the intensive driver programs (lmddrvc, lmfdrvc, lmsdrvc, hyjdrvc, hybdrvc, chkdrvc) as smoke tests and run the FORTRAN example tests in CI

• version 1.3.12 (16/07/2026):

  • Fix non-convergence/NaN in lmder, lmdif and lmstr on problems whose solution is the zero vector, by guarding a 0/0 division in lmpar #76
  • Compare test output using a numeric tolerance instead of an exact text match, so the tests pass across compilers and math libraries #37 #77
  • Fix Windows linking: correct the DLL export/import macros and document that static-library users define CMINPACK_NO_DLL #18
  • Detect and link the CBLAS interface when USE_BLAS is enabled #12
  • Add cminpackcpp.hpp, a header-only C++ wrapper so the solvers accept lambdas, functors and std::function #74
  • Make CMake the standard build system (the Makefile is kept for backward compatibility) and remove the unmaintained Xcode, Visual Studio and Eclipse project files (cminpack.xcodeproj, cminpack*.vcproj/.vcxproj, cminpack.sln, .cproject, .project)
  • Move CI to GitHub Actions (keeping the CMake-based AppVeyor build) and remove the obsolete Travis CI, Coveralls and Coverity configuration and badges

• version 1.3.11 (13/09/2024):

  • Bump installed version number to 1.3.11 #75

• version 1.3.10 (11/09/2024):

  • Disable BLAS by default #66 #70
  • Fix BLAS usage (broken by #58) #68
  • Fix testing on Windows #63
  • Fix pkg-config files #71
  • Bump minimum CMake version #69

• version 1.3.9 (28/05/2024):

  • CMake portability fixes #50 #53 #56 #57 #58 #62
  • MKL-related fixes #51 #52
  • Support more CI build configurations #59 #61

• version 1.3.8 (09/02/2021):

  • CMake now builds by default the single-, double-, and extended-precision versions #45 #48
  • Avoid promoting to doubles in all operations for the single-precision version #47

• version 1.3.7 (09/12/2020):

  • Makefile cleanups #11
  • Cmake-related fixes #20 #21 #23 #27 #28
  • Add Appveyor CI #24
  • Add support for single-precision CBLAS and LAPACK #40

• version 1.3.6 (24/02/2017):

  • Fix FreeBSD build #6
  • CMake: install CMinpackConfig.cmake rather than FindCMinpack.cmake #8
  • CMake: add option USE_BLAS to compile with blas #9

• version 1.3.5 (28/05/2016):

  • Add support for compiling a long double version (Makefile only).
  • CMake: static libraries now have the suffix _s.

• version 1.3.4 (28/05/2014):

  • Add FindCMinpack.cmake cmake module. If you use the cmake install, finding CMinpack from your CMakeLists.txt is as easy as find_package(CMinpack).

• version 1.3.3 (04/02/2014):

  • Add documentation and examples abouts how to add box constraints to the variables.
  • continuous integration https://travis-ci.org/devernay/cminpack

• version 1.3.2 (27/10/2013):

  • Minor change in the CMake build: also set SOVERSION.

• version 1.3.1 (02/10/2013):

  • Fix CUDA examples compilation, and remove non-free files.

• version 1.3.0 (09/06/2012):

  • Optionally use LAPACK and CBLAS in lmpar, qrfac, and qrsolv. Added make lapack to build the LAPACK-based cminpack and "make checklapack" to test it (results of the test may depend on the underlying LAPACK and BLAS implementations). On 64-bits architectures, the preprocessor symbol __LP64__ must be defined (see cminpackP.h) if the LAPACK library uses the LP64 interface (i.e. 32-bits integer, vhereas the ILP interface uses 64 bits integers).

• version 1.2.2 (16/05/2012):

  • Update Makefiles and documentation (see "Using CMinpack" above) for easier building and testing.

• version 1.2.1 (15/05/2012):

  • The library can now be built as double, float or half versions. Standard tests in the "examples" directory can now be lauched using make check (to run common tests, including against the float version), make checkhalf (to test the half version) and make checkfail (to run all the tests, even those that fail).

• version 1.2.0 (14/05/2012):

  • Added original FORTRAN sources for better testing (type make -C fortran, then make -C examples and follow the instructions). Added driver tests lmsdrv, chkdrv, hyjdrv, hybdrv. make -C examples alltest will run all possible test combinations (make sure you have gfortran installed).

• version 1.1.5 (04/05/2012):

  • cminpack now works in CUDA, thanks to Jordi Bataller Mascarell, type make -C cuda (be careful, though: this is a straightforward port from C, and each problem is solved using a single thread). cminpack can now also be compiled with single-precision floating point computation (define __cminpack_real__ to float when compiling and using the library). Fix cmake support for CMINPACK_LIB_INSTALL_DIR. Update the reference files for tests.

• version 1.1.4 (30/10/2011):

  • Translated all the Levenberg-Marquardt code (lmder, lmdif, lmstr, lmder1, lmdif1, lmstr1, lmpar, qrfac, qrsolv, fdjac2, chkder`) to use C-style indices.

• version 1.1.3 (16/03/2011):

  • Minor fix: Change non-standard strnstr() to strstr() in genf77tests.c.

• version 1.1.2 (07/01/2011):

  • Fix Windows DLL building (David Graeff) and document covar in cminpack.h.

• version 1.1.1 (04/12/2010):

  • Complete rewrite of the C functions (without trailing underscore in the function name). Using the original FORTRAN code, the original algorithms structure was recovered, and many goto's were converted to if...then...else. The code should now be both more readable and easier to optimize, both for humans and for compilers. Added lmddrv and lmfdrv test drivers, which test a lot of difficult functions (these functions are explained in Testing Unconstrained Optimization Software by Moré et al.). Also added the pkg-config files to the cmake build, as well as an "uninstall" target, contributed by Geoffrey Biggs.

• version 1.0.4 (18/10/2010):

  • Support for shared library building using CMake, thanks to Goeffrey Biggs and Radu Bogdan Rusu from Willow Garage. Shared libraries can be enabled using cmake options, as in: cmake -DUSE_FPIC=ON -DBUILD_SHARED_LIBS=ON -DBUILD_EXAMPLES=OFF path_to_sources

• version 1.0.3 (18/03/2010):

  • Added CMake support.
  • XCode build is now Universal.
  • Added tfdjac2_ and tfdjac2c examples, which test the accuracy of a finite-differences approximation of the Jacobian.
  • Bug fix in tlmstr1 (signaled by Thomas Capricelli).

• version 1.0.2 (27/02/2009):

  • Added Xcode and Visual Studio project files

• version 1.0.1 (17/12/2007):

  • bug fix in covar() and covar_(), the computation of tolr caused a segfault (signaled by Timo Hartmann).

• version 1.0.0 (24/04/2007):

  • Added FORTRAN and C examples
  • Added documentation from Debian man pages
  • Wrote pure C version
  • Added covar() and covar_(), and use it in tlmdef/tlmdif