1
/* Written by Charles Harris charles.harris@sdl.usu.edu */
3
/* Modified to not depend on Python everywhere by Travis Oliphant.
10
#define ZEROS_PARAM_HEAD int funcalls; int iterations; int error_num
16
static double dminarg1,dminarg2;
17
#define DMIN(a,b) (dminarg1=(a),dminarg2=(b),(dminarg1) < (dminarg2) ?\
18
(dminarg1) : (dminarg2))
20
#define SIGN(a) ((a) > 0.0 ? 1.0 : -1.0)
21
#define ERROR(params,num,val) (params)->error_num=(num); return (val)
25
typedef double (*callback_type)(double,void*);
26
typedef double (*solver_type)(callback_type, double, double, double, double, int,default_parameters*);
28
extern double bisect(callback_type f, double xa, double xb, double xtol, double rtol, int iter, default_parameters *params);
29
extern double ridder(callback_type f, double xa, double xb, double xtol, double rtol, int iter, default_parameters *params);
30
extern double brenth(callback_type f, double xa, double xb, double xtol, double rtol, int iter, default_parameters *params);
31
extern double brentq(callback_type f, double xa, double xb, double xtol, double rtol, int iter, default_parameters *params);
34
extern double fabs(double);
35
extern double sqrt(double);