Overview
Comment:Save old version using out files
Downloads: Tarball | ZIP archive
Timelines: family | ancestors | descendants | both | origin/master | trunk
Files: files | file ages | folders
SHA3-256: b9d71146321e03796ecd7af1af4299ee434b3d3ceebe4b1c6b6e691b6d4e1d9f
User & Date: gawthrop@users.sourceforge.net on 2000-04-11 11:04:08.000
Other Links: branch diff | manifest | tags
Context
2000-04-11
12:22:45
*** empty log message *** check-in: 9b227e22df user: gawthrop@users.sourceforge.net tags: origin/master, trunk
11:04:08
Save old version using out files check-in: b9d7114632 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
2000-04-10
12:08:08
removed echos check-in: 24bf94ae83 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
Changes
Side-by-Side Diff Ignore Whitespace Patch
Modified mttroot/mtt/bin/trans/m/mtt_stime.m from [3cb31be700] to [57b5993ee0]. 
1
2


3
4
5
6
7
8
9



1
2
3
4
5
6
7
8
9

-
-
+
+








function [t,y,y_theta,x] = mtt_stime(system_name,theta,free);
  ## usage: [t,y,y_theta] = mtt_stime(system_name,theta);
function [t,y,y_theta,x,x_last] = mtt_stime(system_name,theta,free);
  ## usage: [t,y,y_theta,x,x_last] = mtt_stime(system_name,theta);
  ##
  ## Simulate system with name system_name and parameter vector theta
  ## The order of components in theta is determined in system_numpar.txt:
  ## y_theta contains the corresponding sensitivity functions
  ## Assumes system generated by the sBG approach
  ## Copyright (C) 1999 by Peter J. Gawthrop
︙
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

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








+
+
+
+
+
-
-
+
+
+
+


-
-
-
+
+
+


-
-
-
-
-
-
+
+
+
+
+
+





  for i=free
    args=sprintf("%i",i);
    for j=1:length(theta)
      args = sprintf("%s %g",args, theta(j));
    endfor

    ## Run system and replace NaN by 1e30 -- easier to handle
    file_name = sprintf("%s_input.dat", system_name);
    if exist(file_name)==2	# Then use data from this file ...
      command = sprintf("./%s_ode2odes.out %s < %s | sed \'s/NAN/Inf/g\' >mtt_out_data.dat\n", \
			system_name, args, file_name);
    else
    command = sprintf("./%s_ode2odes.out %s | sed \'s/NAN/1e30/g\' >mtt_data.dat\n", \
    system_name, args);
      command = sprintf("./%s_ode2odes.out %s | sed \'s/NAN/Inf/g\' >mtt_out_data.dat\n", \
			system_name, args);
    endif
      
   system(command);

    ## Retrieve data
    load -force mtt_data.dat
    y_theta = [y_theta mtt_data(:,3:2:1+ny)];
    ## Retrieve out_data
    load -force mtt_out_data.dat
    y_theta = [y_theta mtt_out_data(:,3:2:1+ny)];
  endfor

  ## System data
  [n,m]=size(mtt_data);
  t = mtt_data(:,1);
  y = mtt_data(:,2:2:ny);
  x = mtt_data(:,3+ny:m);

  ## System out_data
  [n,m]=size(mtt_out_data);
  t = mtt_out_data(:,1);
  y = mtt_out_data(:,2:2:ny);
  x = mtt_out_data(:,3+ny:m);
  x_last = mtt_out_data(n,3+ny:m);
endfunction
MTT: Model Transformation Tools
[ GitHub | SourceHut | Sourceforge | Fossil RSS ]