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

private double deviation ( double [ ] delta , double [ ] [ ] beta ) { final double a = squareSum ( delta ) ; final double b = squareSum ( transposeTimes ( beta , delta ) ) ; return ( a > b ) ? Math . sqrt ( a - b ) : 0. ; }
deviation from a manifold described by beta .
