9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
+
+
+
|
# Copyright (C) 2000 by Peter J. Gawthrop
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.4 2000/05/16 18:57:15 peterg
## Still debugging
##
## Revision 1.3 2000/05/16 11:59:34 peterg
## Stard new version with data files not argv.
##
## Revision 1.2 2000/05/11 19:32:29 peterg
## Put in c version + sensitivity computation
##
## Revision 1.1 2000/04/08 10:43:26 peterg
|
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
|
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
|
-
+
|
if doing_state
i_2 = i_1 + n_y - 1;
i_3 = i_1 + n_y + 1;
i_4 = i_1 + n_y + $Nx;
else
if doing_sensitivities
i_2 = 2 + n_y;
i_2 = 1 + 2*n_y;
ipars = [ipar; [sensitivities(i_sim) 1]];
save -ascii $1_numpar.dat ipars
else
i_2 = i_1 + n_y -1;
endif
endif;
|