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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
/*
* preprocessor macros and structure declarations for ode.
* Copyright Nicholas B. Tufillaro, 1982-1994. All rights reserved.
* GNU enhancements copyright (C) 1996-1997 Free Software Foundation, Inc.
*/
#include <setjmp.h>
#define NAMMAX 32 /* size of identifiers */
#define KMAX 7 /* size of k vector */
#define PASTMAX 7 /* size of history vector */
/*
* symbol table entry
*/
struct sym
{
char sy_name[NAMMAX];
double sy_value;
double sy_val[PASTMAX];
double sy_prime;
double sy_pri[PASTMAX];
double sy_predi; /* predictor value */
double sy_sserr; /* relative single step error */
double sy_aberr; /* absolute single step error */
double sy_acerr; /* accumulated error */
double sy_k[KMAX];
int sy_flags;
struct expr *sy_expr;
struct sym *sy_link;
};
/*
* bits in the flags word
* of a symbol table entry
*/
#define SF_INIT 1 /* set when the initial value is given */
#define SF_ISEQN 2 /* set when an ODE is given */
#define SF_DEPV (SF_INIT|SF_ISEQN)
/*
* enumeration of printable cells in
* a symbol table entry
*/
typedef enum
{
P_VALUE, P_PRIME, P_SSERR, P_ACERR, P_ABERR
} ent_cell;
/*
* the print queue is made
* of these
*/
struct prt
{
ent_cell pr_which; /* which cell to print */
struct prt *pr_link;
struct sym *pr_sym;
};
/*
* yylex returns a pointer to this.
* The receiving yacc action must
* free the space when it's done with it.
*/
struct lex
{
union
{
double lxu_value;
char lxu_name[NAMMAX];
} lx_u;
int lx_lino;
};
/*
* operation types
*/
typedef enum
{
O_NOOP, O_PLUS, O_MINUS, O_MULT, O_DIV, O_POWER, O_SQRT, O_EXP, O_LOG,
O_LOG10, O_SIN, O_COS, O_TAN, O_ASIN, O_ACOS, O_ATAN, O_IDENT, O_CONST,
O_NEG, O_ABS, O_SINH, O_COSH, O_TANH, O_ASINH, O_ACOSH, O_ATANH, O_SQAR,
O_CUBE, O_INV, O_FLOOR, O_CEIL, O_J0, O_J1, O_Y0, O_Y1, O_ERF, O_ERFC,
O_INVERF, O_LGAMMA, O_GAMMA, O_NORM, O_INVNORM, O_IGAMMA, O_IBETA
} op_type;
/*
* an operation in an expression list
*/
struct expr
{
op_type ex_oper;
double ex_value;
struct sym *ex_sym;
struct expr *ex_next;
};
/* integration algorithm type */
typedef enum
{
A_EULER, A_ADAMS_MOULTON, A_RUNGE_KUTTA_FEHLBERG
} integration_type;
/*
* misc. consts
*/
#define TWO (2.0)
#define HALF (0.5)
#define SCALE (1e-4)
#define HMIN (1e-36) /* Minimum step size */
#define HMAX (0.5) /* Maximum step size */
#define LONGMAX (2147483647) /* 2^31-1 */
#define MESH (2)
#define STOPR (tstep>0 ? \
t-0.0625*tstep>=tstop : t-0.0625*tstep<=tstop )
#define STOPA (tstep>0 ? \
t+0.9375*tstep>=tstop : t+0.9375*tstep<=tstop )
|
