/******************************************************************
- The standard deviation calculation is based on:
+- CPUID: processor ID
+
+- IrqId: IRQ ID
+
+- Frequency (Hz): the number of interruptions per second (Hz)
+
+- Total Duration (nsec): the sum of each interrupt duration in nsec
+
+- Duration Standard_deviation = sqrt(1/N Sum ((xi -Xa)^2)) where
+ N: number of interrupts
+ xi: duration of an interrupt (nsec)
+ Xa: average duration (nsec)
+
+- Max IRQ handler duration (nsec) [time interval]: the longest IRQ handler duration in nsec.
+
+-Average period (nsec): 1/Frequency(in HZ)
+
+-Period Standard_deviation = sqrt(1/N Sum ((xi -Xa)^2)) where
+N: number of interruptions
+xi: duration of an interrupt
+Xa: 1/Frequency (in Hz)
+
+-Frequency Standard_deviation = sqrt(1/N Sum ((xi -Xa)^2))
+N: number of interruptions
+xi: duration of an interrupt
+Xa: Frequency (Hz)
+
+
+
+The standard deviation calculation is based on:
http://en.wikipedia.org/wiki/Standard_deviation
Standard_deviation = sqrt(1/N Sum ((xi -Xa)^2))
+
+
To compute the standard deviation, we need to make EventRequests to LTTV. In
the first EventRequest, we compute the average duration (Xa) and the
- frequency (N) of each IrqID. We store the information calculated in the first
+ Number of interruptions (N) of each IrqID. We store the information calculated in the first
EventRequest in an array called FirstRequestIrqExit.
In the second EventRequest, we compute the Sum ((xi -Xa)^2) and store this information
- in a array called SumArray. The function CalculateStandardDeviation() uses FirstRequestIrqExit
+ in a array called SumArray. The function CalculateDurationStandardDeviation() uses FirstRequestIrqExit
and SumArray arrays to calculate the standard deviation.
+
+
*******************************************************************/
LttTime duration;
LttTime start_time;
LttTime end_time;
-}MaxDuration;
+}IrqDuration;
typedef struct {
guint cpu_id;
guint id;
- guint frequency;
+ guint NumerofInterruptions;
LttTime total_duration;
guint average_duration;
- MaxDuration max_irq_handler;
-
+ IrqDuration max_irq_handler;
+ IrqDuration min_irq_handler;
}Irq;
typedef struct {
typedef struct
{
guint irqId;
- guint frequency;
- guint64 sumOfDurations;
+ guint NumerofInterruptions;//frequency;//
+ guint64 sumOfDurations; // to store the Sum ((xi -Xa)^2) of the duration Standard deviation
+ guint64 sumOfPeriods; // to store the Sum ((xi -Xa)^2) of the period Standard deviation
+ guint64 sumOfFrequencies;// to store the Sum ((xi -Xa)^2) of the frequency Standard deviation
}SumId;
static gboolean trace_header(void *hook_data, void *call_data);
static gboolean DisplayViewer (void *hook_data, void *call_data);
static void CalculateData(LttTime time_exit, guint cpu_id, InterruptEventData *event_data);
-static void TotalDurationMaxIrqDuration(irq_entry *e, LttTime time_exit, GArray *FirstRequestIrqExit);
+static void CalculateTotalDurationAndMaxIrqDurationAndMinIrqDuration(irq_entry *e, LttTime time_exit, GArray *FirstRequestIrqExit);
static gboolean FirstRequestIrqEntryCallback(void *hook_data, void *call_data);
static gboolean FirstRequestIrqExitCallback(void *hook_data, void *call_data);
static gboolean SecondRequest(void *hook_data, void *call_data);
static gboolean SecondRequestIrqExitCallback(void *hook_data, void *call_data);
static void CalculateXi(LttEvent *event, InterruptEventData *event_data);
static void SumItems(gint irq_id, LttTime Xi, InterruptEventData *event_data);
-static int CalculateStandardDeviation(gint id, InterruptEventData *event_data);
-static int FrequencyInHZ(gint frequency, TimeWindow time_window);
+static int CalculateDurationStandardDeviation(gint id, InterruptEventData *event_data);
+static int CalculatePeriodStandardDeviation(gint id, InterruptEventData *event_data);
+static int FrequencyInHZ(gint NumerofInterruptions, TimeWindow time_window);
+static guint64 CalculatePeriodInnerPart(guint Xi, guint FrequencyHZ);
+static guint64 CalculateFrequencyInnerPart(guint Xi_in_ns, guint FrequencyHZ);
+static void InterruptFree(InterruptEventData *event_viewer_data);
+static int CalculateFrequencyStandardDeviation(gint id, InterruptEventData *event_data);
/* Enumeration of the columns */
enum{
DURATION_STANDARD_DEV_COLUMN,
MAX_IRQ_HANDLER_COLUMN,
AVERAGE_PERIOD,
- PERIOD_STANDARD_DEV_COLUMN,
+ PERIOD_STANDARD_DEV_COLUMN,
+ FREQUENCY_STANDARD_DEV_COLUMN,
N_COLUMNS
};
G_TYPE_UINT64, /* Duration */
G_TYPE_INT, /* standard deviation */
G_TYPE_STRING, /* Max IRQ handler */
- G_TYPE_DOUBLE, /* Average period */
- G_TYPE_DOUBLE /* period standard deviation */
+ G_TYPE_INT, /* Average period */
+ G_TYPE_INT, /* period standard deviation */
+ G_TYPE_INT /* frequency standard deviation */
+
);
event_viewer_data->TreeView = gtk_tree_view_new_with_model (GTK_TREE_MODEL (event_viewer_data->ListStore));
g_object_unref (G_OBJECT (event_viewer_data->ListStore));
renderer = gtk_cell_renderer_text_new ();
- column = gtk_tree_view_column_new_with_attributes ("CPUID",
+ column = gtk_tree_view_column_new_with_attributes ("CPU ID",
renderer,
"text", CPUID_COLUMN,
NULL);
renderer = gtk_cell_renderer_text_new ();
- column = gtk_tree_view_column_new_with_attributes ("IrqId",
+ column = gtk_tree_view_column_new_with_attributes ("IRQ ID",
renderer,
"text", IRQ_ID_COLUMN,
NULL);
gtk_tree_view_column_set_fixed_width (column, 200);
gtk_tree_view_append_column (GTK_TREE_VIEW (event_viewer_data->TreeView), column);
+ renderer = gtk_cell_renderer_text_new ();
+ column = gtk_tree_view_column_new_with_attributes ("Frequency standard deviation (Hz)",
+ renderer,
+ "text", FREQUENCY_STANDARD_DEV_COLUMN,
+ NULL);
+ gtk_tree_view_column_set_alignment (column, 0.0);
+ gtk_tree_view_column_set_fixed_width (column, 200);
+ gtk_tree_view_append_column (GTK_TREE_VIEW (event_viewer_data->TreeView), column);
event_viewer_data->SelectionTree = gtk_tree_view_get_selection (GTK_TREE_VIEW (event_viewer_data->TreeView));
interrupt_update_time_window,
event_viewer_data);
+ g_object_set_data_full(G_OBJECT(event_viewer_data->Hbox),
+ "event_data",
+ event_viewer_data,
+ (GDestroyNotify) InterruptFree);
FirstRequest(event_viewer_data );
return event_viewer_data;
}
+
/**
*
* For each trace in the traceset, this function:
/**
* This function gets the id of the interrupt. The id is stored in a dynamic structure.
- * Refer to the print.c file for howto extract data from a dynamic structure.
+ * Refer to the print.c file for how to extract data from a dynamic structure.
*/
static guint64 get_interrupt_id(LttEvent *e)
{
element = &g_array_index(FirstRequestIrqEntry,irq_entry,i);
if(element->cpu_id == cpu_id)
{
- TotalDurationMaxIrqDuration(element,time_exit, FirstRequestIrqExit);
+ CalculateTotalDurationAndMaxIrqDurationAndMinIrqDuration(element,time_exit, FirstRequestIrqExit);
g_array_remove_index(FirstRequestIrqEntry, i);
break;
}
}
}
+
/**
- * This function calculates the total duration of an interrupt and the longest Irq handler.
+ * This function calculates the total duration of an interrupt and the longest & shortest Irq handlers.
*
*/
-static void TotalDurationMaxIrqDuration(irq_entry *e, LttTime time_exit, GArray *FirstRequestIrqExit){
+static void CalculateTotalDurationAndMaxIrqDurationAndMinIrqDuration(irq_entry *e, LttTime time_exit, GArray *FirstRequestIrqExit)
+{
Irq irq;
Irq *element;
guint i;
{
irq.cpu_id = e->cpu_id;
irq.id = e->id;
- irq.frequency++;
+ irq.NumerofInterruptions++;
irq.total_duration = ltt_time_sub(time_exit, e->event_time);
irq.max_irq_handler.start_time = e->event_time;
irq.max_irq_handler.end_time = time_exit;
irq.max_irq_handler.duration = ltt_time_sub(time_exit, e->event_time);
-
+ irq.min_irq_handler.start_time = e->event_time;
+ irq.min_irq_handler.end_time = time_exit;
+ irq.min_irq_handler.duration = ltt_time_sub(time_exit, e->event_time);
+
g_array_append_val (FirstRequestIrqExit, irq);
}
else
notFound = TRUE;
duration = ltt_time_sub(time_exit, e->event_time);
element->total_duration = ltt_time_add(element->total_duration, duration);
- element->frequency++;
+ element->NumerofInterruptions++;
+ // Max irq handler
if(ltt_time_compare(duration,element->max_irq_handler.duration) > 0)
{
element->max_irq_handler.duration = duration;
element->max_irq_handler.start_time = e->event_time;
element->max_irq_handler.end_time = time_exit;
}
+ // Min irq handler
+ if(ltt_time_compare(duration,element->min_irq_handler.duration) < 0)
+ {
+ element->min_irq_handler.duration = duration;
+ element->min_irq_handler.start_time = e->event_time;
+ element->min_irq_handler.end_time = time_exit;
+ }
}
}
if(!notFound)
{
irq.cpu_id = e->cpu_id;
irq.id = e->id;
- irq.frequency++;
+ irq.NumerofInterruptions++;
irq.total_duration = ltt_time_sub(time_exit, e->event_time);
-
+ // Max irq handler
irq.max_irq_handler.start_time = e->event_time;
irq.max_irq_handler.end_time = time_exit;
irq.max_irq_handler.duration = ltt_time_sub(time_exit, e->event_time);
+ // Min irq handler
+ irq.min_irq_handler.start_time = e->event_time;
+ irq.min_irq_handler.end_time = time_exit;
+ irq.min_irq_handler.duration = ltt_time_sub(time_exit, e->event_time);
g_array_append_val (FirstRequestIrqExit, irq);
}
return FALSE;
}
+/**
+ * This function calculates the average duration for each Irq Id
+ *
+ */
static void CalculateAverageDurationForEachIrqId(InterruptEventData *event_data)
{
guint64 real_data;
real_data = element->total_duration.tv_sec;
real_data *= NANOSECONDS_PER_SECOND;
real_data += element->total_duration.tv_nsec;
- element->average_duration = real_data / element->frequency;
+ if(element->NumerofInterruptions != 0)
+ element->average_duration = real_data / element->NumerofInterruptions;
+ else
+ element->average_duration = 0;
}
}
static void SumItems(gint irq_id, LttTime Xi, InterruptEventData *event_data)
{
gint i;
- guint time_in_ns;
+ guint Xi_in_ns;
- gint temp;
+ gint duration_inner_part;
+ guint64 period_inner_part;
+ guint64 frequency_inner_part;
+
Irq *average;
SumId *sumItem;
SumId sum;
+ int FrequencyHZ = 0;
gboolean notFound = FALSE;
GArray *FirstRequestIrqExit = event_data->FirstRequestIrqExit;
GArray *SumArray = event_data->SumArray;
- time_in_ns = Xi.tv_sec;
- time_in_ns *= NANOSECONDS_PER_SECOND;
- time_in_ns += Xi.tv_nsec;
+ Xi_in_ns = Xi.tv_sec;
+ Xi_in_ns *= NANOSECONDS_PER_SECOND;
+ Xi_in_ns += Xi.tv_nsec;
for(i = 0; i < FirstRequestIrqExit->len; i++)
{
average = &g_array_index(FirstRequestIrqExit,Irq,i);
if(irq_id == average->id)
{
- temp = time_in_ns - average->average_duration;
- sum.sumOfDurations = pow (temp , 2);
- //printf("one : %d\n", sum.sumOfDurations);
+ duration_inner_part = Xi_in_ns - average->average_duration;
+ FrequencyHZ = FrequencyInHZ(average->NumerofInterruptions, event_data->time_window);
sum.irqId = irq_id;
- sum.frequency = average->frequency;
+ // compute (xi -Xa)^2 of the duration Standard deviation
+ sum.NumerofInterruptions = average->NumerofInterruptions;
+ sum.sumOfDurations = pow (duration_inner_part , 2);
+
+ // compute (xi -Xa)^2 of the period Standard deviation
+ period_inner_part = CalculatePeriodInnerPart(Xi_in_ns, FrequencyHZ);
+
+ // compute (xi -Xa)^2 of the frequency Standard deviation
+ frequency_inner_part = CalculateFrequencyInnerPart(Xi_in_ns, FrequencyHZ);
+
+ sum.sumOfPeriods = period_inner_part;
+
+ sum.sumOfFrequencies = frequency_inner_part;
+
if(event_data->SumArray->len == NO_ITEMS)
{
g_array_append_val (SumArray, sum);
}
else
- {
-
+ {
for(i = 0; i < SumArray->len; i++)
{
sumItem = &g_array_index(SumArray, SumId, i);
{
notFound = TRUE;
sumItem->sumOfDurations += sum.sumOfDurations;
-
- }
+ sumItem->sumOfPeriods += sum.sumOfPeriods;
+ sumItem->sumOfFrequencies += sum.sumOfFrequencies;
+ }
}
if(!notFound)
{
g_array_append_val (SumArray, sum);
}
-
-
+
}
}
}
}
+/**
+ * This function computes the inner part of the period standard deviation = sqrt(1/N Sum ((xi -Xa)^2))
+ * The inner part is: (xi -Xa)^2
+ */
+static guint64 CalculatePeriodInnerPart(guint Xi, guint FrequencyHZ)
+{
+
+ double periodInSec; /*period in sec*/
+ int periodInNSec;
+ gint difference;
+ guint64 result;
+ periodInSec = (double)1/FrequencyHZ;
+ periodInSec *= NANOSECONDS_PER_SECOND;
+ periodInNSec = (int)periodInSec;
+
+ difference = Xi - periodInNSec;
+ result = pow (difference , 2);
+ return result;
+}
+
+/**
+ * This function computes the inner part of the frequency standard deviation = sqrt(1/N Sum ((xi -Xa)^2))
+ * The inner part is: (xi -Xa)^2
+ */
+static guint64 CalculateFrequencyInnerPart(guint Xi_in_ns, guint FrequencyHZ)
+{
+ guint64 result;
+ gint difference;
+
+ difference = Xi_in_ns - FrequencyHZ;
+ result = pow (difference , 2);
+ return result;
+}
/**
* This function displays the result on the viewer
*
GtkTreeIter iter;
guint64 real_data;
guint maxIRQduration;
- double period;
+ double periodInSec;
+ int periodInNsec;
char maxIrqHandler[80];
InterruptEventData *event_data = (InterruptEventData *)hook_data;
GArray *FirstRequestIrqExit = event_data->FirstRequestIrqExit;
int FrequencyHZ = 0;
- period = 0.0;
+ periodInSec = 0;
gtk_list_store_clear(event_data->ListStore);
for(i = 0; i < FirstRequestIrqExit->len; i++)
{
maxIRQduration *= NANOSECONDS_PER_SECOND;
maxIRQduration += element.max_irq_handler.duration.tv_nsec;
- sprintf(maxIrqHandler, "%d [%d.%d - %d.%d]",maxIRQduration, element.max_irq_handler.start_time.tv_sec, \
- element.max_irq_handler.start_time.tv_nsec, element.max_irq_handler.end_time.tv_sec, \
- element.max_irq_handler.end_time.tv_nsec) ;
- FrequencyHZ = FrequencyInHZ(element.frequency,event_data->time_window);
+ sprintf(maxIrqHandler, "%d [%d.%d - %d.%d]",maxIRQduration, element.max_irq_handler.start_time.tv_sec, \
+ element.max_irq_handler.start_time.tv_nsec, element.max_irq_handler.end_time.tv_sec, \
+ element.max_irq_handler.end_time.tv_nsec) ;
+ FrequencyHZ = FrequencyInHZ(element.NumerofInterruptions,event_data->time_window);
+
if(FrequencyHZ != 0)
{
- period =(double)1/FrequencyHZ;
+ periodInSec = (double)1/FrequencyHZ;
+ periodInSec *= NANOSECONDS_PER_SECOND;
+ periodInNsec = (int)periodInSec;
+
}
-
+
gtk_list_store_append (event_data->ListStore, &iter);
gtk_list_store_set (event_data->ListStore, &iter,
CPUID_COLUMN, element.cpu_id,
IRQ_ID_COLUMN, element.id,
FREQUENCY_COLUMN, FrequencyHZ,
DURATION_COLUMN, real_data,
- DURATION_STANDARD_DEV_COLUMN, CalculateStandardDeviation(element.id, event_data),
+ DURATION_STANDARD_DEV_COLUMN, CalculateDurationStandardDeviation(element.id, event_data),
MAX_IRQ_HANDLER_COLUMN, maxIrqHandler,
- AVERAGE_PERIOD , period,
+ AVERAGE_PERIOD , periodInNsec,
+ PERIOD_STANDARD_DEV_COLUMN, CalculatePeriodStandardDeviation(element.id, event_data),
+ FREQUENCY_STANDARD_DEV_COLUMN, CalculateFrequencyStandardDeviation(element.id, event_data),
-1);
+
+ printf("%d %d %lld %d %s %d %d %d\n\n",element.id, FrequencyHZ,real_data,CalculateDurationStandardDeviation(element.id, event_data), maxIrqHandler, periodInNsec, CalculatePeriodStandardDeviation(element.id, event_data), CalculateFrequencyStandardDeviation(element.id, event_data));
}
}
-static int FrequencyInHZ(gint frequency, TimeWindow time_window)
+/**
+ * This function converts the number of interrupts over a time window to
+ * frequency in HZ
+ */
+static int FrequencyInHZ(gint NumerofInterruptions, TimeWindow time_window)
{
guint64 frequencyHz = 0;
- double timeSec;
-
- timeSec = ltt_time_to_double(time_window.time_width);
- double result = (timeSec/NANOSECONDS_PER_SECOND);
- frequencyHz = frequency / result;
+ double timeSec; // time in second
+ double result;
+ result = ltt_time_to_double(time_window.time_width);
+ timeSec = (result/NANOSECONDS_PER_SECOND); //time in second
+ frequencyHz = NumerofInterruptions / timeSec;
return frequencyHz;
}
/**
- * This function calculatees the standard deviation
- *
+ * This function calculates the duration standard deviation
+ * Duration standard deviation = sqrt(1/N Sum ((xi -Xa)^2))
+ * Where:
+ * sumId.sumOfDurations -> Sum ((xi -Xa)^2)
+ * inner_component -> 1/N Sum ((xi -Xa)^2)
+ * deviation-> sqrt(1/N Sum ((xi -Xa)^2))
*/
-static int CalculateStandardDeviation(gint id, InterruptEventData *event_data)
+static int CalculateDurationStandardDeviation(gint id, InterruptEventData *event_data)
{
int i;
SumId sumId;
sumId = g_array_index(event_data->SumArray, SumId, i);
if(id == sumId.irqId)
{
- inner_component = sumId.sumOfDurations/ sumId.frequency;
+ if(sumId.NumerofInterruptions != 0)
+ inner_component = sumId.sumOfDurations/ sumId.NumerofInterruptions;
+ else
+ inner_component = 0.0;
deviation = sqrt(inner_component);
return deviation;
}
}
return deviation;
}
+
+
+/**
+ * This function calculates the period standard deviation
+ * Period standard deviation = sqrt(1/N Sum ((xi -Xa)^2))
+ * Where:
+ * sumId.sumOfPeriods -> Sum ((xi -Xa)^2)
+ * inner_component -> 1/N Sum ((xi -Xa)^2)
+ * period_standard_deviation-> sqrt(1/N Sum ((xi -Xa)^2))
+
+ *
+ */
+static int CalculatePeriodStandardDeviation(gint id, InterruptEventData *event_data)
+{
+ int i;
+ SumId sumId;
+ guint64 inner_component;
+ guint64 period_standard_deviation = 0;
+
+ for(i = 0; i < event_data->SumArray->len; i++)
+ {
+ sumId = g_array_index(event_data->SumArray, SumId, i);
+ if(id == sumId.irqId)
+ {
+ if(sumId.NumerofInterruptions != 0)
+ inner_component = sumId.sumOfPeriods / sumId.NumerofInterruptions;
+ else
+ inner_component = 0;
+
+ period_standard_deviation = sqrt(inner_component);
+ }
+ }
+
+ return period_standard_deviation;
+}
+
+/**
+ * This function calculates the frequency standard deviation
+ * Frequency standard deviation = sqrt(1/N Sum ((xi -Xa)^2))
+ * Where:
+ * sumId.sumOfFrequencies -> Sum ((xi -Xa)^2)
+ * inner_component -> 1/N Sum ((xi -Xa)^2)
+ * frequency_standard_deviation-> sqrt(1/N Sum ((xi -Xa)^2))
+ *
+ */
+static int CalculateFrequencyStandardDeviation(gint id, InterruptEventData *event_data)
+{
+ int i;
+ SumId sumId;
+ guint64 inner_component;
+ guint64 frequency_standard_deviation = 0;
+ for(i = 0; i < event_data->SumArray->len; i++)
+ {
+ sumId = g_array_index(event_data->SumArray, SumId, i);
+ if(id == sumId.irqId)
+ {
+ if(sumId.NumerofInterruptions != 0)
+ inner_component = sumId.sumOfFrequencies / sumId.NumerofInterruptions;
+ else
+ inner_component = 0;
+
+ frequency_standard_deviation = sqrt(inner_component);
+ }
+ }
+ return frequency_standard_deviation;
+}
+
/*
* This function is called by the main window
* when the time interval needs to be updated.
void interrupt_destroy_walk(gpointer data, gpointer user_data)
{
g_info("interrupt_destroy_walk");
+ InterruptEventData *event_data = (InterruptEventData*) data;
interrupt_destructor((InterruptEventData*)data);
-
}
+
void interrupt_destructor(InterruptEventData *event_viewer_data)
{
/* May already been done by GTK window closing */
}
}
+/**
+ This function is called when the viewer is destroyed to free hooks and memory
+*/
+static void InterruptFree(InterruptEventData *event_viewer_data)
+{
+ Tab *tab = event_viewer_data->tab;
+ if(tab != NULL)
+ {
+
+ g_array_free(event_viewer_data->FirstRequestIrqExit, TRUE);
+ g_array_free(event_viewer_data->FirstRequestIrqEntry, TRUE);
+ g_array_free(event_viewer_data->SecondRequestIrqEntry, TRUE);
+ g_array_free(event_viewer_data->SecondRequestIrqExit, TRUE);
+ g_array_free(event_viewer_data->SumArray, TRUE);
+
+ lttvwindow_unregister_time_window_notify(tab, interrupt_update_time_window, event_viewer_data);
+
+ lttvwindow_events_request_remove_all(event_viewer_data->tab,
+ event_viewer_data);
+
+ interrupt_data_list = g_slist_remove(interrupt_data_list, event_viewer_data);
+
+ }
+
+}
+
/**
* plugin's destroy function
*
projects / lttv.git / blobdiff