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
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2004-2006 - INRIA - Farid BELAHCENE
//
// 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_besseli(tree)
// M2SCI function
// Conversion function for Matlab besseli()
// Input: tree = Matlab funcall tree
// Ouput: tree = Scilab equivalent for tree
if size(tree.lhs)==2 then // [y,ierr]=besseli(alpha,x)
error(gettext("Bessel function is not implemented for 2 outputs arguments."))
end
if rhs==2 then
[A,B]=getrhs(tree) // y=besseli(alpha,x)
elseif rhs==3 then //besseli(alpha,x,)
[A,B,C]=getrhs(tree)
end
A=convert2double(A)
B=convert2double(B)
if is_real(B) then
tree.lhs(1).type=Type(Double,Real)
else
tree.lhs(1).type=Type(Double,Unknown)
end
if is_empty(A) | is_empty(B) then // A or B is an empty
tree.lhs(1).dims=list(0,0)
elseif is_a_scalar(A) then // A is a scalar
tree.lhs(1).dims=B.dims
elseif is_a_scalar(B) then// B is a scalar
tree.lhs(1).dims=A.dims
else
Aknowndim=%t // boolean false if A dimension is unknown
Bknowndim=%t // boolean false if B dimension is unknown
for i=1:size(A.dims)
if A.dims(i)==-1 then
Aknowndim=%f
break
end
end
for i=1:size(B.dims)
if B.dims(i)==-1 then
Bknowndim=%f
break
end
end
if is_a_vector(A) then // A is a vector
if is_a_vector(B) then // B is a vector
if A.dims==B.dims then // both A and B are column vectors or both A and B are line vectors
tree.lhs(1).dims=B.dims
else // A is a line vector and B is a column vector, or A is a column vector and B is a line vector
tree.lhs(1).dims(1)=B.dims(1)*B.dims(2)
tree.lhs(1).dims(2)=A.dims(1)*A.dims(2)
end
end
elseif ~(Aknowndim & Bknowndim) then // A or B dimension(s) is unknown
tree.lhs(1).dims==list(-1,-1)
elseif A.dims==B.dims then
tree.lhs(1).dims=A.dims
else
tree.lhs(1).dims==list(-1,-1)
end
end
endfunction
|
