sub ediag {
  abs(shift) * abs(shift);
}