Goncalo at 18:39, 25 September 2013

2013-09-25T18:39:02Z

← Older revision		Revision as of 14:39, 25 September 2013
Line 7:		Line 7:
	else		else
	{		{
	~~double flipRate = r * empiricalFlipRate;~~
	double sum = 0.0;		double sum = 0.0;

Line 13:		Line 12:
	sum += from[i];		sum += from[i];

	sum = r (1.0 - empiricalFlipRate) / states;		sum *= r / states;

	double complement = 1. - r;		double complement = 1. - r;
Line 21:		Line 20:
	{		{
	sum *= 1e15;		sum *= 1e15;
	~~flipRate *= 1e15;~~
	complement *= 1e15;		complement *= 1e15;
	}		}

	~~// printf("r = %g, SUM = %g, COMPLEMENT = %g\n", r, sum, complement);~~

	for (int i = 0; i < states; i++)		for (int i = 0; i < states; i++)
	to[i] = from[i] * complement ~~+ from[i^1] * flipRate~~ + sum;		to[i] = from[i] * complement + sum;
	}		}
	}		}
	</source>		</source>

Goncalo: Created page with " void MarkovModel::Transpose(float * from, float * to, double r) { if (r == 0) for (int i = 0; i < states; i++) to[i] = from[i]; els..."

2013-09-25T18:37:07Z

Created page with "<source lang="cpp"> void MarkovModel::Transpose(float * from, float * to, double r) { if (r == 0) for (int i = 0; i < states; i++) to[i] = from[i]; els..."

New page

<source lang="cpp">
void MarkovModel::Transpose(float * from, float * to, double r)
{
if (r == 0)
for (int i = 0; i < states; i++)
to[i] = from[i];
else
{
double flipRate = r * empiricalFlipRate;
double sum = 0.0;

for (int i = 0; i < states; i++)
sum += from[i];

sum *= r * (1.0 - empiricalFlipRate) / states;

double complement = 1. - r;

// avoid underflows
if (sum < 1e-10)
{
sum *= 1e15;
flipRate *= 1e15;
complement *= 1e15;
}

// printf("r = %g, SUM = %g, COMPLEMENT = %g\n", r, sum, complement);

for (int i = 0; i < states; i++)
to[i] = from[i] * complement + from[i^1] * flipRate + sum;
}
}
</source>

830 - MarkovModel::Transpose() - Revision history

Goncalo at 18:39, 25 September 2013

Goncalo: Created page with " void MarkovModel::Transpose(float * from, float * to, double r) { if (r == 0) for (int i = 0; i < states; i++) to[i] = from[i]; els..."