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
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) INRIA - S. Steer, F. Delebecque, V. 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-en.txt
function M=createstruct(index,N)
// Create a struct with field pointed by index equal to N
// If index is not a list, creation is done by routines
if type(index)==15 then
// When only two indexes
if size(index)==2 then
// Second index is a fieldname
if type(index(2))==10 then
// First index is a numerical value or a vector of numerical values
if type(index(1))==1 then
// X(p).f = y -> index=list(p,f)
I=index(1)
if mtlb_mode() then
Dims=[1 max(I)]
else
Dims=[max(I) 1]
end
nmax=prod(Dims);
if nmax==1 then
M=mlist(["st","dims",index(2)],int32(Dims),N);
else
M=mlist(["st","dims",index(2)],int32(Dims));
//initialize the structure with empty fields
Li=list();for kl=1:nmax, Li(kl)=[];end
//set fields pointed to by index(1) to N
for kl=1:size(I,'*'), Li(I(kl))=N;end
setfield(3,Li,M);
end
return;
// First index is a list of numerical values
elseif type(index(1))==15
// X(p,q[,...]).f=z -> index=list(list(p,q[,...]),f)
I=index(1)
Ndims=lstsize(I) //the number of dimensions
//compute the struct dimensions from max indices values
Dims=[];for kDims=1:Ndims, Dims=[Dims,max(I(kDims)) ];end
//initialize the structure with empty fields
Li=list();for kl=1:prod(Dims),Li(kl)=[];end
//set fields pointed to by I to N
I1=sub2ind(Dims,I(:)) // build one dimensionnal index equivalent to I
for kl=1:size(I1,'*'), Li(I1(kl))=N;end
M=mlist(["st","dims",index(2)],int32(Dims),Li);
return;
// First index is also a fieldname
else
// X.f.g=y -> index=list(f,g)
M=mlist(["st","dims",index(1)],int32([1,1]),...
mlist(["st","dims",index(2)],int32([1,1]),N));
return;
end
// Second index is a fieldname
else
// X.f(p[,q,...])=y : index=list(f,list(p[,q,...]))
if typeof(N)=="st" then // When recursive call of createstruct
if type(index(2))==15 then // index=list(f,list(p,q[,...]))
Dims=list2vec(index(2))'
else // index=list(f,p)
if mtlb_mode() then
Dims=[1 index(2)]
else
Dims=[index(2) 1]
end
end
kmax=prod(Dims)
z=list()
for kl=1:kmax
z(kl)=[]
end
z(kmax)=getfield(3,N)
z=mlist(getfield(1,N),int32(Dims),z);
else
z(index(2)(:))=N;
end
if type(index(1))<>10 then pause,end
M=mlist(["st","dims",index(1)],int32([1,1]),z);
return;
end
// Any number of indexes <> 2
else
// Last index is a fieldname -> init struct
if type(index($))==10 then
M=mlist(["st","dims",index($)],int32([1,1]),N);
index($)=null();
// Last index is a numerical value or a list of numerical values
elseif type(index($))==1 | type(index($))==15 then
// If only one index value (happens when recursive calls to createstruct)
if size(index)==1 then
if typeof(N)=="st" then
if type(index($))==15 then
Dims=list2vec(index($))'
else
if mtlb_mode() then
Dims=[1 index($)]
else
Dims=[index($) 1]
end
end
kmax=prod(Dims)
if kmax==1 then
M=N
else
z=list()
for kl=1:kmax
z(kl)=[]
end
z(kmax)=getfield(3,N)
M=mlist(getfield(1,N),int32(Dims),z);
end
else
M(index($)(:))=N;
end
index($)=null()
// More than one index value
else
if typeof(N)=="st" then
if type(index($))==15 then
Dims=list2vec(index($))'
else
if mtlb_mode() then
Dims=[1 index($)]
else
Dims=[index($) 1]
end
end
kmax=prod(Dims)
z=list()
for kl=1:kmax
z(kl)=[]
end
z(kmax)=getfield(3,N)
z=mlist(getfield(1,N),int32(Dims),z);
else
z(index($)(:))=N;
end
M=mlist(["st","dims",index($-1)],int32([1,1]),z);
index($)=null()
index($)=null()
end
else
// This case should not happen (unknown kind of index)
error("Not implemented");
end
if index==list() then
return
end
M=createstruct(index,M);
return;
end
else
if type(index)==10 then
M=mlist(["st","dims",index($)],int32([1,1]),N);
else
error("Not implemented");
end
end
endfunction
|
