-/*
- 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>
-
- */
+/*
+ This file is part of LilyPond, the GNU music typesetter.
+
+ Copyright (C) 1996--2015 Han-Wen Nienhuys <hanwen@xs4all.nl>
+
+ LilyPond is free software: you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ LilyPond is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with LilyPond. If not, see <http://www.gnu.org/licenses/>.
+*/
#include "least-squares.hh"
#include "warn.hh"
void
-minimise_least_squares (Real * coef, Real * offset,
- Array<Offset> const &input)
+minimise_least_squares (Real *coef, Real *offset,
+ vector<Offset> const &input)
{
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 (vsize i = 0; i < input.size (); i++)
{
Real x = input[i][X_AXIS];
Real y = input[i][Y_AXIS];
sx += x;
sy += y;
sqx += sqr (x);
- sxy += x*y;
+ sxy += x * y;
}
- int N = input.size ();
+
+ int count = input.size ();
*coef = 0.0;
*offset = 0.;
-
- Real den = (N*sqx - sqr (sx));
- if (!N || !den)
+
+ Real den = (count * sqx - sqr (sx));
+ if (!count || !den)
{
- programming_error ("minimise_least_squares (): Nothing to minimise");
+ programming_error ("minimise_least_squares (): Nothing to minimise\n"
+ "This means that vertical spacing is triggered\n"
+ "before line breaking\n");
*coef = 0.0;
- *offset = N ? sy/N : 0.0;
+ *offset = count ? sy / count : 0.0;
}
else
{
- *coef = (N * sxy - sx*sy)/den;
- *offset = (sy - (*coef) * sx)/N;
+ *coef = (count * sxy - sx * sy) / den;
+ *offset = (sy - (*coef) * sx) / count;
}
}