--- /dev/null
+///////////////////////////////////////////////////////////////////////////
+//
+// Copyright (c) 2002, Industrial Light & Magic, a division of Lucas
+// Digital Ltd. LLC
+//
+// 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 Industrial Light & Magic nor the names of
+// its contributors 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.
+//
+///////////////////////////////////////////////////////////////////////////
+
+
+
+#ifndef INCLUDED_IMATHFRAME_H
+#define INCLUDED_IMATHFRAME_H
+
+namespace Imath {
+
+template<class T> class Vec3;
+template<class T> class Matrix44;
+
+//
+// These methods compute a set of reference frames, defined by their
+// transformation matrix, along a curve. It is designed so that the
+// array of points and the array of matrices used to fetch these routines
+// don't need to be ordered as the curve.
+//
+// A typical usage would be :
+//
+// m[0] = Imath::firstFrame( p[0], p[1], p[2] );
+// for( int i = 1; i < n - 1; i++ )
+// {
+// m[i] = Imath::nextFrame( m[i-1], p[i-1], p[i], t[i-1], t[i] );
+// }
+// m[n-1] = Imath::lastFrame( m[n-2], p[n-2], p[n-1] );
+//
+// See Graphics Gems I for the underlying algorithm.
+//
+
+template<class T> Matrix44<T> firstFrame( const Vec3<T>&, // First point
+ const Vec3<T>&, // Second point
+ const Vec3<T>& ); // Third point
+
+template<class T> Matrix44<T> nextFrame( const Matrix44<T>&, // Previous matrix
+ const Vec3<T>&, // Previous point
+ const Vec3<T>&, // Current point
+ Vec3<T>&, // Previous tangent
+ Vec3<T>& ); // Current tangent
+
+template<class T> Matrix44<T> lastFrame( const Matrix44<T>&, // Previous matrix
+ const Vec3<T>&, // Previous point
+ const Vec3<T>& ); // Last point
+
+//
+// firstFrame - Compute the first reference frame along a curve.
+//
+// This function returns the transformation matrix to the reference frame
+// defined by the three points 'pi', 'pj' and 'pk'. Note that if the two
+// vectors <pi,pj> and <pi,pk> are colinears, an arbitrary twist value will
+// be choosen.
+//
+// Throw 'NullVecExc' if 'pi' and 'pj' are equals.
+//
+
+template<class T> Matrix44<T> firstFrame
+(
+ const Vec3<T>& pi, // First point
+ const Vec3<T>& pj, // Second point
+ const Vec3<T>& pk ) // Third point
+{
+ Vec3<T> t = pj - pi; t.normalizeExc();
+
+ Vec3<T> n = t.cross( pk - pi ); n.normalize();
+ if( n.length() == 0.0f )
+ {
+ int i = fabs( t[0] ) < fabs( t[1] ) ? 0 : 1;
+ if( fabs( t[2] ) < fabs( t[i] )) i = 2;
+
+ Vec3<T> v( 0.0, 0.0, 0.0 ); v[i] = 1.0;
+ n = t.cross( v ); n.normalize();
+ }
+
+ Vec3<T> b = t.cross( n );
+
+ Matrix44<T> M;
+
+ M[0][0] = t[0]; M[0][1] = t[1]; M[0][2] = t[2]; M[0][3] = 0.0,
+ M[1][0] = n[0]; M[1][1] = n[1]; M[1][2] = n[2]; M[1][3] = 0.0,
+ M[2][0] = b[0]; M[2][1] = b[1]; M[2][2] = b[2]; M[2][3] = 0.0,
+ M[3][0] = pi[0]; M[3][1] = pi[1]; M[3][2] = pi[2]; M[3][3] = 1.0;
+
+ return M;
+}
+
+//
+// nextFrame - Compute the next reference frame along a curve.
+//
+// This function returns the transformation matrix to the next reference
+// frame defined by the previously computed transformation matrix and the
+// new point and tangent vector along the curve.
+//
+
+template<class T> Matrix44<T> nextFrame
+(
+ const Matrix44<T>& Mi, // Previous matrix
+ const Vec3<T>& pi, // Previous point
+ const Vec3<T>& pj, // Current point
+ Vec3<T>& ti, // Previous tangent vector
+ Vec3<T>& tj ) // Current tangent vector
+{
+ Vec3<T> a(0.0, 0.0, 0.0); // Rotation axis.
+ T r = 0.0; // Rotation angle.
+
+ if( ti.length() != 0.0 && tj.length() != 0.0 )
+ {
+ ti.normalize(); tj.normalize();
+ T dot = ti.dot( tj );
+
+ //
+ // This is *really* necessary :
+ //
+
+ if( dot > 1.0 ) dot = 1.0;
+ else if( dot < -1.0 ) dot = -1.0;
+
+ r = acosf( dot );
+ a = ti.cross( tj );
+ }
+
+ if( a.length() != 0.0 && r != 0.0 )
+ {
+ Matrix44<T> R; R.setAxisAngle( a, r );
+ Matrix44<T> Tj; Tj.translate( pj );
+ Matrix44<T> Ti; Ti.translate( -pi );
+
+ return Mi * Ti * R * Tj;
+ }
+ else
+ {
+ Matrix44<T> Tr; Tr.translate( pj - pi );
+
+ return Mi * Tr;
+ }
+}
+
+//
+// lastFrame - Compute the last reference frame along a curve.
+//
+// This function returns the transformation matrix to the last reference
+// frame defined by the previously computed transformation matrix and the
+// last point along the curve.
+//
+
+template<class T> Matrix44<T> lastFrame
+(
+ const Matrix44<T>& Mi, // Previous matrix
+ const Vec3<T>& pi, // Previous point
+ const Vec3<T>& pj ) // Last point
+{
+ Matrix44<T> Tr; Tr.translate( pj - pi );
+
+ return Mi * Tr;
+}
+
+} // namespace Imath
+
+#endif