diff options
Diffstat (limited to 'macros/Misc/scifunc_block.sci')
| -rw-r--r-- | macros/Misc/scifunc_block.sci | 161 |
1 files changed, 161 insertions, 0 deletions
diff --git a/macros/Misc/scifunc_block.sci b/macros/Misc/scifunc_block.sci new file mode 100644 index 00000000..6407be95 --- /dev/null +++ b/macros/Misc/scifunc_block.sci @@ -0,0 +1,161 @@ +// Scicos +// +// Copyright (C) INRIA - METALAU Project <scicos@inria.fr> +// +// This program is free software; you can redistribute it and/or modify +// it under the terms of the GNU General Public License as published by +// the Free Software Foundation; either version 2 of the License, or +// (at your option) any later version. +// +// This program is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +// GNU General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this program; if not, write to the Free Software +// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. +// +// See the file ../license.txt +// + +function [x,y,typ]=scifunc_block(job,arg1,arg2) + //%Description + // job=='set' : block parameters acquisition + // arg1 is block data structure + // x is returned block data structure + // job=='define' : corresponding block data structure initialisation + // arg1: name of block parameters acquisition macro + // x : block data structure + //%Block data-structure definition + // bl=list('Block',graphics,model,init,'standard_block') + // graphics=list([xo,yo],[l,h],orient,label) + // xo - x coordinate of block origin + // yo - y coordinate of block origin + // l - block width + // h - block height + // orient - boolean, specifies if block is tilded + // label - string block label + // model=list(eqns,#input,#output,#clk_input,#clk_output,state,.. + // rpar,ipar,typ [,firing]) + // eqns - function name (in string form if fortran routine) + // #input - vector of input port sizes + // #output - vector of output port sizes + // #clk_input - vector of clock inputs port sizes + // #clk_output - vector of clock output port sizes + // state - vector (column) of initial condition + // rpar - vector (column) of real parameters + // ipar - vector (column) of integer parameters + // typ - string: 'c' if block is continuous, 'd' if discrete + // 'z' if zero-crossing. + // firing - vector of initial output event firing times + // + x=[]; + y=[]; + typ=[]; + select job + case "set" then + needcompile=0 + x=arg1 + model=arg1.model; + graphics=arg1.graphics; + exprs=graphics.exprs + + if size(exprs(1),"*")==8 then + exprs(1)(9)="0"; + end + while %t do + [ok,i,o,ci,co,xx,z,rpar,auto0,deptime,lab]=scicos_getvalue(.. + ["Set scifunc_block parameters";"only regular blocks supported"],.. + ["input ports sizes"; + "output port sizes"; + "input event ports sizes"; + "output events ports sizes"; + "initial continuous state"; + "initial discrete state"; + "System parameters vector"; + "initial firing vector (<0 for no firing)"; + "is block always active (0:no, 1:yes)" ],.. + list("vec",-1,"vec",-1,"vec",-1,"vec",-1,"vec",-1,"vec",-1,.. + "vec",-1,"vec","sum(%4)","vec",1),exprs(1)) + if ~ok then + break, + end + exprs(1)=lab + xx=xx(:); + z=z(:); + rpar=rpar(:) + nrp=prod(size(rpar)) + // create simulator + i=int(i(:)); + ni=size(i,1); + o=int(o(:)); + no=size(o,1); + ci=int(ci(:)); + nci=size(ci,1); + co=int(co(:)); + nco=size(co,1); + [ok,tt,dep_ut]=genfunc1(exprs(2),i,o,nci,nco,size(xx,1),size(z,1),.. + nrp,"c") + dep_ut(2)=(1==deptime) + if ~ok then + break, + end + [model,graphics,ok]=check_io(model,graphics,i,o,ci,co) + if ok then + auto=auto0 + model.state=xx + model.dstate=z + model.rpar=rpar + if model.ipar <> 0 then + model.opar=model.ipar; + model.ipar=0; + end + if or(model.opar<>tt) then + needcompile=4, + end + model.opar=tt + model.firing=auto + model.dep_ut=dep_ut + x.model=model + exprs(2)=tt + graphics.exprs=exprs + x.graphics=graphics + break + end + end + needcompile=resume(needcompile) + case "define" then + in=1 + out=1 + clkin=[] + clkout=[] + x0=[] + z0=[] + typ="c" + auto=[] + rpar=[] + + model=scicos_model() + model.sim=list("scifunc",3) + model.in=in + model.out=out + model.evtin=clkin + model.evtout=clkout + model.state=x0 + model.dstate=z0 + model.rpar=rpar + model.ipar=0 + model.opar=list() + model.blocktype=typ + model.firing=auto + model.dep_ut=[%t %f] + + exprs=list([sci2exp(in);sci2exp(out);sci2exp(clkin);sci2exp(clkout); + strcat(sci2exp(x0));strcat(sci2exp(z0)); + strcat(sci2exp(rpar));sci2exp(auto)],.. + list("y1=sin(u1)"," "," ","y1=sin(u1)"," "," "," ")) + gr_i=[] + x=standard_define([2 2],model,exprs,gr_i) + end +endfunction |