test-fast-dtoa.cc 10.1 KB
Newer Older
1
// Copyright 2006-2008 the V8 project authors. All rights reserved.
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
// Redistribution and use in source and binary forms, with or without
// modification, are permitted provided that the following conditions are
// met:
//
//     * Redistributions of source code must retain the above copyright
//       notice, this list of conditions and the following disclaimer.
//     * Redistributions in binary form must reproduce the above
//       copyright notice, this list of conditions and the following
//       disclaimer in the documentation and/or other materials provided
//       with the distribution.
//     * Neither the name of Google Inc. nor the names of its
//       contributors may be used to endorse or promote products derived
//       from this software without specific prior written permission.
//
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
// "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
// LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
// A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
// OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
// SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
// LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
// DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
// THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
// (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
// OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
27 28 29

#include <stdlib.h>

30
#include "src/init/v8.h"
31

32
#include "src/base/platform/platform.h"
33 34 35
#include "src/numbers/diy-fp.h"
#include "src/numbers/double.h"
#include "src/numbers/fast-dtoa.h"
36 37 38
#include "test/cctest/cctest.h"
#include "test/cctest/gay-precision.h"
#include "test/cctest/gay-shortest.h"
39

40 41
namespace v8 {
namespace internal {
42
namespace test_fast_dtoa {
43 44 45

static const int kBufferSize = 100;

46 47 48 49 50 51
// Removes trailing '0' digits (modifies {representation}). Can create an empty
// string if all digits are 0.
static void TrimRepresentation(char* representation) {
  size_t len = strlen(representation);
  while (len > 0 && representation[len - 1] == '0') --len;
  representation[len] = '\0';
52 53 54
}

TEST(FastDtoaShortestVariousDoubles) {
55 56
  char buffer_container[kBufferSize];
  Vector<char> buffer(buffer_container, kBufferSize);
57 58 59 60 61
  int length;
  int point;
  int status;

  double min_double = 5e-324;
62 63
  status = FastDtoa(min_double, FAST_DTOA_SHORTEST, 0,
                    buffer, &length, &point);
64
  CHECK(status);
65
  CHECK_EQ(0, strcmp("5", buffer.begin()));
66 67 68
  CHECK_EQ(-323, point);

  double max_double = 1.7976931348623157e308;
69 70
  status = FastDtoa(max_double, FAST_DTOA_SHORTEST, 0,
                    buffer, &length, &point);
71
  CHECK(status);
72
  CHECK_EQ(0, strcmp("17976931348623157", buffer.begin()));
73 74
  CHECK_EQ(309, point);

75 76
  status = FastDtoa(4294967272.0, FAST_DTOA_SHORTEST, 0,
                    buffer, &length, &point);
77
  CHECK(status);
78
  CHECK_EQ(0, strcmp("4294967272", buffer.begin()));
79 80
  CHECK_EQ(10, point);

81 82
  status = FastDtoa(4.1855804968213567e298, FAST_DTOA_SHORTEST, 0,
                    buffer, &length, &point);
83
  CHECK(status);
84
  CHECK_EQ(0, strcmp("4185580496821357", buffer.begin()));
85 86
  CHECK_EQ(299, point);

87 88
  status = FastDtoa(5.5626846462680035e-309, FAST_DTOA_SHORTEST, 0,
                    buffer, &length, &point);
89
  CHECK(status);
90
  CHECK_EQ(0, strcmp("5562684646268003", buffer.begin()));
91 92
  CHECK_EQ(-308, point);

93 94
  status = FastDtoa(2147483648.0, FAST_DTOA_SHORTEST, 0,
                    buffer, &length, &point);
95
  CHECK(status);
96
  CHECK_EQ(0, strcmp("2147483648", buffer.begin()));
97 98
  CHECK_EQ(10, point);

99 100
  status = FastDtoa(3.5844466002796428e+298, FAST_DTOA_SHORTEST, 0,
                    buffer, &length, &point);
101
  if (status) {  // Not all FastDtoa variants manage to compute this number.
102
    CHECK_EQ(0, strcmp("35844466002796428", buffer.begin()));
103 104 105
    CHECK_EQ(299, point);
  }

106
  uint64_t smallest_normal64 = 0x0010'0000'0000'0000;
107
  double v = Double(smallest_normal64).value();
108
  status = FastDtoa(v, FAST_DTOA_SHORTEST, 0, buffer, &length, &point);
109
  if (status) {
110
    CHECK_EQ(0, strcmp("22250738585072014", buffer.begin()));
111 112 113
    CHECK_EQ(-307, point);
  }

114
  uint64_t largest_denormal64 = 0x000F'FFFF'FFFF'FFFF;
115
  v = Double(largest_denormal64).value();
116
  status = FastDtoa(v, FAST_DTOA_SHORTEST, 0, buffer, &length, &point);
117
  if (status) {
118
    CHECK_EQ(0, strcmp("2225073858507201", buffer.begin()));
119 120 121 122 123
    CHECK_EQ(-307, point);
  }
}


124 125 126 127 128 129 130 131 132 133
TEST(FastDtoaPrecisionVariousDoubles) {
  char buffer_container[kBufferSize];
  Vector<char> buffer(buffer_container, kBufferSize);
  int length;
  int point;
  int status;

  status = FastDtoa(1.0, FAST_DTOA_PRECISION, 3, buffer, &length, &point);
  CHECK(status);
  CHECK_GE(3, length);
134
  TrimRepresentation(buffer.begin());
135
  CHECK_EQ(0, strcmp("1", buffer.begin()));
136 137 138 139 140
  CHECK_EQ(1, point);

  status = FastDtoa(1.5, FAST_DTOA_PRECISION, 10, buffer, &length, &point);
  if (status) {
    CHECK_GE(10, length);
141
    TrimRepresentation(buffer.begin());
142
    CHECK_EQ(0, strcmp("15", buffer.begin()));
143 144 145 146 147 148 149
    CHECK_EQ(1, point);
  }

  double min_double = 5e-324;
  status = FastDtoa(min_double, FAST_DTOA_PRECISION, 5,
                    buffer, &length, &point);
  CHECK(status);
150
  CHECK_EQ(0, strcmp("49407", buffer.begin()));
151 152 153 154 155 156
  CHECK_EQ(-323, point);

  double max_double = 1.7976931348623157e308;
  status = FastDtoa(max_double, FAST_DTOA_PRECISION, 7,
                    buffer, &length, &point);
  CHECK(status);
157
  CHECK_EQ(0, strcmp("1797693", buffer.begin()));
158 159 160 161 162 163
  CHECK_EQ(309, point);

  status = FastDtoa(4294967272.0, FAST_DTOA_PRECISION, 14,
                    buffer, &length, &point);
  if (status) {
    CHECK_GE(14, length);
164
    TrimRepresentation(buffer.begin());
165
    CHECK_EQ(0, strcmp("4294967272", buffer.begin()));
166 167 168 169 170 171
    CHECK_EQ(10, point);
  }

  status = FastDtoa(4.1855804968213567e298, FAST_DTOA_PRECISION, 17,
                    buffer, &length, &point);
  CHECK(status);
172
  CHECK_EQ(0, strcmp("41855804968213567", buffer.begin()));
173 174 175 176 177
  CHECK_EQ(299, point);

  status = FastDtoa(5.5626846462680035e-309, FAST_DTOA_PRECISION, 1,
                    buffer, &length, &point);
  CHECK(status);
178
  CHECK_EQ(0, strcmp("6", buffer.begin()));
179 180 181 182 183
  CHECK_EQ(-308, point);

  status = FastDtoa(2147483648.0, FAST_DTOA_PRECISION, 5,
                    buffer, &length, &point);
  CHECK(status);
184
  CHECK_EQ(0, strcmp("21475", buffer.begin()));
185 186 187 188 189 190
  CHECK_EQ(10, point);

  status = FastDtoa(3.5844466002796428e+298, FAST_DTOA_PRECISION, 10,
                    buffer, &length, &point);
  CHECK(status);
  CHECK_GE(10, length);
191
  TrimRepresentation(buffer.begin());
192
  CHECK_EQ(0, strcmp("35844466", buffer.begin()));
193 194
  CHECK_EQ(299, point);

195
  uint64_t smallest_normal64 = 0x0010'0000'0000'0000;
196 197 198
  double v = Double(smallest_normal64).value();
  status = FastDtoa(v, FAST_DTOA_PRECISION, 17, buffer, &length, &point);
  CHECK(status);
199
  CHECK_EQ(0, strcmp("22250738585072014", buffer.begin()));
200 201
  CHECK_EQ(-307, point);

202
  uint64_t largest_denormal64 = 0x000F'FFFF'FFFF'FFFF;
203 204 205 206
  v = Double(largest_denormal64).value();
  status = FastDtoa(v, FAST_DTOA_PRECISION, 17, buffer, &length, &point);
  CHECK(status);
  CHECK_GE(20, length);
207
  TrimRepresentation(buffer.begin());
208
  CHECK_EQ(0, strcmp("22250738585072009", buffer.begin()));
209 210 211 212 213
  CHECK_EQ(-307, point);

  v = 3.3161339052167390562200598e-237;
  status = FastDtoa(v, FAST_DTOA_PRECISION, 18, buffer, &length, &point);
  CHECK(status);
214
  CHECK_EQ(0, strcmp("331613390521673906", buffer.begin()));
215 216 217 218 219
  CHECK_EQ(-236, point);

  v = 7.9885183916008099497815232e+191;
  status = FastDtoa(v, FAST_DTOA_PRECISION, 4, buffer, &length, &point);
  CHECK(status);
220
  CHECK_EQ(0, strcmp("7989", buffer.begin()));
221 222 223 224
  CHECK_EQ(192, point);
}


225
TEST(FastDtoaGayShortest) {
226 227
  char buffer_container[kBufferSize];
  Vector<char> buffer(buffer_container, kBufferSize);
228 229 230 231 232 233 234
  bool status;
  int length;
  int point;
  int succeeded = 0;
  int total = 0;
  bool needed_max_length = false;

235 236
  Vector<const PrecomputedShortest> precomputed =
      PrecomputedShortestRepresentations();
237
  for (int i = 0; i < precomputed.length(); ++i) {
238
    const PrecomputedShortest current_test = precomputed[i];
239 240
    total++;
    double v = current_test.v;
241
    status = FastDtoa(v, FAST_DTOA_SHORTEST, 0, buffer, &length, &point);
242
    CHECK_GE(kFastDtoaMaximalLength, length);
243
    if (!status) continue;
244
    if (length == kFastDtoaMaximalLength) needed_max_length = true;
245 246
    succeeded++;
    CHECK_EQ(current_test.decimal_point, point);
247
    CHECK_EQ(0, strcmp(current_test.representation, buffer.begin()));
248 249 250 251
  }
  CHECK_GT(succeeded*1.0/total, 0.99);
  CHECK(needed_max_length);
}
252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279


TEST(FastDtoaGayPrecision) {
  char buffer_container[kBufferSize];
  Vector<char> buffer(buffer_container, kBufferSize);
  bool status;
  int length;
  int point;
  int succeeded = 0;
  int total = 0;
  // Count separately for entries with less than 15 requested digits.
  int succeeded_15 = 0;
  int total_15 = 0;

  Vector<const PrecomputedPrecision> precomputed =
      PrecomputedPrecisionRepresentations();
  for (int i = 0; i < precomputed.length(); ++i) {
    const PrecomputedPrecision current_test = precomputed[i];
    double v = current_test.v;
    int number_digits = current_test.number_digits;
    total++;
    if (number_digits <= 15) total_15++;
    status = FastDtoa(v, FAST_DTOA_PRECISION, number_digits,
                      buffer, &length, &point);
    CHECK_GE(number_digits, length);
    if (!status) continue;
    succeeded++;
    if (number_digits <= 15) succeeded_15++;
280
    TrimRepresentation(buffer.begin());
281
    CHECK_EQ(current_test.decimal_point, point);
282
    CHECK_EQ(0, strcmp(current_test.representation, buffer.begin()));
283 284 285 286 287 288 289 290 291
  }
  // The precomputed numbers contain many entries with many requested
  // digits. These have a high failure rate and we therefore expect a lower
  // success rate than for the shortest representation.
  CHECK_GT(succeeded*1.0/total, 0.85);
  // However with less than 15 digits almost the algorithm should almost always
  // succeed.
  CHECK_GT(succeeded_15*1.0/total_15, 0.9999);
}
292

293
}  // namespace test_fast_dtoa
294 295
}  // namespace internal
}  // namespace v8