diff options
Diffstat (limited to 'macros/Linear/SUMMATION.sci')
| -rw-r--r-- | macros/Linear/SUMMATION.sci | 159 |
1 files changed, 159 insertions, 0 deletions
diff --git a/macros/Linear/SUMMATION.sci b/macros/Linear/SUMMATION.sci new file mode 100644 index 00000000..f70c7882 --- /dev/null +++ b/macros/Linear/SUMMATION.sci @@ -0,0 +1,159 @@ +// 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]=SUMMATION(job,arg1,arg2) + x=[]; + y=[]; + typ=[]; + select job + case "set" then + x=arg1; + graphics=arg1.graphics + model=arg1.model + exprs=graphics.exprs + if size(exprs,1)==1 then + exprs=[sci2exp(1);exprs;sci2exp(0)]; + elseif size(exprs,1)==2 then + exprs=[exprs;sci2exp(0)]; + end + while %t do + [ok,Datatype,sgn,satur,exprs]=scicos_getvalue("Set sum block parameters",.. + ["Datatype (1=real double 2=complex 3=int32 ...)";.. + "Number of inputs or sign vector (of +1, -1)";.. + "Do on Overflow(0=Nothing 1=Saturate 2=Error)"],.. + list("vec",1,"vec",-1,"vec",1),exprs) + if ~ok then + break, + end + sgn=sgn(:); + if (satur~=0&satur~=1&satur~=2) then + message("Do on overflow must be 0,1,2"); + ok=%f; + end + if size(sgn,1)==1 then + if sgn<1 then + message("Number of inputs must be > 0") + ok=%f + elseif sgn==1 then + in=-1;in2=-2 + sgn=[] + nout=1;nout2=1 + else + in=-ones(sgn,1); + in2=2*in + sgn=ones(sgn,1) + nout=-1;nout2=-2 + end + else + if ~and(abs(sgn)==1) then + message("Signs can only be +1 or -1") + ok=%f + else + in=-ones(size(sgn,1),1); + in2=2*in + nout=-1;nout2=-2 + end + end + it=Datatype*ones(1,size(in,1)); + ot=Datatype; + if Datatype==1 then + model.sim=list("summation",4) + elseif Datatype==2 then + model.sim=list("summation_z",4) + elseif ((Datatype<1) |(Datatype>8)) then + message("Datatype is not supported");ok=%f; + else + if satur==0 then + if Datatype==3 then + model.sim=list("summation_i32n",4) + elseif Datatype==4 then + model.sim=list("summation_i16n",4) + elseif Datatype==5 then + model.sim=list("summation_i8n",4) + elseif Datatype==6 then + model.sim=list("summation_ui32n",4) + elseif Datatype==7 then + model.sim=list("summation_ui16n",4) + elseif Datatype==8 then + model.sim=list("summation_ui8n",4) + end + elseif satur==1 then + if Datatype==3 then + model.sim=list("summation_i32s",4) + elseif Datatype==4 then + model.sim=list("summation_i16s",4) + elseif Datatype==5 then + model.sim=list("summation_i8s",4) + elseif Datatype==6 then + model.sim=list("summation_ui32s",4) + elseif Datatype==7 then + model.sim=list("summation_ui16s",4) + elseif Datatype==8 then + model.sim=list("summation_ui8s",4) + end + elseif satur==2 then + if Datatype==3 then + model.sim=list("summation_i32e",4) + elseif Datatype==4 then + model.sim=list("summation_i16e",4) + elseif Datatype==5 then + model.sim=list("summation_i8e",4) + elseif Datatype==6 then + model.sim=list("summation_ui32e",4) + elseif Datatype==7 then + model.sim=list("summation_ui16e",4) + elseif Datatype==8 then + model.sim=list("summation_ui8e",4) + end + end + end + if ok then + [model,graphics,ok]=set_io(model,graphics,... + list([in,in2],it),... + list([nout,nout2],ot),[],[]) + end + if ok then + model.rpar=satur; + model.ipar=sgn + graphics.exprs=exprs + x.graphics=graphics; + x.model=model + break + end + end + case "define" then + sgn=[1;-1] + model=scicos_model() + model.sim=list("summation",4) + model.in=[-1;-1] + model.out=-1 + model.in2=[-2;-2] + model.out2=-2 + model.ipar=sgn + model.blocktype="c" + model.dep_ut=[%t %f] + + + exprs=sci2exp(sgn) + gr_i=[] + x=standard_define([2 3],model, exprs,gr_i) + end +endfunction |