summaryrefslogtreecommitdiff
path: root/js/PDE/PDE.js
blob: c9339a905890a714f76a6b55361fd5e436d72a93 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
/* autogenerated from "macros/PDE/PDE.sci" */
function PDE() {
    PDE.prototype.define = function PDE() {
        this.model = scicos_model();
        this.model.state = new ScilabDouble([zeros(10,1)]);
        this.model.sim = list(new ScilabString(["PDE"]), new ScilabDouble([0]));
        this.model.in = new ScilabDouble([1],[1],[1],[1],[1]);
        this.model.out = new ScilabDouble([10],[0]);
        this.model.blocktype = new ScilabString(["c"]);
        this.model.dep_ut = new ScilabBoolean([false,true]);
        var params_pde = tlist([["paramspde"],["a"],["b"],["txt_exp"],["check_op1"],["a1"],["b1"],["check_op2"],["a2"],["b2"],["check_op3"],["a3"],["b3"],["check_op4"],["a4"],["b4"],["check_op5"],["a5"],["b5"],["check_op6"],["a6"],["b6"],["check_op7"],["a7"],["b7"],["discr_cst"],["discr_non_cst"],["signe"],["rad_automatique"],["rad_manuel"],["methode"],["ord1"],["ord2"],["ord3"],["degre"],["nnode"],["txt_pas"],["CI"],["dCI"],["CLa"],["CLa_exp"],["CLb"],["CLb_exp"],["points"]],"","","","0","","IN_EDP1(t)","0","","IN_EDP2(t)","0","","IN_EDP3(t)","0","","IN_EDP4(t)","0","","IN_EDP5(t)","0","","IN_EDP6(t)","0","","IN_EDP7(t)","0","0","0","0","0","0","","","","","","","","","0","IN_CL1(t)","0","IN_CL2(t)","");
        var label = list(params_pde,[],"");
        this.gr_i = new ScilabString(["xstringb(orig(1),orig(2),\"PDE\",sz(1),sz(2));"]);
        this.x = new standard_define(new ScilabDouble([3,3]),this.model,label,this.gr_i);
        return new BasicBlock(this.x);
    }
    PDE.prototype.details = function PDE() {
        return this.x;
    }
    PDE.prototype.get = function PDE() {
        var options = {
                okk:["New block\'s name :",this.okk],
        }
        return options;
    }
    PDE.prototype.set = function PDE() {
        var label = this.graphics.exprs;
        var params_pde = label[1-1];
        while (true) {
            [ln,fun]=where()
            if (!or(fun=="do_eval")) {
                var tmpvar0 = IHM_EDP(params_pde);
                var ok = tmpvar0[0];
                var a_domaine = tmpvar0[1];
                var b_domaine = tmpvar0[2];
                var discr = tmpvar0[3];
                var signe = tmpvar0[4];
                var choix = tmpvar0[5];
                var type_meth = tmpvar0[6];
                var degre = tmpvar0[7];
                var Nbr_maillage = tmpvar0[8];
                var CI = tmpvar0[9];
                var CI1 = tmpvar0[10];
                var CLa_type = tmpvar0[11];
                var CLa_exp = tmpvar0[12];
                var CLb_type = tmpvar0[13];
                var CLb_exp = tmpvar0[14];
                var oper = tmpvar0[15];
                var a1 = tmpvar0[16];
                var b1 = tmpvar0[17];
                var a2 = tmpvar0[18];
                var b2 = tmpvar0[19];
                var a3 = tmpvar0[20];
                var b3 = tmpvar0[21];
                var a4 = tmpvar0[22];
                var b4 = tmpvar0[23];
                var a5 = tmpvar0[24];
                var b5 = tmpvar0[25];
                var a6 = tmpvar0[26];
                var b6 = tmpvar0[27];
                var a7 = tmpvar0[28];
                var b7 = tmpvar0[29];
                var k = tmpvar0[30];
                var mesures = tmpvar0[31];
                var params_pde = tmpvar0[32];
                if (ok) {
                    return;
                }
            } else {
                if (exists("%scicos_context")) {
                    var tmpvar1 = setvalue_IHM_EDP(params_pde);
                    var ok = tmpvar1[0];
                    var a_domaine = tmpvar1[1];
                    var b_domaine = tmpvar1[2];
                    var discr = tmpvar1[3];
                    var signe = tmpvar1[4];
                    var choix = tmpvar1[5];
                    var type_meth = tmpvar1[6];
                    var degre = tmpvar1[7];
                    var Nbr_maillage = tmpvar1[8];
                    var CI = tmpvar1[9];
                    var CI1 = tmpvar1[10];
                    var CLa_type = tmpvar1[11];
                    var CLa_exp = tmpvar1[12];
                    var CLb_type = tmpvar1[13];
                    var CLb_exp = tmpvar1[14];
                    var oper = tmpvar1[15];
                    var a1 = tmpvar1[16];
                    var b1 = tmpvar1[17];
                    var a2 = tmpvar1[18];
                    var b2 = tmpvar1[19];
                    var a3 = tmpvar1[20];
                    var b3 = tmpvar1[21];
                    var a4 = tmpvar1[22];
                    var b4 = tmpvar1[23];
                    var a5 = tmpvar1[24];
                    var b5 = tmpvar1[25];
                    var a6 = tmpvar1[26];
                    var b6 = tmpvar1[27];
                    var a7 = tmpvar1[28];
                    var b7 = tmpvar1[29];
                    var k = tmpvar1[30];
                    var mesures = tmpvar1[31];
                    var params_pde = tmpvar1[32];
                }
            }
            this.okk = false;
            this.rdnom = "PDE";
            var ok1 = true;
            while (true) {
                var ok = true;
                this.okk = parseBoolean(arguments[0]["okk"]);
                this.rdnom = arguments[0]["rdnom"];
                this.lab = arguments[0]["lab"];
                if (this.okk==false) {
                    var ok1 = false;
                    return;
                }
                label[3-1] = this.lab;
                this.rdnom = stripblanks(this.rdnom);
                if (this.rdnom==emptystr()) {
                    var ok1 = false;
                    x_message("sorry C file name not defined");
                }
                if (ok1) {
                    break;
                }
            }
            if ((choix==0)) {
                var ind4 = strindex(a4,"x");
                var ind1 = strindex(a1,"x");
                var ind2 = strindex(a2,"x");
                if ((ind4.length!=0||ind1.length!=0||ind2.length!=0)) {
                    if ((signe==1)) {
                        var delta = 1;
                    } else if ((signe==2)) {
                        var delta = -1;
                    } else if ((signe==0)) {
                        var delta = 0;
                    } else {
                        x_message([["le discriminant n\'est pas constant,"],["Vous devez choisir son signe dans l\'IHM"]]);
                        return;
                    }
                } else {
                    var delta = evstr(a4)^2-4*evstr(a1)*evstr(a2);
                }
                if ((delta.length==0)) {
                    var delta = 0;
                }
                var type_meth = arbre_decision(delta);
            }
            var tmpvar2 = translate(CI,CI1,CLa_type,CLa_exp,CLb_type,CLb_exp,oper,type_meth,degre,a_domaine,b_domaine,Nbr_maillage,a1,b1,a2,b2,a3,b3,a4,b4,a5,b5,a6,b6,a7,b7,this.rdnom,mesures);
            var flag_type = tmpvar2[0];
            this.rdnom = tmpvar2[1];
            var DF_type = tmpvar2[2];
            var tt = tmpvar2[3];
            var Nbr = Nbr_maillage;
            if (((CLa_type==1)&&(DF_type==0||DF_type==1))||((CLb_type==1)&&(DF_type==0||DF_type==2))) {
                var Nbr = Nbr+1;
            }
            if ((mesures.length==0)) {
                var out = Nbr_maillage;
            } else {
                var out = [[Nbr_maillage],[size(mesures,"*")]];
            }
            if ((flag_type==1)) {
                this.model.sim = list(new ScilabDouble([this.rdnom]), new ScilabDouble([2004]));
                if ((find(oper==1).length!=0)) {
                    this.model.state = new ScilabDouble([zeros(2*Nbr_maillage,1)]);
                } else {
                    this.model.state = new ScilabDouble([zeros(Nbr_maillage,1)]);
                }
            } else if ((flag_type==2)) {
                this.model.sim = list(new ScilabDouble([this.rdnom]), new ScilabDouble([12004]));
                if ((find(oper==1).length!=0)) {
                    if ((type_meth==3&&(find(oper==2).length!=0||find(oper==4).length!=0))) {
                        this.model.state = new ScilabDouble([zeros(6*Nbr_maillage,1)]);
                    } else if ((type_meth==1)) {
                        this.model.state = new ScilabDouble([zeros(4*Nbr,1)]);
                    } else {
                        this.model.state = new ScilabDouble([zeros(4*Nbr_maillage,1)]);
                    }
                } else {
                    if ((type_meth==3&&(find(oper==2).length!=0||find(oper==4).length!=0))) {
                        this.model.state = new ScilabDouble([zeros(4*Nbr_maillage,1)]);
                    } else if ((type_meth==1)) {
                        this.model.state = new ScilabDouble([zeros(2*Nbr,1)]);
                    } else {
                        this.model.state = new ScilabDouble([zeros(2*Nbr_maillage,1)]);
                    }
                }
            }
            var tmpvar3 = CFORTREDP(this.rdnom,tt);
            var ok1 = tmpvar3[0];
            if (!ok1) {
                break;
            }
            if (!ok) {
                var tmpvar4 = check_io(this.model,this.graphics,ones(k,1),out.slice(),[],[]);
                this.model = tmpvar4[0];
                this.graphics = tmpvar4[1];
                var ok = tmpvar4[2];
            }
            label[1-1] = params_pde;
            label[2-1] = tt;
            this.graphics.exprs = new ScilabDouble([label]);
            this.x.graphics = this.graphics;
            this.x.model = this.model;
            break;
        }
        return new BasicBlock(this.x);
    }
}