ksolve/Gsolve.h

Go to the documentation of this file.

/**********************************************************************
** This program is part of 'MOOSE', the
** Messaging Object Oriented Simulation Environment.
**           Copyright (C) 2003-2014 Upinder S. Bhalla. and NCBS
** It is made available under the terms of the
** GNU Lesser General Public License version 2.1
** See the file COPYING.LIB for the full notice.
**********************************************************************/

#ifndef _GSOLVE_H
#define _GSOLVE_H

class Stoich;
class Gsolve: public ZombiePoolInterface
{
    public: 
        Gsolve();
        ~Gsolve();

        // Field assignment stuff
        Id getStoich() const;
        void setStoich( Id stoich ); 
        Id getCompartment() const;
        void setCompartment( Id compt );

        unsigned int getNumLocalVoxels() const;
        unsigned int getNumAllVoxels() const;
        void setNumAllVoxels( unsigned int num );

        void setNumPools( unsigned int num ); 
        unsigned int getNumPools() const; 
        VoxelPoolsBase* pools( unsigned int i ); 
        double volume( unsigned int i ) const;

        vector< double > getNvec( unsigned int voxel) const;
        void setNvec( unsigned int voxel, vector< double > vec );
        // Dest Finfos
        void process( const Eref& e, ProcPtr p );
        void reinit( const Eref& e, ProcPtr p );
        void initProc( const Eref& e, ProcPtr p );
        void initReinit( const Eref& e, ProcPtr p );

        void updateVoxelVol( vector< double > vols );

        // Solver setup functions
        void rebuildGssaSystem();
        void fillMmEnzDep();
        void fillMathDep();
        void insertMathDepReacs( unsigned int mathDepIndex,
            unsigned int firedReac );
        void makeReacDepsUnique();

        // Solver interface functions
        unsigned int getPoolIndex( const Eref& e ) const;
        unsigned int getVoxelIndex( const Eref& e ) const;

        void setDsolve( Id dsolve );
        
        // ZombiePoolInterface inherited functions

        void setN( const Eref& e, double v );
        double getN( const Eref& e ) const;
        void setNinit( const Eref& e, double v );
        double getNinit( const Eref& e ) const;
        void setDiffConst( const Eref& e, double v );
        double getDiffConst( const Eref& e ) const;

        void getBlock( vector< double >& values ) const;
        void setBlock( const vector< double >& values );

        void updateRateTerms( unsigned int index );


        bool getRandInit() const;
        void setRandInit( bool val );

        static SrcFinfo2< Id, vector< double > >* xComptOut();
        static const Cinfo* initCinfo();
    private:
        GssaSystem sys_;
        vector< GssaVoxelPools > pools_;

        unsigned int startVoxel_;

        Stoich* stoichPtr_;

        Id dsolve_;

        ZombiePoolInterface* dsolvePtr_;
};

#endif  // _GSOLVE_H

---