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
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2002-2004 - INRIA - Vincent COUVERT
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
function [tree]=sci_pow2(tree)
// M2SCI function
// Conversion function for Matlab pow2()
// Input: tree = Matlab funcall tree
// Ouput: tree = Scilab equivalent for tree
// X = pow2(Y)
if rhs==1 then
// %c_pow2 and %b_pow2 are not defined in Scilab
Y = getrhs(tree)
Y = convert2double(Y)
tree=Operation(".^",list(Cste(2),Y),list(Variable(tree.lhs(1).name,Y.infer)))
// X = pow2(E,F)
else
// %c_pow2 and %b_pow2 are not defined in Scilab
[E,F] = getrhs(tree)
E = convert2double(E)
F = convert2double(F)
// Imaginary parts are ignored by Matlab
if is_complex(E) then
E = Funcall("real",1,list(E),list(Variable("",E.infer)))
elseif ~is_real(E) then
newE = Funcall("real",1,list(E),list(Variable("",E.infer)))
repl_poss(newE,E,E,gettext("is Real."));
E=newE
end
if is_complex(F) then
F = Funcall("real",1,list(F),list(Variable("",F.infer)))
elseif ~is_real(F) then
newF = Funcall("real",1,list(F),list(Variable("",F.infer)))
repl_poss(newF,F,F,gettext("is Real."));
F=newF
end
tmp=Operation(".^",list(Cste(2),F),list(Variable("",F.infer)))
tree=Operation(".*",list(E,tmp),list(Variable(tree.lhs(1).name,Infer())))
if is_a_scalar(E) then
tree.out(1).dims=F.dims
elseif is_a_scalar(F) then
tree.out(1).dims=E.dims
elseif E.dims==F.dims then
tree.out(1).dims=E.dims
else
tree.out(1).dims=allunknown(E.dims)
end
tree.out(1).type=Type(Double,Real)
end
endfunction
|
