std::pair<double, double> adjust_for_overtones(std::pair<double, double> base, double *in, unsigned num_samples)
{
double mu = base.first, var = 1.0 / (base.second * base.second);
std::pair<double, double> adjust_for_overtones(std::pair<double, double> base, double *in, unsigned num_samples)
{
double mu = base.first, var = 1.0 / (base.second * base.second);