1 files changed, 96 insertions, 0 deletions
diff --git a/macros/NonLinear/DLRADAPT_f.sci b/macros/NonLinear/DLRADAPT_f.sci
new file mode 100644
index 00000000..b26c662e
--- /dev/null
+++ b/macros/NonLinear/DLRADAPT_f.sci
@@ -0,0 +1,96 @@
+//  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]=DLRADAPT_f(job,arg1,arg2)
+    x=[];
+    y=[];
+    typ=[];
+    select job
+    case "set" then
+        x=arg1;
+        graphics=arg1.graphics;
+        exprs=graphics.exprs
+        model=arg1.model;
+        while %t do
+            [ok,p,rn,rd,g,last_u,last_y,exprs]=scicos_getvalue("Set block parameters",..
+            ["Vector of p mesh points";
+            "Numerator roots (one line for each mesh)";
+            "Denominator roots (one line for each mesh)";
+            "Vector of gain at mesh points" ;
+            "past inputs (Num degree values)";
+            "past outputs (Den degree values)"],..
+            list("vec",-1,"mat",[-1,-1],"mat",..
+            ["size(%1,''*'')","-1"],"vec","size(%1,''*'')",..
+            "vec","size(%2,2)","vec","size(%3,2)"),exprs)
+
+
+            if ~ok then
+                break,
+            end
+            m=size(rn,2)
+            [npt,n]=size(rd)
+            if m>=n then
+                message("Transfer must be strictly proper"),
+            elseif size(rn,1)<>0&size(rn,1)<>size(p,"*") then
+                message("Numerator roots matrix row size''s is incorrect")
+            else
+                rpar=[p(:);real(rn(:));imag(rn(:));real(rd(:));imag(rd(:));g(:)]
+                ipar=[m;n;npt]
+                model.dstate=[last_u(:);last_y(:)]
+                model.rpar=rpar
+                model.ipar=ipar
+                graphics.exprs=exprs
+                x.graphics=graphics;
+                x.model=model
+                break;
+            end
+        end
+    case "define" then
+        p=[0;1]
+        rn=[]
+        rd=[0.2+0.8*%i,0.2-0.8*%i;0.3+0.7*%i,0.3-0.7*%i]
+        g=[1;1]
+        last_u=[]
+        last_y=[0;0]
+
+        model=scicos_model()
+        model.sim="dlradp"
+        model.in=[1;1]
+        model.out=1
+        model.evtin=1
+
+        model.dstate=[last_u;last_y]
+        model.rpar=[p(:);real(rn(:));imag(rn(:));real(rd(:));imag(rd(:));g(:)]
+        model.ipar=[0;2;2]
+        model.blocktype="d"
+        model.firing=[]
+        model.dep_ut=[%t %f]
+
+        exprs=[sci2exp(p);
+        sci2exp(rn);
+        sci2exp(rd,0);
+        sci2exp(g);
+        sci2exp(last_u);
+        sci2exp(last_y)]
+        gr_i=[]
+        x=standard_define([2 2],model,exprs,gr_i)
+    end
+endfunction