Multidimensional Array Indexing to Avoid Nested for loops

I have to optimize this code avoiding for loops:

for i=1:length(PFVFId.poles)
  [T,D]=eig(PFVFId.residue(:,:,i));
  for k=1:M
      temp = zeros(M,M);
      temp(k,k) = D(k,k);
      PFVFId.residuediag(:,:,i,k) = T*temp*T';
  end
end

How can I make this possible?