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
|
// 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_eig(tree)
// M2SCI function
// Conversion function for Matlab eig()
// Input: tree = Matlab funcall tree
// Ouput: tree = Scilab equivalent for tree
// eig(A)
if rhs==1 then
// Because %c_spec and %b_spec are not defined
A = getrhs(tree)
A = convert2double(A)
tree.rhs=Rhs_tlist(A)
tree.name="spec"
// d = eig(A)
if lhs==1 then
tree.lhs(1).dims=list(A.dims(1),1)
tree.lhs(1).type=Type(Double,Unknown)
// [V,D] = eig(A)
else
tree.lhs(1).dims=list(A.dims(1),A.dims(1))
tree.lhs(1).type=Type(Double,Unknown)
tree.lhs(2).dims=list(A.dims(1),A.dims(1))
tree.lhs(2).type=Type(Double,Unknown)
end
// eig(A,B) or eig(A,'nobalance')
elseif rhs==2 then
[A,B] = getrhs(tree)
// 'nobalance'
if B.vtype==String then
no_equiv(gettext("''nobalance'' option, IGNORED."))
tree.rhs=Rhs_tlist(tree.rhs(1))
if lhs==1 then
tree.name="spec"
tree.lhs(1).dims=list(A.dims(1),1)
tree.lhs(1).type=Type(Double,Unknown)
else
tree.rhs(1)=Operation("+",list(tree.rhs(1),Variable("%i",Infer()),list()))
rhs2=Operation("+",list(Cste(1),Variable("%eps",Infer()),list()))
tree.rhs=Rhs_tlist(tree.rhs(1),rhs2)
tree.name="bdiag"
tree.lhs(1).dims=list(A.dims(1),A.dims(1))
tree.lhs(1).type=Type(Double,Unknown)
tree.lhs(2).dims=list(A.dims(1),A.dims(1))
tree.lhs(2).type=Type(Double,Unknown)
end
return
end
if B.vtype<>Double then
set_infos(gettext("eig() with 2 inputs: consider generalized eigen. Check."),2)
end
m=A.dims(1)
n=A.dims(2)
if lhs==1 then
al=gettempvar(1)
be=gettempvar(2)
insert(Equal(list(al,be),Funcall("spec",1,tree.rhs,list(al,be))))
tree=Operation("./",list(al,be),tree.lhs)
tree.out(1).dims=list(m,1)
tree.out(1).type=Type(Double,Unknown)
else
al=gettempvar(1)
be=gettempvar(2)
insert(Equal(list(al,be,tree.lhs(1)),Funcall("eig",1,tree.rhs,list(al,be,tree.lhs(1)))))
tree=Operation("./",list(al,be),list(tree.lhs(1)))
tree=Funcall("diag",1,list(tree),tree.out)
tree.lhs(1).dims=list(A.dims(1),A.dims(1))
tree.lhs(1).type=Type(Double,Unknown)
tree.lhs(2).dims=list(A.dims(1),A.dims(1))
tree.lhs(2).type=Type(Double,Unknown)
end
// eig(A,B,flag)
else
no_equiv(msprintf(gettext("%s, flag IGNORED."),expression2code(tree)))
tree.rhs(3)=null()
tree.lhs(1).dims=list(A.dims(1),A.dims(1))
tree.lhs(1).type=Type(Double,Unknown)
if lhs==2 then
tree.lhs(2).dims=list(A.dims(1),A.dims(1))
tree.lhs(2).type=Type(Double,Unknown)
end
end
endfunction
|
