source file of the GNU LilyPond music typesetter
- (c) 1996--2005 Han-Wen Nienhuys <hanwen@xs4all.nl>
+ (c) 1996--2006 Han-Wen Nienhuys <hanwen@xs4all.nl>
*/
#include "least-squares.hh"
void
minimise_least_squares (Real *coef, Real *offset,
- Array<Offset> const &input)
+ 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];
sqx += sqr (x);
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");
*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;
}
}