[Avida-SVN] r3192 - branches/hjg-dev/source/cpu

[hjg at myxo.css.msu.edu]

Author: hjg
Date: 2009-03-19 09:52:03 -0400 (Thu, 19 Mar 2009)
New Revision: 3192

Modified:
   branches/hjg-dev/source/cpu/cHardwareCPU.cc
Log:


Modified: branches/hjg-dev/source/cpu/cHardwareCPU.cc
===================================================================
--- branches/hjg-dev/source/cpu/cHardwareCPU.cc	2009-03-17 12:24:18 UTC (rev 3191)
+++ branches/hjg-dev/source/cpu/cHardwareCPU.cc	2009-03-19 13:52:03 UTC (rev 3192)
@@ -3460,6 +3460,30 @@
   if (neighbor != NULL) {
     DoDonate(neighbor);
     
+		// Code to track the edit distance between kin donors and recipients
+		const int edit_dist = cGenomeUtil::FindEditDistance(organism->GetGenome(),neighbor->GetGenome());
+		
+		static ofstream rand_file("rand_dists.dat");
+		static int num_rand_donates = 0;
+		static int num_rand_donates_15_dist = 0;
+		static int tot_dist_rand_donate = 0;
+		
+		num_rand_donates++;
+		if (edit_dist > 15) num_rand_donates_15_dist++;
+		tot_dist_rand_donate += edit_dist;
+		
+		if (num_rand_donates == 1000) {
+			rand_file << num_rand_donates << " "
+			<< (double) num_rand_donates_15_dist / (double) num_rand_donates << " "
+			<< (double) tot_dist_rand_donate / (double) num_rand_donates << endl;
+			
+			num_rand_donates = 0;
+			num_rand_donates_15_dist = 0;
+			tot_dist_rand_donate = 0;
+		}
+		
+		
+		
     //print out how often random donations go to kin
     /*
     static ofstream kinDistanceFile("kinDistance.dat");