result of Gamma is never more precise than double 20 largest ulp errors (stored in a double) 1007.51430 ulp for Gamma(164.040279658791491) = 2.46119858953288017e+291) 801.98953 ulp for Gamma(171.624376956302655) = 1.79769313486186445e+308) 742.49998 ulp for Gamma(164.040279658475555) = 2.46119858556998948e+291) 660.48660 ulp for Gamma(121.763991321951437) = 2.60795713456616592e+200) 654.51322 ulp for Gamma(171.62437695630237) = 1.79769313485920759e+308) 574.49999 ulp for Gamma(149.995246887021466) = 3.71933986949636738e+260) 530.69793 ulp for Gamma(134.063437917116232) = 2.02876964214805249e+226) 509.53716 ulp for Gamma(171.624376956302683) = 1.79769313486206882e+308) 492.49996 ulp for Gamma(135.625077550049639) = 4.27308755772562989e+229) 480.42868 ulp for Gamma(135.625077549163308) = 4.27308753914401685e+229) 472.49971 ulp for Gamma(95.832847159702311) = 4.82169199453699892e+147) 466.50069 ulp for Gamma(135.625077550027441) = 4.27308755725975388e+229) 455.53404 ulp for Gamma(134.063437915652969) = 2.02876962761814689e+226) 454.49983 ulp for Gamma(107.262197047543694) = 3.89938829917800549e+170) 452.50053 ulp for Gamma(134.063437915655896) = 2.02876962764697743e+226) 444.50086 ulp for Gamma(164.040279658439005) = 2.46119858511138761e+291) 437.51096 ulp for Gamma(87.6707341875724353) = 4.83784949626082489e+131) 427.49852 ulp for Gamma(164.040279658439317) = 2.46119858511530492e+291) 387.49874 ulp for Gamma(84.5999920667431269) = 5.62388999127852771e+125) 371.50633 ulp for Gamma(107.262197047546593) = 3.89938829923075932e+170) Gamma15751007.51430164.0402796587914912.46119858953288017e+291