R/log.space.R

Defines functions com.log.sum

Documented in com.log.sum

com.log.sum = function(x,y)		# log.sum(x,y) = log( exp(x) + exp(y) )
{
	if (x == -Inf)
		{ return (y); }
	else if (y == -Inf)
		{ return (x); }
	else if (x > y)
		{ return (x + log1p( exp(y - x) ) ); }
	else
		{ return (y + log1p( exp(x - y) ) ); }
}
jarad/compoisson documentation built on May 18, 2019, 3:45 p.m.