DigitizedReflection.h

 
#pragma once
 
#if defined(DigitizedReflection_RECURSES)
#error Recursive header files inclusion detected in DigitizedReflection.h
#else // defined(DigitizedReflection_RECURSES)
#define DigitizedReflection_RECURSES
 
#if !defined DigitizedReflection_h
#define DigitizedReflection_h
 
// Inclusions
#include "GAVector.h"
 
namespace DGtal {
  template<typename TSpace, typename TInputValue = typename TSpace::RealPoint, typename TOutputValue = typename TSpace::Point>
  struct Reflection
  {
    GAVector<TSpace> normalVector;
 
    explicit Reflection ( const GAVector<TSpace> & m=GAVector<TSpace>())
        :normalVector(m){}
 
    inline
    TOutputValue operator()( const TInputValue & aInput ) const
    {
      DGtal::functors::VectorRounding < TInputValue, TOutputValue > roundingOpe;
      Z2i::RealPoint m_r = Z2i::RealPoint(normalVector.my_gavec[0],normalVector.my_gavec[1]);
      Z2i::RealPoint x_r = Z2i::RealPoint(aInput[0],aInput[1]);
      Z2i::RealPoint p=x_r - 2.0*((x_r[0]*m_r[0] + x_r[1]*m_r[1])/(m_r[0]*m_r[0] + m_r[1]*m_r[1]))*m_r;
      return  roundingOpe(p);
    }
 
  };
}
#endif //DigitizedReflection
 
#undef DigitizedReflection_RECURSES
#endif // else defined(DigitizedReflection_RECURSES)