#pragma region CPL License
/*
Nuclex Native Framework
Copyright (C) 2002-2015 Nuclex Development Labs

This library is free software; you can redistribute it and/or
modify it under the terms of the IBM Common Public License as
published by the IBM Corporation; either version 1.0 of the
License, or (at your option) any later version.

This library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
IBM Common Public License for more details.

You should have received a copy of the IBM Common Public
License along with this library
*/
#pragma endregion // CPL License

#ifndef NUCLEX_GEOMETRY_LINES_CURVE1_H
#define NUCLEX_GEOMETRY_LINES_CURVE1_H

#include "Nuclex/Geometry/Config.h"
#include "Nuclex/Geometry/Math.h"

namespace Nuclex { namespace Geometry { namespace Lines {

  // ------------------------------------------------------------------------------------------- //

  /// <summary>Cubic bezier curve on one axis</summary>
  template <typename TScalar>
  struct Curve1 {

    /// <summary>Initializes the cubic bezier curve using the specified coordinates</summary>
    /// <param name="p0">Point at which the cubic bezier curve begins</param>
    /// <param name="p1">Point giving direction and magnitude at the curve start</param>
    /// <param name="p2">Point giving direction and magnitude at the curve end</param>
    /// <param name="p3">Point at which the cubic bezier curve ends</param>
    public: Curve1(TScalar p0, TScalar p1, TScalar p2, TScalar p3) :
      P0(p0), P1(p1), P2(p2), P3(p3) {}

    /// <summary>Offsets the entire curve by the specified amount</summary>
    /// <param name="offset">Amount the curve will be offset</param>
    public: void ShiftBy(TScalar offset) {
      this.P0 += offset;
      this.P1 += offset;
      this.P2 += offset;
      this.P3 += offset;
    }

    /// <summary>Interpolates a point along the bezier curve</summary>
    /// <param name="t">Time (0.0 .. 1.0) at which the point will be interpolated</param>
    /// <returns>The point at the specified time</returns>
    public: TScalar Interpolate(TScalar t) const {
      return
        Math<TScalar>::Power(1 - t, 3) * this->P0 +
        3 * Math<TScalar>::Power(1 - t, 2) * t * this->P1 +
        3 * (1 - t) * Math<TScalar>::Power(t, 2) * this->P2 +
        Math<TScalar>::Power(t, 3) * this->P3;
    }

    /// <summary>Position of the curve's starting point</summary>
    public: TScalar P0;
    /// <summary>Position towards which the curve's starting point extends</summary>
    public: TScalar P1;
    /// <summary>Position from which the curve enters its ending point</summary>
    public: TScalar P2;
    /// <summary>Position at which the curve ends</summary>
    public: TScalar P3;
/*
    /// <summary>The polyline's vertices</summary>
    public: std::vector< Point2<TScalar> > Vertices;
*/
  };

  // ------------------------------------------------------------------------------------------- //

}}} // namespace Nuclex::Geometry::Lines

#endif // NUCLEX_GEOMETRY_LINES_CURVE1_H