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
|
// 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_lu(tree)
// M2SCI function
// Conversion function for Matlab lu()
// Input: tree = Matlab funcall tree
// Ouput: tree = Scilab equivalent for tree
// lu(X)
if rhs==1 then
X = getrhs(tree)
X = convert2double(X)
tree.rhs=Rhs_tlist(X)
// Y = lu(X)
if lhs==1 then
if tree.lhs(1).name=="ans" then
no_equiv(expression2code(tree))
else
no_equiv(tree.lhs(1).name+" = "+expression2code(tree))
end
// [L,U] = lu(X)
elseif lhs==2 then
tree.lhs(1).dims=list(Unknown,Unknown)
tree.lhs(1).type=Type(X.vtype,Real)
tree.lhs(2).dims=list(Unknown,Unknown)
tree.lhs(2).type=X.type
// [L,U,P] = lu(X)
elseif lhs==3 then
tree.lhs(1).dims=list(Unknown,Unknown)
tree.lhs(1).type=Type(X.vtype,Real)
tree.lhs(2).dims=list(Unknown,Unknown)
tree.lhs(2).type=X.type
tree.lhs(3).dims=list(Unknown,Unknown)
tree.lhs(3).type=Type(X.vtype,Real)
// [L,U,P,Q] = lu(X)
else
no_equiv("["+tree.lhs(1).name+","+tree.lhs(2).name+","+tree.lhs(3).name+"] = "+expression2code(tree))
end
// lu(X,thresh)
else
no_equiv(expression2code(tree))
end
endfunction
|
