X-Git-Url: http://git.lttng.org/?a=blobdiff_plain;f=lttv%2Flttv%2Fsync%2Fevent_analysis_chull.c;h=d0dd0a70c88839d5b4564966d982358859e8a788;hb=2f07659468ff3427bb3923cdc5a7ec61b6e8bbce;hp=0bea74a879f325d34f9b2d5dc2cabc33e7488d80;hpb=76be6fc24daf61767bf7f0c2e64f4691fbb56c63;p=lttv.git diff --git a/lttv/lttv/sync/event_analysis_chull.c b/lttv/lttv/sync/event_analysis_chull.c index 0bea74a8..d0dd0a70 100644 --- a/lttv/lttv/sync/event_analysis_chull.c +++ b/lttv/lttv/sync/event_analysis_chull.c @@ -26,6 +26,7 @@ #include #include #include +#include #include #include "sync_chain.h" @@ -33,11 +34,6 @@ #include "event_analysis_chull.h" -#ifndef g_info -#define g_info(format...) g_log (G_LOG_DOMAIN, G_LOG_LEVEL_INFO, format) -#endif - - typedef enum { LOWER, @@ -60,8 +56,8 @@ static void analyzeMessageCHull(SyncState* const syncState, Message* const message); static GArray* finalizeAnalysisCHull(SyncState* const syncState); static void printAnalysisStatsCHull(SyncState* const syncState); -static void writeAnalysisGraphsPlotsCHull(FILE* stream, SyncState* const - syncState, const unsigned int i, const unsigned int j); +static void writeAnalysisGraphsPlotsCHull(SyncState* const syncState, const + unsigned int i, const unsigned int j); // Functions specific to this module static void registerAnalysisCHull() __attribute__((constructor (101))); @@ -77,10 +73,8 @@ static int jointCmp(const Point* const p1, const Point* const p2, const Point* const p3) __attribute__((pure)); static double crossProductK(const Point const* p1, const Point const* p2, const Point const* p3, const Point const* p4) __attribute__((pure)); -static FactorsCHull** calculateAllFactors(SyncState* const syncState); static Factors* calculateFactorsExact(GQueue* const cu, GQueue* const cl, const LineType lineType) __attribute__((pure)); -static void calculateFactorsMiddle(FactorsCHull* factors); static void calculateFactorsFallback(GQueue* const cr, GQueue* const cs, FactorsCHull* const result); static double slope(const Point* const p1, const Point* const p2) @@ -89,8 +83,6 @@ static double intercept(const Point* const p1, const Point* const p2) __attribute__((pure)); static GArray* reduceFactors(SyncState* const syncState, FactorsCHull** allFactors); -static void freeAllFactors(const SyncState* const syncState, FactorsCHull** - const allFactors); static double verticalDistance(Point* p1, Point* p2, Point* const point) __attribute__((pure)); static void floydWarshall(SyncState* const syncState, FactorsCHull** const @@ -108,12 +100,20 @@ static AnalysisModule analysisModuleCHull= { .initAnalysis= &initAnalysisCHull, .destroyAnalysis= &destroyAnalysisCHull, .analyzeMessage= &analyzeMessageCHull, - .analyzeExchange= NULL, - .analyzeBroadcast= NULL, .finalizeAnalysis= &finalizeAnalysisCHull, .printAnalysisStats= &printAnalysisStatsCHull, - .writeAnalysisGraphsPlots= &writeAnalysisGraphsPlotsCHull, - .writeAnalysisGraphsOptions= NULL, + .graphFunctions= { + .writeTraceTraceForePlots= &writeAnalysisGraphsPlotsCHull, + } +}; + +const char* const approxNames[]= { + [EXACT]= "Exact", + [MIDDLE]= "Middle", + [FALLBACK]= "Fallback", + [INCOMPLETE]= "Incomplete", + [ABSENT]= "Absent", + [SCREWED]= "Screwed", }; @@ -167,7 +167,7 @@ static void initAnalysisCHull(SyncState* const syncState) analysisData->stats->allFactors= NULL; } - if (syncState->graphs) + if (syncState->graphsStream) { analysisData->graphsData= malloc(sizeof(AnalysisGraphsDataCHull)); openGraphFiles(syncState); @@ -193,7 +193,7 @@ static void openGraphFiles(SyncState* const syncState) analysisData= (AnalysisDataCHull*) syncState->analysisData; - cwd= changeToGraphDir(syncState->graphs); + cwd= changeToGraphsDir(syncState->graphsDir); analysisData->graphsData->hullPoints= malloc(syncState->traceNb * sizeof(FILE**)); @@ -351,13 +351,13 @@ static void destroyAnalysisCHull(SyncState* const syncState) { if (analysisData->stats->allFactors != NULL) { - freeAllFactors(syncState, analysisData->stats->allFactors); + freeAllFactors(syncState->traceNb, analysisData->stats->allFactors); } free(analysisData->stats); } - if (syncState->graphs) + if (syncState->graphsStream) { if (analysisData->graphsData->hullPoints != NULL) { @@ -366,7 +366,7 @@ static void destroyAnalysisCHull(SyncState* const syncState) if (!syncState->stats && analysisData->graphsData->allFactors != NULL) { - freeAllFactors(syncState, analysisData->graphsData->allFactors); + freeAllFactors(syncState->traceNb, analysisData->graphsData->allFactors); } free(analysisData->graphsData); @@ -510,7 +510,7 @@ static GArray* finalizeAnalysisCHull(SyncState* const syncState) analysisData= (AnalysisDataCHull*) syncState->analysisData; - if (syncState->graphs) + if (syncState->graphsStream && analysisData->graphsData->hullPoints != NULL) { writeGraphFiles(syncState); closeGraphFiles(syncState); @@ -520,21 +520,21 @@ static GArray* finalizeAnalysisCHull(SyncState* const syncState) factors= reduceFactors(syncState, allFactors); - if (syncState->stats || syncState->graphs) + if (syncState->stats || syncState->graphsStream) { if (syncState->stats) { analysisData->stats->allFactors= allFactors; } - if (syncState->graphs) + if (syncState->graphsStream) { analysisData->graphsData->allFactors= allFactors; } } else { - freeAllFactors(syncState, allFactors); + freeAllFactors(syncState->traceNb, allFactors); } return factors; @@ -585,18 +585,19 @@ static void printAnalysisStatsCHull(SyncState* const syncState) FactorsCHull* factorsCHull; factorsCHull= &analysisData->stats->allFactors[j][i]; - printf("\t\t%3d - %-3d: ", i, j); + printf("\t\t%3d - %-3d: %s", i, j, + approxNames[factorsCHull->type]); if (factorsCHull->type == EXACT) { - printf("Exact a0= % 7g a1= 1 %c %7g\n", + printf(" a0= % 7g a1= 1 %c %7g\n", factorsCHull->approx->offset, factorsCHull->approx->drift < 0. ? '-' : '+', fabs(factorsCHull->approx->drift)); } else if (factorsCHull->type == MIDDLE) { - printf("Middle a0= % 7g a1= 1 %c %7g accuracy %7g\n", + printf(" a0= % 7g a1= 1 %c %7g accuracy %7g\n", factorsCHull->approx->offset, factorsCHull->approx->drift - 1. < 0. ? '-' : '+', fabs(factorsCHull->approx->drift - 1.), factorsCHull->accuracy); @@ -609,14 +610,14 @@ static void printAnalysisStatsCHull(SyncState* const syncState) } else if (factorsCHull->type == FALLBACK) { - printf("Fallback a0= % 7g a1= 1 %c %7g error= %7g\n", + printf(" a0= % 7g a1= 1 %c %7g error= %7g\n", factorsCHull->approx->offset, factorsCHull->approx->drift - 1. < 0. ? '-' : '+', fabs(factorsCHull->approx->drift - 1.), factorsCHull->accuracy); } else if (factorsCHull->type == INCOMPLETE) { - printf("Incomplete\n"); + printf("\n"); if (factorsCHull->min->drift != -INFINITY) { @@ -635,7 +636,7 @@ static void printAnalysisStatsCHull(SyncState* const syncState) } else if (factorsCHull->type == SCREWED) { - printf("Screwed\n"); + printf("\n"); if (factorsCHull->min != NULL && factorsCHull->min->drift != -INFINITY) { @@ -654,7 +655,7 @@ static void printAnalysisStatsCHull(SyncState* const syncState) } else if (factorsCHull->type == ABSENT) { - printf("Absent\n"); + printf("\n"); } else { @@ -742,44 +743,19 @@ static double crossProductK(const Point const* p1, const Point const* p2, * Free a container of FactorsCHull * * Args: - * syncState: container for synchronization data. - * allFactors: container of Factors + * traceNb: number of traces + * allFactors: container of FactorsCHull */ -static void freeAllFactors(const SyncState* const syncState, FactorsCHull** - const allFactors) +void freeAllFactors(const unsigned int traceNb, FactorsCHull** const + allFactors) { unsigned int i, j; - for (i= 0; i < syncState->traceNb; i++) + for (i= 0; i < traceNb; i++) { for (j= 0; j <= i; j++) { - FactorsCHull* factorsCHull; - - factorsCHull= &allFactors[i][j]; - if (factorsCHull->type == MIDDLE || factorsCHull->type == - INCOMPLETE || factorsCHull->type == ABSENT) - { - free(factorsCHull->min); - free(factorsCHull->max); - } - else if (factorsCHull->type == SCREWED) - { - if (factorsCHull->min != NULL) - { - free(factorsCHull->min); - } - if (factorsCHull->max != NULL) - { - free(factorsCHull->max); - } - } - - if (factorsCHull->type == EXACT || factorsCHull->type == MIDDLE || - factorsCHull->type == FALLBACK) - { - free(factorsCHull->approx); - } + destroyFactorsCHull(&allFactors[i][j]); } free(allFactors[i]); } @@ -787,6 +763,40 @@ static void freeAllFactors(const SyncState* const syncState, FactorsCHull** } +/* + * Free a FactorsCHull + * + * Args: + * factorsCHull: container of Factors + */ +void destroyFactorsCHull(FactorsCHull* factorsCHull) +{ + if (factorsCHull->type == MIDDLE || factorsCHull->type == + INCOMPLETE || factorsCHull->type == ABSENT) + { + free(factorsCHull->min); + free(factorsCHull->max); + } + else if (factorsCHull->type == SCREWED) + { + if (factorsCHull->min != NULL) + { + free(factorsCHull->min); + } + if (factorsCHull->max != NULL) + { + free(factorsCHull->max); + } + } + + if (factorsCHull->type == EXACT || factorsCHull->type == MIDDLE || + factorsCHull->type == FALLBACK) + { + free(factorsCHull->approx); + } +} + + /* * Analyze the convex hulls to determine the synchronization factors between * each pair of trace. @@ -798,7 +808,7 @@ static void freeAllFactors(const SyncState* const syncState, FactorsCHull** * FactorsCHull*[TraceNum][TraceNum] array. See the documentation for the * member allFactors of AnalysisStatsCHull. */ -static FactorsCHull** calculateAllFactors(SyncState* const syncState) +FactorsCHull** calculateAllFactors(SyncState* const syncState) { unsigned int traceNumA, traceNumB; FactorsCHull** allFactors; @@ -908,7 +918,7 @@ static FactorsCHull** calculateAllFactors(SyncState* const syncState) * Args: * factors: contains the min and max limits, used to store the result */ -static void calculateFactorsMiddle(FactorsCHull* factors) +void calculateFactorsMiddle(FactorsCHull* const factors) { double amin, amax, bmin, bmax, bhat; @@ -1522,20 +1532,19 @@ static void getFactors(FactorsCHull** const allFactors, unsigned int** const * Write the analysis-specific graph lines in the gnuplot script. * * Args: - * stream: stream where to write the data * syncState: container for synchronization data * i: first trace number * j: second trace number, garanteed to be larger than i */ -void writeAnalysisGraphsPlotsCHull(FILE* stream, SyncState* const syncState, - const unsigned int i, const unsigned int j) +void writeAnalysisGraphsPlotsCHull(SyncState* const syncState, const unsigned + int i, const unsigned int j) { AnalysisDataCHull* analysisData; FactorsCHull* factorsCHull; analysisData= (AnalysisDataCHull*) syncState->analysisData; - fprintf(stream, + fprintf(syncState->graphsStream, "\t\"analysis_chull-%1$03d_to_%2$03d.data\" " "title \"Lower half-hull\" with linespoints " "linecolor rgb \"#015a01\" linetype 4 pointtype 8 pointsize 0.8, \\\n" @@ -1547,7 +1556,7 @@ void writeAnalysisGraphsPlotsCHull(FILE* stream, SyncState* const syncState, factorsCHull= &analysisData->graphsData->allFactors[j][i]; if (factorsCHull->type == EXACT) { - fprintf(stream, + fprintf(syncState->graphsStream, "\t%7g + %7g * x " "title \"Exact conversion\" with lines " "linecolor rgb \"black\" linetype 1, \\\n", @@ -1555,25 +1564,25 @@ void writeAnalysisGraphsPlotsCHull(FILE* stream, SyncState* const syncState, } else if (factorsCHull->type == MIDDLE) { - fprintf(stream, + fprintf(syncState->graphsStream, "\t%.2f + %.10f * x " "title \"Min conversion\" with lines " "linecolor rgb \"black\" linetype 5, \\\n", factorsCHull->min->offset, factorsCHull->min->drift); - fprintf(stream, + fprintf(syncState->graphsStream, "\t%.2f + %.10f * x " "title \"Max conversion\" with lines " "linecolor rgb \"black\" linetype 8, \\\n", factorsCHull->max->offset, factorsCHull->max->drift); - fprintf(stream, + fprintf(syncState->graphsStream, "\t%.2f + %.10f * x " "title \"Middle conversion\" with lines " - "linecolor rgb \"gray60\" linetype 1, \\\n", + "linecolor rgb \"black\" linetype 1, \\\n", factorsCHull->approx->offset, factorsCHull->approx->drift); } else if (factorsCHull->type == FALLBACK) { - fprintf(stream, + fprintf(syncState->graphsStream, "\t%.2f + %.10f * x " "title \"Fallback conversion\" with lines " "linecolor rgb \"gray60\" linetype 1, \\\n", @@ -1583,7 +1592,7 @@ void writeAnalysisGraphsPlotsCHull(FILE* stream, SyncState* const syncState, { if (factorsCHull->min->drift != -INFINITY) { - fprintf(stream, + fprintf(syncState->graphsStream, "\t%.2f + %.10f * x " "title \"Min conversion\" with lines " "linecolor rgb \"black\" linetype 5, \\\n", @@ -1592,7 +1601,7 @@ void writeAnalysisGraphsPlotsCHull(FILE* stream, SyncState* const syncState, if (factorsCHull->max->drift != INFINITY) { - fprintf(stream, + fprintf(syncState->graphsStream, "\t%.2f + %.10f * x " "title \"Max conversion\" with lines " "linecolor rgb \"black\" linetype 8, \\\n", @@ -1603,7 +1612,7 @@ void writeAnalysisGraphsPlotsCHull(FILE* stream, SyncState* const syncState, { if (factorsCHull->min != NULL && factorsCHull->min->drift != -INFINITY) { - fprintf(stream, + fprintf(syncState->graphsStream, "\t%.2f + %.10f * x " "title \"Min conversion\" with lines " "linecolor rgb \"black\" linetype 5, \\\n", @@ -1612,7 +1621,7 @@ void writeAnalysisGraphsPlotsCHull(FILE* stream, SyncState* const syncState, if (factorsCHull->max != NULL && factorsCHull->max->drift != INFINITY) { - fprintf(stream, + fprintf(syncState->graphsStream, "\t%.2f + %.10f * x " "title \"Max conversion\" with lines " "linecolor rgb \"black\" linetype 8, \\\n",