diff options
Diffstat (limited to 'macros/Misc/EXPRESSION.sci')
| -rw-r--r-- | macros/Misc/EXPRESSION.sci | 232 |
1 files changed, 232 insertions, 0 deletions
diff --git a/macros/Misc/EXPRESSION.sci b/macros/Misc/EXPRESSION.sci new file mode 100644 index 00000000..aaf1fd24 --- /dev/null +++ b/macros/Misc/EXPRESSION.sci @@ -0,0 +1,232 @@ +// 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]=EXPRESSION(job,arg1,arg2) + x=[]; + y=[]; + typ=[]; + select job + case "set" then + x=arg1 + model=arg1.model; + graphics=arg1.graphics; + exprs=graphics.exprs + %scicos_context=%scicos_context; + for ii=1:8, + execstr("%scicos_context.u"+string(ii)+"=0"), + end + ieee(2) + while %t do + [ok,%nin,%exx,%usenz,exprs]=scicos_getvalue(.. + ["Give a scalar scilab expression using inputs u1, u2,..."; + "If only one input, input is vector [u1,u2,...] (max 8)"; + "ex: (dd*u1+sin(u2)>0)*u3"; + "Note that here dd must be defined in context"],.. + ["number of inputs";"scilab expression";"use zero-crossing (0: no, 1 yes)"],.. + list("vec",1,"vec",1,"vec",1),exprs) + ieee(0) + clear %scicos_context + + if ~ok then + break, + end + %exx=strsubst(exprs(2)," ","") + if %exx==emptystr() then + %exx="0", + end //avoid empty + //expression + + if %nin==1 then + %nini=8, + else + %nini=%nin + end + %head="%foo(" + for %jji=1:%nini-1, + %head=%head+"u"+string(%jji)+",", + end + %head=%head+"u"+string(%nini)+")" + ok=execstr("deff(%head,%exx)","errcatch")==0 + if ~ok then + message(["Erroneous expression";lasterror()]) + else + if %nin>1 then + [model,graphics,ok]=check_io(model,graphics,ones(1,%nin),1,[],... + []) + else + [model,graphics,ok]=check_io(model,graphics,-1,1,[],... + []) + end + if ok then + + [ok,%ok1,ipar,rpar,%nz]=compiler_expression(%foo) + + if ~ok then + message(["Erroneous expression";lasterror()]) + else + if %ok1 then + model.rpar=rpar + model.ipar=ipar + if %usenz then + model.nzcross=%nz + model.nmode=%nz + else + model.nzcross=0 + model.nmode=0 + end + graphics.exprs=exprs + x.graphics=graphics + x.model=model + break + end + end + end + end + end + case "define" then + in=[1;1] + out=1 + txt="(u1>0)*sin(u2)^2" + deff("%foo(u1,u2)",txt) + [%ok1,ipar,rpar,nz]=compile_expr(%foo) + model=scicos_model() + model.sim=list("evaluate_expr",4) + model.in=in + model.out=out + model.rpar=rpar + model.ipar=ipar + model.nzcross=nz + model.nmode=nz + model.dep_ut=[%t %f] + exprs=[string(size(in,"*"));txt;"1"] + gr_i=[] + x=standard_define([5 2],model,exprs,gr_i) + end +endfunction + +function [ok,%ok1,ipar,rpar,%nz]=compiler_expression(%foo) + ok=%t,%ok1=%f,ipar=[],rpar=[],%nz=[] + if exists("%scicos_context") then + %mm=getfield(1,%scicos_context) + for %mi=%mm(3:$) + if execstr(%mi+"=%scicos_context(%mi)","errcatch")<>0 then + ok=%f + end + end + end + if ok then + ok=execstr("[%ok1,ipar,rpar,%nz]=compile_expr(%foo)","errcatch")==0 + end +endfunction + + +function [ok,%ipar,%rpar,%nz]=compile_expr(%foo) + ok=%t + %lst=macr2lst(%foo); + %mm=macrovar(%foo); + %MM=%mm(3); + %FF=["sin";"cos";"tan";"exp";"log"; + "sinh";"cosh";"tanh"; + "int";"round";"ceil";"floor"; + "sign";"abs";"max";"min"; + "asin";"acos";"atan";"asinh";"acosh";"atanh"; + "atan2"; + "log10"; + ]; // ops above 100 + %num_arg=[1;1;1;1;1; + 1;1;1 + 1;1;1;1; + 1;1;2;2; + 1;1;1;1;1;1; + 2; + 1 + ]; //number of arguments + %ZCR=[16,17,18,19,20,21,28,29,30,109,110,111,112,113,114,115,116]; // ops with zero-crossing + %UU=%mm(1) + %ipar=[] + %rpar=[] + %nrpar=0 + %nz=0 + %ijk=4 + while %ijk<length(%lst) + %ijk=%ijk+1 + select evstr(%lst(%ijk)(1)) + case 2 + %jjk=find(%lst(%ijk)(2)==%FF) + if %jjk<> [] then + if evstr(%lst(%ijk)(4))<>%num_arg(%jjk) then + message(%lst(%ijk)(2)+" must have "+string(%num_arg(%jjk))+" arguments") + ok=%f + return + else + %ipar=[%ipar;5;100+%jjk] + if or(100+%jjk==%ZCR) then + %nz=%nz+1, + end + %ijk=%ijk+1 + end + else + %jjk=find(%lst(%ijk)(2)==%MM) + if %jjk<> [] then + if ~exists(%MM(%jjk)) then + message("Variable "+%MM(%jjk)+" is not defined.") + ok=%f + return + end + %var=evstr(%MM(%jjk)) + if size(%var,"*")<>1 then + message("Variable "+%MM(%jjk)+" is not scalar.") + ok=%f + return + else + %nrpar=%nrpar+1 + %rpar(%nrpar)=%var + %ipar=[%ipar;6;%nrpar] + end + else + %jjk=find(%lst(%ijk)(2)==%UU) + if %jjk<> [] then + %ipar=[%ipar;2;%jjk] + else + message("Unknown variable "+%lst(%ijk)(2)) + ok=%f + end + //%ipar=[%ipar;2;evstr(strsubst(%lst(%ijk)(2),'u',''))] + end + end + case 5 + // case of - with one operand (-u1) + if (evstr(%lst(%ijk)(2))==2)&(evstr(%lst(%ijk)(3))==1) then + %ipar=[%ipar;5;99] + else + %ipar=[%ipar;5;evstr(%lst(%ijk)(2))] + if or(evstr(%lst(%ijk)(2))==%ZCR) then + %nz=%nz+1, + end + end + case 6 + // %ipar=[%ipar;6;evstr(%lst(%ijk)(2))] + %nrpar=%nrpar+1 + %rpar(%nrpar)=evstr(%lst(%ijk)(2)) + %ipar=[%ipar;6;%nrpar] + end + end +endfunction |