P_STRAIGHT_ULOOP(16, {
  pd = (ps1 + ps2) >> 1;
}, {
  pd = (ps1 - ps2) >> 1;
})