dist/tests/tgamma.c

45   mpfr_gamma (y, x, MPFR_RNDN);
53   mpfr_gamma (y, x, MPFR_RNDN);
61   mpfr_gamma (y, x, MPFR_RNDN);
68   mpfr_set_ui (x, 0, MPFR_RNDN);
69   mpfr_gamma (y, x, MPFR_RNDN);
76   mpfr_set_ui (x, 0, MPFR_RNDN);
77   mpfr_neg (x, x, MPFR_RNDN);
78   mpfr_gamma (y, x, MPFR_RNDN);
85   mpfr_set_ui (x, 1, MPFR_RNDN);
86   mpfr_gamma (y, x, MPFR_RNDN);
93   mpfr_set_si (x, -1, MPFR_RNDN);
94   mpfr_gamma (y, x, MPFR_RNDN);
107   mpfr_set_str (x, CHECK_X1, 10, MPFR_RNDN);
108   mpfr_gamma (y, x, MPFR_RNDN);
109   mpfr_set_str (x, CHECK_Y1, 10, MPFR_RNDN);
122   mpfr_set_str (x, CHECK_X2, 10, MPFR_RNDN);
123   mpfr_gamma (y, x, MPFR_RNDN);
124   mpfr_set_str (x, CHECK_Y2, 10, MPFR_RNDN);
137   mpfr_set_ui (x, 33, MPFR_RNDN);
149   mpfr_set_ui (y, 120, MPFR_RNDN);
164   inex = mpfr_gamma (y, x, MPFR_RNDN);
181   mpfr_const_pi (x, MPFR_RNDN);
182   mpfr_gamma (y, x, MPFR_RNDN);
187   mpfr_gamma (y, x, MPFR_RNDN);
215   mpfr_gamma (y, x, MPFR_RNDN);
226   mpfr_set_str (x, "-200000000.5", 10, MPFR_RNDN); /* exact */
227   mpfr_gamma (y, x, MPFR_RNDN);
239   mpfr_set_str (x, "-200000000.1", 10, MPFR_RNDN);
240   mpfr_gamma (y, x, MPFR_RNDN);
253   mpfr_set_str (x, "-2.0000000000000000000000005", 10, MPFR_RNDN);
254   mpfr_gamma (y, x, MPFR_RNDN);
271   mpfr_gamma (y, x, MPFR_RNDN);
288   mpfr_gamma (y, x, MPFR_RNDN);
304   mpfr_gamma (y, x, MPFR_RNDN);
319                 MPFR_RNDN);
320   mpfr_div_2exp (x, x, 170, MPFR_RNDN);
321   mpfr_gamma (y, x, MPFR_RNDN);
323   mpfr_set_str (x, "201948391737", 10, MPFR_RNDN);
324   mpfr_mul_2exp (x, x, 92, MPFR_RNDN);
337                 MPFR_RNDN);
338   mpfr_gamma (y, x, MPFR_RNDN);
340   mpfr_set_str (x, "-3.623795987425E-121243", 10, MPFR_RNDN);
353   mpfr_set_str (x, "-200000000.1", 10, MPFR_RNDN);
354   mpfr_gamma (y, x, MPFR_RNDN);
366   mpfr_set_str (x, "423786866", 10, MPFR_RNDN);
367   mpfr_gamma (y, x, MPFR_RNDN);
376   mpfr_set_ui (x, 3, MPFR_RNDN);
377   inex = mpfr_gamma (x, x, MPFR_RNDN);
402   mpfr_set_str (x, "1e19", 10, MPFR_RNDN);
403   inex = mpfr_gamma (x, x, MPFR_RNDN);
425       mpfr_set_ui (x, i+1, MPFR_RNDN);
426       mpfr_gamma (y, x, MPFR_RNDN);
427       mpfr_set_z (x, n, MPFR_RNDN);
453   mpfr_set_str (x, "-1000001.5", 10, MPFR_RNDN);
454   inex = mpfr_gamma (x, x, MPFR_RNDN);
461   mpfr_set_str (x, "-1000000001.5", 10, MPFR_RNDN);
462   inex = mpfr_gamma (x, x, MPFR_RNDN);
475   mpfr_set_str (x, "-4.6308260837372266e+07", 10, MPFR_RNDN);
476   inex = mpfr_gamma (x, x, MPFR_RNDN);
495   mpfr_set_ui_2exp (x, 5, -1, MPFR_RNDN);
497   inex1 = mpfr_gamma (y, x, MPFR_RNDN);
502   inex2 = mpfr_gamma (z, x, MPFR_RNDN);
518   mpfr_set_ui_2exp (x, 32769, -60, MPFR_RNDN);
563   mpfr_set_ui_2exp (x, 1, -60, MPFR_RNDN);
589   mpfr_set_ui (x, 0, MPFR_RNDN);
611   /* MPFR_RNDN: overflow, infinity since 1/x = 2^(emax + 1) */
612   inex2 = mpfr_gamma (z, x, MPFR_RNDN);
644   mpfr_mul_2ui (x, x, 1, MPFR_RNDN);  /* x = 2^emin */
665   /* MPFR_RNDN: overflow, infinity since 1/x = 2^emax */
666   inex2 = mpfr_gamma (z, x, MPFR_RNDN);
718   mpfr_set_ui_2exp (x, 1, e, MPFR_RNDN);
747               if (overflow && rr == MPFR_RNDN && s == 1)
776       mpfr_neg (x, x, MPFR_RNDN);