CBDRFastSolver.h

 
#pragma once
 
#if defined(CBDRFASTSOLVER_RECURSES)
#error Recursive header files inclusion detected in CBDRFastSolver.h
#else // defined(CBDRFASTSOLVER_RECURSES)
#define CBDRFASTSOLVER_RECURSES
 
#if !defined CBDRFASTSOLVER_h
#define CBDRFASTSOLVER_h
 
// Inclusions
#include "NBijectiveReflectionGenerator.h"
namespace DGtal {
  template<typename TSpace,typename TDomain>
  class CBDRFastSolver {
    typedef CBDR_naiverotation<TSpace> BijectiveReflections;
    typedef Reflection<TSpace> DigitizedReflection;
    public:
    CBDRFastSolver(const std::vector<BijectiveReflections>& vecOptimisedReflections,const double rotAngle,
                    const typename TDomain::Point center,
                    const int km):my_vecOptimisedReflections(vecOptimisedReflections),my_angle(rotAngle),my_center(center),kmax(km){}
 
    BijectiveReflections solve() {
      double angleInDegrees = my_angle * 180.0 / M_PI;
      angleInDegrees = fmod(angleInDegrees, 360.0);
      if (angleInDegrees < 0) {
        angleInDegrees += 360.0;
      }
      return my_vecOptimisedReflections[static_cast<int>(std::round(angleInDegrees))];
    }
    protected:
    std::vector<BijectiveReflections> my_vecOptimisedReflections;
    double my_angle;
    typename TDomain::Point my_center;
    int kmax;
  };
}
 
 
 
#endif //CBDRFASTSOLVER
#undef CBDRFASTSOLVER_RECURSES
#endif // else defined(CBDRFASTSOLVER_RECURSES)