1
2
3
4
5
6
7
8
9
10
11
12
|
#! /bin/sh
######################################
##### Model Transformation Tools #####
######################################
# Bourne shell script: make_ode2odes
# Copyright (c) P.J.Gawthrop July 1998.
Sys=$1
|
>
>
>
>
>
>
>
>
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
#! /bin/sh
######################################
##### Model Transformation Tools #####
######################################
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
###############################################################
# Bourne shell script: make_ode2odes
# Copyright (c) P.J.Gawthrop July 1998.
Sys=$1
|
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
|
MTTu = $1_input(MTTx,t); # Evaluate initial input
MTTy = $1_odeo(MTTx,MTTu,t); # Evaluate initial output
mtt_write(t,MTTx,MTTy,$Nx,$Ny); # And write them
#Integration loop
for it = 1:iLast
for i = 1:STEPFACTOR
MTTu = $1_input(MTTx,t);# Input
MTTx = $1_switch(MTTx); # Switches
MTTdx = $1_ode(MTTx,MTTu,t);
for j = 1:$Nx
MTTx(j) = MTTx(j) + MTTdx(j)*DDT;
end;
MTTy = $1_odeo(MTTx,MTTu,t);
|
|
|
|
|
<
<
|
|
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
MTTu = $1_input(MTTx,t); # Evaluate initial input
MTTy = $1_odeo(MTTx,MTTu,t); # Evaluate initial output
mtt_write(t,MTTx,MTTy,$Nx,$Ny); # And write them
#Integration loop
for it = 1:iLast
for i = 1:STEPFACTOR
[MTTu] = $1_input(MTTx,t);# Input
[MTTx] = $1_switch(MTTx); # Switches
[MTTdx] = $1_ode(MTTx,MTTu,t);
[MTTx] = mtt_update(MTTdx,MTTx,DDT,$Nx,METHOD);
[MTTy] = $1_odeo(MTTx,MTTu,t);
t = t + DDT;
end;
mtt_write(t,MTTx,MTTy,$Nx,$Ny);
end;
EOF
|