BaseCurve.hxx

/*--------------------------------------------------------------------------*\
 |                                                                          |
 |  Copyright (C) 2017                                                      |
 |                                                                          |
 |         , __                 , __                                        |
 |        /|/  \               /|/  \                                       |
 |         | __/ _   ,_         | __/ _   ,_                                |
 |         |   \|/  /  |  |   | |   \|/  /  |  |   |                        |
 |         |(__/|__/   |_/ \_/|/|(__/|__/   |_/ \_/|/                       |
 |                           /|                   /|                        |
 |                           \|                   \|                        |
 |                                                                          |
 |      Enrico Bertolazzi                                                   |
 |      Dipartimento di Ingegneria Industriale                              |
 |      Universita` degli Studi di Trento                                   |
 |      email: enrico.bertolazzi@unitn.it                                   |
 |                                                                          |
\*--------------------------------------------------------------------------*/
 
 
namespace G2lib {
 
  // - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
  /*\
   |  ____                  ____
   | | __ )  __ _ ___  ___ / ___|   _ _ ____   _____
   | |  _ \ / _` / __|/ _ \ |  | | | | '__\ \ / / _ \
   | | |_) | (_| \__ \  __/ |__| |_| | |   \ V /  __/
   | |____/ \__,_|___/\___|\____\__,_|_|    \_/ \___|
  \*/
 
  using Ipair         = std::pair<real_type,real_type>; 
  using IntersectList = std::vector<Ipair>;             
 
  /*\
   |   _       _                          _
   |  (_)_ __ | |_ ___ _ __ ___  ___  ___| |_
   |  | | '_ \| __/ _ \ '__/ __|/ _ \/ __| __|
   |  | | | | | ||  __/ |  \__ \  __/ (__| |_
   |  |_|_| |_|\__\___|_|  |___/\___|\___|\__|
  \*/
 
  class BaseCurve;
 
  bool
  collision( BaseCurve const * pC1, BaseCurve const * pC2 );
 
  bool
  collision_ISO(
    BaseCurve const * pC1,
    real_type         offs_C1,
    BaseCurve const * pC2,
    real_type         offs_C2
  );
 
  inline
  bool
  collision_SAE(
    BaseCurve const * pC1,
    real_type         offs_C1,
    BaseCurve const * pC2,
    real_type         offs_C2
  ) {
    return collision_ISO( pC1, -offs_C1, pC2, -offs_C2 );
  }
 
  void
  intersect(
    BaseCurve const * pC1,
    BaseCurve const * pC2,
    IntersectList   & ilist
  );
 
  void
  intersect_ISO(
    BaseCurve const * pC1, real_type offs_C1,
    BaseCurve const * pC2, real_type offs_C2,
    IntersectList   & ilist
  );
 
  inline
  void
  intersect_SAE(
    BaseCurve const * pC1, real_type offs_C1,
    BaseCurve const * pC2, real_type offs_C2,
    IntersectList   & ilist
  ) {
    G2lib::intersect_ISO( pC1, -offs_C1, pC2, -offs_C2, ilist );
  }
 
  #ifdef CLOTHOIDS_BACK_COMPATIBILITY
 
  inline
  bool
  collision(
    BaseCurve const * pC1,
    real_type         offs_C1,
    BaseCurve const * pC2,
    real_type         offs_C2
  ) {
    if ( G2lib::use_ISO ) return collision_ISO( pC1, offs_C1, pC2, offs_C2 );
    else                  return collision_SAE( pC1, offs_C1, pC2, offs_C2 );
  }
 
  inline
  void
  intersect(
    BaseCurve const * pC1, real_type offs_C1,
    BaseCurve const * pC2, real_type offs_C2,
    IntersectList   & ilist
  ) {
    if ( G2lib::use_ISO ) G2lib::intersect_ISO( pC1,  offs_C1, pC2,  offs_C2, ilist );
    else                  G2lib::intersect_SAE( pC1, -offs_C1, pC2, -offs_C2, ilist );
  }
 
  #endif
 
  class BaseCurve {
 
    string const m_name;
 
  public:
 
    // block default constructor
    BaseCurve( ) = delete;
    BaseCurve( BaseCurve const & ) = delete;
    BaseCurve const & operator = ( BaseCurve const & ) = delete;
 
    BaseCurve( string const & name ) : m_name(name) {}
 
    virtual
    ~BaseCurve() = default;
 
    virtual void setup( GenericContainer const & gc ) = 0;
 
    void build( GenericContainer const & gc ) { this->setup( gc ); }
 
    virtual CurveType type() const = 0;
 
    string name()      const { return m_name; }
    string type_name() const { return to_string(type()); }
 
    // . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
 
    virtual real_type length() const = 0;
 
    virtual real_type length_ISO( real_type offs ) const = 0;
 
    real_type
    length_SAE( real_type offs ) const
    { return this->length_ISO(-offs); }
 
    /*\
     |   _     _
     |  | |__ | |__   _____  __
     |  | '_ \| '_ \ / _ \ \/ /
     |  | |_) | |_) | (_) >  <
     |  |_.__/|_.__/ \___/_/\_\
    \*/
 
    virtual
    void
    bbox(
      real_type & xmin,
      real_type & ymin,
      real_type & xmax,
      real_type & ymax
    ) const = 0;
 
    virtual
    void
    bbox_ISO(
      real_type   offs,
      real_type & xmin,
      real_type & ymin,
      real_type & xmax,
      real_type & ymax
    ) const = 0;
 
    void
    bbox_SAE(
      real_type   offs,
      real_type & xmin,
      real_type & ymin,
      real_type & xmax,
      real_type & ymax
    ) const {
      this->bbox_ISO( -offs, xmin, ymin, xmax, ymax );
    }
 
    /*\
     |   _    _   _____    _                _
     |  | |__| |_|_   _| _(_)__ _ _ _  __ _| |___
     |  | '_ \ '_ \| || '_| / _` | ' \/ _` | / -_)
     |  |_.__/_.__/|_||_| |_\__,_|_||_\__, |_\___|
     |                                |___/
    \*/
 
    virtual
    void
    bb_triangles(
      vector<Triangle2D> & tvec,
      real_type            max_angle = Utils::m_pi/18,
      real_type            max_size  = 1e100,
      integer              icurve    = 0
    ) const = 0;
 
    virtual
    void
    bb_triangles_ISO(
      real_type            offs,
      vector<Triangle2D> & tvec,
      real_type            max_angle = Utils::m_pi/18,
      real_type            max_size  = 1e100,
      integer              icurve    = 0
    ) const = 0;
 
    virtual
    void
    bb_triangles_SAE(
      real_type            offs,
      vector<Triangle2D> & tvec,
      real_type            max_angle = Utils::m_pi/18,
      real_type            max_size  = 1e100,
      integer              icurve    = 0
    ) const = 0;
 
    /*\
     |   ____             _          _______           _
     |  | __ )  ___  __ _(_)_ __    / / ____|_ __   __| |
     |  |  _ \ / _ \/ _` | | '_ \  / /|  _| | '_ \ / _` |
     |  | |_) |  __/ (_| | | | | |/ / | |___| | | | (_| |
     |  |____/ \___|\__, |_|_| |_/_/  |_____|_| |_|\__,_|
     |              |___/
    \*/
 
    virtual real_type theta_begin() const { return this->theta(0); }
 
    virtual real_type theta_end() const { return this->theta(this->length()); }
 
    virtual real_type kappa_begin() const { return this->kappa(0); }
 
    virtual real_type kappa_end() const { return this->kappa(this->length()); }
 
    virtual real_type x_begin() const { return this->X(0); }
 
    virtual real_type y_begin() const { return this->Y(0); }
 
    virtual real_type x_end() const { return this->X(this->length()); }
 
    virtual real_type y_end() const { return this->Y(this->length()); }
 
    virtual real_type x_begin_ISO( real_type offs ) const { return this->X_ISO(0,offs); }
 
    virtual real_type y_begin_ISO( real_type offs ) const { return this->Y_ISO(0,offs); }
 
    virtual real_type x_end_ISO( real_type offs ) const { return this->X_ISO(this->length(),offs); }
 
    virtual real_type y_end_ISO( real_type offs ) const { return this->Y_ISO(this->length(),offs); }
 
    real_type x_begin_SAE( real_type offs ) const { return this->x_begin_ISO(-offs); }
 
    real_type y_begin_SAE( real_type offs ) const { return this->y_begin_ISO(-offs); }
 
    real_type x_end_SAE( real_type offs ) const { return this->x_end_ISO(-offs); }
 
    real_type y_end_SAE( real_type offs ) const { return this->y_end_ISO(-offs); }
 
    virtual real_type tx_begin() const { return this->tx(0); }
 
    virtual real_type ty_begin() const { return this->ty(0); }
 
    virtual real_type tx_end() const { return this->tx(this->length()); }
 
    virtual real_type ty_end() const { return this->ty(this->length()); }
 
    virtual real_type nx_begin_ISO() const { return -this->ty(0); }
 
    virtual real_type ny_begin_ISO() const { return this->tx(0); }
 
    virtual real_type nx_end_ISO() const { return -this->ty(this->length()); }
 
    virtual real_type ny_end_ISO() const { return this->tx(this->length()); }
 
    real_type nx_begin_SAE() const { return -nx_begin_ISO(); }
 
    real_type ny_begin_SAE() const { return -ny_begin_ISO(); }
 
    real_type nx_end_SAE() const { return -nx_end_ISO(); }
 
    real_type ny_end_SAE() const { return -ny_end_ISO(); }
 
    /*\
     |  _   _          _
     | | |_| |__   ___| |_ __ _
     | | __| '_ \ / _ \ __/ _` |
     | | |_| | | |  __/ || (_| |
     |  \__|_| |_|\___|\__\__,_|
    \*/
 
    virtual real_type theta( real_type s ) const = 0;
 
    virtual real_type theta_D( real_type s ) const = 0;
 
    virtual real_type theta_DD( real_type s ) const = 0;
 
    virtual real_type theta_DDD( real_type s ) const = 0;
 
    /*\
     |   _
     |  | | ____ _ _ __  _ __   __ _
     |  | |/ / _` | '_ \| '_ \ / _` |
     |  |   < (_| | |_) | |_) | (_| |
     |  |_|\_\__,_| .__/| .__/ \__,_|
     |            |_|   |_|
    \*/
 
    real_type kappa( real_type s ) const { return theta_D(s); }
 
    real_type kappa_D( real_type s ) const { return theta_DD(s); }
 
    real_type kappa_DD( real_type s ) const { return theta_DDD(s); }
 
    /*\
     |  _____                   _   _   _
     | |_   _|   __ _ _ __   __| | | \ | |
     |   | |    / _` | '_ \ / _` | |  \| |
     |   | |   | (_| | | | | (_| | | |\  |
     |   |_|    \__,_|_| |_|\__,_| |_| \_|
    \*/
 
    virtual real_type tx( real_type s ) const;
 
    virtual real_type ty( real_type s ) const;
 
    virtual real_type tx_D( real_type s ) const;
 
    virtual real_type ty_D( real_type s ) const;
 
    virtual real_type tx_DD( real_type s ) const;
 
    virtual real_type ty_DD( real_type s ) const;
 
    virtual real_type tx_DDD( real_type s ) const;
 
    virtual real_type ty_DDD( real_type s ) const;
 
    // . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
    real_type nx_ISO( real_type s ) const { return -ty(s); }
 
    real_type nx_ISO_D( real_type s ) const { return -ty_D(s); }
 
    real_type nx_ISO_DD( real_type s ) const { return -ty_DD(s); }
 
    real_type nx_ISO_DDD( real_type s ) const { return -ty_DDD(s); }
 
    real_type ny_ISO( real_type s ) const { return tx(s); }
 
    real_type ny_ISO_D( real_type s ) const { return tx_D(s); }
 
    real_type ny_ISO_DD( real_type s ) const { return tx_DD(s); }
 
    real_type ny_ISO_DDD( real_type s ) const { return tx_DDD(s); }
 
    real_type nx_SAE( real_type s ) const { return ty(s); }
 
    real_type nx_SAE_D( real_type s ) const { return ty_D(s); }
 
    real_type nx_SAE_DD( real_type s ) const { return ty_DD(s); }
 
    real_type nx_SAE_DDD( real_type s ) const { return ty_DDD(s); }
 
    real_type ny_SAE( real_type s ) const { return -tx(s); }
 
    real_type ny_SAE_D( real_type s ) const { return -tx_D(s); }
 
    real_type ny_SAE_DD ( real_type s ) const { return -tx_DD(s); }
 
    real_type ny_SAE_DDD( real_type s ) const { return -tx_DDD(s); }
 
    // . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
 
    virtual
    void
    tg( real_type s, real_type & tg_x, real_type & tg_y ) const {
      tg_x = this->tx(s);
      tg_y = this->ty(s);
    }
 
    virtual
    void
    tg_D( real_type s, real_type & tg_x_D, real_type & tg_y_D ) const {
      tg_x_D = this->tx_D(s);
      tg_y_D = this->ty_D(s);
    }
 
    virtual
    void
    tg_DD( real_type s, real_type & tg_x_DD, real_type & tg_y_DD ) const {
      tg_x_DD = this->tx_DD(s);
      tg_y_DD = this->ty_DD(s);
    }
 
    virtual
    void
    tg_DDD( real_type s, real_type & tg_x_DDD, real_type & tg_y_DDD ) const {
      tg_x_DDD = this->tx_DDD(s);
      tg_y_DDD = this->ty_DDD(s);
    }
 
    // . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
 
    void
    nor_ISO( real_type s, real_type & nx, real_type & ny ) const
    { tg( s, ny, nx ); nx = -nx; }
 
    void
    nor_ISO_D( real_type s, real_type & nx_D, real_type & ny_D ) const
    { tg_D( s, ny_D, nx_D ); nx_D = -nx_D; }
 
    void
    nor_ISO_DD( real_type s, real_type & nx_DD, real_type & ny_DD ) const
    { tg_DD( s, ny_DD, nx_DD ); nx_DD = -nx_DD; }
 
    void
    nor_ISO_DDD( real_type s, real_type & nx_DDD, real_type & ny_DDD ) const
    { tg_DDD( s, ny_DDD, nx_DDD ); nx_DDD = -nx_DDD; }
 
    void
    nor_SAE( real_type s, real_type & nx, real_type & ny ) const
    { tg( s, ny, nx ); ny = -ny; }
 
    void
    nor_SAE_D( real_type s, real_type & nx_D, real_type & ny_D ) const
    { tg_D( s, ny_D, nx_D ); ny_D = -ny_D; }
 
    void
    nor_SAE_DD( real_type s, real_type & nx_DD, real_type & ny_DD ) const
    { tg_DD( s, ny_DD, nx_DD ); ny_DD = -ny_DD; }
 
    void
    nor_SAE_DDD( real_type s, real_type & nx_DDD, real_type & ny_DDD ) const
    { tg_DDD( s, ny_DDD, nx_DDD ); ny_DDD = -ny_DDD; }
 
    // . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
 
    virtual
    void
    evaluate(
      real_type   s,
      real_type & th,
      real_type & k,
      real_type & x,
      real_type & y
    ) const {
      eval( s, x, y );
      th = theta( s );
      k  = theta_D( s );
    }
 
    virtual
    void
    evaluate_ISO(
      real_type   s,
      real_type   offs,
      real_type & th,
      real_type & k,
      real_type & x,
      real_type & y
    ) const {
      eval_ISO( s, offs, x, y );
      th = theta( s );
      k  = theta_D( s );
      k /= 1+offs*k; // scale curvature
    }
 
    virtual
    void
    evaluate_SAE(
      real_type   s,
      real_type   offs,
      real_type & th,
      real_type & k,
      real_type & x,
      real_type & y
    ) const {
      eval_SAE( s, offs, x, y );
      th = theta( s );
      k  = theta_D( s );
      k /= 1-offs*k; // scale curvature
    }
 
    // . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
 
    virtual real_type X( real_type s ) const = 0;
 
    virtual real_type Y( real_type s ) const = 0;
 
    virtual real_type X_D( real_type s ) const = 0;
 
    virtual real_type Y_D( real_type s ) const = 0;
 
    virtual real_type X_DD( real_type s ) const = 0;
 
    virtual real_type Y_DD( real_type s ) const = 0;
 
    virtual real_type X_DDD( real_type s ) const = 0;
 
    virtual real_type Y_DDD( real_type s ) const = 0;
 
    // . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
 
    virtual
    void
    eval( real_type s, real_type & x, real_type & y ) const = 0;
 
    virtual
    void
    eval_D( real_type s, real_type & x_D, real_type & y_D ) const = 0;
 
    virtual
    void
    eval_DD( real_type s, real_type & x_DD, real_type & y_DD ) const = 0;
 
    virtual
    void
    eval_DDD( real_type s, real_type & x_DDD, real_type & y_DDD ) const = 0;
 
    /*\
     |         __  __          _
     |   ___  / _|/ _|___  ___| |_
     |  / _ \| |_| |_/ __|/ _ \ __|
     | | (_) |  _|  _\__ \  __/ |_
     |  \___/|_| |_| |___/\___|\__|
    \*/
 
    virtual real_type X_ISO( real_type s, real_type offs ) const;
 
    virtual real_type Y_ISO( real_type s, real_type offs ) const;
 
    virtual real_type X_ISO_D( real_type s, real_type offs ) const;
 
    virtual real_type Y_ISO_D( real_type s, real_type offs ) const;
 
    virtual real_type X_ISO_DD( real_type s, real_type offs ) const;
 
    virtual real_type Y_ISO_DD( real_type s, real_type offs ) const;
 
    virtual real_type X_ISO_DDD( real_type s, real_type offs ) const;
 
    virtual real_type Y_ISO_DDD( real_type s, real_type offs ) const;
 
    real_type X_SAE( real_type s, real_type offs ) const { return this->X_ISO(s,-offs); }
 
    real_type Y_SAE( real_type s, real_type offs ) const { return this->Y_ISO(s,-offs); }
 
    real_type X_SAE_D( real_type s, real_type offs ) const { return this->X_ISO_D(s,-offs); }
 
    real_type Y_SAE_D( real_type s, real_type offs ) const { return this->Y_ISO_D(s,-offs); }
 
    real_type X_SAE_DD( real_type s, real_type offs ) const { return this->X_ISO_DD(s,-offs); }
 
    real_type Y_SAE_DD( real_type s, real_type offs ) const { return this->Y_ISO_DD(s,-offs); }
 
    real_type X_SAE_DDD( real_type s, real_type offs ) const { return this->X_ISO_DDD(s,-offs); }
 
    real_type Y_SAE_DDD( real_type s, real_type offs ) const { return this->Y_ISO_DDD(s,-offs); }
 
    // . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
 
    virtual
    void
    eval_ISO(
      real_type   s,
      real_type   offs,
      real_type & x,
      real_type & y
    ) const;
 
    void
    eval_SAE(
      real_type   s,
      real_type   offs,
      real_type & x,
      real_type & y
    ) const {
      this->eval_ISO( s, -offs, x, y );
    }
 
    virtual
    void
    eval_ISO_D(
      real_type   s,
      real_type   offs,
      real_type & x_D,
      real_type & y_D
    ) const;
 
    void
    eval_SAE_D(
      real_type   s,
      real_type   offs,
      real_type & x_D,
      real_type & y_D
    ) const {
      this->eval_ISO_D( s, -offs, x_D, y_D );
    }
 
    virtual
    void
    eval_ISO_DD(
      real_type   s,
      real_type   offs,
      real_type & x_DD,
      real_type & y_DD
    ) const;
 
    void
    eval_SAE_DD(
      real_type   s,
      real_type   offs,
      real_type & x_DD,
      real_type & y_DD
    ) const {
      this->eval_ISO_DD( s, -offs, x_DD, y_DD );
    }
 
    virtual
    void
    eval_ISO_DDD(
      real_type   s,
      real_type   offs,
      real_type & x_DDD,
      real_type & y_DDD
    ) const;
 
    void
    eval_SAE_DDD(
      real_type   s,
      real_type   offs,
      real_type & x_DDD,
      real_type & y_DDD
    ) const {
      this->eval_ISO_DDD( s, -offs, x_DDD, y_DDD );
    }
 
    /*\
     |  _                        __
     | | |_ _ __ __ _ _ __  ___ / _| ___  _ __ _ __ ___
     | | __| '__/ _` | '_ \/ __| |_ / _ \| '__| '_ ` _ \
     | | |_| | | (_| | | | \__ \  _| (_) | |  | | | | | |
     |  \__|_|  \__,_|_| |_|___/_|  \___/|_|  |_| |_| |_|
    \*/
 
    virtual
    void
    translate( real_type tx, real_type ty ) = 0;
 
    virtual
    void
    rotate( real_type angle, real_type cx, real_type cy ) = 0;
 
    virtual
    void
    scale( real_type sc ) = 0;
 
    virtual
    void
    reverse() = 0;
 
    virtual
    void
    change_origin( real_type newx0, real_type newy0 ) = 0;
 
    virtual
    void
    trim( real_type s_begin, real_type s_end ) = 0;
 
    /*\
     |   _       _                          _
     |  (_)_ __ | |_ ___ _ __ ___  ___  ___| |_
     |  | | '_ \| __/ _ \ '__/ __|/ _ \/ __| __|
     |  | | | | | ||  __/ |  \__ \  __/ (__| |_
     |  |_|_| |_|\__\___|_|  |___/\___|\___|\__|
    \*/
 
    virtual bool collision( BaseCurve const * pC ) const = 0;
 
    virtual
    bool
    collision_ISO(
      real_type         offs,
      BaseCurve const * pC,
      real_type         offs_C
    ) const = 0;
 
    bool
    collision_SAE(
      real_type         offs,
      BaseCurve const * pC,
      real_type         offs_C
    ) const {
      return this->collision_ISO( -offs, pC, -offs_C );
    }
 
    virtual
    void
    intersect(
      BaseCurve const * pC,
      IntersectList   & ilist
    ) const = 0;
 
    virtual
    void
    intersect_ISO(
      real_type         offs,
      BaseCurve const * pC,
      real_type         offs_C,
      IntersectList   & ilist
    ) const = 0;
 
    void
    intersect_SAE(
      real_type         offs,
      BaseCurve const * pC,
      real_type         offs_C,
      IntersectList   & ilist
    ) const {
      this->intersect_ISO( -offs, pC, -offs_C, ilist );
    }
 
    /*\
     |      _ _     _
     |   __| (_)___| |_ __ _ _ __   ___ ___
     |  / _` | / __| __/ _` | '_ \ / __/ _ \
     | | (_| | \__ \ || (_| | | | | (_|  __/
     |  \__,_|_|___/\__\__,_|_| |_|\___\___|
    \*/
 
    virtual
    integer
    closest_point_ISO(
      real_type   qx,
      real_type   qy,
      real_type & x,
      real_type & y,
      real_type & s,
      real_type & t,
      real_type & dst
    ) const = 0;
 
    integer
    closest_point_SAE(
      real_type   qx,
      real_type   qy,
      real_type & x,
      real_type & y,
      real_type & s,
      real_type & t,
      real_type & dst
    ) const {
      integer res = this->closest_point_ISO( qx, qy, x, y, s, t, dst );
      t = -t;
      return res;
    }
 
    virtual
    integer // true if projection is unique and orthogonal
    closest_point_ISO(
      real_type   qx,
      real_type   qy,
      real_type   offs,
      real_type & x,
      real_type & y,
      real_type & s,
      real_type & t,
      real_type & dst
    ) const = 0;
 
    integer
    closest_point_SAE(
      real_type   qx,
      real_type   qy,
      real_type   offs,
      real_type & x,
      real_type & y,
      real_type & s,
      real_type & t,
      real_type & dst
    ) const {
      integer res = this->closest_point_ISO( qx, qy, -offs, x, y, s, t, dst );
      t = -t;
      return res;
    }
 
    virtual
    real_type
    distance( real_type qx, real_type qy ) const {
      real_type x, y, s, t, dst;
      closest_point_ISO( qx, qy, x, y, s, t, dst );
      return dst;
    }
 
    real_type
    distance_ISO(
      real_type qx,
      real_type qy,
      real_type offs
    ) const {
      real_type x, y, s, t, dst;
      this->closest_point_ISO( qx, qy, offs, x, y, s, t, dst );
      return dst;
    }
 
    real_type
    distance_SAE(
      real_type qx,
      real_type qy,
      real_type offs
    ) const {
      real_type x, y, s, t, dst;
      this->closest_point_SAE( qx, qy, offs, x, y, s, t, dst );
      return dst;
    }
 
    /*\
     |    __ _           _ ____ _____
     |   / _(_)_ __   __| / ___|_   _|
     |  | |_| | '_ \ / _` \___ \ | |
     |  |  _| | | | | (_| |___) || |
     |  |_| |_|_| |_|\__,_|____/ |_|
    \*/
 
    bool
    findST_ISO(
      real_type   x,
      real_type   y,
      real_type & s,
      real_type & t
    ) const {
      real_type X, Y, dst;
      integer icode = this->closest_point_ISO( x, y, X, Y, s, t, dst );
      return icode >= 0;
    }
 
    bool
    findST_SAE(
      real_type   x,
      real_type   y,
      real_type & s,
      real_type & t
    ) const {
      real_type X, Y, dst;
      integer icode = this->closest_point_SAE( x, y, X, Y, s, t, dst );
      return icode >= 0;
    }
 
    // - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
 
    virtual
    void
    info( ostream_type & stream ) const = 0;
 
#ifdef CLOTHOIDS_BACK_COMPATIBILITY
#include "BaseCurve_compatibility.hxx"
#endif
 
  };
 
}