docs: private double deviation(double[] delta, double[][] beta) {
      final double a = squareSum(delta);
      final double b = squareSum(transposeTimes(beta, delta));
  
      return (a > b) ? FastMath.sqrt(a - b) : 0.;
    }