830 - Haplotyper::ScoreLeftConditional()

From Genome Analysis Wiki
Jump to navigationJump to search
void Haplotyper::ScoreLeftConditional()
   {
   ResetMemoryPool();
   GetMemoryBlock(0);

   SetupPrior(leftMatrices[0]);
   ConditionOnData(leftMatrices[0], 0, current_genotypes[0]);

   double theta = 0.0;
   for (int i = 1; i < markers; i++)
      {
      GetMemoryBlock(i);

      Transpose(leftMatrices[i-1], leftMatrices[i], thetas[i-1]);
      ConditionOnData(leftMatrices[i], i, current_genotypes[i]);
      }
   }