commented out debugging routines in container.cpp

scripts/veltest.nsl

@@ -1,8 +1,8 @@
 
 add edmexperiment "Velocity Test Experiment" {
-	set phase_steps 30 #phase average loops
+	set phase_steps 20 #phase average loops
 #	set bounces 10000 #number of bounces to run the simulation over
-	set lifetime 120	#The lifetime of a particle (seconds)
+	set lifetime 20	#The lifetime of a particle (seconds)
 
 	#Set the distance to 'lift' the particle off of the bottles
 	#surface after a collision has happened. This is in to fix
@@ -35,7 +35,7 @@ add edmexperiment "Velocity Test Experiment" {
 	add particle neutron {
 		set gamma -29.1646870E+6  # (Hertz per tesla)
 
-		#vary velocity from 10 to 200 in 20 steps
+		#vary velocity from 1 to 200 in 50 steps
 		#set velocity 3 # (Meters per second)
 		#vary velocity from 100 to 500 in 3 steps
 		set velocity 500

src/container.cpp

@@ -282,6 +282,7 @@ void container::reflect ( vector3& velocity_vec, const vector3& normal, const lo
 			if (normal.z < 0)
 				velocity_vec.z *= -1.;
 		}
+	/*
 		// Dump out the dot product between this and the normal
 		static std::ofstream dotps("dotps.txt");
 		static std::ofstream dotwrong("dotwrong.txt");
@@ -300,7 +301,7 @@ void container::reflect ( vector3& velocity_vec, const vector3& normal, const lo
 		if (dotp > 1)
 			dotwrong << dotp << "\t(" << velocity_vec << ")\t(" << normal*velocity << ")" << endl;
 		reflects << velocity_vec << endl;
-	}
+	}*/
 	}
 	else
 		throw runtime_error("Unrecognised reflection type requested for reflection");