+++ /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