AccFlower2D.ih

#include <cstdlib>

// IMPLEMENTATION of inline methods.

// ----------------------- Standard services ------------------------------

template <typename T>
inline
DGtal::AccFlower2D<T>::~AccFlower2D()
{
}

template <typename T>
inline
DGtal::AccFlower2D<T>::AccFlower2D(const double x0, const double y0,
           const double radius, const double smallRadius,
           const unsigned int k, const double phi)
  :  myCenter(x0,y0), myRadius(radius), myVarRadius(smallRadius), 
     myK(k), myPhi(phi)
{
  myKp = 2 * myK/ ( M_PI * M_PI );
}


template <typename T>
inline
DGtal::AccFlower2D<T>::AccFlower2D(const RealPoint2D &aPoint, const double radius,
           const double smallRadius,
           const unsigned int k, const double phi)
  : myCenter(aPoint), myRadius(radius), myVarRadius(smallRadius), 
    myK(k), myPhi(phi)
{
  myKp = 2 * myK/ ( M_PI * M_PI );
}

template <typename T>
inline
DGtal::AccFlower2D<T>::AccFlower2D(const Point &aPoint, const double radius,
           const double smallRadius,
           const unsigned int k, const double phi)
  : myRadius(radius), myVarRadius(smallRadius), myK(k), myPhi(phi)
{
  myCenter = aPoint;
  myKp = 2 * myK/ ( M_PI * M_PI );
}

// ------------- Implementation of 'StarShaped' services ------------------

template <typename T>
inline
double
DGtal::AccFlower2D<T>::parameter( const RealPoint2D & pp ) const
{
  RealPoint2D p( pp );
  p -= myCenter;
  
  double angle = 0.0;
  if ( ( p.at( 0 ) == 0.0 ) && ( p.at( 1 ) == 0.0 ) )
    return angle;

  if ( p.at( 0 ) >= p.at( 1 ) )
    {
      if ( p.at( 0 ) >= -p.at( 1 ) )
  angle = atan( p.at( 1 ) / p.at( 0 ) );
      else
  angle = 1.5* M_PI + atan( - p.at( 0 ) / p.at( 1 ) );
    }
  else // ( p.at( 0 ) >= p.at( 1 ) )
    {
      if ( p.at( 0 ) >= -p.at( 1 ) )
  angle = 0.5*M_PI - atan( p.at( 0 ) / p.at( 1 ) );
      else
  angle = M_PI + atan( p.at( 1 ) / p.at( 0 ) );
    }
  angle = ( angle < 0.0 ) ? angle + 2*M_PI : angle;
  return angle;
}

template <typename T>
inline
typename DGtal::AccFlower2D<T>::RealPoint2D 
DGtal::AccFlower2D<T>::x(const double t ) const
{
  double tt =t;
  while ( tt >= M_PI ) tt -= 2 * M_PI; 
  while ( tt < -M_PI ) tt += 2 * M_PI;
  
  double ktn = myKp * tt * tt * tt;
  double r = myRadius + myVarRadius * cos( ktn );
  RealPoint2D c( r * cos( tt ), r * sin( tt ) );
  c += myCenter;
  return c;
}


template <typename T>
inline
typename DGtal::AccFlower2D<T>::RealVector2D 
DGtal::AccFlower2D<T>::xp( const double tt ) const
{
  double t= tt;
  while ( t >= M_PI ) t -= 2 * M_PI; 
  while ( t < -M_PI ) t += 2 * M_PI;
  
  double ktn = myKp * t * t * t;
  double ktnp = 3 * myKp * t * t ;

  double r = myRadius + myVarRadius * cos( ktn );
  double rp = - myVarRadius * sin( ktn ) * ktnp;
  RealPoint2D c( rp * cos( t ) - r * sin( t ),
     rp * sin( t ) + r * cos( t ) );
  return c;
}

template <typename T>
inline
typename DGtal::AccFlower2D<T>::RealVector2D
DGtal::AccFlower2D<T>::xpp( const double tt ) const
{
  double t=tt;
  while ( t >= M_PI ) t -= 2 * M_PI; 
  while ( t < -M_PI ) t += 2 * M_PI;
  
  double ktn = myKp * t * t * t;
  double ktnp = 3 * myKp * t * t;
  double ktnpp = 6 * myKp * t;
  
  double r = myRadius + myVarRadius * cos( ktn );
  double rp = - myVarRadius * sin( ktn ) * ktnp;
  double rpp = - myVarRadius * cos( ktn ) * ktnp * ktnp - 
    myVarRadius * sin( ktn ) * ktnpp;
  
  RealPoint2D c( rpp * cos( t ) - 2 * rp * sin( t ) - r * cos( t ),
     rpp * sin( t ) + 2 * rp * cos( t ) - r * sin( t ) );
  return c;
}


// Interface - public :

template <typename T>
inline
void
DGtal::AccFlower2D<T>::selfDisplay ( std::ostream & out ) const
{
  out << "[AccFlower2D] center= "<<myCenter<<" radius="<<myRadius<<" smallradius="<<myVarRadius
      <<" myK="<<myK<<" phase-shift="<<myPhi;
}

template <typename T>
inline
bool
DGtal::AccFlower2D<T>::isValid() const
{
    return true;
}



// Implementation of inline functions                                        //

template <typename T>
inline
std::ostream&
DGtal::operator<< ( std::ostream & out, 
      const AccFlower2D<T> & object )
{
  object.selfDisplay( out );
  return out;
}

//                                                                           //