LookupTable Class Reference

#include <RateLookup.h>

Collaboration diagram for LookupTable:

Public Member Functions
void	addColumns (int species, const vector< double > &C1, const vector< double > &C2)
	Adds the columns for a given species. Columns supplied are C1 and C2. More...
void	column (unsigned int species, LookupColumn &column)
void	lookup (const LookupColumn &column, const LookupRow &row, double &C1, double &C2)
	Actually performs the lookup and the linear interpolation. More...
	LookupTable ()
	LookupTable (double min, double max, unsigned int nDivs, unsigned int nSpecies)
	number of species (no. of columns / 2) More...
void	row (double x, LookupRow &row)

Private Attributes
double	dx_
double	max_
	max of the voltage / caConc range More...
double	min_
	min of the voltage / caConc range More...
unsigned int	nColumns_
	(# columns) = 2 * (# species) More...
unsigned int	nPts_
vector< double >	table_
	Flattened table. More...

Detailed Description

Definition at line 27 of file RateLookup.h.

Constructor & Destructor Documentation

LookupTable::LookupTable ( )

inline

Definition at line 30 of file RateLookup.h.

{ ; }

LookupTable::LookupTable	(	double	min,
		double	max,
		unsigned int	nDivs,
		unsigned int	nSpecies
	)

number of species (no. of columns / 2)

Parameters

min	min of range
max	max of range
nDivs	number of divisions (~ no. of rows)

Definition at line 15 of file RateLookup.cpp.

{
    min_ = min;
    max_ = max;
    // Number of points is 1 more than number of divisions.
    // Then add one more since we may interpolate at the last point in the table.
    nPts_ = nDivs + 1 + 1;
    dx_ = ( max - min ) / nDivs;
    // Every row has 2 entries for each type of gate
    nColumns_ = 2 * nSpecies;

    //~ interpolate_.resize( nSpecies );
    table_.resize( nPts_ * nColumns_ );
}

Member Function Documentation

void LookupTable::addColumns	(	int	species,
		const vector< double > &	C1,
		const vector< double > &	C2
	)

Adds the columns for a given species. Columns supplied are C1 and C2.

Definition at line 31 of file RateLookup.cpp.

Referenced by HSolveActive::createLookupTables().

{
    vector< double >::const_iterator ic1 = C1.begin();
    vector< double >::const_iterator ic2 = C2.begin();
    vector< double >::iterator iTable = table_.begin() + 2 * species;
    // Loop until last but one point
    for ( unsigned int igrid = 0; igrid < nPts_ - 1 ; ++igrid ) {
        *( iTable )     = *ic1;
        *( iTable + 1 ) = *ic2;

        iTable += nColumns_;
        ++ic1, ++ic2;
    }
    // Then duplicate the last point
    *( iTable )     = C1.back();
    *( iTable + 1 ) = C2.back();

    //~ interpolate_[ species ] = interpolate;
}

Here is the caller graph for this function:

void LookupTable::column	(	unsigned int	species,
		LookupColumn &	column
	)

Definition at line 56 of file RateLookup.cpp.

References LookupColumn::column.

Referenced by HSolveActive::createLookupTables().

{
    column.column = 2 * species;
    //~ column.interpolate = interpolate_[ species ];
}

Here is the caller graph for this function:

void LookupTable::lookup	(	const LookupColumn &	column,
		const LookupRow &	row,
		double &	C1,
		double &	C2
	)

Actually performs the lookup and the linear interpolation.

Definition at line 76 of file RateLookup.cpp.

References LookupColumn::column, LookupRow::fraction, and LookupRow::row.

Referenced by HSolveActive::reinitChannels().

{
    double a, b;
    double *ap, *bp;

    ap = row.row + column.column;

    //~ if ( ! column.interpolate ) {
        //~ C1 = *ap;
        //~ C2 = *( ap + 1 );
        //~
        //~ return;
    //~ }

    bp = ap + nColumns_;

    a = *ap;
    b = *bp;
    C1 = a + ( b - a ) * row.fraction;

    a = *( ap + 1 );
    b = *( bp + 1 );
    C2 = a + ( b - a ) * row.fraction;
}

Here is the caller graph for this function:

void LookupTable::row	(	double	x,
		LookupRow &	row
	)

Returns the row corresponding to x in the "row" parameter. i.e., returns the leftover fraction and the row's start address.

Definition at line 62 of file RateLookup.cpp.

References LookupRow::fraction, and LookupRow::row.

Referenced by HSolveActive::reinitChannels().

{
    if ( x < min_ )
        x = min_;
    else if ( x > max_ )
        x = max_;

    double div = ( x - min_ ) / dx_;
    unsigned int integer = ( unsigned int )( div );

    row.fraction = div - integer;
    row.row = &( table_.front() ) + integer * nColumns_;
}

Here is the caller graph for this function:

Member Data Documentation

double LookupTable::dx_

private

This is the smallest difference: (max - min) / nDivs

Definition at line 73 of file RateLookup.h.

double LookupTable::max_

private

max of the voltage / caConc range

Definition at line 69 of file RateLookup.h.

double LookupTable::min_

private

min of the voltage / caConc range

Definition at line 68 of file RateLookup.h.

unsigned int LookupTable::nColumns_

private

(# columns) = 2 * (# species)

Definition at line 75 of file RateLookup.h.

unsigned int LookupTable::nPts_

private

Number of rows in the table. Equal to nDivs + 2, so that interpol. is safe at either end.

Definition at line 70 of file RateLookup.h.

vector< double > LookupTable::table_

private

Flattened table.

Definition at line 67 of file RateLookup.h.

---

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

• moose-core/hsolve/RateLookup.h
• moose-core/hsolve/RateLookup.cpp