Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkCrossHelper_h
00018 #define __itkCrossHelper_h
00019
00020 #include "itkNumericTraits.h"
00021
00022 namespace itk
00023 {
00033 template< typename TVector >
00034 class CrossHelper
00035 {
00036 public:
00037 typedef TVector VectorType;
00038 typedef typename VectorType::ValueType ValueType;
00039
00040 itkStaticConstMacro ( Dimension, unsigned int, VectorType::Dimension );
00041
00047 VectorType operator ( ) ( const VectorType& iU,
00048 const VectorType& iV ) const
00049 {
00050 VectorType oCross;
00051
00052 if( Dimension > 2 )
00053 {
00054 oCross[0] = iU[1] * iV[2] - iV[1] * iU[2];
00055 oCross[1] = iV[0] * iU[2] - iU[0] * iV[2];
00056 oCross[2] = iU[0] * iV[1] - iV[0] * iU[1];
00057
00058 for( unsigned int dim = 3; dim < Dimension; dim++ )
00059 {
00060 oCross[dim] = 0.0;
00061 }
00062 }
00063 else
00064 {
00065 oCross.Fill( 0. );
00066 }
00067
00068 return oCross;
00069 }
00070 };
00071 }
00072
00073 #endif
00074