VMS/C_Libraries/Histogram: Histogram.c comparison

comparison Histogram.c @ 10:7a39408f9ea3

Version 0

author	Merten Sach <msach@mailbox.tu-berlin.de>
date	Mon, 11 Jul 2011 18:10:52 +0200
parents	060d63cb5d34
children	31a248920429

comparison

equal deleted inserted replaced

-:d68258caa53d
+:e15db7e192f0
 * Author: seanhalle@yahoo.com
 *
 */
 #include <stdio.h>
 #include "Histogram.h"
+#include "../VMS.h"
 #include "../vutilities.h"
 #include "../vmalloc.h"
+//External variables for saving of the histogram
+//These have to be defined by to plugins in order to enable VMS to print this
+//information to the histogram file
+extern char __ProgrammName[];   //Defined in main.c
+extern char __Scheduler[];      //Defined in VPThread_PluginFns.c
+extern char __DataSet[255];
 /*This Histogram Abstract Data Type has a number of bins plus a range of
 * values that the bins span, both chosen at creation.
 *
 *One creates a Histogram instance using the makeHistogram function, then
 /*Inline because use with RDTSC in innermost code so need ultra-fast
 */
 void inline
-addIntervalToHist( int32 startIntvl, int32 endIntvl, Histogram *hist )
+addIntervalToHist( uint32 startIntvl, uint32 endIntvl, Histogram *hist )
 {
 int32 value;
 value = endIntvl - startIntvl;
 if( value < 0 || value > 10000000 ) return; //sanity check
 value = endIntvl - startIntvl;
 if( value < 0 || value > 10000000 ) return; //sanity check
 subFromHist( value, hist );
 }
+void
+saveHistToFile(Histogram *hist)
+{
+FILE *output;
+int32 binIdx, binStart, binEnd, centerValue, width;
+int32 maxHeight, i,n;
+float32 total, total2, binPercent, expectedValue1, expectedValue2;
+if(hist == NULL || hist->name == NULL)
+return;
+//Calculate the average
+//do all except the top bin
+maxHeight = 0; total = 0.0; expectedValue1 = 0.0;
+for( i = 0; i < hist->numBins -1; i++ )
+{
+if( maxHeight < hist->bins[ i ] ) maxHeight = hist->bins[ i ];
+total += hist->bins[ i ];
+binStart = hist->startOfRange + hist->binWidth * i;
+expectedValue1 += hist->bins[ i ] * (binStart + hist->binWidth/2.0);
+}
+//copy and calc expected value minus the top bin
+expectedValue2  = expectedValue1;
+expectedValue2 /= total;
+total2          = total;
+//now do last iteration, to add the top bin
+if(maxHeight < hist->bins[ i ])
+maxHeight = hist->bins[ i ];
+total += hist->bins[ i ];
+binStart = hist->startOfRange + hist->binWidth * i;
+expectedValue1 += hist->bins[ i ] * (binStart + hist->binWidth/2.0);
+expectedValue1 /= total;
+//If a histogram directory does not exist, do not save to file.
+//TODO change to argument
+char filename[255];
+for(n=0;n<255;n++)
+{
+sprintf(filename, "./histograms/%s.%d.dat", hist->name,n);
+output = fopen(filename,"r");
+if(output)
+{
+fclose(output);
+}else{
+break;
+}
+}
+printf("Saving Hist to File: %s ...\n", filename);
+output = fopen(filename,"w+");
+if(output == NULL){
+printf("[!]No histogram was saved. To save histograms create folder 'histograms'.\n");
+return;
+}
+/*
+* Write the header of the measurement file.
+*/
+//--------------------------
+//Build Environment
+fprintf(output, "# >> Build Environment <<\n");
+fprintf(output, "# Hardware Architecture: ");
+#ifdef __x86_64
+fprintf(output, "x86_64");
+#endif __
+#ifdef __i386
+fprintf(output, "x86");
+#endif
+fprintf(output, "\n");
+fprintf(output, "# GCC VERSION: %d.%d.%d\n",__GNUC__,__GNUC_MINOR__,__GNUC_PATCHLEVEL__);
+fprintf(output, "# Build Date: %s %s\n", __DATE__, __TIME__);
+fprintf(output, "# Number of Cores: %d\n", NUM_CORES);
+//--------------------------
+//--------------------------
+//VMS Plugins
+fprintf(output, "#\n# >> VMS Plugins <<\n");
+fprintf(output, "# Language : ");
+#ifdef VPTHREAD
+fprintf(output, "VPThread");
+#endif
+#ifdef VCILK
+fprintf(output, "VCilk");
+#endif
+#ifdef SSR
+fprintf(output, "SSR");
+#endif
+fprintf(output, "\n");
+fprintf(output, "# Scheduler: %s\n", __Scheduler);
+//--------------------------
+//Application
+fprintf(output, "#\n# >> Application <<\n");
+fprintf(output, "# Name: %s\n", __ProgrammName);
+fprintf(output, "# Data Set:\n%s\n",__DataSet);
+//--------------------------
+//Histogram
+fprintf(output, "#\n# Histogram Name: %s\n", hist->name);
+fprintf(output, "# Expected Values\n");
+fprintf(output, "#\tnum samples: %d | expected value: %3.2f \n",
+(int)total, expectedValue1 );
+fprintf(output, "#\tminus top bin, num samples: %d | expected value: %3.2f \n",
+(int)total2, expectedValue2 );
+fprintf(output, "#\n# [Interval] [Center Value] [Count] [relative Count] [Width]\n");
+for( binIdx = 0; binIdx < hist->numBins; binIdx++ )
+{
+binStart = hist->startOfRange + hist->binWidth * binIdx;
+binEnd = binStart + hist->binWidth - 1;
+centerValue = (binStart+binEnd)/2;
+width = (binEnd-binStart)+1;
+fprintf(output, "%d-%d\t%d\t%d\t%.4f\t%d\n", binStart, binEnd, centerValue,
+hist->bins[ binIdx ],
+hist->bins[ binIdx ]/total, width);
+}
+fclose(output);
+fflush(stdout);
+}
 void
 printHist( Histogram *hist )
 {
 int32   binIdx, i, numBars, maxHeight, barValue, binStart, binEnd;
 binStart = hist->startOfRange + hist->binWidth * i;
 expectedValue1 += hist->bins[ i ] * (binStart + hist->binWidth/2.0);
 expectedValue1 /= total;
-barValue = maxHeight / 60;  //60 spaces across page for tallest bin
 printf( "histogram: " );
 if( hist->name != NULL ) printf( "%s\n", hist->name );
 else printf( "\n" );
 printf( "               num samples: %d | expected value: %3.2f \n",
 (int)total, expectedValue1 );
 printf( "minus top bin, num samples: %d | expected value: %3.2f \n",
 (int)total2, expectedValue2 );
+if(maxHeight < 60){
+barValue = 1;
+printf("Single Bar Value: %i\n", barValue);
+}else{
+barValue = maxHeight / 60;  //60 spaces across page for tallest bin
+printf("Single Bar Value: %0.3f\n", (float)maxHeight /60);
+}
 if( barValue == 0 ) { printf("error: bar val zero\n"); return; }
 for( binIdx = 0; binIdx < hist->numBins; binIdx++ )
 {
 binStart = hist->startOfRange + hist->binWidth * binIdx;
 binEnd = binStart + hist->binWidth - 1;

Mercurial > cgi-bin > hgwebdir.cgi > VMS > C_Libraries > Histogram

comparison Histogram.c @ 10:7a39408f9ea3