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
|
function [tree]=sci_repmat(tree)
// Copyright INRIA
// M2SCI function
// Conversion function for Matlab repmat()
// Input: tree = Matlab funcall tree
// Ouput: tree = Scilab equivalent for tree
// Emulation function: mtlb_repmat()
// V.C.
// repmat(A,m)
if rhs==2 then
[A,m] = getrhs(tree)
if is_complex(m) then
m=Funcall("real",1,list(m),list(Variable("",m.infer)))
elseif ~is_real(m) then
newm=Funcall("real",1,list(m),list(Variable("",m.infer)))
repl_poss(newm,m,m,"is Real");
m=newm
end
vtype=A.vtype
if A.vtype==Unknown then
tree.name="mtlb_repmat"
if typeof(m)=="cste" then
for kd=1:lstsize(A.dims)
tree.lhs(1).dims(kd)=A.dims(kd)*m.value
if tree.lhs(1).dims(kd)<0 then
tree.lhs(1).dims(kd)=Unknown
end
end
end
tree.lhs(1).vtype=A.vtype
return
elseif A.vtype==Boolean then
A=convert2double(A)
elseif A.vtype==String then
A=convert2double(A)
end
// repmat(A,m) equivalent to repmat(A,m,m)
if is_a_scalar(m) then
tree.rhs=list(m,m);
tree.name="ones";
tree=Operation(".*.",list(tree,A),tree.lhs)
if typeof(m)=="cste" then
tree.out(1).dims=A.dims
for kd=1:2
tree.out(1).dims(kd)=A.dims(kd)*m.value
if tree.out(1).dims(kd)<0 then
tree.out(1).dims(kd)=Unknown
end
end
end
tree.out(1).type=Type(A.vtype,A.property)
// repmat(A,[m,n,...])
elseif not_a_scalar(m)
tmp=m
if typeof(m)<>"variable" then
tmp=gettempvar()
insert(Equal(list(tmp),m))
end
newrhs=list()
for k=1:size(m.dims)
newrhs(k)=Operation("ext",list(tmp,Cste(k)),list())
end
tree.rhs=newrhs;
tree.name="ones";
tree=Operation(".*.",list(tree,A),tree.lhs)
tree.out(1).type=Type(Double,A.property)
// repmat(A,m) or repmat(A,[m,n,...])
else
tree.name="mtlb_repmat"
tree.lhs(1).type=Type(A.vtype,A.property)
if or(vtype==[String,Boolean]) then
tree.lhs(1).type=Type(vtype,A.property)
end
return
end
// repmat(A,m,n)
else
[A,m,n] = getrhs(tree)
if is_complex(m) then
m=Funcall("real",1,list(m),list(Variable("",m.infer)))
elseif ~is_real(m) then
newm=Funcall("real",1,list(m),list(Variable("",m.infer)))
repl_poss(newm,m,m,"is Real");
m=newm
end
if is_complex(n) then
n=Funcall("real",1,list(n),list(Variable("",n.infer)))
elseif ~is_real(n) then
newn=Funcall("real",1,list(n),list(Variable("",n.infer)))
repl_poss(newn,n,n,"is Real");
n=newn
end
vtype=A.vtype
if A.vtype==Unknown then
tree.name="mtlb_repmat"
if is_a_scalar(A) & typeof(m)=="cste" then
tree.lhs(1).dims(1)=m.value
end
if is_a_scalar(A) & typeof(n)=="cste" then
tree.lhs(1).dims(2)=n.value
end
tree.lhs(1).vtype=A.vtype
return
elseif A.vtype==Boolean then
A=convert2double(A)
elseif A.vtype==String then
A=convert2double(A)
end
tree.rhs=Rhs(m,n)
tree.name="ones";
tree=Operation(".*.",list(tree,A),tree.lhs)
if typeof(m)=="cste" & typeof(n)=="cste" then
tree.out(1).dims=A.dims
tree.out(1).dims(1)=A.dims(1)*m.value
if tree.out(1).dims(1)<0 then
tree.out(1).dims(1)=Unknown
end
tree.out(1).dims(2)=A.dims(2)*n.value
if tree.out(1).dims(2)<0 then
tree.out(1).dims(2)=Unknown
end
end
end
if vtype==String then
tree=Funcall("asciimat",1,list(tree),tree.out);
tree.lhs(1).type=Type(String,Real)
elseif vtype==Boolean then
tree.out(1).type=Type(Boolean,Real)
end
endfunction
|
