NeuroProject/neurone.cpp
36fc94483816master

neurone.cpp
View Options

			#include "neurone.hpp"
			#include "constant.h"
			#include <vector>
			#include <cmath>

			Neurone::Neurone() :
			membrane_potential_(Standard_potential),rest_time_(0.0) {}


			void Neurone::update(double const Input_current, double time_i) {

			double new_potential(0.0);

			if(not isInRest()) {
			new_potential = exp(-DeltaTime/Tau)*(membrane_potential_);
			new_potential += Input_current(Resistance)(1-exp(-DeltaTime/Tau));
			/// 1) calculate the new potential

			membrane_potential_=new_potential;
			/// 2) update the membrane potential

			Spike(time_i);
			/// 3) make the neurone spike
			}
			}

			double Neurone::getMembranePotential () const {

			return membrane_potential_-0.070;
			/// -0.070 because we use a 0 origin for the membrane potential
			}

			bool Neurone::isInRest() {

			if (not (rest_time_ <= 0.0)) {
			rest_time_-=DeltaTime;
			/// retires an unit of time
			return true;
			}
			return false;
			}

			void Neurone::Spike(double time_i) {

			if(membrane_potential_ >= Spike_Treshold) {
			spikes_times_.push_back(time_i);
			/// stores the spike time
			membrane_potential_=Standard_potential;
			/// reinitialize the membrane potential
			rest_time_=TauR;
			/// set the rest time
			};
			}

			std::vector<double> Neurone::getSpikeTimes() const {

			return spikes_times_;
			}