release: 0.1.57

[lilypond.git] / lily / leastsquares.cc

diff --git a/lily/leastsquares.cc b/lily/leastsquares.cc
index 04dd26bdcc2d6cbf873d3b8d259575852c1c6490..0f1568f29f82f44f2ac714b2e1d064666a685dbd 100644 (file)
--- a/lily/leastsquares.cc
+++ b/lily/leastsquares.cc
@@ -1,35 +1,38 @@
 #include "leastsquares.hh"
+
 void
 Least_squares::OK() const
 {
-    assert( input.size() > 1 );
-    Real dx = 0.0;
-    for (int i=1; i < input.size(); i++)
-       dx += abs (input[i-1].x() - input[i].x());
+  assert (input.size() > 1);
+  Real dx = 0.0;
+  for (int i=1; i < input.size(); i++)
+       dx += abs (input[i-1].x() - input[i].x ());
 
-    assert(dx);
+  assert (dx);
 }
+
 void
-Least_squares::minimise(Real &coef, Real &offset)
+Least_squares::minimise (Real &coef, Real &offset)
 {
-    OK();
-    Real sx = 0.0;
-    Real sy = 0.0;
-    Real sqx =0.0;
-    Real sxy = 0.0;
+  OK();
+  Real sx = 0.0;
+  Real sy = 0.0;
+  Real sqx =0.0;
+  Real sxy = 0.0;
 
-    for (int i=0; i < input.size();i++) {
+  for (int i=0; i < input.size();i++) 
+    {
        Real x=input[i].x();
        Real y = input[i].y();
        sx += x;
        sy += y;
-       sqx += sqr(x);
+       sqx += sqr (x);
        sxy += x*y;
     }
-    int N = input.size();
-    
+  int N = input.size();
+  
 
-    coef = (N * sxy - sx*sy )/(N*sqx - sqr(sx));
-    offset = (sy - coef * sx)/N;
+  coef = (N * sxy - sx*sy)/(N*sqx - sqr (sx));
+  offset = (sy - coef * sx)/N;
        
 }