1#include <math.h>
 2
 3float logbf(float x)
 4{
 5	if (!isfinite(x))
 6		return x * x;
 7	if (x == 0)
 8		return -1/(x*x);
 9	return ilogbf(x);
10}