[<prev] [next>] [<thread-prev] [thread-next>] [day] [month] [year] [list]
Message-ID: <1371006559.9844.250.camel@gandalf.local.home>
Date: Tue, 11 Jun 2013 23:09:19 -0400
From: Steven Rostedt <rostedt@...dmis.org>
To: Juri Lelli <juri.lelli@...il.com>
Cc: fweisbec@...il.com, mingo@...hat.com, linux-kernel@...r.kernel.org,
Chase Douglas <chase.douglas@...onical.com>
Subject: Re: [PATCH 3/3] ftrace: fix stddev calculation
On Tue, 2013-06-11 at 11:08 +0200, Juri Lelli wrote:
> When FUNCTION_GRAPH_TRACER is enabled, ftrace can profile kernel functions
> and print basic statistics about them. Unfortunately, running stddev
> calculation is wrong. This patch corrects it implementing Welford’s method:
>
> s^2 = 1 / (n * (n-1)) * (n * \Sum (x_i)^2 - (\Sum x_i)^2) .
Looking at this further, we only need this fix. We don't need the other
two patches, as that's just verifying the algorithm, and not something
we need to do for run time tests. The run time tests is to test
functionality, not calculations that can be done out of the kernel.
Can you resubmit with just this change. And add the above line as a
comment below.
>
> Signed-off-by: Juri Lelli <juri.lelli@...il.com>
> Cc: Steven Rostedt <rostedt@...dmis.org>
> Cc: Frederic Weisbecker <fweisbec@...il.com>
> Cc: Ingo Molnar <mingo@...hat.com>
> ---
> kernel/trace/ftrace.c | 8 ++++++--
> 1 file changed, 6 insertions(+), 2 deletions(-)
>
> diff --git a/kernel/trace/ftrace.c b/kernel/trace/ftrace.c
> index 6caaa0e..073a328 100644
> --- a/kernel/trace/ftrace.c
> +++ b/kernel/trace/ftrace.c
> @@ -593,13 +593,17 @@ void function_stat_calc(struct ftrace_profile *rec,
> if (rec->counter <= 1)
> *stddev = 0;
> else {
> - *stddev = rec->time_squared - rec->counter * (*avg) * (*avg);
> + /*
> + * Apply Welford's method.
Welford's method is not well known. Please add:
* s^2 = 1 / (n * (n-1)) * (n * \Sum (x_i)^2 - (\Sum x_i)^2)
Thanks!
-- Steve
> + */
> + *stddev = rec->counter * rec->time_squared -
> + rec->time * rec->time;
>
> /*
> * Divide only 1000 for ns^2 -> us^2 conversion.
> * trace_print_graph_duration will divide 1000 again.
> */
> - do_div(*stddev, (rec->counter - 1) * 1000);
> + do_div(*stddev, rec->counter * (rec->counter - 1) * 1000);
> }
> }
>
--
To unsubscribe from this list: send the line "unsubscribe linux-kernel" in
the body of a message to majordomo@...r.kernel.org
More majordomo info at http://vger.kernel.org/majordomo-info.html
Please read the FAQ at http://www.tux.org/lkml/
Powered by blists - more mailing lists