home *** CD-ROM | disk | FTP | other *** search
/ Hackers Magazine 57 / CdHackersMagazineNr57.iso / Software / Multimedia / k3d-setup-0.7.11.0.exe / include / k3d / k3dsdk / bezier.h < prev    next >
Encoding:
C/C++ Source or Header  |  2008-01-23  |  2.9 KB  |  90 lines
  1. #ifndef K3DSDK_BEZIER_H
  2. #define K3DSDK_BEZIER_H
  3.  
  4. // K-3D
  5. // Copyright (c) 1995-2004, Timothy M. Shead
  6. //
  7. // Contact: tshead@k-3d.com
  8. //
  9. // This program is free software; you can redistribute it and/or
  10. // modify it under the terms of the GNU General Public
  11. // License as published by the Free Software Foundation; either
  12. // version 2 of the License, or (at your option) any later version.
  13. //
  14. // This program is distributed in the hope that it will be useful,
  15. // but WITHOUT ANY WARRANTY; without even the implied warranty of
  16. // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  17. // General Public License for more details.
  18. //
  19. // You should have received a copy of the GNU General Public
  20. // License along with this program; if not, write to the Free Software
  21. // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  22.  
  23. /** \file
  24.         \brief Implements basic Bezier functions
  25.         \author Tim Shead (tshead@k-3d.com)
  26. */
  27.  
  28. #include "bezier_private.h"
  29. #include "vectors.h"
  30.  
  31. #include <cassert>
  32. #include <iterator>
  33. #include <vector>
  34.  
  35. namespace k3d
  36. {
  37.  
  38. /// Computes a Bezier curve value with given degree, control points, and parameter value
  39. template<int n, typename value_t, typename iterator_t>
  40. value_t bezier(const iterator_t First, const iterator_t Last, const double t)
  41. {
  42.     // Sanity checks ...
  43.     assert(std::distance(First, Last) == n+1);
  44.     
  45.     return detail::bezier_implementation<n, n, value_t, iterator_t>().value(t, Last-1);
  46. }
  47.  
  48. /// Treats a Bezier curve as a two-dimensional function, returning the curve value Y for a specific X (instead of using the curve parameter t, which is nonlinear)
  49. template<int n, typename value_t, typename iterator_t>
  50. double bezier_function(const iterator_t First, const iterator_t Last, const double X, const double Hint, const double MaxError, const unsigned long MaxIterations, double& Error, unsigned long& Iterations)
  51. {
  52.     assert(MaxIterations);
  53.  
  54.     double t = Hint;
  55.     value_t result;
  56.     for(Iterations = 1; Iterations <= MaxIterations; ++Iterations)
  57.         {
  58.             result = bezier<n, value_t, iterator_t>(First, Last, t);
  59.             Error = X - result[0];
  60.  
  61.             if(std::abs(Error) < MaxError)
  62.                 return result[1];
  63.  
  64.             t += 0.4 * Error;
  65.         }
  66.  
  67.     return result[1];
  68. }
  69.  
  70. /// Returns a Bezier / Bernstein basis for the given order, control point number, and parameter value
  71. double BernsteinBasis(const unsigned long Order, const unsigned long ControlPoint, const double Parameter);
  72.  
  73. /// Computes a Bezier curve value with given order,  control points, and parameter value
  74. template<class Type>
  75. Type Bezier(const std::vector<Type>& ControlPoints, const double Parameter)
  76. {
  77.     // Sanity checks ...
  78.     assert(ControlPoints.size() > 1);
  79.  
  80.     Type result = BernsteinBasis(ControlPoints.size(), 0, Parameter) * ControlPoints[0];
  81.     for(unsigned long i = 1; i < ControlPoints.size(); i++)
  82.         result += to_vector(BernsteinBasis(ControlPoints.size(), i, Parameter) * ControlPoints[i]);
  83.  
  84.     return result;
  85. }
  86.  
  87. } // namespace k3d
  88.  
  89. #endif // !K3DSDK_BEZIER_H
  90.