HSolveActive.h

Go to the documentation of this file.

/**********************************************************************
** This program is part of 'MOOSE', the
** Messaging Object Oriented Simulation Environment.
**   copyright (C) 2003-2007 Upinder S. Bhalla, Niraj Dudani 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 _HSOLVE_ACTIVE_H
#define _HSOLVE_ACTIVE_H

#include "../basecode/header.h"
#include <set>
#include <limits> // Max and min 'double' values needed for lookup table init.
#include "../biophysics/CaConcBase.h"
#include "../biophysics/HHGate.h"
#include "../biophysics/ChanBase.h"
#include "../biophysics/ChanCommon.h"
#include "../biophysics/HHChannelBase.h"
#include "../biophysics/HHChannel.h"
#include "../biophysics/SpikeGen.h"
#include "HSolveUtils.h"
#include "HSolveStruct.h"
#include "HinesMatrix.h"
#include "HSolvePassive.h"
#include "RateLookup.h"

class HSolveActive: public HSolvePassive
{
    typedef vector< CurrentStruct >::iterator currentVecIter;

public:
    HSolveActive();

    void setup( Id seed, double dt );
    void step( ProcPtr info );          
    void reinit( ProcPtr info );

protected:
    int                       caAdvance_;

    double                    vMin_;
    double                    vMax_;
    int                       vDiv_;
    double                    caMin_;
    double                    caMax_;
    int                       caDiv_;

    vector< CurrentStruct >   current_;         
    vector< double >          state_;           
    //~ vector< int >             instant_;
    vector< ChannelStruct >   channel_;         
    vector< SpikeGenStruct >  spikegen_;
    vector< SynChanStruct >   synchan_;
    vector< CaConcStruct >    caConc_;          
    vector< double >          ca_;              
    vector< double >          caActivation_;    
    vector< double* >         caTarget_;        
    LookupTable               vTable_;
    LookupTable               caTable_;
    vector< bool >            gCaDepend_;       
    vector< unsigned int >    caCount_;         
    vector< int >             caDependIndex_;   
    vector< LookupColumn >    column_;          
    vector< LookupRow >       caRowCompt_;      
    vector< LookupRow* >      caRow_;           
    vector< int >             channelCount_;    
    vector< currentVecIter >  currentBoundary_; 
    vector< unsigned int >    chan2compt_;      
    vector< unsigned int >    chan2state_;      
    vector< double >          externalCurrent_; 
    vector< double >          externalCalcium_; 
    vector< Id >              caConcId_;        
    vector< Id >              channelId_;       
    vector< Id >              gateId_;          
    //~ vector< vector< Id > >    externalChannelId_;
    vector< unsigned int >    outVm_;           
    vector< unsigned int >    outCa_;           
     vector< unsigned int >    outIk_;

private:
    void readHHChannels();
    void readGates();
    void readCalcium();
    void readSynapses();
    void readExternalChannels();
    void createLookupTables();
    void manageOutgoingMessages();

    void cleanup();

    void reinitSpikeGens( ProcPtr info );
    void reinitCompartments();
    void reinitCalcium();
    void reinitChannels();

    void calculateChannelCurrents();
    void updateMatrix();
    void forwardEliminate();
    void backwardSubstitute();
    void advanceCalcium();
    void advanceChannels( double dt );
    void advanceSynChans( ProcPtr info );
    void sendSpikes( ProcPtr info );
    void sendValues( ProcPtr info );

    static const int INSTANT_X;
    static const int INSTANT_Y;
    static const int INSTANT_Z;
};

#endif // _HSOLVE_ACTIVE_H