1,555
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search830 - Haplotyper::SampleHaplotypes() (view source)
Revision as of 14:10, 2 October 2013
, 14:10, 2 October 2013no edit summary
<source lang="cpp">
<source lang="cpp">
float Haplotyper::SampleHaplotypes(float * vector, int & first, int & second)
float Haplotyper::SampleHaplotypes(float * vector, int & first, int & second, Random & rand)
{
{
float sum = 0.0;
float sum = 0.0;
float *probability = vector;
// Calculate sum over all states
// Calculate sum over all states
sum = 0.0;
sum = 0.0;
for (probability = leftMatrices[markers - 1]; first < states; first++)
for (probability = vector, first = 0; first < states; first++)
{
{
for (second = 0; second <= first; second++)
for (second = 0; second <= first; second++)
