/usr/include/InsightToolkit/Common/itkElasticBodyReciprocalSplineKernelTransform.txx 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 | /*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkElasticBodyReciprocalSplineKernelTransform.txx
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 __itkElasticBodyReciprocalSplineKernelTransform_txx
#define __itkElasticBodyReciprocalSplineKernelTransform_txx
#include "itkElasticBodyReciprocalSplineKernelTransform.h"
namespace itk
{
template <class TScalarType, unsigned int NDimensions>
ElasticBodyReciprocalSplineKernelTransform<TScalarType, NDimensions>::
ElasticBodyReciprocalSplineKernelTransform()
{
// Alpha = 8 ( 1 - \nu ) - 1
m_Alpha = 8.0 * ( 1.0 - .25 ) - 1;
}
template <class TScalarType, unsigned int NDimensions>
ElasticBodyReciprocalSplineKernelTransform<TScalarType, NDimensions>::
~ElasticBodyReciprocalSplineKernelTransform()
{
}
/**
* This method has been deprecated as of ITK 3.6.
* Please use the method: void ComputeG(vector,gmatrix) instead.
*/
#if !defined(ITK_LEGACY_REMOVE)
template <class TScalarType, unsigned int NDimensions>
const typename ElasticBodyReciprocalSplineKernelTransform<TScalarType, NDimensions>::GMatrixType &
ElasticBodyReciprocalSplineKernelTransform<TScalarType, NDimensions>::
ComputeG( const InputVectorType & ) const
{
itkLegacyReplaceBodyMacro(itkElasticBodyReciprocalSplineKernelTransform::ComputeG_vector,
3.6,itkElasticBodyReciprocalSplineKernelTransform::ComputeG_vector_gmatrix);
return this->m_GMatrix;
}
#endif
template <class TScalarType, unsigned int NDimensions>
void
ElasticBodyReciprocalSplineKernelTransform<TScalarType, NDimensions>
::ComputeG(const InputVectorType & x, GMatrixType & gmatrix ) const
{
const TScalarType r = x.GetNorm();
const TScalarType factor =
( r > 1e-8 ) ? ( -1.0 / r ): NumericTraits<TScalarType>::Zero;
const TScalarType radial = m_Alpha * r;
for(unsigned int i=0; i<NDimensions; i++)
{
const typename InputVectorType::ValueType xi = x[i] * factor;
// G is symmetric
for(unsigned int j=0; j<i; j++)
{
const TScalarType value = xi * x[j];
gmatrix[i][j] = value;
gmatrix[j][i] = value;
}
gmatrix[i][i] = radial + xi * x[i];
}
}
template <class TScalarType, unsigned int NDimensions>
void
ElasticBodyReciprocalSplineKernelTransform<TScalarType, NDimensions>
::PrintSelf(std::ostream& os, Indent indent) const
{
Superclass::PrintSelf(os,indent);
os << indent << "m_Alpha: " << m_Alpha << std::endl;
}
} // namespace itk
#endif
|