MATLAB: Preallocation of matrix with undetermined dimensions.

odometerpreallocationunknown dimensions

I need to generate a matrix with a number of dimensions which will be determined by user inputs.

So zeros(3,3,3) gives me a 3x3x3 matrix. My problem is that I need zeros(3,3,3,…..,3) where there are N number of dimensions all of size 3. N is determined by user inputs. Is there a way to efficiently accomplish this, or am I stuck creating a matrix of large dimensions and just under filling it?

Final note: as it is a user-input determined value, I'd really like to avoid using eval() if possible.

Thank you for any assistance you can offer!

Best Answer

zeros(3 * ones(1,N))

Related Solutions

MATLAB: How to avoid for loops when generating index arrays

Here's another method, less memory consuming than NDGRID

    mm=sparse(-n:n);
    ll=sparse(0:n);
    map=bsxfun(@le,abs(mm.'),ll);
    idx=nonzeros(bsxfun(@times,map,1:length(ll)  ));
       l=full(ll(idx));
    idx=nonzeros(bsxfun(@times,map,(1:length(mm)).'))  ;
       m=full(mm(idx));

MATLAB: How to loop through an unknown number of matrix dimensions

You can use linear indexing, if you want to operate on all elements:

X = rand(2,3,4);
Y = zeros(size(X));
for k = 1:numel(X)
  Y(k) = X(k) ^ 2;
end

Of course Y = X .^ 2 would be easier without a loop, but this is thought as demonstration only.

If you want to operate on submatrices, you can move them to the first dimensions temporarily:

X = rand(2,3,4,5,6);
operateOn = [3,5];  % Operate on submatrix along 3rd and 5th dimension
v   = 1:ndims(X);
XX  = permute(X, [operateOn, setdiff(v, operateOn)]);
sXX = size(XX);
XX  = reshape(XX, [sXX(1:2), prod(sXX(3:end)]);
for k = 1:size(XX, 3)
   submatrix = XX(:, :, k);
   ... do what you want
end

Another example is using a cell vector for indexing. If you want a subarray with indexing the the last two dimensions - independent from knowing how many dimensions the input has:

X = rand(2,3,4,5,6);
n = ndims(X);
index        = cell(1, n);
index(:)     = {':'};
index{end-1} = 3;
index{end}   = 4;
Y = X(index{:});  % Equivalent to: X(:, :, :, 3, 4)

You can replace a bunch of nested loops by a single loop also using an index vector instead of scalar indices, see Answers: 333926-recursive-function-for-replacing-multiple-for-loops :

X = rand(2,3,4,5,6);
nv   = ndims(X);
v    = ones(1, nv);
vLim = size(X);
ready = false;
while ~ready
  % Do what you need with X and the index vector v

  ...
    % Update the index vector:
    ready = true;       % Assume that the WHILE loop is ready
    for k = 1:nv
      v(k) = v(k) + 1;
      if v(k) <= vLim(k)
        ready = false;  % No, WHILE loop is not ready now
        break;          % v(k) increased successfully, leave "for k" loop
      end
      v(k) = 1;         % Reset v(k), proceed to next k
    end
  end

This is equivalent to the following without the need to know the number of dimensions before:

for i1 = 1:size(X, 1)
  for i2 = 1:size(X, 2)
    for i3 = 1:size(X, 3)
      for i4 = 1:size(X, 4)
        for i5 = 1:size(X, 5)
          v = [i1, i2, i3, i4, i5];
          % Do what you need with X and the index vector v
          ...
        end
      end
    end
  end
end

See an example: https://www.mathworks.com/matlabcentral/answers/447469-looping-through-a-matrix-of-unknown-dimensions

Best Answer

Related Solutions

MATLAB: How to avoid for loops when generating index arrays

MATLAB: How to loop through an unknown number of matrix dimensions

Related Question