Overview
| Comment: | No change |
|---|---|
| Downloads: | Tarball | ZIP archive |
| Timelines: | family | ancestors | descendants | both | origin/master | trunk |
| Files: | files | file ages | folders |
| SHA3-256: |
f271eec1cf5a652b210d006aae57fc21 |
| User & Date: | gawthrop@users.sourceforge.net on 1998-06-27 14:49:16.000 |
| Other Links: | branch diff | manifest | tags |
Context
|
1998-06-27
| ||
| 15:06:40 | Initial revision check-in: d740da803b user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
| 14:49:16 | No change check-in: f271eec1cf user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
| 13:24:04 |
Causality now set correctly for: multi-port C and I C and I with arrows pointing in check-in: 074b696b35 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
Changes
Modified mttroot/mtt/bin/trans/ode2odes_m
from [0c76dc5319]
to [e2fc3e22e9].
| ︙ | ︙ | |||
11 12 13 14 15 16 17 18 19 20 21 22 23 24 | # Copyright (c) P.J.Gawthrop, 1996. ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ ## Revision 1.15 1998/05/21 16:20:27 peterg ## Modified to include explicit algebraic loop solution ## ## Revision 1.14 1998/05/19 19:48:02 peterg ## Read the simpar file now. ## ## Revision 1.13 1998/05/14 08:05:10 peterg | > > > | 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 | # Copyright (c) P.J.Gawthrop, 1996. ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ ## Revision 1.16 1998/06/25 08:47:23 peterg ## Put correct arguments for _input -- (x,t) ## ## Revision 1.15 1998/05/21 16:20:27 peterg ## Modified to include explicit algebraic loop solution ## ## Revision 1.14 1998/05/19 19:48:02 peterg ## Read the simpar file now. ## ## Revision 1.13 1998/05/14 08:05:10 peterg |
| ︙ | ︙ | |||
79 80 81 82 83 84 85 | #else # PARAMS=$2; #fi # PARAMS="$PARAMS ;" | | > | | 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 | #else # PARAMS=$2; #fi # PARAMS="$PARAMS ;" $MATRIX --verbose << EOF > ode2odes_m.log 2>mtt_error.txt %System structure [nx,ny,nu,nz,nyz] = $1_def; %Read in parameters $1_numpar; %Read in state x = $1_state; %Set the initial output %if ny>0 % y = $1_odeo(x,0); %end; |
| ︙ | ︙ | |||
129 130 131 132 133 134 135 136 137 138 139 140 141 142 |
method = tolower(METHOD)
if nx>0
if strcmp(method,'lsode')
X = lsode('$1_ode', x, T);
elseif strcmp(method,'euler')
%Euler integration
X=[];
dt = (T(2)-T(1))/STEPFACTOR;
for t=T'
X = [X; xx'];
ts = t;
for i=1:STEPFACTOR
x = xx(1:nx);
| > | 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 |
method = tolower(METHOD)
if nx>0
if strcmp(method,'lsode')
X = lsode('$1_ode', x, T);
elseif strcmp(method,'euler')
%Euler integration
disp("Euler")
X=[];
dt = (T(2)-T(1))/STEPFACTOR;
for t=T'
X = [X; xx'];
ts = t;
for i=1:STEPFACTOR
x = xx(1:nx);
|
| ︙ | ︙ | |||
193 194 195 196 197 198 199 | else X = zeros(size(T)); end; if ny>0 % compute y and print it i = 0; Y=[]; for t=T' | | | | 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 |
else
X = zeros(size(T));
end;
if ny>0 % compute y and print it
i = 0; Y=[];
for t=T'
i = i+1; X(i,:);
y = $1_odeo(X(i,:)',t);
Y = [Y; y'];
end;
write_matrix([T,Y], '$1_odeso');
end;
EOF
|
| ︙ | ︙ |