AngDiff.m

function [d, e] = AngDiff(x, y)
%ANGDIFF  Compute angle difference accurately
%
%   [d, e] = ANGDIFF(x, y) computes z = y - x, reduced to (-180,180].  d =
%   round(z) and e = z - round(z).  x and y can be any compatible shapes.

  [d, t] = sumx(AngNormalize(-x), AngNormalize(y));
  d = AngNormalize(d);
  d(d == 180 & t > 0) = -180;
  [d, e] = sumx(d, t);
end