]> git.donarmstrong.com Git - lilypond.git/blobdiff - lily/least-squares.cc
*** empty log message ***
[lilypond.git] / lily / least-squares.cc
index d589ea06fcf348f2b74f98f49fbdc18f95b52a79..e7f6b3058499bac455b5d39559a11e885d2ad591 100644 (file)
@@ -1,18 +1,17 @@
-/*   
+/*
   least-squares.cc -- implement minimise_least_squares
-  
+
   source file of the GNU LilyPond music typesetter
-  
+
   (c) 1996--2005 Han-Wen Nienhuys <hanwen@cs.uu.nl>
-  
- */
+*/
 
 #include "least-squares.hh"
 
 #include "warn.hh"
 
 void
-minimise_least_squares (Real * coef, Real * offset,
+minimise_least_squares (Real *coef, Real *offset,
                        Array<Offset> const &input)
 {
   Real sx = 0.0;
@@ -20,7 +19,7 @@ minimise_least_squares (Real * coef, Real * offset,
   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_AXIS];
       Real y = input[i][Y_AXIS];
@@ -33,18 +32,18 @@ minimise_least_squares (Real * coef, Real * offset,
 
   *coef = 0.0;
   *offset = 0.;
-  
-  Real den = (N*sqx - sqr (sx));
+
+  Real den = (N *sqx - sqr (sx));
   if (!N || !den)
     {
       programming_error ("minimise_least_squares ():  Nothing to minimise");
       *coef = 0.0;
-      *offset = N ? sy/N : 0.0;
+      *offset = N ? sy / N : 0.0;
     }
   else
     {
-      *coef = (N * sxy - sx*sy)/den;
-      *offset = (sy - (*coef) * sx)/N;
+      *coef = (N *sxy - sx*sy) / den;
+      *offset = (sy - (*coef) * sx) / N;
     }
 }