How to replace use of cell storage with use of arrays?

Question

POLLY 2018-9-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/420567-how-to-replace-use-of-cell-storage-with-use-of-arrays

评论： Stephen23 2018-9-27

Thanks for taking the time to read this. This is a part of the script I'm working on:

resultX0=cell(length(n),length(q), ntrials);
resultY0=cell(length(n),length(q), ntrials);
resultG=cell(length(n),length(q), ntrials);
for i=1:length(n)
  for j=1:length(q)
      for k=1:1:ntrials
          [G,X0,Y0]=matrix_plantsubm(M,N,c,n(i),p,q(j));
          resultX0{i,j,k}=X0;
          resultY0{i,j,k}=Y0;
          resultG{i,j,k}=G;
      end
  end
end

I would like to replace use of cell storage with use of arrays (e.g. resultG{i,j,k} = G). I want to run the function 10 times (in this case ntrials) for each value q and n. How can I do that?

Thanks

5 个评论
显示 3更早的评论隐藏 3更早的评论

POLLY 2018-9-24

They are the matrices 500*500 that were randomly generated for each i,j,k

James Tursa 2018-9-24

What code do you currently use to save and access the cells?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Nicole Peltier 2018-9-24

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/420567-how-to-replace-use-of-cell-storage-with-use-of-arrays#answer_338230

编辑：Nicole Peltier 2018-9-24

在 MATLAB Online 中打开

According to a previous time you posted about matrix_plantsubm.m on Matlab Answers (thanks Google!), it looks like G, X0, and Y0 are all MxN matrices. Therefore, you could declare resultX0, resultY0, and resultG as follows:

resultX0 = zeros(length(n), length(q), ntrials, M, N);
resultY0 = zeros(length(n), length(q), ntrials, M, N);
resultG = zeros(length(n), length(q), ntrials, M, N);

Then your for-loops should look like this:

for i=1:length(n)
  for j=1:length(q)
      for k=1:ntrials
          [G,X0,Y0]=matrix_plantsubm(M,N,c,n(i),p,q(j));
          resultX0(i,j,k,:,:)=X0;
          resultY0(i,j,k,:,:)=Y0;
          resultG(i,j,k,:,:)=G;
      end
  end
end

Hope this helps!

5 个评论
显示 3更早的评论隐藏 3更早的评论

Nicole Peltier 2018-9-26

在 MATLAB Online 中打开

Is there a reason why this loop must be separate from the previous one? If not, I'd recommend combining the loops like this:

% Declare all variables here
resultX0 = zeros(length(n), length(q), ntrials, M, N);
resultY0 = zeros(length(n), length(q), ntrials, M, N);
resultG = zeros(length(n), length(q), ntrials, M, N);
resultX=zeros(length(n),length(q), ntrials, M, N);
resultY=zeros(length(n),length(q), ntrials, M, N);
for i=1:length(n)
  for j=1:length(q)
    gamma = 6/((q(j) - p)*n(i));
      for k=1:ntrials
        % The following code is from the first loop
        [G,X0,Y0]=matrix_plantsubm(M,N,c,n(i),p,q(j));
        resultX0(i,j,k,:,:)=X0;
        resultY0(i,j,k,:,:)=Y0;
        resultG(i,j,k,:,:)=G;
        % Consolidate code from second loop into first one
        [X,Y,Q, iter] = ADMM(G, c, n(i), gamma, tau, opt_tol, verbose);
        resultX(i,j,k, :, :)=X;
        resultY(i,j,k, :, :)=Y;
      end
    end
  end

This way, you can use the variable G which you've already calculated instead of having to select the correct index out of resultG. (Otherwise you'll have to use resultG(i,j,k,:,:) as your input argument to ADMM.)