function getdctfreturn = getdctf(f)
  for u = 0:7
    for v = 0:7
      F(u + 1,v + 1) = dctf(u, v, f);
      cgret = F(u + 1,v + 1);
    end
  end
  getdctfreturn = F;
