Hi,
I have a problem in speeding my code. In my code, there is a step as below,
a=repmat(1:6,3,1)'; index=[1:3;2:4]; result=zeros(3,3,size(index,1)); for i=1:size(index,1) result(:,:,i)=a(index(i,:),:) end
In practice, the matrix a and index are very large n*3 matrix. a is the node coordinates and the index is the triangle face index of the node. The surface is too big. So I really need to speed this loop. I have an idea that vectorization could make the code faster. But I could not get the ideal output result, even with some matrix "resize" or matrix rotation functions, like "resize" or "reshape". I am very thankful for any suggestion.
Best Answer