How can I preallocate it?

I want preallocate GGP{F} in my program. How can I do it?

    % 
for F=1:PE*AE*RE
    GGP{F}=zeros(2,12);
    for i=1:12
        n=0;
        [GGP{F}(2,i),GGP{F}(1,i)]=cart2pol(GG{F}(1,i),GG{F}(2,i));
        for j=1:NodesTETA*NodesR*NodesZ
            if abs((GGP{F}(1,i)-p{j}(1))/Radiale)<1 && abs((GGP{F}(2,i)-p{j}(2))/peripherale)<1 && abs((GG{F}(3,i)-p{j}(3))/Axiale)<1
                n=n+1;
                NODZ{F,i}(n)=j;
            end
        end
    end
end

[James Tursa]

What's going on with your cart2pol call? Looks like you are always passing it 0's.

Answer 1

[James Tursa]

GGP = cell(1,PE*AE*RE);

You could move the GGP{F}=zeros(2,12) stuff outside the loop as well, but it wouldn't save you anything to do so, so I would leave that where it is.