#include <linterp.h>

Inheritance diagram for NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >:

Collaboration diagram for NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >:

Public Types
typedef T	value_type

typedef ArrayRefCountT	array_ref_count_type

typedef GridRefCountT	grid_ref_count_type

typedef boost::numeric::ublas::array_adaptor< T >	grid_type

typedef boost::const_multi_array_ref< T, N >	array_type

typedef std::unique_ptr< array_type >	array_type_ptr

Public Member Functions
template<class IterT1 , class IterT2 , class IterT3 >
	NDInterpolator (IterT1 grids_begin, IterT2 grids_len_begin, IterT3 f_begin, IterT3 f_end)

template<class IterT1 , class IterT2 , class IterT3 , class RefCountIterT >
	NDInterpolator (IterT1 grids_begin, IterT2 grids_len_begin, IterT3 f_begin, IterT3 f_end, ArrayRefCountT &refF, RefCountIterT grid_refs_begin)

template<class IterT1 , class IterT2 , class IterT3 >
void	init (IterT1 grids_begin, IterT2 grids_len_begin, IterT3 f_begin, IterT3 f_end)

template<class IterT1 , class IterT2 , class IterT3 , class RefCountIterT >
void	init_refcount (IterT1 grids_begin, IterT2 grids_len_begin, IterT3 f_begin, IterT3 f_end, ArrayRefCountT &refF, RefCountIterT grid_refs_begin)

template<class IterT1 , class IterT2 >
void	set_grids (IterT1 grids_begin, IterT2 grids_len_begin, bool bCopy)

template<class IterT1 , class RefCountIterT >
void	set_grids_refcount (RefCountIterT refs_begin, RefCountIterT refs_end)

template<class IterT >
void	set_f_array (IterT f_begin, IterT f_end, bool bCopy)

void	set_f_refcount (ArrayRefCountT &refF)

int	find_cell (int dim, T x) const

T	get_f_val (array< int, N > const &cell_index, array< int, N > const &v_index) const

T	get_f_val (array< int, N > const &cell_index, int v) const

Public Attributes
array_type_ptr	m_pF

ArrayRefCountT	m_ref_F

vector< T >	m_F_copy

vector< grid_type >	m_grid_list

vector< GridRefCountT >	m_grid_ref_list

vector< vector< T > >	m_grid_copy_list

Static Public Attributes
static const int	m_N = N

static const bool	m_bCopyData = CopyData

static const bool	m_bContinuous = Continuous

Detailed Description

template<int N, class T, bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>
class NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >

Definition at line 71 of file linterp.h.

Member Typedef Documentation

◆ array_ref_count_type

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

typedef ArrayRefCountT NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::array_ref_count_type

Definition at line 74 of file linterp.h.

◆ array_type

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

typedef boost::const_multi_array_ref<T, N> NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::array_type

Definition at line 82 of file linterp.h.

◆ array_type_ptr

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

typedef std::unique_ptr<array_type> NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::array_type_ptr

Definition at line 83 of file linterp.h.

◆ grid_ref_count_type

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

typedef GridRefCountT NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::grid_ref_count_type

Definition at line 75 of file linterp.h.

◆ grid_type

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

typedef boost::numeric::ublas::array_adaptor<T> NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::grid_type

Definition at line 81 of file linterp.h.

◆ value_type

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

typedef T NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::value_type

Definition at line 73 of file linterp.h.

Constructor & Destructor Documentation

◆ NDInterpolator() [1/2]

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

template<class IterT1 , class IterT2 , class IterT3 >

NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::NDInterpolator	(	IterT1	grids_begin,
		IterT2	grids_len_begin,
		IterT3	f_begin,
		IterT3	f_end
	)

inline

Definition at line 96 of file linterp.h.

                                                                                            {
     init(grids_begin, grids_len_begin, f_begin, f_end);  
   }

◆ NDInterpolator() [2/2]

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

template<class IterT1 , class IterT2 , class IterT3 , class RefCountIterT >

NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::NDInterpolator	(	IterT1	grids_begin,
		IterT2	grids_len_begin,
		IterT3	f_begin,
		IterT3	f_end,
		ArrayRefCountT &	refF,
		RefCountIterT	grid_refs_begin
	)

inline

Definition at line 102 of file linterp.h.

                                                                                                                                                 {
     init_refcount(grids_begin, grids_len_begin, f_begin, f_end, refF, grid_refs_begin);
   } 

Member Function Documentation

◆ find_cell()

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

int NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::find_cell	(	int	dim,
		T	x
	)		const

inline

Definition at line 174 of file linterp.h.

                                     {  
     grid_type const &grid(m_grid_list[dim]);
     if (x < *(grid.begin())) return -1;
     else if (x >= *(grid.end()-1)) return grid.size()-1;
     else {
       auto i_upper = std::upper_bound(grid.begin(), grid.end(), x);
       return i_upper - grid.begin() - 1;
     }   
   }

◆ get_f_val() [1/2]

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

T NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::get_f_val	(	array< int, N > const &	cell_index,
		array< int, N > const &	v_index
	)		const

inline

Definition at line 185 of file linterp.h.

                                                                                  {
     array<int,N> f_index;
     
     if (m_bContinuous) {      
       for (int i=0; i<N; i++) {
         if (cell_index[i] < 0) {
           f_index[i] = 0;         
         } else if (cell_index[i] >= m_grid_list[i].size()-1) {
           f_index[i] = m_grid_list[i].size()-1;       
         } else {
           f_index[i] = cell_index[i] + v_index[i];        
         }
       }
     } else {
       for (int i=0; i<N; i++) {
         if (cell_index[i] < 0) {
           f_index[i] = 0;
         } else if (cell_index[i] >= m_grid_list[i].size()-1) {
           f_index[i] = (2*m_grid_list[i].size())-1;
         } else {
           f_index[i] = 1 + (2*cell_index[i]) + v_index[i];
         }
       }
     }
     return (*m_pF)(f_index);
   }

◆ get_f_val() [2/2]

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

T NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::get_f_val	(	array< int, N > const &	cell_index,
		int	v
	)		const

inline

Definition at line 212 of file linterp.h.

                                                            {
     array<int,N> v_index;
     for (int dim=0; dim<N; dim++) {
       v_index[dim] = (v >> (N-dim-1)) & 1;                      // test if the i-th bit is set
     }
     return get_f_val(cell_index, v_index);
   } 

◆ init()

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

template<class IterT1 , class IterT2 , class IterT3 >

void NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::init	(	IterT1	grids_begin,
		IterT2	grids_len_begin,
		IterT3	f_begin,
		IterT3	f_end
	)

inline

Definition at line 107 of file linterp.h.

                                                                                       {    
     set_grids(grids_begin, grids_len_begin, m_bCopyData);
     set_f_array(f_begin, f_end, m_bCopyData);
   }  

◆ init_refcount()

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

template<class IterT1 , class IterT2 , class IterT3 , class RefCountIterT >

void NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::init_refcount	(	IterT1	grids_begin,
		IterT2	grids_len_begin,
		IterT3	f_begin,
		IterT3	f_end,
		ArrayRefCountT &	refF,
		RefCountIterT	grid_refs_begin
	)

inline

Definition at line 112 of file linterp.h.

                                                                                                                                                     {       
     set_grids(grids_begin, grids_len_begin, m_bCopyData);
     set_grids_refcount(grid_refs_begin, grid_refs_begin + N);
     set_f_array(f_begin, f_end, m_bCopyData);
     set_f_refcount(refF);
   } 

◆ set_f_array()

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

template<class IterT >

void NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::set_f_array	(	IterT	f_begin,
		IterT	f_end,
		bool	bCopy
	)

inline

Definition at line 144 of file linterp.h.

                                                            {
     unsigned int nGridPoints = 1;
     array<int,N> sizes;
     for (unsigned int i=0; i<m_grid_list.size(); i++) {
       sizes[i] = m_grid_list[i].size();
       nGridPoints *= sizes[i];
     }
 
     int f_len = f_end - f_begin;
     if ( (m_bContinuous && f_len != nGridPoints) || (!m_bContinuous && f_len != 2 * nGridPoints) ) {
       throw std::invalid_argument("f has wrong size");
     }
     for (unsigned int i=0; i<m_grid_list.size(); i++) {
       if (!m_bContinuous) { sizes[i] *= 2; }      
     }
 
     m_F_copy.clear();
     if (bCopy == false) {
       m_pF.reset(new array_type(f_begin, sizes));
     } else {
       m_F_copy = vector<T>(f_begin, f_end);
       m_pF.reset(new array_type(&m_F_copy[0], sizes));
     }
   }  

◆ set_f_refcount()

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

void NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::set_f_refcount ( ArrayRefCountT & refF )

inline

Definition at line 168 of file linterp.h.

                                             {    
     m_ref_F = refF;
   }

◆ set_grids()

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

template<class IterT1 , class IterT2 >

void NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::set_grids	(	IterT1	grids_begin,
		IterT2	grids_len_begin,
		bool	bCopy
	)

inline

Definition at line 120 of file linterp.h.

                                                                          {
     m_grid_list.clear();
     m_grid_ref_list.clear();
     m_grid_copy_list.clear();
     for (int i=0; i<N; i++) {
       int gridLength = grids_len_begin[i];
       if (bCopy == false) {     
         T const *grid_ptr = &(*grids_begin[i]);
         m_grid_list.push_back(grid_type(gridLength, (T*) grid_ptr ));                   // use the given pointer
       } else {
         m_grid_copy_list.push_back( vector<T>(grids_begin[i], grids_begin[i] + grids_len_begin[i]) );   // make our own copy of the grid
         T *begin = &(m_grid_copy_list[i][0]);
         m_grid_list.push_back(grid_type(gridLength, begin));                            // use our copy
       }
     }
   }    

◆ set_grids_refcount()

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

template<class IterT1 , class RefCountIterT >

void NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::set_grids_refcount	(	RefCountIterT	refs_begin,
		RefCountIterT	refs_end
	)

inline

Definition at line 137 of file linterp.h.

                                                                             {
     assert(refs_end - refs_begin == N); 
     m_grid_ref_list.assign(refs_begin, refs_begin + N);
   } 

Member Data Documentation

◆ m_bContinuous

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

const bool NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::m_bContinuous = Continuous

static

Definition at line 79 of file linterp.h.

◆ m_bCopyData

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

const bool NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::m_bCopyData = CopyData

static

Definition at line 78 of file linterp.h.

◆ m_F_copy

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

vector<T> NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::m_F_copy

Definition at line 87 of file linterp.h.

◆ m_grid_copy_list

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

vector<vector<T> > NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::m_grid_copy_list

Definition at line 91 of file linterp.h.

◆ m_grid_list

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

vector<grid_type> NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::m_grid_list

Definition at line 89 of file linterp.h.

◆ m_grid_ref_list

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

vector<GridRefCountT> NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::m_grid_ref_list

Definition at line 90 of file linterp.h.

◆ m_N

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

const int NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::m_N = N

static

Definition at line 77 of file linterp.h.

◆ m_pF

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

array_type_ptr NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::m_pF

Definition at line 85 of file linterp.h.

◆ m_ref_F

template<int N, class T , bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>

ArrayRefCountT NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >::m_ref_F

Definition at line 86 of file linterp.h.

The documentation for this class was generated from the following file:

lib/private/linterp.h

Public Types

Public Member Functions

Public Attributes

Static Public Attributes

Detailed Description

template<int N, class T, bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass> class NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >

Member Typedef Documentation

◆ array_ref_count_type

◆ array_type

◆ array_type_ptr

◆ grid_ref_count_type

◆ grid_type

◆ value_type

Constructor & Destructor Documentation

◆ NDInterpolator() [1/2]

◆ NDInterpolator() [2/2]

Member Function Documentation

◆ find_cell()

◆ get_f_val() [1/2]

◆ get_f_val() [2/2]

◆ init()

◆ init_refcount()

◆ set_f_array()

◆ set_f_refcount()

◆ set_grids()

◆ set_grids_refcount()

Member Data Documentation

◆ m_bContinuous

◆ m_bCopyData

◆ m_F_copy

◆ m_grid_copy_list

◆ m_grid_list

◆ m_grid_ref_list

◆ m_N

◆ m_pF

◆ m_ref_F

template<int N, class T, bool CopyData = true, bool Continuous = true, class ArrayRefCountT = EmptyClass, class GridRefCountT = EmptyClass>
class NDInterpolator< N, T, CopyData, Continuous, ArrayRefCountT, GridRefCountT >