ITK
4.1.0
Insight Segmentation and Registration Toolkit
|
00001 /*========================================================================= 00002 * 00003 * Copyright Insight Software Consortium 00004 * 00005 * Licensed under the Apache License, Version 2.0 (the "License"); 00006 * you may not use this file except in compliance with the License. 00007 * You may obtain a copy of the License at 00008 * 00009 * http://www.apache.org/licenses/LICENSE-2.0.txt 00010 * 00011 * Unless required by applicable law or agreed to in writing, software 00012 * distributed under the License is distributed on an "AS IS" BASIS, 00013 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 00014 * See the License for the specific language governing permissions and 00015 * limitations under the License. 00016 * 00017 *=========================================================================*/ 00018 #ifndef __itkCrossHelper_h 00019 #define __itkCrossHelper_h 00020 00021 #include "itkNumericTraits.h" 00022 00023 namespace itk 00024 { 00035 template< typename TVector > 00036 class CrossHelper 00037 { 00038 public: 00039 typedef TVector VectorType; 00040 typedef typename VectorType::ValueType ValueType; 00041 00042 itkStaticConstMacro (Dimension, unsigned int, VectorType::Dimension); 00043 00049 VectorType operator()(const VectorType & iU, 00050 const VectorType & iV) const 00051 { 00052 VectorType oCross; 00053 00054 if ( Dimension > 2 ) 00055 { 00056 oCross[0] = iU[1] * iV[2] - iV[1] * iU[2]; 00057 oCross[1] = iV[0] * iU[2] - iU[0] * iV[2]; 00058 oCross[2] = iU[0] * iV[1] - iV[0] * iU[1]; 00059 00060 for ( unsigned int dim = 3; dim < Dimension; dim++ ) 00061 { 00062 oCross[dim] = 0.0; 00063 } 00064 } 00065 else 00066 { 00067 oCross.Fill(0.); 00068 } 00069 00070 return oCross; 00071 } 00072 }; 00073 } 00074 00075 #endif 00076