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
|
function [tree]=sci_MFUN(tree)
// Prototype for max and min
// Copyright INRIA
// M2SCI function
// Conversion function for Matlab MFUN()
// Input: tree = Matlab funcall tree
// Ouput: tree = Scilab equivalent for tree
// Emulation function: mtlb_MFUN()
// V.C.
// C = MFUN(A) or [C,I] = MFUN(A)
if rhs==1 then
A = getrhs(tree)
vtype=Double
if or(A.vtype==[Boolean,Unknown]) then
vtype=Unknown // If A is a scalar then Matlab return Double type value else a Boolean type value
end
A = convert2double(A)
tree.rhs = Rhs(A)
dim = first_non_singleton(A)
if dim==-1 then
// All output dims are unknown
tree.lhs(1).dims=allunknown(A.dims)
tmp=gettempvar()
insert(Equal(list(tmp),A))
// First non singleton dimension will be computed at execution
tree.rhs=Rhs(tmp,Funcall("firstnonsingleton",1,list(tmp),list()))
else
tree.lhs(1).dims=A.dims
if dim==0 then
tree.rhs=Rhs(A)
tree.lhs(1).dims=list(1,1)
elseif dim==1 then
tree.rhs=Rhs(A,"r")
tree.lhs(1).dims(dim)=1
elseif dim==2 then
tree.rhs=Rhs(A,"c")
tree.lhs(1).dims(dim)=1
else
tree.rhs=Rhs(A,dim)
tree.lhs(1).dims(dim)=1
end
end
// C = MFUN(A) or [C,I] = MFUN(A)
if is_real(A) then
tree.lhs(1).type=Type(vtype,Real)
elseif is_complex(A) then
// Scilab MFUN() does not work with complexes so mtlb_MFUN() is called
tree.name="mtlb_MFUN"
tree.lhs(1).type=Type(Double,Unknown)
else
// Scilab MFUN() does not work with complexes so mtlb_MFUN() is called but can be replaced
tree.name="mtlb_MFUN"
tree.lhs(1).type=Type(vtype,Unknown)
end
// [C,I] = MFUN(A)
if lhs==2 then
// All dims Unknown because if A is [] then I is also [] else I is a row vector
tree.lhs(2).dims=allunknown(A.dims)
tree.lhs(2).type=Type(Double,Real)
end
// C = MFUN(A,B)
elseif rhs==2 then
[A,B]=getrhs(tree)
vtype=Double
if A.vtype==Boolean & B.vtype==Boolean | A.vtype==Unknown & B.vtype==Unknown then
vtype=Unknown
end
A=convert2double(A)
B=convert2double(B)
tree.rhs=Rhs(A,B)
if is_real(A) & is_real(B) then
if not_empty(A) & not_empty(B) then
tree.lhs(1).dims=A.dims
tree.lhs(1).type=Type(vtype,Real)
else
// Perhaps an input is an empty matrix
tree.name="mtlb_MFUN"
tree.lhs(1).dims=A.dims
tree.lhs(1).type=Type(vtype,Real)
end
else
// Inputs can be complexes and/or empty matrices
tree.name="mtlb_MFUN"
tree.lhs(1).dims=A.dims
tree.lhs(1).type=Type(vtype,Unknown)
end
// C = MFUN(A,[],dim) or [C,I] = MFUN(A,[],dim)
else
[A,tmp,dim]=getrhs(tree)
vtype=Double
if or(A.vtype==[Boolean,Unknown]) then
vtype=Unknown
end
A=convert2double(A)
tree.rhs=Rhs(A)
// C = MFUN(A,[],dim) or [C,I] = MFUN(A,[],dim)
if or(lhs==[1,2]) then
if is_real(A) then
tree.lhs(1).type=Type(Double,Real)
if typeof(dim)=="cste" then
if dim.value==1 then
tree.rhs=Rhs(A,"r")
tree.lhs(1).dims=A.dims
tree.lhs(1).dims(1)=Unknown // 0 or 1
elseif dim.value==2 then
tree.rhs=Rhs(A,"c")
tree.lhs(1).dims=A.dims
tree.lhs(1).dims(2)=Unknown // 0 or 1
elseif dim.value<=size(A.dims) then
tree.lhs(1).dims=A.dims
tree.lhs(1).dims(dim.value)=1
else
// Scilab MFUN() does not work when dim is greater than number of dims of A
tree.name="mtlb_MFUN"
tree.rhs=Rhs(A,tmp,dim)
tree.lhs(1).dims=A.dims
end
else
// If dim is 1 it can be replaced by 'r'
// If dim is 2 it can be replaced by 'c'
tree.name="mtlb_MFUN"
tree.rhs=Rhs(A,tmp,dim)
tree.lhs(1).dims=allunknown(A.dims)
end
else
// A can be complex....
tree.name="mtlb_MFUN"
tree.rhs=Rhs(A,tmp,dim)
tree.lhs(1).dims=allunknown(A.dims)
tree.lhs(1).type=Type(Double,Unknown)
end
end
// [C,I] = MFUN(A,[],dim)
if lhs==2 then
tree.lhs(2).type=Type(Double,Real)
if is_real(A) then
if typeof(dim)=="cste" then
if dim.value<=size(A.dims) then
tree.lhs(2).dims=allunknown(A.dims)
else
tree.lhs(2).dims=A.dims
end
else
tree.lhs(2).dims=allunknown(A.dims)
end
else
tree.lhs(2).dims=allunknown(A.dims)
end
end
end
endfunction
|
