stat_progress_archive.hpp

//| This file is a part of the ERC ResiBots project.
//| Copyright 2015, ISIR / Universite Pierre et Marie Curie (UPMC)
//| Main contributor(s): Jean-Baptiste Mouret, mouret@isir.fr
//|                      Antoine Cully, cully@isir.upmc.fr
//|
//| This software is governed by the CeCILL license under French law
//| and abiding by the rules of distribution of free software.  You
//| can use, modify and/ or redistribute the software under the terms
//| of the CeCILL license as circulated by CEA, CNRS and INRIA at the
//| following URL "http://www.cecill.info".
//|
//| As a counterpart to the access to the source code and rights to
//| copy, modify and redistribute granted by the license, users are
//| provided only with a limited warranty and the software's author,
//| the holder of the economic rights, and the successive licensors
//| have only limited liability.
//|
//| In this respect, the user's attention is drawn to the risks
//| associated with loading, using, modifying and/or developing or
//| reproducing the software by the user in light of its specific
//| status of free software, that may mean that it is complicated to
//| manipulate, and that also therefore means that it is reserved for
//| developers and experienced professionals having in-depth computer
//| knowledge. Users are therefore encouraged to load and test the
//| software's suitability as regards their requirements in conditions
//| enabling the security of their systems and/or data to be ensured
//| and, more generally, to use and operate it in the same conditions
//| as regards security.
//|
//| The fact that you are presently reading this means that you have
//| had knowledge of the CeCILL license and that you accept its terms.

#ifndef STAT_BEST_HPP
#define STAT_BEST_HPP
#include <sferes/stat/stat.hpp>
namespace sferes
{
  namespace stat
  {

    SFERES_STAT(ProgressArchive, Stat)
    {
    public:
      template<typename E>
	void refresh(const E& ea)
      {
	if(ea.gen()%50!=0)
	  return;

	this->_create_log_file(ea, "progress_archive.dat");	

	int archive_size=0;
	float archive_mean=0;
	float archive_max=0;
	float mean_dist_center=0;

	for (size_t i = 0; i < ea.array().size(); ++i)
	  for (size_t j = 0; j < ea.array().size(); ++j)
	    for (size_t k = 0; k < ea.array().size(); ++k)
	      for (size_t l = 0; l < ea.array().size(); ++l)
		for (size_t m = 0; m < ea.array().size(); ++m)
		  for (size_t n = 0; n < ea.array().size(); ++n)
		    if (ea.array()[i][j][k][l][m][n])
		      {
			typename E::indiv_t indiv=ea.array()[i][j][k][l][m][n];
			archive_size++;
			archive_mean+=indiv->fit().behavior().covered_distance;
			if(archive_max<indiv->fit().behavior().covered_distance)
			  archive_max=indiv->fit().behavior().covered_distance;
			float dist=0;
			for(int dc=0;dc<indiv->fit().behavior().duty_cycle.size();dc++)
			  {
			    float diff=(indiv->fit().behavior().duty_cycle[dc]-round(indiv->fit().behavior().duty_cycle[dc]*4)/4);
			    dist+=diff*diff;
			  }
			mean_dist_center+=sqrtf(dist);
		      }
	archive_mean/=archive_size;
	mean_dist_center/=archive_size;


	(*this->_log_file)<<ea.gen()<< " "<<archive_size<<" "<<archive_mean<<" "<<archive_max<<" "<<mean_dist_center<<std::endl;
	  
      } 
    };
  }
}
#endif