staticmaillage.hpp

#ifndef _STATIC_MAILLAGE_HPP_
#define _STATIC_MAILLAGE_HPP_

#include <map>

#include "kernel/jafarMacro.hpp"

namespace jafar {
  namespace facetsmap {
    class MaillagePoint {
      template<class _T1, class _T2>
      friend struct std::pair;
      private:
        MaillagePoint() : i(0), j(0) {}
      public:
        MaillagePoint(int ni, int nj) : i(ni), j(nj) { }
        int i,j;
        inline bool operator<(const MaillagePoint& m ) const
        {
          if( i < m.i) return true;
          if( i == m.i) return j < m.j;
          return false;
        }
    };
    template<typename T>
    class StaticMaillage {
      protected:
        typedef typename std::map<MaillagePoint, T> maillage;
        typedef typename maillage::const_iterator maillage_cit;
      public:
        StaticMaillage() { };
        inline T get(int i, int j) {
          typename maillage::iterator it;
          if( (it = m_maillage.find(MaillagePoint(i,j))) != m_maillage.end() )
          {
            return it->second;
          } else {
            return 0.;
          }
        }
        inline void set(int i, int j, T v)
        {
          m_maillage[MaillagePoint(i,j)] = v;
        }
        T totalValue()
        {
          T tv = 0;
          for(maillage_cit it = begin(); it != end(); it++)
          {
              tv += it->second;
          }
          return tv;
        }
      protected:
        inline maillage_cit begin() { return m_maillage.begin(); }
        inline maillage_cit end() { return m_maillage.end(); }
      private:
        maillage m_maillage;
    };
  };
};

#endif