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
|
function [tree]=sci_besselk(tree)
// Copyright INRIA
// M2SCI function
// Conversion function for Matlab besselk()
// Input: tree = Matlab funcall tree
// Ouput: tree = Scilab equivalent for tree
// F.B
if size(tree.lhs)==2 then // [y,ierr]=besselk(alpha,x)
error("Bessel function is not implemented for 2 outputs arguments")
end
if rhs==2 then
[A,B]=getrhs(tree) // y=besselk(alpha,x)
elseif rhs==3 then // y=besselk(alpha,x,1)
[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
|
