prev_no[ 1 ] = _TScalar( 1 );
for( unsigned int i = 0; i < N; ++i )
{
+ auto ntg = tg[ i ].GetNorm( );
+ if( double( ntg ) > double( 0 ) )
+ tg[ i ] /= ntg;
+
_TScalar ct = prev_tg * tg[ i ];
TVector a = itk::CrossProduct( prev_tg, tg[ i ] );
_TScalar st = a.GetNorm( );
( prev_no * ct ) +
( itk::CrossProduct( a, prev_no ) * st ) +
( a * ( ( a * prev_no ) * ( _TScalar( 1 ) - ct ) ) );
+ auto nno = no[ i ].GetNorm( );
+ if( double( nno ) > double( 0 ) )
+ no[ i ] /= nno;
+
bn[ i ] = itk::CrossProduct( tg[ i ], no[ i ] );
prev_tg = tg[ i ];
} // rof
this->m_Frames.push_back( m );
-
} // rof
this->m_FramesUpdated = true;