Sigbuilder_intfunc

Fonction d'interfaçage
fr

##
Sigbuilder

### Contenu du fichier

    function [x,y,typ]=Sigbuilder(job,arg1,arg2)
    x=[];y=[],typ=[]
    select job
    case 'plot' then
      standard_draw(arg1)
    case 'getinputs' then
      [x,y,typ]=standard_inputs(arg1)
    case 'getoutputs' then
      [x,y,typ]=standard_outputs(arg1)
    case 'getorigin' then
      [x,y]=standard_origin(arg1)
    case 'set' then
      //paths to updatable parameters or states
      ppath = list(1,3)
      newpar=list();
      y=0;
      for path=ppath do
        np=size(path,'*')
        spath=list()
        for k=1:np
          spath($+1)='model'
          spath($+1)='rpar'
          spath($+1)='objs'
          spath($+1)=path(k)
        end
        xx=arg1(spath)// get the block
        execstr('xxn='+xx.gui+'(''set'',xx)')
        if diffobjs(xxn,xx) then
          model=xx.model
          model_n=xxn.model
          if ~is_modelica_block(xx) then
            modified=or(model.sim<>model_n.sim)|..
                     ~isequal(model.state,model_n.state)|..
                     ~isequal(model.dstate,model_n.dstate)|..
                     ~isequal(model.odstate,model_n.odstate)|..
                     ~isequal(model.rpar,model_n.rpar)|..
                     ~isequal(model.ipar,model_n.ipar)|..
                     ~isequal(model.opar,model_n.opar)|..
                     ~isequal(model.label,model_n.label)
            if or(model.in<>model_n.in)|or(model.out<>model_n.out)|..
               or(model.in2<>model_n.in2)|or(model.out2<>model_n.out2)|..
               or(model.outtyp<>model_n.outtyp)|or(model.intyp<>model_n.intyp) then
              needcompile=1
            end
            if or(model.firing<>model_n.firing) then
              needcompile=2
            end
            if (size(model.in,'*')<>size(model_n.in,'*'))|..
              (size(model.out,'*')<>size(model_n.out,'*')) then
              needcompile=4
            end
            if model.sim=='input'|model.sim=='output' then
              if model.ipar<>model_n.ipar then
                needcompile=4
              end
            end
            if or(model.blocktype<>model_n.blocktype)|..
               or(model.dep_ut<>model_n.dep_ut) then
              needcompile=4
            end
            if (model.nzcross<>model_n.nzcross)|(model.nmode<>model_n.nmode) then
              needcompile=4
            end
            if prod(size(model_n.sim))>1 then
              if model_n.sim(2)>1000 then
                if model.sim(1)<>model_n.sim(1) then
                  needcompile=4
                end
              end
            end
          else
            modified=or(model_n<>model)
            eq=model.equations;eqn=model_n.equations;
            if or(eq.model<>eqn.model)|or(eq.inputs<>eqn.inputs)|..
               or(eq.outputs<>eqn.outputs) then
              needcompile=4
            end
          end
         //parameter or states changed
          arg1(spath)=xxn// Update
          newpar(size(newpar)+1)=path// Notify modification
          y=max(y,needcompile)
        end
      end
      x=arg1
      typ=newpar
    case 'define' then
    scs_m_1=scicos_diagram(..
            version="scicos4.2",..
            props=scicos_params(..
                  wpar=[600,450,0,0,600,450],..
                  Title=["Sigbuilder","./"],..
                  tol=[0.0001;0.000001;1.000D-10;100001;0;0;0],..
                  tf=100,..
                  context=" ",..
                  void1=[],..
                  options=tlist(["scsopt","3D","Background","Link","ID","Cmap"],list(%t,33),[8,1],[1,5],..
                  list([5,1],[4,1]),[0.8,0.8,0.8]),..
                  void2=[],..
                  void3=[],..
                  doc=list()))
    scs_m_1.objs(1)=scicos_block(..
                    gui="CURVE_c",..
                    graphics=scicos_graphics(..
                             orig=[329.63473,606.18517],..
                             sz=[40,40],..
                             flip=%t,..
                             theta=0,..
                             exprs=["3";"[0,1,2]";"[10,20,-30]";"y";"n"],..
                             pin=[],..
                             pout=6,..
                             pein=4,..
                             peout=2,..
                             gr_i=list(..
                             ["rpar=arg1.model.rpar;n=model.ipar(1);order=model.ipar(2);";
                             "xx=rpar(1:n);yy=rpar(n+1:2*n);";
                             "[XX,YY,rpardummy]=Do_Spline(n,order,xx,yy)";
                             "xmx=maxi(XX);xmn=mini(XX);";
                             "ymx=maxi(YY);ymn=mini(YY);";
                             "dx=xmx-xmn;if dx==0 then dx=maxi(xmx/2,1);end";
                             "xmn=xmn-dx/20;xmx=xmx+dx/20;";
                             "dy=ymx-ymn;if dy==0 then dy=maxi(ymx/2,1);end;";
                             "ymn=ymn-dy/20;ymx=ymx+dy/20;";
                             "xx2=orig(1)+sz(1)*((XX-xmn)/(xmx-xmn));";
                             "yy2=orig(2)+sz(2)*((YY-ymn)/(ymx-ymn));";
                             "xset(''color'',2)";
                             "xpoly(xx2,yy2,''lines'');"],8),..
                             id="",..
                             in_implicit=[],..
                             out_implicit="E"),..
                    model=scicos_model(..
                             sim=list("curve_c",4),..
                             in=[],..
                             in2=[],..
                             intyp=1,..
                             out=1,..
                             out2=[],..
                             outtyp=1,..
                             evtin=1,..
                             evtout=1,..
                             state=[],..
                             dstate=[],..
                             odstate=list(),..
                             rpar=[0;1;2;10;20;-30],..
                             ipar=[3;3;1],..
                             opar=list(),..
                             blocktype="c",..
                             firing=0,..
                             dep_ut=[%f,%t],..
                             label="",..
                             nzcross=0,..
                             nmode=0,..
                             equations=list()),..
                    doc=list())
    scs_m_1.objs(2)=scicos_link(..
                      xx=[349.63473;349.49528],..
                      yy=[600.47089;565.10704],..
                      id="drawlink",..
                      thick=[0,0],..
                      ct=[5,-1],..
                      from=[1,1,0],..
                      to=[3,1,1])
    scs_m_1.objs(3)=scicos_block(..
                    gui="CLKSPLIT_f",..
                    graphics=scicos_graphics(..
                             orig=[349.49528;565.10704],..
                             sz=[0.3333333,0.3333333],..
                             flip=%t,..
                             theta=0,..
                             exprs=[],..
                             pin=[],..
                             pout=[],..
                             pein=2,..
                             peout=[8;4],..
                             gr_i=list([],8),..
                             id="",..
                             in_implicit=[],..
                             out_implicit=[]),..
                    model=scicos_model(..
                             sim="split",..
                             in=[],..
                             in2=[],..
                             intyp=1,..
                             out=[],..
                             out2=[],..
                             outtyp=1,..
                             evtin=1,..
                             evtout=[1;1],..
                             state=[],..
                             dstate=[],..
                             odstate=list(),..
                             rpar=[],..
                             ipar=[],..
                             opar=list(),..
                             blocktype="d",..
                             firing=[%f,%f,%f],..
                             dep_ut=[%f,%f],..
                             label="",..
                             nzcross=0,..
                             nmode=0,..
                             equations=list()),..
                    doc=list())
    scs_m_1.objs(4)=scicos_link(..
                      xx=[349.49528;266.69602;266.69602;270.35525;342.80795;342.80795;349.63473],..
                      yy=[565.10704;565.10704;680.99483;680.99483;680.99483;651.89946;651.89946],..
                      id="drawlink",..
                      thick=[0,0],..
                      ct=[5,-1],..
                      from=[3,2,0],..
                      to=[1,1,1])
    scs_m_1.objs(5)=scicos_block(..
                    gui="OUT_f",..
                    graphics=scicos_graphics(..
                             orig=[398.20616,616.18517],..
                             sz=[20,20],..
                             flip=%t,..
                             theta=0,..
                             exprs="1",..
                             pin=6,..
                             pout=[],..
                             pein=[],..
                             peout=[],..
                             gr_i=list(" ",8),..
                             id="",..
                             in_implicit="E",..
                             out_implicit=[]),..
                    model=scicos_model(..
                             sim="output",..
                             in=-1,..
                             in2=-2,..
                             intyp=-1,..
                             out=[],..
                             out2=[],..
                             outtyp=1,..
                             evtin=[],..
                             evtout=[],..
                             state=[],..
                             dstate=[],..
                             odstate=list(),..
                             rpar=[],..
                             ipar=1,..
                             opar=list(),..
                             blocktype="c",..
                             firing=[],..
                             dep_ut=[%f,%f],..
                             label="",..
                             nzcross=0,..
                             nmode=0,..
                             equations=list()),..
                    doc=list())
    scs_m_1.objs(6)=scicos_link(..
                      xx=[378.20616;398.20616],..
                      yy=[626.18517;626.18517],..
                      id="drawlink",..
                      thick=[0,0],..
                      ct=[1,1],..
                      from=[1,1,0],..
                      to=[5,1,1])
    scs_m_1.objs(7)=scicos_block(..
                    gui="CLKOUTV_f",..
                    graphics=scicos_graphics(..
                             orig=[339.49528,505.10704],..
                             sz=[20,30],..
                             flip=%t,..
                             theta=0,..
                             exprs="1",..
                             pin=[],..
                             pout=[],..
                             pein=8,..
                             peout=[],..
                             gr_i=list(" ",8),..
                             id="",..
                             in_implicit=[],..
                             out_implicit=[]),..
                    model=scicos_model(..
                             sim="output",..
                             in=[],..
                             in2=[],..
                             intyp=1,..
                             out=[],..
                             out2=[],..
                             outtyp=1,..
                             evtin=1,..
                             evtout=[],..
                             state=[],..
                             dstate=[],..
                             odstate=list(),..
                             rpar=[],..
                             ipar=1,..
                             opar=list(),..
                             blocktype="d",..
                             firing=[],..
                             dep_ut=[%f,%f],..
                             label="",..
                             nzcross=0,..
                             nmode=0,..
                             equations=list()),..
                    doc=list())
    scs_m_1.objs(8)=scicos_link(..
                      xx=[349.49528;349.49528],..
                      yy=[565.10704;535.10704],..
                      id="drawlink",..
                      thick=[0,0],..
                      ct=[5,-1],..
                      from=[3,1,0],..
                      to=[7,1,1])
    model=scicos_model(..
             sim="csuper",..
             in=[],..
             in2=[],..
             intyp=1,..
             out=-1,..
             out2=[],..
             outtyp=1,..
             evtin=[],..
             evtout=1,..
             state=[],..
             dstate=[],..
             odstate=list(),..
             rpar=scs_m_1,..
             ipar=[],..
             opar=list(),..
             blocktype="h",..
             firing=[],..
             dep_ut=[%f,%f],..
             label="",..
             nzcross=0,..
             nmode=0,..
             equations=list())
    //  gr_i='xstringb(orig(1),orig(2),''Sigbuilder'',sz(1),sz(2),''fill'')';
        gr_i=['ipar=arg1.model.rpar.objs(1).model.ipar';
    	'rpar=arg1.model.rpar.objs(1).model.rpar';
    	'n=ipar(1);order=ipar(2);';
    	'xx=rpar(1:n);yy=rpar(n+1:2*n);';
    	'[XX,YY,rpardummy]=Do_Spline2(n,order,xx,yy)';
    	'xmx=maxi(XX);xmn=mini(XX);';
    	'ymx=maxi(YY);ymn=mini(YY);';
    	'dx=xmx-xmn;if dx==0 then dx=maxi(xmx/2,1);end';
    	'xmn=xmn-dx/20;xmx=xmx+dx/20;';
    	'dy=ymx-ymn;if dy==0 then dy=maxi(ymx/2,1);end;';
    	'ymn=ymn-dy/20;ymx=ymx+dy/20;';
    	'rect=[xmn,ymn;xmx,ymx];';
    	'xx2=orig(1)+sz(1)*((XX-xmn)/(xmx-xmn));';
    	'yy2=orig(2)+sz(2)*((YY-ymn)/(ymx-ymn));';
    	'xset(''color'',2)';
    	'xpoly(xx2,yy2,''lines'');']

      x=standard_define([2 2],model,[],gr_i)
    end
    endfunction

    //=========================================================
    function [X,Y,orpar]=Do_Spline2(N,order,x,y)

      X=[];Y=[];orpar=[];
      METHOD=getmethod(order);

      if (METHOD=='zero order') then
        X=x(1);Y=y(1);
        for i=1:N-1
          X=[X;x(i);x(i+1);x(i+1)];
          Y=[Y;y(i);y(i);y(i+1)];
        end
        return;
      end
      //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      if (METHOD=='linear') then
        X=[];
        for i=1:N
          X=[X;x(i)];
          Y=[Y;y(i)];
        end
        return;
      end
      //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
       if (N<20) then NP=4;else
       if (N<40) then NP=2;else
        if (N<100) then NP=1;else
            NP=0;end;end;
      end
      for i=1:N-1
        X=[X;linspace(x(i),x(i+1),NP+2)']; // pour tous sauf "linear" et "zero order"
      end
      //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      if (N>2) & (METHOD=='order 2') then
        Z=ORDER2(x,y);
        A=Z(1:N-1);
        B=Z(N:2*N-2);
        C=Z(2*N-1:3*N-3);

        for j=1:size(X,'*')
          for i=N-1:-1:1
    	if X(j)>=x(i) then,break;end
          end
          Y(j)=A(i)*(X(j)-x(i))^2+B(i)*(X(j)-x(i))+C(i);
        end
        orpar=matrix(Z,-1,1)
      end
       //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      if (METHOD=='not_a_knot') then
        try
        d = splin(x, y, METHOD);
        Y = interp(X, x, y, d);
        orpar=d(:);
        catch
         xinfo('ERROR in SPLINE: '+METHOD)
        end

      end
       //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      if (METHOD=='periodic') then
        if y(1)<>y(N) then
          y(N)=y(1)
        end
        try
          d = splin(x, y,METHOD);
          Y = interp(X, x, y, d);
          orpar=d(:);
        catch
        xinfo('ERROR in SPLINE: '+METHOD)
        end
      end
        //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      if (METHOD=='monotone' ) then
        try
          d = splin(x, y, METHOD);
          Y = interp(X, x, y, d);
          orpar=d(:);
        catch
        xinfo('ERROR in SPLINE: '+METHOD)
        end

      end
       //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      if (METHOD=='fast') then
        try
          d = splin(x, y, METHOD);
          Y = interp(X, x, y, d);
          orpar=d(:);
        catch
          xinfo('ERROR in SPLINE:  '+METHOD)
        end
      end
       //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      if (METHOD=='clamped') then
        try
          d = splin(x, y, METHOD,[0;0]);
          Y = interp(X, x, y, d);
          orpar=d(:);
        catch
          xinfo('ERROR in SPLINE: '+METHOD)
        end
      end

    endfunction

    function METHOD=getmethod(order)
      select order
       case 0 then, METHOD='zero order'
       case 1 then, METHOD='linear'
       case 2 then, METHOD='order 2'
       case 3 then, METHOD='not_a_knot'
       case 4 then, METHOD='periodic'
       case 5 then, METHOD='monotone'
       case 6 then, METHOD='fast'
       case 7 then, METHOD='clamped'
      end
    endfunction

    function [Z]=ORDER2(x,y)
    N=size(x,'*')-1;
    A=zeros(3*N-1,N*3);
    B=zeros(3*N-1,1);
    for i=1:N
       j=3*(i-1)+1;
       A(j,i+2*N)=1;
       B(j)=y(i);
       A(j+1,i)=(x(i+1)-x(i))^2;
       A(j+1,i+N)=x(i+1)-x(i);
       A(j+1,i+2*N)=1;
       B(j+1)=y(i+1);
    end

    for i=1:N-1
       j=3*(i-1)+1;
       A(j+2,i)=2*(x(i+1)-x(i));
       A(j+2,i+N)=1;
       A(j+2,i+N+1)=-1;
    end

    Q=zeros(3*N,3*N);
    for i=1:N
      Q(i,i)=4*(x(i+1)-x(i))^2
      Q(i,i+N)=2*(x(i+1)-x(i))
      Q(i+N,i)=2*(x(i+1)-x(i))
      Q(i+N,i+N)=1;
    end

    At=[Q,A';A,zeros(3*N-1,3*N-1)]
    Bt=[zeros(3*N,1);B]
    Zt=At\Bt;
    Z=Zt(1:3*N,1)
    endfunction
    //===================================================