#include <cern_mroot_root.h>
This one-dimensional root-finding routine, based on cern_mroot, is probably slower than the more typical 1-d routines, but also tends to converge for a larger class of functions than cern_root, gsl_root_brent, or gsl_root_stef. It has been modified from cern_mroot and slightly optimized, but has the same basic behavior.
If denotes the current iteration, and
denotes the previous iteration, then the calculation is terminated if either (or both) of the following tests is successful
Definition at line 63 of file cern_mroot_root.h.
Public Member Functions | |
cern_mroot_root () | |
virtual | ~cern_mroot_root () |
int | get_info () |
Get the value of INFO from the last call to solve() (default 0). | |
virtual const char * | type () |
Return the type, "cern_mroot_root" . | |
virtual int | solve (double &ux, param_t &pa, func_t &func) |
Solve func using x as an initial guess, returning x . | |
Data Fields | |
int | maxf |
Maximum number of function evaluations. | |
double | scale |
The original scale parameter from CERNLIB (default 10.0). | |
double | eps |
The smallest floating point number (default ![]() | |
Protected Attributes | |
int | info |
Internal storage for the value of info . |
int get_info | ( | ) | [inline] |
Get the value of INFO
from the last call to solve() (default 0).
The value of info is assigned according to the following list. The values 1-8 are the standard behavior from CERNLIB. 0 - The function solve() has not been called. 1 - Test 1 was successful.
2 - Test 2 was successful.
3 - Both tests were successful.
4 - Number of iterations is greater than cern_mroot_root::maxf.
5 - Approximate (finite difference) Jacobian matrix is singular.
6 - Iterations are not making good progress.
7 - Iterations are diverging.
8 - Iterations are converging, but either cern_mroot_root::tolx is too small or the Jacobian is nearly singular or the variables are badly scaled.
9 - Either root::tolf or root::tolx is not greater than zero.
Definition at line 96 of file cern_mroot_root.h.
int maxf |
Maximum number of function evaluations.
If , then 200 (which is the CERNLIB default) is used. The default value of
maxf
is zero which then implies the default from CERNLIB.
Definition at line 105 of file cern_mroot_root.h.
double eps |
The smallest floating point number (default ).
The original prescription from CERNLIB for eps
is given below:
#if !defined(CERNLIB_DOUBLE) PARAMETER (EPS = 0.84293 69702 17878 97282 52636 392E-07) #endif #if defined(CERNLIB_IBM) PARAMETER (EPS = 0.14901 16119 38476 562D-07) #endif #if defined(CERNLIB_VAX) PARAMETER (EPS = 0.37252 90298 46191 40625D-08) #endif #if (defined(CERNLIB_UNIX))&&(defined(CERNLIB_DOUBLE)) PARAMETER (EPS = 0.14901 16119 38476 600D-07) #endif
Definition at line 134 of file cern_mroot_root.h.
Documentation generated with Doxygen and provided under the GNU Free Documentation License. See License Information for details.
Project hosting provided by
,
O2scl Sourceforge Project Page