biophysics/HHChannel2D.h

Go to the documentation of this file.

#ifndef _HHChannel2D_h
#define _HHChannel2D_h
/**********************************************************************
** This program is part of 'MOOSE', the
** Messaging Object Oriented Simulation Environment,
** also known as GENESIS 3 base code.
**           copyright (C) 2003-2005 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.
*********************************************************************
*/

// Ported to asyn13 on 2014-05-30 by Subhasis Ray

typedef double ( *PFDD )( double, double );
class HHChannel2D: public ChanCommon
{
#ifdef DO_UNIT_TESTS
    friend void testHHChannel2D();
#endif // DO_UNIT_TESTS

    public:
        HHChannel2D();

        // Value field access function definitions
        void setUseConcentration( int value );
        int getUseConcentration( );
        void setXindex( string index );
        string getXindex() const;
        void setYindex( string index );
        string getYindex() const;
        void setZindex( string index );
        string getZindex() const;

        void setXpower( const Eref& e, double Xpower );
        double getXpower( const Eref& e) const;
        void setYpower( const Eref& e, double Ypower );
        double getYpower( const Eref& e ) const;
        void setZpower( const Eref& e, double Zpower );
        double getZpower( const Eref& e ) const;
        void setInstant( int Instant );
        int getInstant() const;
        void setX( double X );
        double getX() const;
        void setY( double Y );
        double getY() const;
        void setZ( double Z );
        double getZ() const;

        HHGate2D* getXgate( unsigned int i );
        HHGate2D* getYgate( unsigned int i );
        HHGate2D* getZgate( unsigned int i );
        void setNumGates( unsigned int num );

        unsigned int getNumXgates() const;
        unsigned int getNumYgates() const;
        unsigned int getNumZgates() const;

        // Dest function definitions

        void vProcess( const Eref& e, ProcPtr p );

        void vReinit( const Eref& e, ProcPtr p );
        void conc1( double conc );
        void conc2( double conc );

        double ( *takeXpower_ )( double, double );
        double ( *takeYpower_ )( double, double );
        double ( *takeZpower_ )( double, double );

        void createGate( const Eref& e, string gateType );

        void innerCreateGate(
             const string& gateName,
            HHGate2D** gatePtr, Id chanId, Id gateId );

        bool checkOriginal( Id chanId ) const;

        void destroyGate( const Eref& e, string gateType );

        void innerDestroyGate( const string& gateName,
            HHGate2D** gatePtr, Id chanId );

        bool setGatePower( const Eref& e, double power,
            double* assignee, const string& gateType );

        static PFDD selectPower( double power);

        static const Cinfo* initCinfo();
        
    private:
        int dependency( string index, unsigned int dim );
        double depValue( int dependency );
        double integrate( double state, double dt, double A, double B );

        double Xpower_; 
        double Ypower_; 
        double Zpower_; 

        int instant_;   
        double X_;   
        double Y_;   
        double Z_;   

        bool xInited_, yInited_, zInited_; 

        double g_;  
        
        double conc1_;
        double conc2_;
        
        string Xindex_;
        string Yindex_;
        string Zindex_;

        int Xdep0_;
        int Xdep1_;
        int Ydep0_;
        int Ydep1_;
        int Zdep0_;
        int Zdep1_;

        HHGate2D* xGate_;
        HHGate2D* yGate_; 
        HHGate2D* zGate_; 

        static const double EPSILON;
        static const int INSTANT_X;
        static const int INSTANT_Y;
        static const int INSTANT_Z;

        static double power1( double x, double p ) {
            return x;
        }
        static double power2( double x, double p ) {
            return x * x;
        }
        static double power3( double x, double p ) {
            return x * x * x;
        }
        static double power4( double x, double p ) {
            return power2( x * x, p );
        }
        static double powerN( double x, double p );
};


#endif // _HHChannel2D_h

---