function f = gfeat(gam)
gam = gam(:);
m = round(sqrt(length(gam)));
gam = reshape(gam,[m m]);

totsum = sum(abs(sum(gam')));
maxsum = max(abs(sum(gam')));
abssum = abs(sum(sum(gam)));

sumj = abs(sum(gam'));

[tt,tm,mm,dj,tj] = dish(gam);

%f = [totsum maxsum tt tm mm]';
%f = [tt tm mm]';
%f = [tt mm]';
%f = [tt]';
%f = [tt]';
%f = [mm]';
%f = [tm]';

%f = tt + tm + mm;
%f = totsum + maxsum;
%f = totsum + tt;
%f = totsum;
%f = totsum + maxsum + tt + tm + mm;
%f = totsum  + tm;
%f = totsum + tt;  % this one has potential
%f = abssum + tm;  % looks like this one is it *more likely
%f = abssum + mm; % this one could be it too
%f = abssum + tt;
%f = abssum + tt - tm;

%f = totsum + mm;
%sumj = sumj(:);
%dj = dj(:);
%f = max(sumj + dj);
%f = max(sumj) + max(dj);
%dj
%tj
%f = sum(dj);
%f = abs(sum(diag(gam)));
%f = crossums(gam);
f = abs( gam(1,1) + gam(1,2) + gam(1,3) );
return


function [tt,tm,mm,dj,tj] = dish(gam)
tt = 0;
tm = 0;
mm = 0;
dj = [];
tj = [];

for col = 1:size(gam,2)
  g = gam(:,col);
  dg = 0;
  tj(col) = 0;
  for i=1:length(g)
    for j=i+1:length(g)
      tt = tt + abs(g(i)-g(j));
      tj(col) = tj(col) + abs(g(i)-g(j));
      dg = max(dg,abs(g(i)-g(j)));
    end
  end
  tm = tm + dg;
  mm = max(mm,dg);
  dj(col) = dg;
end

return