Implement Math.expm1 using port from fdlibm.

third_party/fdlibm/fdlibm.cc

Index: third_party/fdlibm/fdlibm.cc
diff --git a/third_party/fdlibm/fdlibm.cc b/third_party/fdlibm/fdlibm.cc
index 2f6eab17e8ac3fc9b9ccb6a14018803ec531500d..a92d5c374da68f17b1a96bd41306526fe2c10689 100644
--- a/third_party/fdlibm/fdlibm.cc
+++ b/third_party/fdlibm/fdlibm.cc
@@ -34,19 +34,19 @@ const double MathConstants::constants[] = {
     2.02226624879595063154e-21,   // pio2_2t   4
     2.02226624871116645580e-21,   // pio2_3    5
     8.47842766036889956997e-32,   // pio2_3t   6
-    -1.66666666666666324348e-01,  // S1        7
+    -1.66666666666666324348e-01,  // S1        7  coefficients for sin
     8.33333333332248946124e-03,   //           8
     -1.98412698298579493134e-04,  //           9
     2.75573137070700676789e-06,   //          10
     -2.50507602534068634195e-08,  //          11
     1.58969099521155010221e-10,   // S6       12
-    4.16666666666666019037e-02,   // C1       13
+    4.16666666666666019037e-02,   // C1       13  coefficients for cos
     -1.38888888888741095749e-03,  //          14
     2.48015872894767294178e-05,   //          15
     -2.75573143513906633035e-07,  //          16
     2.08757232129817482790e-09,   //          17
     -1.13596475577881948265e-11,  // C6       18
-    3.33333333333334091986e-01,   // T0       19
+    3.33333333333334091986e-01,   // T0       19  coefficients for tan
     1.33333333333201242699e-01,   //          20
     5.39682539762260521377e-02,   //          21
     2.18694882948595424599e-02,   //          22
@@ -65,13 +65,20 @@ const double MathConstants::constants[] = {
     1.90821492927058770002e-10,   // ln2_lo   35
     1.80143985094819840000e+16,   // 2^54     36
     6.666666666666666666e-01,     // 2/3      37
-    6.666666666666735130e-01,     // LP1      38
+    6.666666666666735130e-01,     // LP1      38  coefficients for log1p
     3.999999999940941908e-01,     //          39
     2.857142874366239149e-01,     //          40
     2.222219843214978396e-01,     //          41
     1.818357216161805012e-01,     //          42
     1.531383769920937332e-01,     //          43
     1.479819860511658591e-01,     // LP7      44
+    7.09782712893383973096e+02,   //          45  overflow threshold for expm1
+    1.44269504088896338700e+00,   // 1/ln2    46
+    -3.33333333333331316428e-02,  // Q1       47  coefficients for expm1
+    1.58730158725481460165e-03,   //          48
+    -7.93650757867487942473e-05,  //          49
+    4.00821782732936239552e-06,   //          50
+    -2.01099218183624371326e-07   // Q5       51
 };