reverted to Unity 2019. Doing so prevented ISynapse errors. Program runs but voltage is locked at 50m/V. Will look into this later.

Author: Carlos Eckert

Assets/Save.meta

@@ -32,6 +32,9 @@ public class GameManager : MonoBehaviour
         public double[] Npre;
         public double[] Hpre;
 
+        // Dictionary to store saved states
+        private Dictionary<string, double[]> savedStates = new Dictionary<string, double[]>();
+
         // for loading a file
         private bool loading = false;
         public bool Loading
@@ -204,5 +207,32 @@ private void OnApplicationPause(bool pause)
         {
             isRunning = !pause;
         }
+        /// <summary>
+        /// Saves the state data to the dictionary for later retrieval.
+        /// </summary>
+        /// <param name="variableName">The name of the variable being saved.</param>
+        /// <param name="data">The data to be saved.</param>
+        public void SaveState(string variableName, double[] data)
+        {
+            savedStates[variableName] = data;
+        }
+
+        /// <summary>
+        /// Retrieves the saved state data for a specific variable.
+        /// </summary>
+        /// <param name="variableName">The name of the variable to retrieve.</param>
+        /// <returns>The saved state data, or an empty array if the state does not exist.</returns>
+        public double[] GetSavedState(string variableName)
+        {
+            if (savedStates.ContainsKey(variableName))
+            {
+                return savedStates[variableName];
+            }
+            else
+            {
+                Debug.LogError($"No saved state found for {variableName}");
+                return new double[0];
+            }
+        }
     }
 }

Assets/Scripts/C2M2/GameManager.cs

@@ -0,0 +1,60 @@
+using System;
+using System.Collections.Generic;
+using Vector = MathNet.Numerics.LinearAlgebra.Vector<double>;
+namespace C2M2.NeuronalDynamics.Simulation
+{
+
+    public class GatingVariable
+    {
+        // Name of the gating variable (i.e: n, m, h)
+        public string Name { get; }
+
+        // Alpha and Beta functions for the gating variable
+        public Func<Vector, Vector> Alpha { get; }
+        public Func<Vector, Vector> Beta { get; }
+
+        public double Exponent { get; }
+
+        // Probability is the initial state probability for the gating variable
+        public double Probability { get; set; }
+        public Vector CurrentState { get; set; }
+        public Vector PreviousState { get; set; }
+
+        public GatingVariable(string name, Func<Vector, Vector> alpha, Func<Vector, Vector> beta, double exponent, double probability, int nodeCount)
+        {
+            Name = name;
+            Alpha = alpha;
+            Beta = beta;
+            Exponent = exponent;
+            Probability = probability;
+            CurrentState = Vector.Build.Dense(nodeCount, probability);
+            PreviousState = CurrentState.Clone();
+        }
+        
+    }
+
+    public class IonChannel
+    {
+        public string Name { get; set; }                  // Name of the ion channel (e.g., K+ channel)
+        public double Conductance { get; set; }           // g (conductance)
+        public double ReversalPotential { get; set; }     // E (reversal potential)
+        public List<GatingVariable> GatingVariables { get; set; } // List of gating variables (can be dynamic)
+
+        public IonChannel(string name, double conductance, double reversalPotential)
+        {
+            Name = name;
+            Conductance = conductance;
+            ReversalPotential = reversalPotential;
+            GatingVariables = new List<GatingVariable>();
+        }
+        public void AddGatingVariable(GatingVariable gatingVariable)
+        {
+            GatingVariables.Add(gatingVariable);
+        }
+
+        public void RemoveGatingVariable(string name)
+        {
+            GatingVariables.RemoveAll(gv => gv.Name == name);
+        }
+    }
+}