Fixes for RNG

Fig13A.hoc

@@ -175,6 +175,7 @@ assign_synapses(0.01)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"

Fig13B.hoc

@@ -176,6 +176,7 @@ assign_synapses(0.005)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"

Fig5A.hoc

@@ -175,6 +175,7 @@ assign_synapses(1)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"

Fig5B.hoc

@@ -177,6 +177,7 @@ assign_synapses(0.5)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"

Fig6A.hoc

@@ -209,6 +209,7 @@ assign_synapses(0.5)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"

Fig6B.hoc

@@ -209,6 +209,7 @@ assign_synapses(0.333)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"

Fig6C.hoc

@@ -214,6 +214,7 @@ assign_synapses(0.25)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"

Fig6D.hoc

@@ -209,6 +209,7 @@ assign_synapses(0.2)			// assign value (in uS) to each synapse
 objectvar RG,Pulse[ncells]			// create object variables
 
 RG = new Random()				// create random generator
+RG.Random123()
 
 for (i=0; i<ncells; i=i+1) {
 	print "defining current pulse in RE[",i,"]"

README

@@ -122,3 +122,4 @@ ModelDB-Note: This model updated for NEURON version 4.3.1 and higher.
 10-2007: AD, better synaptic mechanisms using counters
 20120106 changed mod files for better solve methods (avoid euler)
 as per web page http://www.neuron.yale.edu/phpbb/viewtopic.php?f=28&t=592
+2024-11: Updated hoc files to use Random123