itkSpecialFunctionGenerator.h

Go to the documentation of this file.

#ifndef __itkSpecialFunctionGenerator_h
#define __itkSpecialFunctionGenerator_h

#include "utlCore.h"
#include "utlNDArray.h"


namespace utl
{
inline double 
Hyperg1F1(double a, double b, double x); 


template < class T >
T
Lagurre ( const int n, const double a, const T x );


inline double 
Gamma(const double x);

inline double 
GammaLower (const double s, const double x);

inline double 
BesselJa(const double a, const double x);

inline double 
BesselJInteger(const int n, const double x);

inline double 
BesselJIntegerPrime(const int n, const double x);

// template<class T> 
// NDArray<T,1>
// GetRotatedSHCoefficients(const NDArray<T,1>& shInput, const NDArray<T,2>& rotationMatrix);

template < class T >
utl_shared_ptr<utl::NDArray<T,1> >
GetSymmetricTensorSHCoef(const T b, const T e1, const T e2, const int lMax, const T theta=0, const T phi=0 );

template < class T >
std::vector< std::vector<T> >
GetSymmetricTensorSHCoefDerivative(const T b, const T e1, const T e2, const int lMax, const T theta=0, const T phi=0 );
 
inline double
GetExpProductLegendreCoef(const double a, const double b, const int l );

template < class T >
utl_shared_ptr<utl::NDArray<T,1> >
ComputeDWISHCoefficientsForGPDCylinder (const T radius, const T diffusivity, const T deltaBig, const T deltaSmall, const T qq, const int lMax , const T theta=0, const T phi=0);


template < class T >
double
ComputeOrientationalOrderFromSHCoefficients(const utl::NDArray<T,1>& shCoef, const utl::NDArray<T,1>& axis);

inline double
ComputeOrientationalOrderFromSymmetricTensor(const double e1, const double e2, const double phi=0);
    
}


#if !defined(__itkSpecialFunctionGenerator_hxx)
#include "itkSpecialFunctionGenerator.hxx"
#endif 

#endif