Commit ffff9da7 authored by niklas.baumgarten's avatar niklas.baumgarten
Browse files

preparing parallel averages

parent 3f203357
......@@ -25,6 +25,14 @@ void Averages::Update(Sums sums, int M) {
Q2 = sums.Q2 / M;
Q3 = sums.Q3 / M;
Q4 = sums.Q4 / M;
// Y = abs(PPM->SumOnCommSplit(sums.Y, 0) / M);
// Y2 = PPM->SumOnCommSplit(sums.Y2, 0) / M;
// Y3 = PPM->SumOnCommSplit(sums.Y3, 0) / M;
// Y4 = PPM->SumOnCommSplit(sums.Y4, 0) / M;
// Q = abs(PPM->SumOnCommSplit(sums.Q, 0) / M);
// Q2 = PPM->SumOnCommSplit(sums.Q2, 0) / M;
// Q3 = PPM->SumOnCommSplit(sums.Q3, 0) / M;
// Q4 = PPM->SumOnCommSplit(sums.Q4, 0) / M;
}
void Variances::Update(Averages avgs) {
......@@ -37,6 +45,5 @@ void Kurtosis::Update(Averages avgs, Variances vars) {
3.0 * avgs.Q * avgs.Q * avgs.Q * avgs.Q) / vars.Q / vars.Q;
Y = (avgs.Y4 - 4.0 * avgs.Y3 * avgs.Y + 6.0 * avgs.Y2 * avgs.Y * avgs.Y -
3.0 * avgs.Y * avgs.Y * avgs.Y * avgs.Y) / vars.Y / vars.Y;
if (Y > 100.0)
Warning("Kurtosis of Y above 100!")
if (Y > 100.0) Warning("Kurtosis of Y above 100!")
}
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment