facetscloseness.hpp

//
// C++ Interface: %{MODULE}
//
// Description: 
//
//
// Author: %{AUTHOR} <%{EMAIL}>, (C) %{YEAR}
//
// Copyright: See COPYING file that comes with this distribution
//
//
#ifndef JAFAR_FACETSMAPFACETSCLOSENESS_H
#define JAFAR_FACETSMAPFACETSCLOSENESS_H

#include "jmath/jblas.hpp"
#include "facetsmap/facet.hpp"

namespace jafar {

namespace facetsmap {

class FacetsCloseness{
  public:
    FacetsCloseness(Facet* f1, Facet* f2);
    ~FacetsCloseness();
  public:
//     void computeCloseness();
    inline double distance() const {
      if(m_distanceNeedUpdate)
      {
        m_distance = ublas::norm_2( m_f1->globalOrigin() - m_f2->globalOrigin()  );
        m_distanceNeedUpdate = false;
      }
      return m_distance;
    }
    inline double parallelism() const
    {
      if(m_parallelismNeedUpdate)
      {
        m_parallelism = ublas::inner_prod(m_f1->globalNormal(), m_f2->globalNormal() );
        m_parallelismNeedUpdate = false;
      }
      return m_parallelism;
    }
    inline double zncc() const {
      if(m_znccNeedUpdate)
      {
        m_zncc = m_f1->compare(m_f2);
        m_znccNeedUpdate = false;
      }
      return m_zncc;
    }
  public:
    Facet * getF1() { return m_f1; };
    Facet * getF2() { return m_f2; };
  private:
    Facet* m_f1;
    Facet* m_f2;
    mutable double m_distance;
    mutable double m_parallelism;
    mutable double m_zncc;
    mutable bool m_distanceNeedUpdate, m_parallelismNeedUpdate, m_znccNeedUpdate;
};

};

};

#endif