---

viongate.h

#ifndef __VIONGATE_H__
#define __VIONGATE_H__

#include <math.h>
#include <stdlib.h>
#include "advanceable.h"
#include "membranepatchsimple.h"
#include "iongate.h"

#define VIONGATE_VM_MIN -0.100     // -100mV

#define VIONGATE_VM_MAX +0.100     // +100mV

#define VIONGATE_VM_INC +0.00002   // 0.05mV

#define VIONGATE_TABLE_SIZE      ((int)((VIONGATE_VM_MAX - VIONGATE_VM_MIN) / VIONGATE_VM_INC + 1))



class VIonGate : public IonGate {

public:
  VIonGate(void) { Vm=0; Vresting=0; VmScale=0; };

  // NOTE: The destructor comes within the macro IONGATE_TABLES

  virtual void reset(void);

  virtual int updateInternal(void);

  virtual int advance(void);

  virtual double alpha(double V) { return (infty(V))/tau(V); }

  virtual double beta(double V) { return (1.0-infty(V))/tau(V); }

  virtual double tau(double V) { return  1.0/(beta(V)+alpha(V)); }

  virtual double infty(double V) { return alpha(V)/(beta(V)+alpha(V)); }

  virtual double pInfty(MembranePatchSimple *m) { return infty(m->Vm);}

  virtual void ConnectToMembrane(MembranePatchSimple *m) { this->Vm = &(m->Vm); this->Vresting = &(m->Vresting); this->VmScale = &(m->VmScale);}

protected:

  IONGATE_TABLES(VIonGate);

  float  *Vresting;

  double *Vm;

  float *VmScale;
};

#endif

---