672bc08caeca3828a8c93a9dd7e00e0f1ec70091
[scilab.git] / scilab / modules / scicos_blocks / macros / NonLinear / INTRP2BLK_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]=INTRP2BLK_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,c,exprs]=scicos_getvalue("Set Interpolation block parameters",..
39             ["X coord.";"Y coord.";"Z values"],list("vec",-1,"vec",-1,"mat",[-1,-1]),exprs)
40             if ~ok then break,end
41             if size(a,"*") <> size(c,"c") | size(b,"*") <> size(c,"r") then
42                 message("incompatible dimension")
43             elseif min(a(2:$)-a(1:$-1)) <=0 | min(b(2:$)-b(1:$-1)) <=0 then
44                 message("X and Y must be strictly increasing")
45             else
46                 if ok then
47                     graphics.exprs=exprs
48                     model.rpar=[a(:);b(:);c(:)]
49                     model.ipar=[size(a,"*");size(b,"*")]
50                     x.graphics=graphics;x.model=model
51                     break
52                 end
53             end
54         end
55     case "define" then
56         a=[0;1];b=[0;1];c=[0 1;1 2]
57         model=scicos_model()
58         model.sim=list("intrp2",1)
59         model.in=[1;1]
60         model.out=1
61         model.rpar=[a;b;c(:)]
62         model.ipar=[2;2]
63         model.blocktype="c"
64         model.dep_ut=[%t %f]
65         exprs=[strcat(sci2exp(a));strcat(sci2exp(b));strcat(sci2exp(c,0))]
66         gr_i=["xstringb(orig(1),orig(2),''interp2'',sz(1),sz(2),''fill'');"]
67         x=standard_define([3 2],model,exprs,gr_i)
68     end
69 endfunction