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
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 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.1-en.txt
function M=%hm_i_hm(varargin)
//insertion of a matrix in an hypermatrix
[lhs,rhs]=argn(0)
M=varargin(rhs)
N=varargin(rhs-1)//inserted matrix
dims=matrix(double(M.dims),-1,1);
v=matrix(M.entries,-1,1);
nd=size(dims,"*")
if rhs-2>nd then dims(nd+1:rhs-2)=1;end
//convert N-dimensional indexes to 1-D
[Ndims,I]=convertindex(list(dims,double(matrix(N.dims,1,-1))),varargin(1:$-2))
if or(Ndims>dims) then
//extend the destination matrix
I1=0
for k=size(Ndims,"*"):-1:1
ik1=(1:dims(k))'
if ik1<>[] then
if Ndims(k)>1 then
if size(I1,"*")>1 then
I1=(Ndims(k)*I1).*.ones(ik1)+ones(I1).*.(ik1-1)
else
I1=Ndims(k)*I1+ik1-1
end
else
I1=Ndims(k)*I1+ik1-1
end
end
end
select type(v)
case 1
v1=zeros(prod(Ndims),1)
case 2 then
v1=zeros(prod(Ndims),1)
case 4 then
v1=(zeros(prod(Ndims),1)==1)
case 8 then
v1=iconvert(zeros(prod(Ndims),1),inttype(v))
case 10 then
v1=emptystr(prod(Ndims),1)
end
v1(I1+1)=v;v=v1
end
v(I)=matrix(N.entries,-1,1)
while Ndims($)==1 then Ndims($)=[],end
select size(Ndims,"*")
case 0
M=v
case 1
M=v
case 2
M=matrix(v,Ndims(1),Ndims(2))
else
M=mlist(["hm","dims","entries"],int32(matrix(Ndims,1,-1)),v)
end
endfunction
|
