1
2
3
4
5
6
7
8
9
10
11
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
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
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
|
function Info=scicos_simulate(scs_m,Info,%scicos_context,flag)
// Function for running scicos simulation in batch mode
// Info=scicos_simulate(scs_m,Info[,%scicos_context][,flag])
//
// scs_m: scicos diagram (obtained by "load file.cos"). Note that
// the version of file.cos must be the current version. If not, load
// into scicos and save.
//
// %scicos_context: a scilab struct containing values of
// symbolic variables used in the context and Scicos blocks. This
// is often used to change a parameter in the diagram context. In that
// case, make sure that in the diagram context the variable is defined such
// that it can be modified. Say a variable "a" is to be defined in the
// context having value 1, and later in batch mode, we want to change
// the value of "a". In that case, in the context of the diagram place:
// if ~exists('a') then a=1,end
// If you want then to run the simulation in batch mode using the value
// a=2, set:
// %scicos_context.a=2
//
// Info: a list. It must be list() at the first call, then use output
// Info as input Info for the next calls. Info contains compilation and
// simulation information and is used to avoid recompilation when not
// needed.
//
// flag: string. If it equals 'nw' (no window), then blocks using
// graphical windows are not executed. Note that the list of such
// blocks must be updated as new blocks are added.
//
//
// list of blocks to ignore (blocks using graphics) in nw mode
Ignoreb=['cscope','cmscope','scope','mscope','scopexy','evscpe','affich']
//
load SCI/macros/scicos/lib
exec(loadpallibs,-1)
//redefine gui functions
prot=funcprot();funcprot(0);
deff('disablemenus()',' ')
deff('enablemenus()',' ')
do_terminate=do_terminate1
funcprot(prot)
if argn(2)==3 then
if type(%scicos_context)==10&(stripblanks(%scicos_context)=='nw') then
Ignore=Ignoreb
%scicos_context=struct()
else
Ignore=[]
end
elseif argn(2)==2 then
Ignore=[]
%scicos_context=struct()
elseif stripblanks(flag)=='nw' then
Ignore=Ignoreb
else
Ignore=[]
end
//
if Info<>list() then
[%tcur,%cpr,alreadyran,needstart,needcompile,%state0]=Info(:)
else
%tcur=0;%cpr=list();alreadyran=%f;needstart=%t;needcompile=4;%state0=list();
end
//
tolerances=scs_m.props.tol
solver=tolerances(6)
[%scicos_context,ierr]=script2var(scs_m.props.context, ...
%scicos_context);
if ierr==0 then
[scs_m,%cpr,needcompile,ok]=do_eval(scs_m,%cpr)
if needcompile<>4&size(%cpr)>0 then %state0=%cpr.state,end
alreadyran=%f
else
error(['Incorrect context definition, ';lasterror()] )
end
if %cpr==list() then need_suppress=%t, else need_suppress=%f,end
[%cpr,%state0_n,needcompile,alreadyran,ok]=..
do_update(%cpr,%state0,needcompile)
if ~ok then error('Error updating parameters.'),end
if or(%state0_n<>%state0) then //initial state has been changed
%state0=%state0_n
[alreadyran,%cpr]=do_terminate1(scs_m,%cpr)
choix=[]
end
if %cpr.sim.xptr($)-1<size(%cpr.state.x,'*') & solver<100 then
warning(['Diagram has been compiled for implicit solver'
'switching to implicit Solver'])
solver=100
tolerances(6)=solver
elseif (%cpr.sim.xptr($)-1==size(%cpr.state.x,'*')) & ..
( solver==100 & size(%cpr.state.x,'*')<>0) then
warning(['Diagram has been compiled for explicit solver'
'switching to explicit Solver'])
solver=0
tolerances(6)=solver
end
if need_suppress then //this is done only once
for i=1:length(%cpr.sim.funs)
if type(%cpr.sim.funs(i))<>13 then
if find(%cpr.sim.funs(i)(1)==Ignore)<>[] then
%cpr.sim.funs(i)(1)='trash';
end
end
end
end
switch_to_old_graphics_style()
if needstart then //scicos initialisation
if alreadyran then
[alreadyran,%cpr]=do_terminate1(scs_m,%cpr)
alreadyran=%f
end
%tcur=0
%cpr.state=%state0
tf=scs_m.props.tf;
if tf*tolerances==[] then
error(['Simulation parameters not set']);
end
ierr=execstr('[state,t]=scicosim(%cpr.state,%tcur,tf,%cpr.sim,'+..
'''start'',tolerances)','errcatch')
%cpr.state=state
if ierr<>0 then
restore_graphics_style()
error(['Initialisation problem:'])
end
end
ierr=execstr('[state,t]=scicosim(%cpr.state,%tcur,tf,%cpr.sim,'+..
'''run'',tolerances)','errcatch')
%cpr.state=state
if ierr==0 then
alreadyran=%t
if tf-t<tolerances(3) then
needstart=%t
[alreadyran,%cpr]=do_terminate1(scs_m,%cpr)
else
%tcur=t
end
else
restore_graphics_style()
error(['Simulation problem:';lasterror()])
end
restore_graphics_style()
Info=list(%tcur,%cpr,alreadyran,needstart,needcompile,%state0)
endfunction
function restore_graphics_style()
global bak
if Ignore==[] then
//gg=xget('window') // for bug in figure_style and winsid
//xset('window',0) // for bug in figure_style and winsid
set('old_style',bak) //set('figure_style',bak)
//xset('window',gg) // for bug in figure_style and winsid
end
endfunction
function switch_to_old_graphics_style()
global bak
if Ignore==[] then
bak=stripblanks(get("old_style")) //bak=get('figure_style')
set("old_style","on")
//set('figure_style','old')
end
endfunction
function [alreadyran,%cpr]=do_terminate1(scs_m,%cpr)
// Copyright INRIA
if prod(size(%cpr))<2 then alreadyran=%f,return,end
par=scs_m.props;
if alreadyran then
alreadyran=%f
//terminate current simulation
ierr=execstr('[state,t]=scicosim(%cpr.state,par.tf,par.tf,'+..
'%cpr.sim,''finish'',par.tol)','errcatch')
%cpr.state=state
if ierr<>0 then
error(['End problem:';lasterror()])
end
end
endfunction
|
