Program.cs

Go to the documentation of this file.

namespace SpikingNeuronNetwork.Console
{
    using Lib;
    using Lib.NeuronModels;
    using System;
    using System.Collections.Generic;
    using System.Linq;

    public class Program
    {
        static void Main()
        {
            const int numInputs = 4;
            const int maxNumOutputSpikes = 3;
            var inputSpikeTimes = new List<Spike>
                {
                    new Spike { Time = 2, NeuronIndex = 0}, 
                    new Spike { Time = 1, NeuronIndex = 2}, 
                    new Spike { Time = 1.5, NeuronIndex = 2}, 
                    new Spike { Time = 4, NeuronIndex = 3}, 
                };
            
            // Construct theta neuron object with randomized weights
            var thetaNeuron = new ThetaNeuron(numInputs, 0.01);

            // Display Neuron Properties
            Console.WriteLine(thetaNeuron);

            // Simulate Theta Neuron Numerically
            thetaNeuron.Method = SimulationMethod.Numerical;
            var outputSpikeTimes = thetaNeuron.RunThetaNeuron(inputSpikeTimes, maxNumOutputSpikes).ToList();

            //Simulate Theta Neuron with Event-Driven Methodology
            thetaNeuron.ResetState();
            thetaNeuron.Method = SimulationMethod.EventDriven;
            var outputSpikeTimes2 = thetaNeuron.RunThetaNeuron(inputSpikeTimes, maxNumOutputSpikes).ToList();

            //Find Mean-Squared in Spike Times Between the Different Methods
            if (outputSpikeTimes.Count != outputSpikeTimes2.Count) 
            {
                Console.WriteLine("The Two Simulation Methods Produced Different Numbers of Output Spikes!");
                Console.WriteLine();
            }
            else if (outputSpikeTimes.Count == 0) 
            {
                Console.WriteLine("The Two Simulation Methods Produced No Output Spikes!");
                Console.WriteLine();
            }
            else 
            {
                double meanSquaredError = 0;
                meanSquaredError += outputSpikeTimes.Select((t, j) => Math.Pow(t.Time - outputSpikeTimes2[j].Time, 2)).Sum();
                meanSquaredError /= outputSpikeTimes.Count;
                Console.WriteLine("Output Spike Mean-Squared Error Between Simulation Methods: ");
                Console.WriteLine("\t" + meanSquaredError);
                Console.WriteLine();
            }

            //Stack two neurons sequentially
            if (outputSpikeTimes.Count > 0)
            {
                Console.WriteLine("A second neuron is simulated in connection to the output of the previous neuron");
                var thetaNeuron2 = new ThetaNeuron(1);
                var outputSpikeTimesList = outputSpikeTimes.Select((t, j) => new Spike{Time = t.Time, NeuronIndex = j}).ToList();
                thetaNeuron2.RunThetaNeuron(outputSpikeTimesList, maxNumOutputSpikes);
            }

            Console.ReadKey();
        }
    }
}