830 - Haplotyper::SampleOneHaplotype()

float Haplotyper::SampleOneHaplotype(float * vector, double first, double & second, float & choice)
   {
   for (second = 0; second < states; second++)
      {
      choice += vector[CalculateIndex(first, second)] * 0.5;
      if (choice >= 0.0) return;
      }

   second = first;
   }