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_private.h < prev    next >
Encoding:
C/C++ Source or Header  |  2008-01-23  |  2.0 KB  |  80 lines
  1. #ifndef K3DSDK_BEZIER_PRIVATE_H
  2. #define K3DSDK_BEZIER_PRIVATE_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 internal Bezier computation functions
  25.         \author Tim Shead (tshead@k-3d.com)
  26. */
  27.  
  28. #include <cmath>
  29.  
  30. namespace k3d
  31. {
  32.  
  33. namespace detail
  34. {
  35.  
  36. /// Compile-time computation of an integer factorial
  37. template<int n>
  38. struct factorial
  39. {
  40.     enum { result = n * factorial<n-1>::result };
  41. };
  42.  
  43. template<>
  44. struct factorial<0>
  45. {
  46.     enum { result = 1 };
  47. };
  48.  
  49. template<int n, int i>
  50. double bernstein_basis(const double t)
  51. {
  52.     return (static_cast<double>(factorial<n>::result) / static_cast<double>(factorial<i>::result * factorial<n-i>::result)) * std::pow(t, i) * std::pow(1-t, n-i);
  53. }
  54.  
  55. template<int n, int i, typename value_t, typename iterator_t>
  56. struct bezier_implementation
  57. {
  58.     value_t value(double t, iterator_t B)
  59.     {
  60.         value_t a(bernstein_basis<n, i>(t) * (*B));
  61.         value_t b(bezier_implementation<n, i-1, value_t, iterator_t>().value(t, --B));
  62.         return a + b;
  63.     }
  64. };
  65.  
  66. template<int n, typename value_t, typename iterator_t>
  67. struct bezier_implementation<n, 0, value_t, iterator_t>
  68. {
  69.     value_t value(double t, iterator_t B)
  70.     {
  71.         return bernstein_basis<n, 0>(t) * (*B);
  72.     }
  73. };
  74.  
  75. } // namespace detail
  76.  
  77. } // namespace k3d
  78.  
  79. #endif // !K3DSDK_BEZIER_PRIVATE_H
  80.