4c27f978ac9230b393d60d4b9c6ce479651724f2
[scilab.git] / scilab / modules / scicos_blocks / macros / NonLinear / INTRPLBLK_f.sci
1 //  Scicos
2 //
3 //  Copyright (C) INRIA - METALAU Project <scicos@inria.fr>
4 //
5 // This program is free software; you can redistribute it and/or modify
6 // it under the terms of the GNU General Public License as published by
7 // the Free Software Foundation; either version 2 of the License, or
8 // (at your option) any later version.
9 //
10 // This program is distributed in the hope that it will be useful,
11 // but WITHOUT ANY WARRANTY; without even the implied warranty of
12 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
13 // GNU General Public License for more details.
14 //
15 // You should have received a copy of the GNU General Public License
16 // along with this program; if not, write to the Free Software
17 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
18 //
19 // See the file ../license.txt
20 //
21
22 function [x,y,typ]=INTRPLBLK_f(job,arg1,arg2)
23     x=[];y=[];typ=[];
24     select job
25     case "plot" then
26         standard_draw(arg1)
27     case "getinputs" then
28         [x,y,typ]=standard_inputs(arg1)
29     case "getoutputs" then
30         [x,y,typ]=standard_outputs(arg1)
31     case "getorigin" then
32         [x,y]=standard_origin(arg1)
33     case "set" then
34         x=arg1;
35         graphics=arg1.graphics;exprs=graphics.exprs
36         model=arg1.model;
37         while %t do
38             [ok,a,b,exprs]=scicos_getvalue("Set Interpolation block parameters",..
39             ["X coord.";"Y coord."],list("vec",-1,"vec",-1),exprs)
40             if ~ok then break,end
41             if size(a,"*") <> size(b,"*") then
42                 message("X and Y must have the same size")
43             elseif min(a(2:$)-a(1:$-1)) <=0 then
44                 message("X must be strictly increasing")
45             else
46                 if ok then
47                     graphics.exprs=exprs
48                     model.rpar=[a(:);b(:)]
49                     x.graphics=graphics;x.model=model
50                     break
51                 end
52             end
53         end
54     case "define" then
55         a=[0;1];b=[0;1]
56         model=scicos_model()
57         model.sim="intrpl"
58         model.in=1
59         model.out=1
60         model.rpar=[a;b]
61         model.blocktype="c"
62         model.dep_ut=[%t %f]
63
64         exprs=[strcat(sci2exp(a));strcat(sci2exp(b))]
65         gr_i=["xstringb(orig(1),orig(2),''interp'',sz(1),sz(2),''fill'');"]
66         x=standard_define([2 2],model,exprs,gr_i)
67     end
68 endfunction