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
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) INRIA
//
// 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.1-en.txt
function y=%cblock_e(varargin)
x=varargin($)
dims=[];
for k=2:length(x)
dims=[dims;size(getfield(k,x))];
end
m=dims(1,1) // common row number
dims=cumsum([1 dims(:,2)']);
n=dims($)-1 // number of columns
y=mlist("cblock")
j=varargin($-1)
if type(j)==2|type(j)==129 then
j=horner(j,n)
elseif type(j)==4 then
j=find(j)
elseif and(size(j)==[-1,-1]) then
j=1:n
end
if length(varargin)==2 then //y=x(i)
if m<>1 then
error(msprintf(_("%s: Syntax x(j) only allowed for row cblock''s.\n"),"cblock_e"));
end
i=1
else //y=x(i,j)
i=varargin(1)
end
first=%t
for k=1:size(j,"*")
jk=j(k)
I=find(jk>=dims(1:$-1)&jk<dims(2:$))
if I>1 then jk=jk-dims(I)+1,end
v=getfield(I+1,x)
if first then
temp=v(i,jk)
first=%f
else
if type(v)==type(temp) then
temp=[temp v(i,jk)]
else
setfield(length(y)+1,temp,y)
temp=v(i,jk)
end
end
end
if length(y)==1 then
y=temp
else
setfield(length(y)+1,temp,y)
end
endfunction
|
