16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
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
|
+
+
+
-
+
|
# Copyright (c) P.J.Gawthrop 1997.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
# Revision 1.2 1997/01/21 22:57:17 peterg
# Various bug fixes.
#
## Revision 1.1 1997/01/21 10:52:23 peterg
## Initial revision
##
###############################################################
# Inform user
echo Creating $1_ode.c
# Remove the old files
rm -f $1_ode.c1 $1_ode.c2 $1_ode.c3 $1_ode.c
rm -f $1_ode.c1 $1_ode.c2 $1_ode.c
# Remove the old log file
rm -f ode_r2c.log
# Use reduce to accomplish the transformation
reduce >ode_r2c.log << EOF
|
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
+
|
ON BigFloat, NumVal;
PRECISION 16; %Compatible with Matlab
%OFF Nat;
ON NERO; % Suppress zero elements
%Generate the Header part
OUT "$1_ode.c1";
write "/*"$
write "Differential algebraic eqns in c form for system $1"$
write "NB Arrays should be defined to be one larger than expected"$
|
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
119
120
121
122
123
124
125
126
127
128
129
130
131
132
|
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
119
120
121
122
|
-
-
-
-
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
-
+
-
+
-
+
+
+
|
write "void $1_ode()"$
write " "$
write "{"$
write "/* Declare standard arrays */"$
write " extern float y[", MTTNy+1, "]; /* $1_ode output */"$
write " extern float dx[", MTTNx+1, "]; /* $1_ode state derivative */"$
write " extern float x[", MTTNx+1, "]; /* $1_ode state */"$
write " extern float u[", MTTNu+1, "]; /* $1_ode input */"$
write " extern double y[", MTTNy+1, "]; /* $1_ode output */"$
write " extern double dx[", MTTNx+1, "]; /* $1_ode state derivative */"$
write " extern double x[", MTTNx+1, "]; /* $1_ode state */"$
write " extern double u[", MTTNu+1, "]; /* $1_ode input */"$
%External (global) variable list
write "/* External (global) variable list */ "$
IF MTTNvar>0 THEN
BEGIN
FOR i := 1:MTTNvar DO
IF numberp(MTTVar(i,1))
THEN
BEGIN
% Do nowt
END
ELSE
BEGIN
write "extern float ", MTTVar(i,1), ";"$
END$
END$
write " float t0;"$
write " double t0;"$
FOR i := 1:9 DO
BEGIN
write " float t", i, ";"$
write " double t", i, ";"$
END$
%Declarations$
write "/* State variable list */ "$
FOR i := 1:MTTNx DO
BEGIN
write " float mttx", i, ";"$
write " double mttx", i, ";"$
END$
write "/* Input variable list */ "$
FOR i := 1:MTTNu DO
BEGIN
write " float mttu", i, ";"$
write " double mttu", i, ";"$
END$
write "/* Counter */ "$
write " int i;"$
write "/* Parameter list */ "$
write "#include ", """$1_sympar.c"""$
write " "$
write " /*====== Set up the state variables ======*/"$
FOR i := 1:MTTNx DO
BEGIN
write " mttx", i, " = x[", i, "];"$
END$
|