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
|
// 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 [M] = %hm_e(varargin)
// version modified by Bruno
// extract an sub_hypermatrix
[lhs,rhs]=argn(0)
M=varargin(rhs)
dims=double(M.dims) // new hypermatrices have an int32 dims field
nbdims = size(dims,"*");
nb_index_vect = rhs-1;
if nb_index_vect > nbdims then
error(msprintf(_("%s: Too many subscripts.\n"),"%hm_e"));
elseif nb_index_vect < nbdims then
nbdims = nb_index_vect;
dims(nbdims) = prod(dims(nbdims:$))
end
I = 0
dims1 = [];
for k=nb_index_vect:-1:1
ik=varargin(k)
type_ik = type(ik)
if type_ik==1 then // usual subscript
if min(size(ik))<0 then
if rhs==2 then // subscript is : in fact
M=matrix(M.entries,-1,1)
return
end
ik=1:dims(k)
else
ik=round(ik)
if min(ik) <= 0 then, error(21), end
if max(ik) > dims(k) then, error(21), end
end
elseif type_ik==2 then // poly subscript
ik=round(horner(ik,dims(k)))
if min(ik) <= 0 then, error(21), end
if max(ik) > dims(k) then, error(21), end
elseif type_ik == 129 then // implicit poly subscript (p.e. 1:$)
ik=round(horner(ik(:),dims(k)))
if min(ik(1),ik(3)) <= 0 then, error(21), end
if max(ik(1),ik(3)) > dims(k) then, error(21), end
ik = ik(1):ik(2):ik(3)
elseif type_ik==4 | type_ik==6 then // boolean and sparse boolean subscript
if size(ik,"*") ~= dims(k) then, error(21), end
ik=find(ik)
elseif typeof(ik) == "hypermat" then // hm boolean subscript
if type(ik.entries) ~= 4 then, error(21), end
if size(ik,"*") ~= dims(k) then, error(21), end
ik=find(ik.entries)
else
error(21)
end
nik = size(ik, "*");
if (nik == 0) then, M=[], return, end
dims1 = [nik, dims1]
if nik > 1 then
ik=ik(:)
if size(I,"*") > 1 then
I=(dims(k)*I).*.ones(ik)+ones(I).*.(ik-1)
else
I=dims(k)*I+ik-1
end
else
I=dims(k)*I+ik-1
end
end
//
dims1(max(find(dims1>1))+1:$)=[] // elimine la fin si les dims sont 1
// exemple dims1 = [2 4 1 6 7 1 1] alors on obtient dims1 = [2 4 1 6 7]
// (pour compatibilite avec Matlab)
select size(dims1,"*")
case 0
M = M.entries(I+1)
case 1
M = M.entries(I+1)
case 2
M = matrix(M.entries(I+1),dims1(1),dims1(2))
else
M = mlist(["hm","dims","entries"],int32(dims1),M.entries(I+1))
end
endfunction
|
