11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
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.9 1998/02/26 15:12:26 peterg
## Removed parameter stuff
##
# Revision 1.8 1996/08/30 14:54:36 peter
# Took back lock.
#
## Revision 1.7 1996/08/24 14:23:00 peter
## Global parameter passing.
## Proper error handling.
#
|
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
|
63
64
65
66
67
68
69
70
71
72
73
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
|
+
+
+
+
-
+
-
-
+
+
+
+
+
+
+
+
+
+
+
+
|
# echo Using default parameter $PARAMS
#else
# PARAMS=$2;
#fi
$MATRIX << EOF > dm2fr_m.log 2>mtt_error.txt
# Data
$1_numpar;
# Simulation parameters
$1_simpar;
[nx,ny,nu,nz,nyz] = $1_def
% $PARAMS
if exist('W')==0
W = logspace(-1,2,100)';
if exist('mttWMIN')==0
mttWMIN = -1;
end;
if exist('mttWMAX')==0
mttWMAX = 2;
end;
if exist('mttWSTEPS')==0
mttWSTEPS = 100;
end;
W = logspace(mttWMIN,mttWMAX,mttWSTEPS)';
if exist('u0')==0
u0 = zeros(nu,1);
u0(1) = 1;
end;
|