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
|
function [tree]=sci_hankel(tree)
// Copyright INRIA
// M2SCI function
// Conversion function for Matlab hankel()
// Input: tree = Matlab funcall tree
// Ouput: tree = Scilab equivalent for tree
// V.C.
tree.name="hank"
// H = hankel(c)
if rhs==1 then
// %c_hank and %b_hank are not defined
if tree.rhs(1).vtype<>Double then
tree.rhs(1)=convert2double(tree.rhs(1))
end
if typeof(tree.rhs(1))<>"variable" then
c=gettempvar();
insert(Equal(list(c),tree.rhs(1)))
else
c=tree.rhs(1);
end
n=gettempvar()
insert(Equal(list(n),Funcall("size",1,Rhs(c,"*"),list())));
if is_a_scalar(tree.rhs(1)) then
v=c
else
// v=[matrix(c,1,-1),zeros(1,n)]
matrix_funcall=Funcall("matrix",1,Rhs(c,1,-1),list())
zeros_funcall=Funcall("zeros",1,Rhs(1,n),list())
v=Operation("rc",list(matrix_funcall,zeros_funcall),list())
end
dim=1
for kd=1:lstsize(tree.rhs(1).dims)
if tree.rhs(1).dims(kd)<>Unknown then
dim=dim*tree.rhs(1).dims(kd)
else
dim=Unknown
break
end
end
tree.lhs(1).dims=list(dim,dim)
tree.lhs(1).type=tree.rhs(1).type
tree.rhs=Rhs(n,n,v);
// H = hankel(c,r)
else
// %c_hank and %b_hank are not defined
if tree.rhs(1).vtype<>Double then
tree.rhs(1)=convert2double(tree.rhs(1))
end
if typeof(tree.rhs(1))<>"variable" then
c=gettempvar();
insert(Equal(list(c),tree.rhs(1)))
else
c=tree.rhs(1);
end
if tree.rhs(2).vtype<>Double then
tree.rhs(2)=convert2double(tree.rhs(2))
end
if typeof(tree.rhs(2))<>"variable" then
r=gettempvar();
insert(Equal(list(r),tree.rhs(2)))
else
r=tree.rhs(2);
end
m=gettempvar();
n=gettempvar();
insert(Equal(list(m),Funcall("size",1,Rhs(c,"*"),list())));
insert(Equal(list(n),Funcall("size",1,Rhs(r,"*"),list())));
// v=[matrix(c,1,-1),matrix(r(2:$),1,-1)]
cmatrix_funcall=Funcall("matrix",1,Rhs(c,1,-1),list())
imp=Operation(":",list(Cste(2),Cste("$")),list())
ext=Operation("ext",list(r,imp),list())
rmatrix_funcall=Funcall("matrix",1,Rhs(ext,1,-1),list())
v=Operation("rc",list(cmatrix_funcall,rmatrix_funcall),list())
dim1=1
for kd=1:lstsize(tree.rhs(1).dims)
if tree.rhs(1).dims(kd)<>Unknown then
dim1=dim1*tree.rhs(1).dims(kd)
else
dim1=Unknown
break
end
end
dim2=1
for kd=1:lstsize(tree.rhs(2).dims)
if tree.rhs(2).dims(kd)<>Unknown then
dim2=dim2*tree.rhs(2).dims(kd)
else
dim2=Unknown
break
end
end
tree.lhs(1).dims=list(dim1,dim2)
tree.lhs(1).type=tree.rhs(1).type
tree.rhs=Rhs(m,n,v);
end
endfunction
|
