diff options
Diffstat (limited to 'macros/Misc/generic_block2.sci')
| -rw-r--r-- | macros/Misc/generic_block2.sci | 148 |
1 files changed, 148 insertions, 0 deletions
diff --git a/macros/Misc/generic_block2.sci b/macros/Misc/generic_block2.sci new file mode 100644 index 00000000..ae05dc2b --- /dev/null +++ b/macros/Misc/generic_block2.sci @@ -0,0 +1,148 @@ +// 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]=generic_block2(job,arg1,arg2) + // + x=[]; + y=[]; + typ=[]; + select job + case "set" then + x=arg1 + model=arg1.model; + graphics=arg1.graphics; + label=graphics.exprs + if size(label,"*")==14 then + label(9)=[], + end //compatiblity + while %t do + [ok,function_name,funtyp,i,o,ci,co,xx,z,rpar,ipar,nmode,nzcr,auto0,depu,dept,lab]=.. + scicos_getvalue("Set GENERIC block parameters",.. + ["simulation function"; + "function type (0,1,2,..)"; + "input ports sizes"; + "output port sizes"; + "input event ports sizes"; + "output events ports sizes"; + "initial continuous state"; + "initial discrete state"; + "Real parameters vector"; + "Integer parameters vector"; + "number of modes"; + "number of zero_crossings"; + "initial firing vector (<0 for no firing)"; + "direct feedthrough (y or n)"; + "time dependence (y or n)"],.. + list("str",1,"vec",1,"vec",-1,"vec",-1,"vec",-1,"vec",-1,.. + "vec",-1,"vec",-1,"vec",-1,"vec",-1,"vec",1,"vec",1,"vec","sum(%6)",.. + "str",1,"str",1),label) + if ~ok then + break, + end + label=lab + function_name=stripblanks(function_name) + xx=xx(:); + z=z(:); + rpar=rpar(:); + ipar=int(ipar(:)); + i=int(i(:)); + o=int(o(:)); + ci=int(ci(:)); + co=int(co(:)); + funtyp=int(funtyp) + if funtyp<0 then + message("function type cannot be negative"); + ok=%f; + end + if [ci;co]<>[] then + if max([ci;co])>1 then + message("vector event links not supported"); + ok=%f; + end + end + depu=stripblanks(depu); + if part(depu,1)=="y" then + depu=%t; + else + depu=%f; + end + dept=stripblanks(dept); + if part(dept,1)=="y" then + dept=%t; + else + dept=%f; + end + dep_ut=[depu dept]; + if ok then + [model,graphics,ok]=check_io(model,graphics,i,o,ci,co) + end + if ok then + // AVERIFIER + if funtyp==3 then + needcompile=4; + end + model.sim=list(function_name,funtyp); + model.state=xx + model.dstate=z + model.rpar=rpar + model.ipar=ipar + // needcompile=4 AVERIFIER CANEMARCHEQUAVECFORTRAN + model.firing=auto0 + model.nzcross=nzcr + model.nmode=nmode + model.dep_ut=dep_ut + arg1.model=model + graphics.exprs=label + arg1.graphics=graphics + x=arg1 + break + end + end + needcompile=resume(needcompile) + case "define" then + model=scicos_model() + function_name="sinblk"; + funtyp=1; + model.sim=list(function_name,funtyp) + + model.in=1 + model.out=1 + model.evtin=[] + model.evtout=[] + model.state=[] + model.dstate=[] + model.rpar=[] + model.ipar=[] + model.blocktype="c" + model.firing=[] + model.dep_ut=[%t %f] + label=[function_name;sci2exp(funtyp); + sci2exp(model.in);sci2exp(model.out); + sci2exp(model.evtin);sci2exp(model.evtout); + sci2exp(model.state);sci2exp(model.dstate); + sci2exp(model.rpar);sci2exp(model.ipar); + sci2exp(model.nmode);sci2exp(model.nzcross); + sci2exp(model.firing);"y";"n"]; + gr_i=[] + x=standard_define([2 2],model,label,gr_i) + end +endfunction + |