summaryrefslogblamecommitdiff
path: root/src/math/logb.c
blob: 7f8bdfae19f8ca0a85b611d38a73695ccf588e43 (plain) (tree)
1
2
3
4
5
6
7
8
9
                 


              
                                             

                          

  


                         
                             
                   
                                
                        
 
#include <math.h>

/*
special cases:
	logb(+-0) = -inf, and raise divbyzero
	logb(+-inf) = +inf
	logb(nan) = nan
*/

double logb(double x)
{
	if (!isfinite(x))
		return x * x;
	if (x == 0)
		return -1/(x*x);
	return ilogb(x);
}