text: ##########################################################################
  #
  #  Copyright (c) 2008, Image Engine Design Inc. All rights reserved.
  #
  #  Redistribution and use in source and binary forms, with or without
  #  modification, are permitted provided that the following conditions are
  #  met:
  #
  #     * Redistributions of source code must retain the above copyright
  #       notice, this list of conditions and the following disclaimer.
  #
  #     * Redistributions in binary form must reproduce the above copyright
  #       notice, this list of conditions and the following disclaimer in the
  #       documentation and/or other materials provided with the distribution.
  #
  #     * Neither the name of Image Engine Design nor the names of any
  #       other contributors to this software may be used to endorse or
  #       promote products derived from this software without specific prior
  #       written permission.
  #
  #  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS
  #  IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO,
  #  THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
  #  PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
  #  CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
  #  EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
  #  PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
  #  PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
  #  LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
  #  NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
  #  SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  #
  ##########################################################################
  
  import unittest
  import imath
  import IECore
  import math
  
  class BezierAlgoTest( unittest.TestCase ) :
  
  	def testCubic( self ) :
  
  		def c( p, v ) :
  
  			p.append( imath.V3f( v[0], v[1], 0 ) )
  
  		p = IECore.V3fVectorData()
  		IECore.bezierSubdivide( imath.V2f( 0 ), imath.V2f( 0, 1 ), imath.V2f( 1, 1 ), imath.V2f( 1, 0 ), 0.01, lambda v : c( p, v ) )
  
  		self.assertEqual( p[0], imath.V3f( 0 ) )
  		self.assertEqual( p[-1], imath.V3f( 1, 0, 0 ) )
  
  
  		b = imath.Box3f( imath.V3f( 0, 0, 0 ), imath.V3f( 1, 1, 0 ) )
  		for i in range( 0, p.size()-1 ) :
  
  			self.assertTrue( b.intersects( p[i] ) )
  			self.assertTrue( (p[i] - p[i+1]).length() < 0.18 )
  
  	def testQuadratic( self ) :
  
  		def c( p, v ) :
  
  			p.append( imath.V3f( v[0], v[1], 0 ) )
  
  		p = IECore.V3fVectorData()
  		IECore.bezierSubdivide( imath.V2f( 0 ), imath.V2f( 1, 1 ), imath.V2f( 2, 0 ), 0.01, lambda v : c( p, v ) )
  
  		self.assertEqual( p[0], imath.V3f( 0 ) )
  		self.assertEqual( p[-1], imath.V3f( 2, 0, 0 ) )
  
  
  		b = imath.Box3f( imath.V3f( 0, 0, 0 ), imath.V3f( 2, 1, 0 ) )
  		for i in range( 0, p.size()-1 ) :
  
  			self.assertTrue( b.intersects( p[i] ) )
  			self.assertTrue( (p[i] - p[i+1]).length() < 0.18 )
  
  		# useful for making a visual check
  		#m = MeshPrimitive( IntVectorData( [ p.size() ] ), IntVectorData( range(p.size() ) ), "linear", p )
  		#Writer.create( m, "/tmp/bezier.cob" ).write()
  
  
  if __name__ == "__main__":
      unittest.main()
  