830 - Haplotyper::SummarizeOptions()

float Haplotyper::SummarizeOptions(float * vector, float ** left, float ** right, float theta, int first, int second)
   {
   right[0][0] = (1. - theta) * (1. - theta);
   right[1][0] = right[0][1] = (1. - theta) * theta / states;
   right[1][1] = theta / states * theta / states;

   left[0][0] = left[1][0] = left[0][1] = left[0][0] = 0.;

   float * probability = vector;

   for (int k = 0; k < states; k++)
      {
      for (int l = 0; l < k; l++, probability++)
         {
         left[1][1] += *probability;
         if (first == k || first == l) left[0][1] += *probability * 0.5;
         if (second == k || second == l) left[1][0] += *probability * 0.5;
         if (first == k && second == l || first == l && second == k) left[0][0] += *probability * 0.5;
         }

      left[1][1] += *probability;
      if (first == k) left[0][1] += *probability;
      if (second == k) left[1][0] += *probability;
      if (first == k && second == l) left[0][0] += *probability;
      }

   return left[0][0] * right[0][0] + left[1][0] * right[1][0] +
          left[0][1] * right[0][1] + left[1][1] * right[1][1];
   }