1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
/* Written by Charles Harris charles.harris@sdl.usu.edu */
/* Modified to not depend on Python everywhere by Travis Oliphant.
*/
#ifndef ZEROS_H
#define ZEROS_H
#define ZEROS_PARAM_HEAD int funcalls; int iterations; int error_num
typedef struct {
ZEROS_PARAM_HEAD;
} default_parameters;
static double dminarg1,dminarg2;
#define DMIN(a,b) (dminarg1=(a),dminarg2=(b),(dminarg1) < (dminarg2) ?\
(dminarg1) : (dminarg2))
#define SIGN(a) ((a) > 0.0 ? 1.0 : -1.0)
#define ERROR(params,num,val) (params)->error_num=(num); return (val)
#define SIGNERR -1
#define CONVERR -2
typedef double (*callback_type)(double,void*);
typedef double (*solver_type)(callback_type, double, double, double, double, int,default_parameters*);
extern double bisect(callback_type f, double xa, double xb, double xtol, double rtol, int iter, default_parameters *params);
extern double ridder(callback_type f, double xa, double xb, double xtol, double rtol, int iter, default_parameters *params);
extern double brenth(callback_type f, double xa, double xb, double xtol, double rtol, int iter, default_parameters *params);
extern double brentq(callback_type f, double xa, double xb, double xtol, double rtol, int iter, default_parameters *params);
extern double fabs(double);
extern double sqrt(double);
#endif
|
