void PrintHistograms(TString fname) { TH1D *h[4]; TH1D *he[4]; double min[] = {-8.0,1.0,0.0,1400.0}; double max[] = {2.0,91.0,1.0,4000.0}; double emin[] = {0.0,0.0,0.0,0.0}; double emax[] = {0.5,10.0,0.05,100.0}; for (int i=0;i<4;i++) { stringstream ss; ss << i; TString name = "c" + TString(ss.str()); h[i] = new TH1D(name,"TAGH timewalk: "+name+";"+name+";counters",100,min[i],max[i]); name = "e_c" + TString(ss.str()); he[i] = new TH1D(name,"TAGH timewalk: "+name+";"+name+";counters"+";"+name+";counters",100,emin[i],emax[i]); } TH1D hChisq("Chisq","TAGH timewalk: chi-square / Ndf;chi-square / Ndf;counters",150,0.0,150.0); ifstream fin(fname); for (int i=0;i<274;i++) { // counters char sep;int n,ndf; double chisq,c[4],dc[4]; fin >> n >> sep >> ndf >> sep >> chisq; for (int j=0;j<4;j++) fin >> sep >> c[j]; for (int j=0;j<4;j++) fin >> sep >> dc[j]; if (chisq==0.0) continue; for (int j=0;j<4;j++) h[j]->Fill(c[j]); for (int j=0;j<4;j++) he[j]->Fill(dc[j]); hChisq.Fill(chisq/ndf); } fin.close(); for (int i=0;i<4;i++) {PrintHistogram(*h[i]);PrintHistogram(*he[i]);delete h[i];delete he[i];} PrintHistogram(hChisq); }
void PrintResults(void) { int ix; #ifdef DEBUG fprintf(ofp, "Debug dump of globals global struct:\n"); fprintf(ofp, " files_are_different = %d\n", globals.files_are_different); fprintf(ofp, " file1_size = %d\n", globals.file1_size); fprintf(ofp, " file2_size = %d\n", globals.file2_size); fprintf(ofp, " compare_size = %d\n", globals.compare_size); fprintf(ofp, " filename1 = '%s'\n", globals.filename1); fprintf(ofp, " filename2 = '%s'\n", globals.filename2); fprintf(ofp, " rng_size = %d\n", globals.rng_size); fprintf(ofp, " rng_numel = %d\n", globals.rng_numel); fprintf(ofp, "Debug dump of range array:\n"); for (ix = 0; ix < globals.rng_numel; ix++) { fprintf(ofp, " %8d %8d\n", globals.rng[ix].start, globals.rng[ix].end); } #endif if (histogram) { PrintHistogram(); } else { PrintHeader(); PrintSeparator(); for (ix = 0; ix < globals.rng_numel; ix++) { PrintRange(globals.rng[ix].start, globals.rng[ix].end); } } }
int main(int argc, char **argv) { long *Histogram ; int P, q ; printf("Test PrintHistogram: "); P = 68 ; /* number of processors >= 1 */ Histogram = (long *) malloc((P+1)* sizeof(long)) ; if ( Histogram == NULL ){ printf("Error\n") ; exit(1); } /* Initialise histogram*/ for (q=0; q<P/2; q++) Histogram[q] = q ; for (q=P/2; q<=P; q++) Histogram[q] = 0 ; printf("\n"); printf("*** Histogram for processor q should be q\n"); PrintHistogram(0, P, Histogram) ; printf("OK\n") ; exit(0); } /* end main */
int main(int argc, char *argv[]) { char *indexFileName=NULL; char *fastaFileName=NULL; int numThreads = 1; int whichStrand = 0; int space = NTSpace; int c; RGBinary rg; RGIndex index; while((c = getopt(argc, argv, "f:i:n:s:A:h")) >= 0) { switch(c) { case 'f': fastaFileName=strdup(optarg); break; case 'h': return PrintUsage(); case 'i': indexFileName=strdup(optarg); break; case 's': whichStrand=atoi(optarg); break; case 'n': numThreads=atoi(optarg); break; case 'A': space=atoi(optarg); break; default: fprintf(stderr, "Unrecognized option: -%c\n", c); return 1; } } if(1 == argc || argc != optind) { return PrintUsage(); } if(NULL == indexFileName) { PrintError(Name, "indexFileName", "Command line option", Exit, InputArguments); } if(NULL == fastaFileName) { PrintError(Name, "fastaFileName", "Command line option", Exit, InputArguments); } assert(whichStrand == BothStrands || whichStrand == ForwardStrand || whichStrand == ReverseStrand); /* Read in the rg binary file */ RGBinaryReadBinary(&rg, space, fastaFileName); /* Read the index */ RGIndexRead(&index, indexFileName); assert(index.space == rg.space); fprintf(stderr, "%s", BREAK_LINE); PrintHistogram(&index, &rg, whichStrand, numThreads); fprintf(stderr, "%s", BREAK_LINE); fprintf(stderr, "%s", BREAK_LINE); fprintf(stderr, "Cleaning up.\n"); /* Delete the index */ RGIndexDelete(&index); /* Delete the rg */ RGBinaryDelete(&rg); fprintf(stderr, "%s", BREAK_LINE); fprintf(stderr, "Terminating successfully!\n"); fprintf(stderr, "%s", BREAK_LINE); return 0; }
static bool PrintRvmResEvent() { int sum=0; SmonViceId Vice; unsigned long Time; unsigned long VolID; FileResStats FileRes; DirResStats DirRes; long lshsize; HistoElem *LogSizeHisto; long lmhsize; HistoElem *LogMaxHisto; ResConflictStats Conflicts; long shhsize; HistoElem *SuccHierHist; long fhhsize; HistoElem *FailHierHist; ResLogStats ResLog; long vlhsize; HistoElem *VarLogHisto; long lssize; HistoElem *LogSize; sum = ReadRvmResEvent(&Vice, &Time, &VolID, &FileRes, &DirRes, &lshsize, &LogSizeHisto, &lmhsize, &LogMaxHisto, &Conflicts, &shhsize, &SuccHierHist, &fhhsize, &FailHierHist, &ResLog, &vlhsize, &VarLogHisto, &lssize, &LogSize); if (sum != 0) LogMsg(0,LogLevel,LogFile,"**** BAD RESOLUTION RECORD ****"); printf("%-12s: ","SRVRVMRES"); PrintViceId(&Vice); PrintDecimal(Time); PrintHex(VolID); printf("\nGeneral:\n"); PrintDecimal(ResLog.NumWraps); PrintDecimal(ResLog.NumAdmGrows); PrintDecimal(ResLog.NumVAllocs); PrintDecimal(ResLog.NumVFrees); PrintDecimal(ResLog.Highest); printf("\nFile:\n"); PrintDecimal(FileRes.Resolves); PrintDecimal(FileRes.NumSucc); PrintDecimal(FileRes.NumConf); PrintDecimal(FileRes.RuntForce); PrintDecimal(FileRes.WeakEq); printf("\n"); PrintDecimal(FileRes.NumReg); PrintDecimal(FileRes.UsrResolver); PrintDecimal(FileRes.SuccUsrResolver); PrintDecimal(FileRes.PartialVSG); printf("\nDir:\n"); PrintDecimal(DirRes.Resolves); PrintDecimal(DirRes.NumSucc); PrintDecimal(DirRes.NumConf); PrintDecimal(DirRes.NumNoWork); PrintDecimal(DirRes.Problems); PrintDecimal(DirRes.PartialVSG); printf("\n"); PrintHistogram(lshsize,LogSizeHisto); PrintHistogram(lmhsize,LogMaxHisto); printf("Conflicts:\n"); PrintDecimal(Conflicts.NameName); PrintDecimal(Conflicts.RemoveUpdate); PrintDecimal(Conflicts.UpdateUpdate); PrintDecimal(Conflicts.Rename); PrintDecimal(Conflicts.LogWrap); PrintDecimal(Conflicts.Other); printf("\n"); PrintHistogram(shhsize,SuccHierHist); PrintHistogram(fhhsize,FailHierHist); PrintHistogram(vlhsize,VarLogHisto); PrintHistogram(lssize,LogSize); if (sum) return mtrue; else return mfalse; }