/usr/include/InsightToolkit/Numerics/itkSymmetricEigenSystem.h is in libinsighttoolkit3-dev 3.20.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 | /*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkSymmetricEigenSystem.h
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#ifndef __itkSymmetricEigenSystem_h
#define __itkSymmetricEigenSystem_h
#include "itkObject.h"
#include "vnl/algo/vnl_symmetric_eigensystem.h"
#include "itkObjectFactory.h"
#include "itkObject.h"
#include "itkFixedArray.h"
#include "itkMatrix.h"
#include "itkObjectFactory.h"
namespace itk
{
/** \class SymmetricEigenSystem
* \brief wrapper of the vnl_symmetric_eigensystem algorithm
*
* This class is not thread-safe. If you are interested in thread-safety please
* use the class SymmetricEigenAnalysis in Insight/Code/Common.
*
* \warning THIS CLASS IS DEPRECATED AND IT IS SCHEDULED FOR BEING REMOVED
* FROM THE TOOLKIT IN RELEASE 2.4
*
* \sa SymmetricEigenAnalysis
*
* \ingroup Numerics Deprecated
*/
template< class TMatrixElement, int VNumberOfRows >
class ITK_EXPORT SymmetricEigenSystem : public Object
{
public:
/** Standard "Self" typedef. */
typedef SymmetricEigenSystem Self;
typedef Object Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Run-time type information (and related methods). */
itkTypeMacro( SymmetricEigenSystem, Object );
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** 1D array typedef */
typedef FixedArray< TMatrixElement, VNumberOfRows > ArrayType;
/** 2D array typedef */
typedef FixedArray< ArrayType, VNumberOfRows > Array2DType;
/** Array for eigen vectors */
typedef Array2DType EigenVectorArrayType;
/** Array type for eigen values */
typedef ArrayType EigenValueArrayType;
/** Matrix Type */
typedef Matrix< TMatrixElement, VNumberOfRows, VNumberOfRows > MatrixType;
/** Internal eigen sytem type. */
typedef vnl_symmetric_eigensystem< TMatrixElement > InternalEigenSystemType;
/** Set/Get the target matrix for the eigen analysis */
itkSetObjectMacro(Matrix, MatrixType);
MatrixType* GetMatrix()
{ return m_Matrix; }
/** Set/Get the absolute order flag.
* By setting this flag true, after the calculation of eigen vectors
* and values, if the absolute eigen value of eigen vector[j > i] is
* greater that of eigen vector[i], reorder the eigen vectors so that
* every absolute eigen values of eigen vector[j < i] is always greater than or
* equal to that of the eigen vectors[i] */
itkSetMacro(UseAbsoluteOrder, bool);
itkGetMacro(UseAbsoluteOrder, bool);
/** returns the eigen vectors in a 2D array */
EigenVectorArrayType* GetEigenVectors()
{ return &m_EigenVectors; }
/** returns the eigen values in an 1D array */
EigenValueArrayType* GetEigenValues()
{ return &m_EigenValues; }
/** dummy method that calls the GenerateData method to
* produce the eigen vectors and values. */
void Update()
{ this->GenerateData(); }
protected:
SymmetricEigenSystem();
virtual ~SymmetricEigenSystem();
void PrintSelf(std::ostream& os, Indent indent) const;
/** Produces the eigen vectors and values. */
void GenerateData();
private:
SymmetricEigenSystem(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
/** the target matrix */
MatrixType* m_Matrix;
/** eigen vectors output */
EigenVectorArrayType m_EigenVectors;
/** eigen values output */
EigenValueArrayType m_EigenValues;
/** flag for absolute ordering of eigen vectors and
* eigen values */
bool m_UseAbsoluteOrder;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkSymmetricEigenSystem.txx"
#endif
#endif
|