23 #ifndef COPASI_CCSPMethod
24 #define COPASI_CCSPMethod
174 virtual void step(
const double & deltaT);
456 #endif // COPASI_CCSPMethod
void CSPParticipationIndex(C_INT &N, C_INT &M, C_FLOAT64 &tauM1, CMatrix< C_FLOAT64 > &B0)
std::vector< CMatrix< C_FLOAT64 > > mVec_mImportanceIndex
CMatrix< C_FLOAT64 > mAmplitudeTab
CCSPMethod(const CCopasiContainer *pParent=NULL)
const CArrayAnnotation * getSlowParticipationIndexAnn() const
CMatrix< C_FLOAT64 > mParticipationIndexNormedColumn
CArrayAnnotation * pParticipationIndexNormedRowAnn
CArrayAnnotation * pSlowParticipationIndexAnn
CMatrix< C_FLOAT64 > mParticipationIndexTab
std::vector< CMatrix< C_FLOAT64 > > mVec_mFastParticipationIndex
void calculateJacobian(C_INT &n, CVector< C_FLOAT64 > &y, CMatrix< C_FLOAT64 > &J)
void findCandidatesNumber(C_INT &n, C_INT &k, CVector< C_FLOAT64 > &tsc, C_INT &info)
const CArrayAnnotation * getImportanceIndexAnn() const
CArrayAnnotation * pTmp5NormedRow
void CSPOutput(C_INT &N, C_INT &M, C_INT &R)
virtual void step(const double &deltaT)
void findTimeScaleSeparation(C_INT &n, C_INT &k, CVector< C_FLOAT64 > &tsc, C_INT &info)
std::vector< CMatrix< C_FLOAT64 > > mVec_mRadicalPointer
void sminverse(C_INT &n, CMatrix< C_FLOAT64 > &A, CMatrix< C_FLOAT64 > &B)
static CTSSAMethod * createMethod(CCopasiMethod::SubType subType=CCopasiMethod::unset)
std::vector< CMatrix< C_FLOAT64 > > mVec_mSlowParticipationIndex
bool modesAreExhausted(C_INT &N, C_INT &M, C_FLOAT64 &tauM, C_FLOAT64 &tauM1, CVector< C_FLOAT64 > &g, CMatrix< C_FLOAT64 > &A, CMatrix< C_FLOAT64 > &B, CMatrix< C_FLOAT64 > &F)
const CArrayAnnotation * getFastParticipationIndexAnn() const
CMatrix< C_FLOAT64 > mFastParticipationIndexTab
void CSPImportanceIndex(C_INT &N, C_FLOAT64 &tauM1, CMatrix< C_FLOAT64 > &Q)
void CSPradicalPointer(C_INT &N, C_INT &M, CMatrix< C_FLOAT64 > &A, CMatrix< C_FLOAT64 > &B)
CMatrix< C_FLOAT64 > mImportanceIndexNormedRow
void smadd(CMatrix< C_FLOAT64 > &A, CMatrix< C_FLOAT64 > &B, CMatrix< C_FLOAT64 > &C, C_INT &n1, C_INT &n2)
CMatrix< C_FLOAT64 > mFastReactionPointerNormedTab
std::vector< CMatrix< C_FLOAT64 > > mVec_mParticipationIndexNormedRow
void yCorrection(C_INT &N, C_INT &M, CVector< C_FLOAT64 > &y, CMatrix< C_FLOAT64 > &TAUM, CMatrix< C_FLOAT64 > &F, CMatrix< C_FLOAT64 > &A)
CMatrix< C_FLOAT64 > mImportanceIndex
const CArrayAnnotation * getFastReactionPointerNormedAnn() const
CVector< C_FLOAT64 > mFastParticipationIndex
CArrayAnnotation * pFastReactionPointerNormedAnn
void modesAmplitude(C_INT &N, C_INT &M, CVector< C_FLOAT64 > &g, CMatrix< C_FLOAT64 > &B, CMatrix< C_FLOAT64 > &F)
std::vector< CMatrix< C_FLOAT64 > > mVec_mAmplitude
CArrayAnnotation * pTmp4NormedColumn
CVector< C_FLOAT64 > mSlowParticipationIndex
CMatrix< C_FLOAT64 > mParticipationIndexNormedRowTab
void emptyOutputData(C_INT &N, C_INT &M, C_INT &R)
CArrayAnnotation * pParticipationIndexAnn
CMatrix< C_FLOAT64 > mRadicalPointer
const CArrayAnnotation * getParticipationIndexNormedRowAnn() const
void createAnnotationsM()
const CArrayAnnotation * getImportanceIndexNormedRowAnn() const
CMatrix< C_FLOAT64 > mParticipationIndexNormedRow
std::vector< CMatrix< C_FLOAT64 > > mVec_mParticipationIndex
CMatrix< C_FLOAT64 > mFastReactionPointerNormed
CArrayAnnotation * pImportanceIndexNormedRowAnn
CArrayAnnotation * pTmp4Fast
CVector< C_FLOAT64 > mYerror
void smmult(CMatrix< C_FLOAT64 > &A, CMatrix< C_FLOAT64 > &B, CMatrix< C_FLOAT64 > &C, C_INT &n1, C_INT &n2, C_INT &n3)
CArrayAnnotation * pTmp4NormedRow
CMatrix< C_FLOAT64 > mSlowParticipationIndexTab
CArrayAnnotation * pAmplitudeAnn
const CArrayAnnotation * getParticipationIndexNormedColumnAnn() const
void perturbateA(C_INT &n, CMatrix< C_FLOAT64 > &A, C_FLOAT64 delta)
void initializeParameter()
CArrayAnnotation * pImportanceIndexAnn
CMatrix< C_FLOAT64 > mParticipationIndex
CArrayAnnotation * pTmp3Normed
CMatrix< C_FLOAT64 > mFastReactionPointerTab
CMatrix< C_FLOAT64 > mImportanceIndexTab
CArrayAnnotation * pTmp4Slow
std::vector< CMatrix< C_FLOAT64 > > mVec_mImportanceIndexNormedRow
CArrayAnnotation * pParticipationIndexNormedColumnAnn
CArrayAnnotation * pFastReactionPointerAnn
C_INT isBlockDiagonal(C_INT &N, C_INT &M, CMatrix< C_FLOAT64 > &ALA, C_FLOAT64 SMALL)
const CArrayAnnotation * getFastReactionPointerAnn() const
void smsubst(CMatrix< C_FLOAT64 > &A, CMatrix< C_FLOAT64 > &B, CMatrix< C_FLOAT64 > &C, C_INT &n1, C_INT &n2)
virtual void predifineAnnotation()
std::vector< CMatrix< C_FLOAT64 > > mVec_mFastReactionPointerNormed
virtual bool setAnnotationM(size_t step)
CArrayAnnotation * pRadicalPointerAnn
CMatrix< C_FLOAT64 > mImportanceIndexNormedRowTab
void cspstep(const double &deltaT, C_INT &n, C_INT &m, CMatrix< C_FLOAT64 > &A, CMatrix< C_FLOAT64 > &B)
void setVectors(int fast)
CMatrix< C_FLOAT64 > mRadicalPointerTab
std::vector< CMatrix< C_FLOAT64 > > mVec_mParticipationIndexNormedColumn
CArrayAnnotation * pFastParticipationIndexAnn
CMatrix< C_FLOAT64 > mFastReactionPointer
void smnorm(C_INT &n, CMatrix< C_FLOAT64 > &A, CMatrix< C_FLOAT64 > &B, C_INT &n1)
const CArrayAnnotation * getAmplitudeAnn() const
const CArrayAnnotation * getParticipationIndexAnn() const
void printResult(std::ostream *ostream) const
void basisRefinement(C_INT &N, C_INT &M, CMatrix< C_FLOAT64 > &ALA, CMatrix< C_FLOAT64 > &TAU, CMatrix< C_FLOAT64 > &A, CMatrix< C_FLOAT64 > &B, CMatrix< C_FLOAT64 > &A0, CMatrix< C_FLOAT64 > &B0)
std::vector< CMatrix< C_FLOAT64 > > mVec_mFastReactionPointer
virtual void start(const CState *initialState)
CVector< C_FLOAT64 > mAmplitude
CMatrix< C_FLOAT64 > mParticipationIndexNormedColumnTab
const CArrayAnnotation * getRadicalPointerAnn() const