83 for (
size_t i = 0; i < record.
mSize; i++)
85 std::cout << record.
mpValues[i] <<
" ";
88 std::cout << std::endl;
102 mpInterpolationState(NULL),
103 mShiftBeginIndex(false)
110 mBeginRecordIndex(0),
112 mShiftBeginIndex(false)
119 for (
size_t i = 0; i < stateNum; i++)
168 for (
size_t i = 1; i < size; i++, y++)
193 std::cout <<
"x " << x <<
" is out of range" << std::endl;
211 tmpNumer *= x -
mpX[i];
217 numer = tmpNumer / (x -
mpX[i]);
226 y +=
mpY[i] * numer / denom;
277 for (
size_t species = 1; species <= speciesNum; species++)
280 mpY[j] = ppArray[j][species];
296 std::cout <<
"State Num: " << interp.
mStateNum << std::endl;
297 std::cout <<
"Record Num: " << interp.
mHaveRecordNum << std::endl;
299 std::cout <<
"Shift or not: ";
302 std::cout <<
"Yes" << std::endl;
304 std::cout <<
"No" << std::endl;
310 std::cout <<
"Record #:" << i << std::endl;
315 std::cout <<
"Finish, aha" << std::endl;
CStateRecord * mpInterpolationState
CStateRecord & operator=(const CStateRecord &rhs)
void printInterpolation(const CInterpolation &interp)
void recordState(const C_FLOAT64 time, const C_FLOAT64 *y)
const C_FLOAT64 & getSlowPropensitySum() const
void printRecord(const CStateRecord &record)
void calculateFiringTime()
C_FLOAT64 interpolation(C_FLOAT64 x)
C_FLOAT64 * getArray() const
const C_FLOAT64 & getTime() const
CStateRecord * getInterpolationState()
const size_t & getArrayLen() const