1 /* tpow -- test file for mpc_pow. 2 3 Copyright (C) 2009, 2011 INRIA 4 5 This file is part of GNU MPC. 6 7 GNU MPC is free software; you can redistribute it and/or modify it under 8 the terms of the GNU Lesser General Public License as published by the 9 Free Software Foundation; either version 3 of the License, or (at your 10 option) any later version. 11 12 GNU MPC is distributed in the hope that it will be useful, but WITHOUT ANY 13 WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS 14 FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for 15 more details. 16 17 You should have received a copy of the GNU Lesser General Public License 18 along with this program. If not, see http://www.gnu.org/licenses/ . 19 */ 20 21 #include "mpc-tests.h" 22 23 static void 24 reuse_bug (void) 25 { 26 /* bug found by the automatic builds on 27 http://hydra.nixos.org/build/1469029/log/raw */ 28 mpc_t x, y, z; 29 mp_prec_t prec = 2; 30 31 for (prec = 2; prec <= 20; prec ++) 32 { 33 mpc_init2 (x, prec); 34 mpc_init2 (y, prec); 35 mpc_init2 (z, prec); 36 37 mpfr_set_ui (mpc_realref (x), 0ul, GMP_RNDN); 38 mpfr_set_ui_2exp (mpc_imagref (x), 3ul, -2, GMP_RNDN); 39 mpc_set_ui (y, 8ul, MPC_RNDNN); 40 41 mpc_pow (z, x, y, MPC_RNDNN); 42 mpc_pow (y, x, y, MPC_RNDNN); 43 if (mpfr_signbit (mpc_imagref (y)) != mpfr_signbit (mpc_imagref (z))) 44 { 45 printf ("Error: regression, reuse_bug reproduced\n"); 46 exit (1); 47 } 48 49 mpc_clear (x); 50 mpc_clear (y); 51 mpc_clear (z); 52 } 53 } 54 55 56 int 57 main (void) 58 { 59 DECL_FUNC (C_CC, f, mpc_pow); 60 61 test_start (); 62 63 reuse_bug (); 64 65 data_check (f, "pow.dat"); 66 tgeneric (f, 2, 1024, 7, 10); 67 68 test_end (); 69 70 return 0; 71 } 72