/usr/include/ableton/link/LinearRegression.hpp is in ableton-link-dev 1.0.0+dfsg-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 | /* Copyright 2016, Ableton AG, Berlin. All rights reserved.
*
* This program 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 2 of the License, or
* (at your option) any later version.
*
* This program 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 this program. If not, see <http://www.gnu.org/licenses/>.
*
* If you would like to incorporate Link into a proprietary software application,
* please contact <link-devs@ableton.com>.
*/
#pragma once
#include <cfloat>
#include <cmath>
#include <numeric>
#include <utility>
namespace ableton
{
namespace link
{
template <typename It>
std::pair<double, double> linearRegression(It begin, It end)
{
using namespace std;
using Point = pair<double, double>;
const double numPoints = static_cast<double>(distance(begin, end));
const double meanX =
accumulate(begin, end, 0.0, [](double a, Point b) { return a + b.first; })
/ numPoints;
const double productXX = accumulate(begin, end, 0.0,
[&meanX](double a, Point b) { return a + pow(b.first - meanX, 2.0); });
const double meanY =
accumulate(begin, end, 0.0, [](double a, Point b) { return a + b.second; })
/ numPoints;
const double productXY =
inner_product(begin, end, begin, 0.0, [](double a, double b) { return a + b; },
[&meanX, &meanY](Point a, Point b) {
return ((a.first - meanX) * (b.second - meanY));
});
const double slope = productXX == 0.0 ? 0.0 : productXY / productXX;
const double intercept = meanY - (slope * meanX);
return make_pair(slope, intercept);
}
} // namespace link
} // namespace ableton
|