Overview
Comment:Parameter passing added.
Downloads: Tarball | ZIP archive
Timelines: family | ancestors | descendants | both | origin/master | trunk
Files: files | file ages | folders
SHA3-256: 1f88275ca2b136ca5d8ae86f7ae4c1e351a74b848213a3cbd3ba99a5a65c8c57
User & Date: gawthrop@users.sourceforge.net on 1996-08-11 19:48:29.000
Other Links: branch diff | manifest | tags
Context
1996-08-11
19:49:27
Parameter passing (via $4) added. check-in: bf28d3a1de user: gawthrop@users.sourceforge.net tags: origin/master, trunk
19:48:29
Parameter passing added. check-in: 1f88275ca2 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
19:47:37
Changed 'exp' to 'expm' whoops! check-in: aa391be669 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
Changes
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
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







+
+
+






+
+
+
+
+
+
+
+
+



+
+
+
-
-
+
+
+

-
+




-
+

-
+





# Copyright (c) P.J.Gawthrop, 1996.

###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.2  1996/08/11 10:40:00  peter
## Now handles e and NaN numbers.
##
## Revision 1.1  1996/08/10 14:15:03  peter
## Initial revision
##
###############################################################

echo Creating $1_ir.dat

if [ "$2" = "" ]; 
then
  PARAMS='t=[0:0.1:10]'
  echo Using default parameter $PARAMS
else
  PARAMS=$2;
fi


$MATRIX << EOF | awk '$1 !~ /[A-MO-Zb-cf-z]/ && NF>0 {print} '  >$1_ir.dat

  $PARAMS;
  [n,m]=size(t);
  if m>n
  tmax=2; tstep=tmax/100;
  T = [0:tstep:tmax]';
    t=t';
  end;

  [A,B,C,D,E] = $1_dm($1_numpar);
  y = dm2ir(A,B,C,D,E,T);
  y = dm2ir(A,B,C,D,E,t);

  [Ny,Nu]=size(D);
  NN = Ny*Nu;

  xy = [T y(:,1)];
  xy = [t y(:,1)];
  for i=2:NN
    xy = [xy T y(:,i)];
    xy = [xy t y(:,i)];
  end;
xy

EOF

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
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







+
+
+









+
+
+
+
+
+
+
+
+
+



+
+
+
+
+
-
-
+
+

-
+




-
+

-
+





# Copyright (c) P.J.Gawthrop, 1996.

###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.3  1996/08/11 10:38:30  peter
## Minor debugging changes.
##
## Revision 1.2  1996/08/11 10:06:50  peter
## Now allows e format and NaN numbers.
##
## Revision 1.1  1996/08/11 09:45:46  peter
## Initial revision
##
###############################################################

echo Creating $1_sr.dat


if [ "$2" = "" ]; 
then
  PARAMS='t=[0:0.1:10]'
  echo Using default parameter $PARAMS
else
  PARAMS=$2;
fi


$MATRIX << EOF  | awk '$1 !~ /[A-MO-Zb-cf-z]/ && NF>0 {print} '  >$1_sr.dat

  $PARAMS;
  [n,m]=size(t);
  if m>n
    t=t';
  end;
  tmax=10.0; tstep=tmax/100;
  T = [0:tstep:tmax]';


  [A,B,C,D,E] = $1_dm($1_numpar);
  y = dm2sr(A,B,C,D,E,T);
  y = dm2sr(A,B,C,D,E,t);

  [Ny,Nu]=size(D);
  NN = Ny*Nu;

  xy = [T y(:,1)];
  xy = [t y(:,1)];
  for i=2:NN
    xy = [xy T y(:,i)];
    xy = [xy t y(:,i)];
  end;
xy

EOF


MTT: Model Transformation Tools
GitHub | SourceHut | Sourceforge | Fossil RSS ]