function S = itersum(K,phi,m,n)
dims = repmat([m],[1 n]);
dims1 = repmat([m],[1 n-1]);

%X = list_all_ind(dims);
X1 = list_all_ind(dims1);

S = 0;
for i=1:size(X1,1)
  y = X1(i,:);
  Xy = [[1:m]' repmat(y,[m 1])];
  Xyi = coord2ind(Xy,dims);
  Fx = phi(Xyi);
  Kx = K(Xyi);
  
  lam = min(Fx);
  %S = S + lam*max(sum(Kx)) + sum(Kx.*(Kx>0));
  S = S + lam*max(0,sum(Kx)) + sum(Kx.*(Kx>0));
end

return