Field Detail

• testMode

  private boolean testMode

• DBL_EPSILON

  private final double DBL_EPSILON

• DBL_MIN

  private final double DBL_MIN

  2**-1022 = D1MACH(1).

• DBL_MAX

  private final double DBL_MAX

  See Also:

  Constant Field Values

• LM_MACHEP

  private final double LM_MACHEP

  resolution of arithmetic

• LM_DWARF

  private final double LM_DWARF

• LM_SQRT_DWARF

  private final double LM_SQRT_DWARF

  Square should not underflow

• LM_SQRT_GIANT

  private final double LM_SQRT_GIANT

  Square should not overflow

• LM_USERTOL

  private final double LM_USERTOL

• ftol

  protected double ftol

  ftol is a nonnegative input variable. Termination occurs when both the actual and predicted relative reductions in the sum of squares are at most ftol. Therefore, ftol measures the relative error desired in the sum of squares.

• xtol

  protected double xtol

  xtol is a nonnegative input variable. Termination occurs when the relative error between two consecutive iterates is at most xtol. Therefore, xtol measures the relative error desired in the approximate solution.

• gtol

  protected double gtol

  gtol is a nonnegative input variable. Termination occurs when the cosine of the angle between fvec and any column of the jacobian is at most gtol in absolute value. Therefore, gtol measures the orthogonality desired between the function vector and the columns of the jacobian.

• epsfcn

  protected double epsfcn

  step used to calculate the jacobian epsfcn is an input variable used in choosing a step length for the forward-difference approximation. The relative errors in the functions are assumed to be of the order of epsfcn.

• factor

  protected double factor

  initial bound to steps in the outer loop factor is a positive input variable used in determining the initial step bound. This bound is set to the product of factor and the euclidean norm of diag*x if nonzero, or else to factor itself. In most cases factor should lie in the interval (0.1,100.0). Generally, the value 100.0 is recommended.

• maxcall

  protected int maxcall

  maximum number of iterations

• maxfev

  private int maxfev

  maxfev is a positive integer input variable. Termination occurs when the number of calls to lm_fcn is at least maxfev by the end of an iteration.

• mode

  protected int mode

  mode is an integer input variable. If mode = 1, the variables will be scaled internally. If mode = 2, the scaling is specified by the input diag.

• printflags

  protected int printflags

  Bits ored to show where to output diagnostics

• lm_infmsg

  private final java.lang.String[] lm_infmsg

  Set message texts (indexed by info)

• lm_shortmsg

  private final java.lang.String[] lm_shortmsg

• debugMessages

  protected boolean debugMessages

• debugMatrix

  protected boolean debugMatrix

• m

  private int m

  m is a positive integer input variable set to the number of functions.

• n

  private int n

  n is a positive integer variable set to the number of variables; n must not exceed m.

• x

  protected double[] x

  x is an array of length n. On input x must contain an initial estimate of the solution vector. On output x contains the final estimate of the solution vector.

• fnorm

  protected double fnorm

  norm of the residue vector fvec

• nfev

  protected int nfev

  actual number of iterations

• info

  protected int[] info

  status of minimization info is an integer OUTPUT variable that indicates the termination status of lm_lmdif as follows: info

  • fvec

    private double[] fvec

    fvec is an output array of length m which contains the functions evaluated at the output x.

  • diag

    protected double[] diag

    diag is an array of length n. If mode = 1, diag is internally set. If mode = 2, diag must contain positive entries that serve as multiplicative scale factors for the variables.

  • qtf

    private double[] qtf

    qtf is an OUTPUT array of length n which contains the first n elements of the vector (q transpose)*fvec.

  • fjac

    private double[] fjac

    fjac is an OUTPUT m by n array. The upper n by n submatrix of fjac contains an upper triangular matrix r with diagonal elements of nonincreasing magnitude such that pT*(jacT*jac)*p = rT*r (NOTE: T stands for matrix transposition), where p is a permutation matrix and jac is the final calculated jacobian. Column j of p is column ipvt(j) (see below) of the identity matrix. The lower trapezoidal part of fjac contains information generated during the computation of r.

  • wa1

    private double[] wa1

    wa1, wa2, and wa3 are work arrays of length n.

  • wa2

    private double[] wa2

  • wa3

    private double[] wa3

  • wa4

    private double[] wa4

    wa4 is a work array of length m, used among others to hold residuals from evaluate.

  • ipvt

    private int[] ipvt

    ipvt is an integer OUTPUT array of length n. It defines a permutation matrix p such that jac*p = q*r, where jac is the final calculated jacobian, q is orthogonal (not stored), and r is upper triangular with diagonal elements of nonincreasing magnitude. Column j of p is column ipvt(j) of the identity matrix.

  • par

    private double[] par

  • testCase

    private int testCase

  • xSeries

    private double[] xSeries

  • ySeries

    private double[] ySeries

  • DRAPER24D

    private final int DRAPER24D

    See Also:

    Constant Field Values

  • ROSENBROCK

    private final int ROSENBROCK

    See Also:

    Constant Field Values

  • FREUDENSTEIN_AND_ROTH

    private final int FREUDENSTEIN_AND_ROTH

    See Also:

    Constant Field Values

  • JENNRICH_AND_SAMPSON

    private final int JENNRICH_AND_SAMPSON

    See Also:

    Constant Field Values

  • HELICAL_VALLEY

    private final int HELICAL_VALLEY

    See Also:

    Constant Field Values

  • BARD

    private final int BARD

    See Also:

    Constant Field Values

  • MEYER

    private final int MEYER

    See Also:

    Constant Field Values

  • BOX_3D

    private final int BOX_3D

    See Also:

    Constant Field Values

  • POWELL_SINGULAR

    private final int POWELL_SINGULAR

    See Also:

    Constant Field Values

  • KOWALIK_AND_OSBORNE

    private final int KOWALIK_AND_OSBORNE

    See Also:

    Constant Field Values

  • BROWN_AND_DENNIS

    private final int BROWN_AND_DENNIS

    See Also:

    Constant Field Values

  • OSBORNE1

    private final int OSBORNE1

    See Also:

    Constant Field Values

  • OSBORNE2

    private final int OSBORNE2

    See Also:

    Constant Field Values

  • WATSON

    private final int WATSON

    See Also:

    Constant Field Values

  • HOCK25

    private final int HOCK25

    See Also:

    Constant Field Values

  • BROWN_ALMOST_LINEAR

    private final int BROWN_ALMOST_LINEAR

    See Also:

    Constant Field Values

  • LINEAR_FULL_RANK

    private final int LINEAR_FULL_RANK

    See Also:

    Constant Field Values

  • LINEAR_RANK1

    private final int LINEAR_RANK1

    See Also:

    Constant Field Values

  • LINEAR_RANK1_WITH_ZERO_COLUMNS_AND_ROWS

    private final int LINEAR_RANK1_WITH_ZERO_COLUMNS_AND_ROWS

    See Also:

    Constant Field Values

  • CHEBYQUAD

    private final int CHEBYQUAD

    See Also:

    Constant Field Values

  • LEVMAR_ROSENBROCK

    private final int LEVMAR_ROSENBROCK

    See Also:

    Constant Field Values

  • MODIFIED_ROSENBROCK

    private final int MODIFIED_ROSENBROCK

    See Also:

    Constant Field Values

  • POWELL_2_PARAMETER

    private final int POWELL_2_PARAMETER

    See Also:

    Constant Field Values

  • WOOD

    private final int WOOD

    See Also:

    Constant Field Values

  • HOCK1

    private final int HOCK1

    See Also:

    Constant Field Values

  • HOCK21_MODIFIED

    private final int HOCK21_MODIFIED

    See Also:

    Constant Field Values

  • HATFLDB

    private final int HATFLDB

    See Also:

    Constant Field Values

  • HATFLDC

    private final int HATFLDC

    See Also:

    Constant Field Values

  • EQUILIBRIUM_COMBUSTION

    private final int EQUILIBRIUM_COMBUSTION

    See Also:

    Constant Field Values

  • Constructor Detail

    • Lmmin

      public Lmmin()

    • Lmmin

      public Lmmin​(int m,
                   int n)

      Parameters:

      m - the number of functions

      n - the number of variables; n must not exceed m.

    • Lmmin

      public Lmmin​(int m,
                   int n,
                   double[] x)

      Parameters:

      m - the number of functions

      n - the number of variables; n must not exceed m.

      x -

  • Method Detail

    • statusMessage

      public void statusMessage​(int status)

    • getParameters

      public double[] getParameters()

    • fitTestModel

      private void fitTestModel()

    • dumpResults

      public void dumpResults()

    • getStatus

      public int getStatus()

    • driver

      public void driver()

      driver.

    • lm_lmdif

      private void lm_lmdif()

    • lm_qrfac

      private void lm_qrfac​(int m,
                            int n,
                            double[] a,
                            int pivot,
                            int[] ipvt,
                            double[] rdiag,
                            double[] acnorm,
                            double[] wa)

    • fitToFunction

      public abstract void fitToFunction​(double[] x,
                                         double[] fvec,
                                         int[] info)

    • lm_printout_std

      private void lm_printout_std​(double[] par,
                                   double[] fvec,
                                   int printflags,
                                   int iflag,
                                   int iter,
                                   int nfev)

    • lm_lmpar

      private void lm_lmpar​(int n,
                            double[] r,
                            int ldr,
                            int[] ipvt,
                            double[] diag,
                            double[] qtb,
                            double delta,
                            double[] par,
                            double[] x,
                            double[] sdiag,
                            double[] aux,
                            double[] xdi)

    • lm_qrsolv

      private void lm_qrsolv​(int n,
                             double[] r,
                             int ldr,
                             int[] ipvt,
                             double[] diag,
                             double[] qtb,
                             double[] x,
                             double[] sdiag,
                             double[] wa)

    • lm_enorm

      double lm_enorm​(int n,
                      double[] x,
                      int offset)

    • fitToTestFunction

      private void fitToTestFunction​(double[] par,
                                     double[] fvec)

      lm_enorm.

    • shiftedChebyshev

      private double shiftedChebyshev​(double x,
                                      int n)