Skip to main content
This is a pre-production deployment of Warehouse. Changes made here affect the production instance of PyPI (
Help us improve Python packaging - Donate today!

A Python module for solving optimization problems with nonlinear least-squares

Project Description

Consult the module API page at

for all information related to this module, including information regarding the latest changes to the code. The page at the URL shown above lists all of the module functionality you can invoke in your own code.

With regard to the basic purpose of this module, it provides a domain agnostic implementation of nonlinear least-squares algorithms (gradient-descent and Levenberg-Marquardt) for fitting a model to observed data. Typically, a model involves several parameters and each observed data element can be expressed as a function of those parameters plus noise. The goal of nonlinear least-squares is to estimate the best values for the parameters given all of the observed data. In order to illustrate how to use the NonlinearLeastSquares class, the module also comes with another class, OptimizeSurfaceFit, whose job is to fit the best surface to noisy height data over an XY-plane. The model in this case would be an analytical expression for the height surface and the goal of nonlinear least-squares would be to estimate the best values for the parameters in the model.

Typical usage syntax for invoking the domain-agnostic NonlinearLeastSquares through your own domain-specific class such as OptimizeSurfaceFit is shown below:

optimizer =  NonlinearLeastSquares(
                 max_iterations = 200,
                 delta_for_jacobian = 0.000001,
                 delta_for_step_size = 0.0001,

surface_fitter = OptimizeSurfaceFit(
                     gen_data_synthetically = True,
                     datagen_functional = "7.8*(x - 0.5)**4 + 2.2*(y - 0.5)**2",
                     data_dimensions = (16,16),
                     how_much_noise_for_synthetic_data = 0.3,
                     model_functional = "a*(x-b)**4 + c*(y-d)**2",
                     initial_param_values = {'a':2.0, 'b':0.4, 'c':0.8, 'd':0.4},
                     display_needed = True,
                     debug = True,


result = surface_fitter.calculate_best_fitting_surface('lm')


result = surface_fitter.calculate_best_fitting_surface('gd')

Release History

This version
History Node


History Node


History Node


Download Files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Filename, Size & Hash SHA256 Hash Help File Type Python Version Upload Date
(56.9 kB) Copy SHA256 Hash SHA256
Source None Dec 14, 2016

Supported By

Elastic Elastic Search Pingdom Pingdom Monitoring Dyn Dyn DNS Sentry Sentry Error Logging CloudAMQP CloudAMQP RabbitMQ Heroku Heroku PaaS Kabu Creative Kabu Creative UX & Design Fastly Fastly CDN DigiCert DigiCert EV Certificate Google Google Cloud Servers DreamHost DreamHost Log Hosting