Scale Average

scale_mean.cpp

MVector scaleAverage(const std::vector<MVector>& values)
{
    if (values.empty()) return MVector::zero;
    const MVector sum = std::accumulate(values.begin(), values.end(), MVector::zero,
                                        [](const MVector& a, const MVector& b)
                                        {
                                            // assumes [Vx, Vy, Vz] > 0
                                            return MVector(a.x + std::log(b.x), a.y + std::log(b.y), a.z + std::log(b.z));
                                        });
    const MVector average = sum / values.length();
    return MVector(std::exp(average.x), std::exp(average.y), std::exp(average.z))
}