X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=lily%2Fleast-squares.cc;h=9c250bea915e7590396ef1bb071ec06ddf5ab413;hb=a6a51abfd0195a3cf7d6ea095cf69808852f21ce;hp=d95279a54e0d801420b28b6348f3b0f76f2583d1;hpb=f9214bac21e9926dc3248416f58190c98c4167a9;p=lilypond.git diff --git a/lily/least-squares.cc b/lily/least-squares.cc index d95279a54e..9c250bea91 100644 --- a/lily/least-squares.cc +++ b/lily/least-squares.cc @@ -1,9 +1,20 @@ /* - least-squares.cc -- implement minimise_least_squares + This file is part of LilyPond, the GNU music typesetter. - source file of the GNU LilyPond music typesetter + Copyright (C) 1996--2015 Han-Wen Nienhuys - (c) 1996--2007 Han-Wen Nienhuys + 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 . */ #include "least-squares.hh" @@ -12,21 +23,21 @@ void minimise_least_squares (Real *coef, Real *offset, - vector const &input) + vector const &input) { Real sx = 0.0; Real sy = 0.0; Real sqx = 0.0; Real sxy = 0.0; - for (vsize 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 count = input.size (); @@ -37,7 +48,9 @@ minimise_least_squares (Real *coef, Real *offset, 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 = count ? sy / count : 0.0; }