...

[FrontAlgorithms.git] / appli / CTArteries / FourierSeries.hxx

diff --git a/appli/CTArteries/FourierSeries.hxx b/appli/CTArteries/FourierSeries.hxx
deleted file mode 100644 (file)
index 6398da8..0000000
--- a/appli/CTArteries/FourierSeries.hxx
+++ /dev/null
@@ -1,102 +0,0 @@
-/* =======================================================================
- * @author: Leonardo Florez-Valencia
- * @email: florez-l@javeriana.edu.co
- * =======================================================================
- */
-#ifndef __FourierSeries__hxx__
-#define __FourierSeries__hxx__
-
-// -------------------------------------------------------------------------
-template< class S, unsigned int D >
-template< class S2, unsigned int D2 >
-FourierSeries< S, D >::
-FourierSeries( const FourierSeries< S2, D2 >& o )
-  : m_RealAxis( D - 2 ),
-    m_ImagAxis( D - 1 )
-{
-  int q = o.GetNumberOfHarmonics( );
-  this->SetNumberOfHarmonics( q );
-  for( int l = -q; l <= q; ++l )
-    ( *this )[ l ] =
-      TComplex( S( std::real( o[ l ] ) ), S( std::imag( o[ l ] ) ) );
-}
-
-// -------------------------------------------------------------------------
-template< class S, unsigned int D >
-template< class _TIt >
-FourierSeries< S, D >::
-FourierSeries( const _TIt& b, const _TIt& e, unsigned int q )
-  : m_RealAxis( D - 2 ),
-    m_ImagAxis( D - 1 )
-{
-  std::vector< TComplex > aux;
-  for( _TIt it = b; it != e; ++it )
-    aux.push_back( TComplex( S( ( *it )[ 0 ] ), S( ( *it )[ 1 ] ) ) );
-  this->_DFT( aux, q );
-}
-
-// -------------------------------------------------------------------------
-template< class S, unsigned int D >
-template< class S2, unsigned int D2 >
-typename FourierSeries< S, D >::
-Self& FourierSeries< S, D >::
-operator=( const FourierSeries< S2, D2 >& o )
-{
-  if( D == D2 )
-  {
-    this->m_RealAxis = o.GetRealAxis( );
-    this->m_ImagAxis = o.GetImagAxis( );
-  }
-  else
-  {
-    this->m_RealAxis = D - 2;
-    this->m_ImagAxis = D - 1;
-
-  } // fi
-  int q = o.GetNumberOfHarmonics( );
-  this->SetNumberOfHarmonics( q );
-  for( int l = -q; l <= q; ++l )
-    ( *this )[ l ] =
-      TComplex( S( std::real( o[ l ] ) ), S( std::imag( o[ l ] ) ) );
-  return( *this );
-}
-
-// -------------------------------------------------------------------------
-template< class S, unsigned int D >
-template< class _TOtherScalar >
-void FourierSeries< S, D >::
-SetEllipses( const std::vector< _TOtherScalar >& ellipses )
-{
-  int Q = ( ellipses.size( ) - 2 ) >> 2;
-  Q = ( Q < 1 )? 1: Q;
-  this->SetNumberOfHarmonics( Q );
-  auto cIt = ellipses.begin( );
-  S cx = ( cIt != ellipses.end( ) )? *cIt: S( 0 );
-  if( cIt != ellipses.end( ) )
-    ++cIt;
-  S cy = ( cIt != ellipses.end( ) )? *cIt: S( 0 );
-  if( cIt != ellipses.end( ) )
-    ++cIt;
-  ( *this )[ 0 ] = TComplex( cx, cy );
-  for( int l = 1; l <= Q; ++l )
-  {
-    S a = ( cIt != ellipses.end( ) )? *cIt: S( 0 );
-    if( cIt != ellipses.end( ) )
-      ++cIt;
-    S b = ( cIt != ellipses.end( ) )? *cIt: S( 0 );
-    if( cIt != ellipses.end( ) )
-      ++cIt;
-    S t = ( cIt != ellipses.end( ) )? *cIt: S( 0 );
-    if( cIt != ellipses.end( ) )
-      ++cIt;
-    S p = ( cIt != ellipses.end( ) )? *cIt: S( 0 );
-    if( cIt != ellipses.end( ) )
-      ++cIt;
-    this->SetEllipse( l, a, b, t, p );
-
-  } // rof
-}
-
-#endif // __FourierSeries__hxx__
-
-// eof - $RCSfile$