summaryrefslogtreecommitdiff
path: root/modules/compatibility_functions/macros/enlarge_shape.sci
blob: c9208d07826975b7b2329da01f706d96bc99ab5b (plain)
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
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2001 - INRIA - Serge STEER
//
// 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 v1=enlarge_shape(dims,dims1,v)
    //Utility function for cell insertion
    //dims :  vector of old dimensions
    //dims1 : vector of new dimensions
    //v     : list of entries
    //v1    : new list of entries
    //
    // undefined elements are set to []
    // and(dims1>=dims) must be true
    ki=find(dims1>dims)
    if ki==nd then // the last dimension -> add empty entries at the end
        for k=prod(dims)+1:prod(dims1),v(k)=[];v1=v;end
    else // enlarge and merge
        v1=list();for k=1:prod(dims1),v1(k)=[],end
        //merge compute the new index of M entries into enlarged cell
        I1=[];
        for k=size(dims1,"*"):-1:1
            ik1=(1:dims(k))'
            if ik1<>[] then
                if dims1(k)>1 then
                    if size(I1,"*")>1 then
                        I1=(dims1(k)*I1).*.ones(ik1)+ones(I1).*.(ik1-1)
                    else
                        I1=dims1(k)*I1+ik1-1
                    end
                else
                    I1=dims1(k)*I1+ik1-1
                end
            end
        end
        for iw=1:size(I1,"*"),v1(I1(iw)+1)=v(iw);end
    end
endfunction