 #include "CUDAErrorCheckHelpers.h"
 #include "RecordingElectrodes.h"
 // Silences the printfs
-//#define QUIETSTART
-
-
-// CUDA Spiking Simulator Run
-//	INPUT (incomplete):
-//			total_number_of_neurons = Total number of neurons in simulation
-//			total_number_of_connections = Total number of Synapses
-//			presynaptic_neuron_indices = vector- entry corresponds to pre-syn neuron ID
-//			postsynaptic_neuron_indices = vector- entry corresponds to post-syn neuron ID
-//			delays = vector- entry corresponds to axonal delay of syn
-//			weight = vector- entry correponds to synaptic weight
-//			spikes = vector- entry corresponds to countdown of spike on syn
-//			stdp = vector- entry corresponding to whether or not syn is stdp
+// #define QUIETSTART
+
+
 //			lastactive = vector- indicating last time synapse emitted current
-//			parameters=
-//			w_max = stdp parameter
-//			a_minus = stdp parameter
-//			a_plus = stdp parameter
-//			tau_minus = stdp parameter
-//			tau_plus = stdp parameter
-//			timestep = timestep of the simulation
-//			total_time_per_epoch = total time for which the sim should run per epoch
 void GPUDeviceComputation (
 					Neurons * neurons,
 					Connections * connections,

				// Current simulation timestep
				current_time_in_seconds = float(timestep_index)*float(timestep);
				// Start by resetting all the things
				CudaSafeCall(cudaMemset(currentinjection, 0.0f, total_number_of_neurons*sizeof(float)));	
				// If there are poisson populations
				if (numPoisson > 0) {
					// First create the set of random numbers of poisson neurons
					randoms<<<vectorblocksPerGrid, threadsPerBlock>>>(states, gpu_randfloats, total_number_of_neurons);
					CudaCheckError();

					// Update Poisson neuron states

			}
			// Reset the variables necessary
			// CAN GO INTO CLASSES EVENTUALLY
			CudaSafeCall(cudaMemcpy(neurons->d_neuron_variables, neurons->neuron_variables, sizeof(float)*neurons->total_number_of_neurons, cudaMemcpyHostToDevice));
			CudaSafeCall(cudaMemset(neurons->d_lastspiketime, -1000.0f, neurons->total_number_of_neurons*sizeof(float)));
			CudaSafeCall(cudaMemset(connections->d_spikes, 0, sizeof(int)*total_number_of_connections));
			CudaSafeCall(cudaMemset(connections->d_lastactive, -1000.0f, sizeof(float)*total_number_of_connections));
			CudaSafeCall(cudaMemset(connections->d_spikebuffer, -1, total_number_of_connections*sizeof(int)));
