densitymap.hpp

#ifndef _DENSITY_MAP_HPP_
#define _DENSITY_MAP_HPP_
#define ABS(x) ( ( x ) < 0 ? - (x) : (x) )

#include "facetsmap/staticmaillage.hpp"
#include "facetsmap/params.hpp"

namespace jafar {
  namespace facetsmap {
    class Facet;
    class DensityMap : private StaticMaillage<double> {
      public:
        DensityMap() : StaticMaillage<double>() { }
      public:
        double computePotentiel(Facet* f);
        void updateDensity(Facet* f);
        double mapTotalValue();
      private:
        inline double computeValue(double v, int i, int j, double oo)
        {
          return MIN(v + (1. - (ABS(i) + ABS(j)) / (2.*oo)), 1.);
        }
    };
  };
};


#endif