private double diffNormF(AMatrix tempA, AMatrix tempB) {
   AMatrix temp = tempA.copy();
   temp.sub(tempB);
   double total = temp.elementSquaredSum();
   temp.abs();
   double scale = temp.elementMax();
   return Math.abs(scale - 0) > 1e-12 ? total / scale : 0;
 }