34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
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/Inf/g\' >mtt_out_data.dat\n", \
system_name, args);
endif
system(command);
## Retrieve out_data
load -force mtt_out_data.dat
|
|
|
|
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
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/Inf/g\' >mtt_out_data.dat\n", ...
system_name, args);
endif
system(command);
## Retrieve out_data
load -force mtt_out_data.dat
|