result of j1 is never more precise than double 1 results were exact to double the precision 20 largest ulp errors (stored in a double) 395722335422678.87500 ulp for j1(-7.01558666981561885) = -3.06922212802401322e-17) 395722335422678.87500 ulp for j1(7.01558666981561885) = 3.06922212802401322e-17) 391978696084206.18750 ulp for j1(3.83170597020751247) = -5.66665631165062675e-17) 391978696084206.18750 ulp for j1(-3.83170597020751247) = 5.66665631165062675e-17) 91419592227950.34375 ulp for j1(13.3236919363142228) = -5.5655522888335781e-17) 91419592227950.34375 ulp for j1(-13.3236919363142228) = 5.5655522888335781e-17) 87931867404032.92188 ulp for j1(7.01558666981561796) = -2.33967564928602544e-16) 87931867404032.92188 ulp for j1(-7.01558666981561796) = 2.33967564928602544e-16) 73299515135831.57812 ulp for j1(29.0468285349168553) = -4.09897693402069815e-17) 73299515135831.57812 ulp for j1(-29.0468285349168553) = 4.09897693402069815e-17) 50946649633552.52344 ulp for j1(-3.83170597020751336) = 4.1670841340517604e-16) 50946649633552.52344 ulp for j1(3.83170597020751336) = -4.1670841340517604e-16) 45425211871697.38281 ulp for j1(-3.83170597020751202) = -1.16243210297381289e-16) 45425211871697.38281 ulp for j1(3.83170597020751202) = 1.16243210297381289e-16) 42352788017203.80469 ulp for j1(60.4694578453474918) = -2.45945428650101575e-17) 34427272864953.46094 ulp for j1(16.4706300508776344) = -3.19778671886064104e-16) 34427272864953.46094 ulp for j1(-16.4706300508776344) = 3.19778671886064104e-16) 32967217250379.86719 ulp for j1(7.01558666981562151) = 8.24671579906767941e-16) 32967217250379.86719 ulp for j1(-7.01558666981562151) = -8.24671579906767941e-16) 25725882141706.30859 ulp for j1(41.6170942128144503) = 7.01792398635892254e-17) j11024395722335422678.87500-7.01558666981561885-3.06922212802401322e-17