1*3e12c5d1SDavid du Colombier #include <errno.h> 2*3e12c5d1SDavid du Colombier #include <float.h> 3*3e12c5d1SDavid du Colombier #include <math.h> 4*3e12c5d1SDavid du Colombier 5*3e12c5d1SDavid du Colombier static long tab[] = 6*3e12c5d1SDavid du Colombier { 7*3e12c5d1SDavid du Colombier 1L, 10L, 100L, 1000L, 10000L, 8*3e12c5d1SDavid du Colombier 100000L, 1000000L, 10000000L 9*3e12c5d1SDavid du Colombier }; 10*3e12c5d1SDavid du Colombier 11*3e12c5d1SDavid du Colombier double pow10(int n)12*3e12c5d1SDavid du Colombierpow10(int n) 13*3e12c5d1SDavid du Colombier { 14*3e12c5d1SDavid du Colombier int m; 15*3e12c5d1SDavid du Colombier 16*3e12c5d1SDavid du Colombier if(n > DBL_MAX_10_EXP){ 17*3e12c5d1SDavid du Colombier errno = ERANGE; 18*3e12c5d1SDavid du Colombier return HUGE_VAL; 19*3e12c5d1SDavid du Colombier } 20*3e12c5d1SDavid du Colombier if(n < DBL_MIN_10_EXP){ 21*3e12c5d1SDavid du Colombier errno = ERANGE; 22*3e12c5d1SDavid du Colombier return 0.0; 23*3e12c5d1SDavid du Colombier } 24*3e12c5d1SDavid du Colombier if(n < 0) 25*3e12c5d1SDavid du Colombier return 1/pow10(-n); 26*3e12c5d1SDavid du Colombier if(n < sizeof(tab)/sizeof(tab[0])) 27*3e12c5d1SDavid du Colombier return tab[n]; 28*3e12c5d1SDavid du Colombier m = n/2; 29*3e12c5d1SDavid du Colombier return pow10(m) * pow10(n-m); 30*3e12c5d1SDavid du Colombier } 31