// ################################################################################//// name:          statistics.cc//// author:        Martin Pelikan//// purpose:       functions that compute and print out the statistics during and//                after the run//// last modified: February 1999//// #################################################################################include <stdio.h>#include <stdlib.h>#include "statistics.h"#include "population.h"#include "fitness.h"#include "boa.h"#include "mymath.h"#include "memalloc.h"#include "graph.h"// ================================================================================//// name:          initializeBasicStatistics//// function:      initializes basic statistics stuff//// parameters:    boaParams....the parameters sent to the BOA//// returns:       (int) 0//// ================================================================================int intializeBasicStatistics(BasicStatistics *statistics, BoaParams *boaParams){  // allocate memory for univariate frequencies  statistics->p1 = (float*) Calloc(boaParams->n,sizeof(float));  // get back  return 0;};// ================================================================================//// name:          doneBasicStatistics//// function:      done method for basic statistics//// parameters:    (none)//// returns:       (int) 0//// ================================================================================int doneBasicStatistics(BasicStatistics *statistics){  // free the memory used by the univariate frequencies  Free(statistics->p1);  // get back  return 0;};// ================================================================================//// name:          computeBasicStatistics//// function:      computes some basic statistics (on fitness and so) and sets some//                variables required for printing this information out in the//                future//// parameters:    t............number of generation//                population...a current population//                boaParams....the parameters sent to the BOA//// returns:       (int) 0//// ================================================================================int computeBasicStatistics(BasicStatistics *statistics, long t, Population *population, BoaParams *boaParams){  long i;  // set some variables  statistics->N                 = population->N;  statistics->n                 = population->n;  statistics->generation        = t;  statistics->guidanceThreshold = boaParams->guidanceThreshold;  // compute the maximal, minimal and average fitness in the population  statistics->max = 0;  statistics->minF = statistics->maxF = statistics->avgF = population->f[0];  if (isBestDefined())    statistics->numOptimal = (isOptimal(population->x[0],statistics->n))? 1:0;  for (i=1; i<statistics->N; i++)    {      statistics->avgF += population->f[i];      if (population->f[i]<statistics->minF)	statistics->minF = population->f[i];      else	if (population->f[i]>statistics->maxF)	  {	    statistics->maxF = population->f[i];	    statistics->max = i;	  }      if (isBestDefined())	if (isOptimal(population->x[i],statistics->n))	  statistics->numOptimal++;    };  statistics->avgF /= statistics->N;  // allocate memory for and compute the univariate frequencies  computeUnivariateFrequencies(population,statistics->p1);  // set the best guy (if defined...)  if (isBestDefined())    statistics->bestX = population->x[statistics->max];  // get back  return 0;};// ================================================================================//// name:          generationStatistics//// function:      prints out various statistics//// parameters:    out..........output stream//// returns:       (int) 0//// ================================================================================int generationStatistics(FILE *out, BasicStatistics *statistics){  // output not null?  if (out==NULL)    return 0;  // print it all out  fprintf(out,"--------------------------------------------------------\n");  fprintf(out,"Generation                   : %lu\n",statistics->generation);  fprintf(out,"Fitness evaluations          : %lu\n",getFitnessCalls());  fprintf(out,"Fitness (max/avg/min)        : (%5f %5f %5f)\n",statistics->maxF,statistics->avgF,statistics->minF);  if (isBestDefined())    fprintf(out,"Percentage of optima in pop. : %1.2f\n",((float)statistics->numOptimal/statistics->N)*100);  fprintf(out,"Population bias              : ");  printGuidance(out,statistics->p1,statistics->n,statistics->guidanceThreshold);  fprintf(out,"\n");  fprintf(out,"Best solution in the pop.    : ");  printIndividual(out,statistics->bestX,statistics->n);  fprintf(out,"\n");  // get back  return 0;};// ================================================================================//// name:          fitnessStatistics//// function:      prints out the statistics on fitness//// parameters:    out..........output stream//// returns:       (int) 0//// ================================================================================int fitnessStatistics(FILE *out, BasicStatistics *statistics){  // output not null?  if (out==NULL)    return 0;  // print it all out  fprintf(out,"%3lu %7lu %10f %10f %10f\n",statistics->generation,getFitnessCalls(),statistics->maxF,statistics->avgF,statistics->minF);  // get back  return 0;};// ================================================================================//// name:          finalStatistics//// function:      prints out statistics (called after the run is over)//// parameters:    out...........output stream//                termination...the reason for terminating the algorithm (string)//// returns:       (int) 0//// ================================================================================int finalStatistics(FILE *out, char *termination, BasicStatistics *statistics){  // is output stream null?  if (out==NULL)    return 0;  // print it all out  fprintf(out, "\n=================================================================\n");  fprintf(out, "FINAL STATISTICS\n");  fprintf(out, "Termination reason           : %s\n",termination);  fprintf(out, "Generations performed        : %lu\n",statistics->generation);  fprintf(out, "Fitness evaluations          : %lu\n",getFitnessCalls());  fprintf(out, "Fitness (max/avg/min)        : (%5f %5f %5f)\n",statistics->maxF,statistics->avgF,statistics->minF);  if (isBestDefined())    fprintf(out, "Percentage of optima in pop. : %1.2f\n",((float)statistics->numOptimal/statistics->N)*100);  fprintf(out,"Population bias              : ");  printGuidance(out,statistics->p1,statistics->n,statistics->guidanceThreshold);  fprintf(out,"\n");  fprintf(out,"Best solution in the pop.    : ");  printIndividual(out,statistics->bestX,statistics->n);  fprintf(out,"\n\nThe End.\n");  // get back  return 0;};// ================================================================================//// name:          printModel//// function:      prints out the constructed model (called each generation)//// parameters:    out..........output stream//                t............number of generation//                G............the constructed network//// returns:       (int) 0//// ================================================================================int printModel(FILE *out, long t, AcyclicOrientedGraph *G){  int i,k;  // is output stream null?  if (out==NULL)    return 0;  // print out the generation number  fprintf(out,"--------------------------------------------------------\n");  fprintf(out,"Generation: %3lu\n\n",t);  // print out the model  for (i=0; i<G->size(); i++)    {      fprintf(out,"%3u",i);      fprintf(out," <- ");      for (k=0; k<G->getNumIn(i); k++)	{	  if (k!=0)	    fprintf(out,", ");	  fprintf(out,"%u",G->getParentList(i)[k]);	};      fprintf(out,"\n");    };  // get back  return 0;};// ================================================================================//// name:          printGuidance//// function:      prints out a population bias (where each bit is biased) according//                to a threshold ("0" if p1<treshold, "1" if p1>1-treshold, else //                ".")//// parameters:    out..........output stream//                p1...........univariate frequencies for 1's//                n............string length//                threshold....threshold for bias (closer than that is biased)//// returns:       (int) 0//// ================================================================================int printGuidance(FILE *out, float *p1, int n, float threshold){  int k;  float threshold1;    // compute upper threshold  threshold1 = 1-threshold;  // print where the frequencies are biased  for (k=0; k<n; k++)    if (p1[k]<threshold)      fprintf(out,"0");    else      if (p1[k]>threshold1)	fprintf(out,"1");      else	fprintf(out,".");  // get back  return 0;};