00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039 module mango.format.Double;
00040
00041 private import mango.format.Int,
00042 mango.format.Number;
00043
00044 extern (C) double frexp (double, int*);
00045
00046 version (Posix)
00047 {
00048 extern (C) int strncasecmp (char *, char*, uint);
00049 alias strncasecmp memicmp;
00050 }
00051 else
00052 extern (C) int memicmp (char*, char*, uint);
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063 class Double : Number
00064 {
00065 alias double tValue;
00066
00067
00068
00069
00070
00071 private static tValue pow10 (uint exp)
00072 in {
00073 assert (exp < 512);
00074 }
00075 body
00076 {
00077 static tValue[] Powers =
00078 [
00079 1.0e1,
00080 1.0e2,
00081 1.0e4,
00082 1.0e8,
00083 1.0e16,
00084 1.0e32,
00085 1.0e64,
00086 1.0e128,
00087 1.0e256,
00088 ];
00089
00090 tValue mult = 1.0;
00091 foreach (tValue power; Powers)
00092 {
00093 if (exp & 1)
00094 mult *= power;
00095 if ((exp >>= 1) == 0)
00096 break;
00097 }
00098 return mult;
00099 }
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110 final static tValue parse (tChar[] src, uint* ate=null)
00111 {
00112 char c;
00113 char* p;
00114 int exp;
00115 bool sign;
00116 tValue value = 0.0;
00117
00118
00119 c = *(p = src.ptr + trim (src, sign));
00120
00121
00122
00123 while (c >= '0' && c <= '9')
00124 {
00125 value = value * 10 + (c - '0');
00126 c = *++p;
00127 }
00128
00129
00130 if (c == '.')
00131 c = *++p;
00132
00133
00134
00135
00136
00137
00138
00139
00140 while (c >= '0' && c <= '9')
00141 {
00142 value = value * 10 + (c - '0');
00143 c = *++p;
00144 --exp;
00145 }
00146
00147
00148 if (value)
00149 {
00150
00151 if (c == 'e' || c == 'E')
00152 {
00153 uint len;
00154 exp += Int.parse (src[(p-src.ptr)+1..length], Radix.Decimal, &len);
00155 p += len;
00156 }
00157
00158
00159
00160 if (exp < 0)
00161 value /= pow10 (-exp);
00162 else
00163 value *= pow10 (exp);
00164 }
00165 else
00166
00167 if (p == src.ptr)
00168 if (memicmp (p, "inf", 3) == 0)
00169 p += 3, value = tValue.infinity;
00170 else
00171 if (memicmp (p, "nan", 3) == 0)
00172 p += 3, value = tValue.nan;
00173
00174
00175 if (ate)
00176 *ate = p - src.ptr;
00177 return sign ? -value : value;
00178 }
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193 static final tChar[] format (tChar[] dst, tValue x, uint decimals = 6, bool scientific = false)
00194 in {
00195 assert (dst.length >= 32);
00196 }
00197 body
00198 {
00199
00200
00201 static int toDigit (inout tValue v, inout int count)
00202 {
00203 int digit;
00204
00205
00206 if (++count > 17)
00207 digit = 0;
00208 else
00209 {
00210
00211 digit = cast(int) v;
00212 v = (v - digit) * 10.0;
00213 }
00214 return digit + '0';
00215 }
00216
00217
00218 if (((cast(ushort*) &x)[3] & 0x7ff0) == 0x7ff0)
00219 if (*(cast(ulong*) &x) & 0x000f_ffff_ffff_ffff)
00220 return "nan";
00221 else
00222 return "inf";
00223
00224 int exp;
00225 bool sign;
00226
00227
00228 if (x < 0.0)
00229 {
00230 x = -x;
00231 sign = true;
00232 }
00233
00234
00235 if (x > 0.0)
00236 {
00237
00238 x += 0.5 / pow10 (decimals);
00239
00240
00241 frexp (x, &exp);
00242 exp = cast(int) (0.301029995664 * exp);
00243
00244
00245 int len = exp;
00246 if (exp < 0)
00247 {
00248 --exp;
00249 x *= pow10 (len = -exp);
00250 }
00251 else
00252 x /= pow10 (exp);
00253
00254
00255 if (len + 32 > dst.length)
00256 scientific = true;
00257 }
00258
00259 char*p = dst;
00260 int count = 0;
00261
00262
00263 if (sign)
00264 *p++ = '-';
00265
00266
00267 if (scientific)
00268 {
00269
00270 *p++ = toDigit (x, count);
00271 *p++ = '.';
00272
00273
00274 while (decimals-- > 0)
00275 *p++ = toDigit (x, count);
00276
00277
00278 if (exp)
00279 {
00280 *p++ = 'e';
00281 *p++ = (exp < 0) ? '-' : '+';
00282 if (exp < 0)
00283 exp = -exp;
00284
00285 if (exp >= 100)
00286 {
00287 *p++ = (exp/100) + '0';
00288 exp %= 100;
00289 }
00290
00291 *p++ = (exp/10) + '0';
00292 *p++ = (exp%10) + '0';
00293 }
00294 }
00295 else
00296 {
00297
00298 if (exp < 0)
00299 *p++ = '0';
00300 else
00301
00302 for (; exp >= 0; --exp)
00303 *p++ = toDigit (x, count);
00304
00305
00306 *p++ = '.';
00307
00308
00309 for (++exp; exp < 0 && decimals > 0; --decimals, ++exp)
00310 *p++ = '0';
00311
00312
00313
00314 while (decimals-- > 0)
00315 *p++ = toDigit (x, count);
00316 }
00317
00318 return dst [0..(p - dst.ptr)];
00319 }
00320
00321
00322 /+
00323
00324
00325
00326
00327
00328
00329 version (ReadLimit)
00330 {
00331 int len;
00332 const int MaxDigits = 17;
00333
00334 p = d.ptr + trimZero (d) - 1;
00335
00336 while ((c = *++p) >= '0' && c <= '9')
00337 {
00338 if (++len < MaxDigits)
00339 value = value * 10 + (c - '0');
00340 else
00341
00342 ++exp;
00343 }
00344
00345 if (c == '.')
00346 {
00347
00348
00349 if (value == 0)
00350 {
00351 while (*++p == '0')
00352 --exp;
00353 --p;
00354 }
00355
00356 while ((c = *++p) >= '0' && c <= '9')
00357 {
00358 if (++len < MaxDigits)
00359 {
00360 --exp;
00361 value = value * 10 + (c - '0');
00362 }
00363 }
00364 }
00365
00366 if (len == 0)
00367 {
00368 if (c != '0' && d.length >= 3)
00369 if (d == "inf")
00370 return tValue.infinity;
00371 else
00372 if (d == "nan")
00373 return tValue.nan;
00374 return 0.0 * sign;
00375 }
00376
00377 if (c == 'e' || c == 'E')
00378 exp += strtoi (d[(p-d.ptr)+1..length]);
00379 }
00380 +/
00381 }
00382
00383