How do I preallocate memory

Question

Chuck Dalby 2024-7-2，16:38

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2133836-how-do-i-preallocate-memory

回答： Rik 2024-7-2，18:23

I get the error message at the end of the first for loop (after code: OutData(ii).hucID = hucID;) that I should consider preallocating memory. How do I do that for the following code:

for ii = 1:length(hucPoly)   
    hucSub      =   polyshape(hucPoly(ii).X, hucPoly(ii).Y);
    hucID       =   hucPoly(ii).Name;
    OutData(ii).hucID   =   hucID;
    
    for i = 1:size(inFiles,1)
        inShapeOld  =   shaperead([dataFolder char(inFiles.OldShape(i))]);
        
        inPolyOld   =   polyshape([inShapeOld.X],[inShapeOld.Y], 'Simplify', false);
        
        inShapeNew  =   shaperead([dataFolder char(inFiles.NewShape(i))]);
        inPolyNew   =   polyshape([inShapeNew.X],[inShapeNew.Y], 'Simplify', false);
        
        inPolyNewsub    =   intersect(inPolyNew, hucSub,'KeepCollinearPoints',true);
        inPolyOldsub    =   intersect(inPolyOld, hucSub);
        [OutData(ii).intsecArea(i).Data, S] = polyIntersect(CoRegErrData, inPolyNewsub, inPolyOldsub);
        OutData(ii).intsecArea(i).ID    =   char(inFiles.NewShape(i));
        if ~isnan(S(1).X(1))
            shapewrite(S, [dataFolder 'outShapes_S\' hucID '_outInt_' char(inFiles.NewShape(i))])
        end
        T   =   table(OutData(ii).intsecArea(i).Data, 'VariableNames', {'outInt_acres1' });
        writetable(T, [dataFolder 'outTables_S\' hucID '_outInt_' char(inFiles.NewShape(i)) '.txt'])
    end
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Voss 2024-7-2，17:04

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2133836-how-do-i-preallocate-memory#answer_1480006

编辑：Voss 2024-7-2，17:05

在 MATLAB Online 中打开

One way is to do this before the outer loop:

clear OutData
OutData(length(hucPoly)) = struct();

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Voss 2024-7-2，17:07

在 MATLAB Online 中打开

Another way is to do this before the outer loop:

OutData = struct('hucID',num2cell(zeros(1,length(hucPoly))),'intsecArea',[]);

请先登录，再进行评论。

Answer 2

Rik 2024-7-2，18:23

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2133836-how-do-i-preallocate-memory#answer_1480061

A third way to do it: loop backwards.

for ii = numel(hucPoly):-1:1

This also replaces the call to the length function. I have yet to see a situation where people actually mean max(size(A)) instead of the number of elements.

But the main point is that for a struct it is much less important to pre-allocate properly.