| ︙ | | | ︙ | |
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
# Copyright (c) P.J.Gawthrop 1997.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
# Revision 1.12 1997/05/15 08:39:56 peterg
# Don't initialise states - now done in numpar file.
#
# Revision 1.11 1997/05/12 16:00:54 peterg
# Removed itime again.,
#
# Revision 1.10 1997/05/10 10:05:15 peterg
|
>
>
>
|
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
# Copyright (c) P.J.Gawthrop 1997.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.13 1998/02/24 13:34:45 peterg
## Back under RCS
##
# Revision 1.12 1997/05/15 08:39:56 peterg
# Don't initialise states - now done in numpar file.
#
# Revision 1.11 1997/05/12 16:00:54 peterg
# Removed itime again.,
#
# Revision 1.10 1997/05/10 10:05:15 peterg
|
| ︙ | | | ︙ | |
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
|
%Read the reduce definitions file
in "$1_def.r";
%Set up the number of argument variables to zero in case the user has forgotten
MTTNVar := 0;
%Read the symbolic parameters file
in "$1_sympar.r";
ON BigFloat, NumVal;
PRECISION 16; %Compatible with Matlab
%OFF Nat;
ON NERO; % Suppress zero elements
|
|
|
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
|
%Read the reduce definitions file
in "$1_def.r";
%Set up the number of argument variables to zero in case the user has forgotten
MTTNVar := 0;
%Read the symbolic parameters file
%%in "$1_sympar.r";
ON BigFloat, NumVal;
PRECISION 16; %Compatible with Matlab
%OFF Nat;
ON NERO; % Suppress zero elements
|
| ︙ | | | ︙ | |
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
write "#define MTTNY1 ", MTTNy+1 $
write "#define MTTNU1 ", MTTNu+1 $
write "#define MTTNX2 ", MTTNx+2 $
write "#define MTTNY2 ", MTTNy+2 $
write "#define MTTNU2 ", MTTNu+2 $
write "#include <stdio.h>"$
write "#include ""$1_odes.h"" "$
write "#include ""$1_ode.c"" "$
write "#include ""$1_input.c"" "$
write "#include ""$1_numpar.c"" "$
%External (global) variable list
write "#include ""$1_sympar.h"" "$
write "/* Declare standard arrays */"$
write "double y[MTTNY1]; /* $1_ode output */"$
|
|
>
|
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
|
write "#define MTTNY1 ", MTTNy+1 $
write "#define MTTNU1 ", MTTNu+1 $
write "#define MTTNX2 ", MTTNx+2 $
write "#define MTTNY2 ", MTTNy+2 $
write "#define MTTNU2 ", MTTNu+2 $
write "#include <stdio.h>"$
write "#include ""$1_simpar.h"" "$
write "#include ""$1_ode.c"" "$
write "#include ""$1_input.c"" "$
write "#include ""$1_numpar.c"" "$
write "#include ""$1_state.c"" "$
%External (global) variable list
write "#include ""$1_sympar.h"" "$
write "/* Declare standard arrays */"$
write "double y[MTTNY1]; /* $1_ode output */"$
|
| ︙ | | | ︙ | |
152
153
154
155
156
157
158
159
160
161
162
163
164
165
|
write "fps = fopen(""$1_odes.m"", ""w""); "$
write "fpso = fopen(""$1_odeso.m"", ""w""); "$
%Set up user-defined constants
write "/* Set up user-defined constants */"$
write " $1_numpar();"$
%Initialise main (Euler) integration loop
write "/* Initialise main (Euler) integration loop */"$
write " time = 0.0;"$
write " dt = DT/STEPFACTOR;"$
%% write " for (i=1; i<=MTTNX; i++)"$
%% write " x[i] = 0.0;"$
|
>
>
>
>
|
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
|
write "fps = fopen(""$1_odes.m"", ""w""); "$
write "fpso = fopen(""$1_odeso.m"", ""w""); "$
%Set up user-defined constants
write "/* Set up user-defined constants */"$
write " $1_numpar();"$
%Set up initial state
write "/* Set up initial state */"$
write " $1_state();"$
%Initialise main (Euler) integration loop
write "/* Initialise main (Euler) integration loop */"$
write " time = 0.0;"$
write " dt = DT/STEPFACTOR;"$
%% write " for (i=1; i<=MTTNX; i++)"$
%% write " x[i] = 0.0;"$
|
| ︙ | | | ︙ | |