Home
last modified time | relevance | path

Searched refs:mpc_realref (Results 1 – 25 of 70) sorted by relevance

123

/dflybsd-src/contrib/mpc/src/
H A Dasin.c32 if (mpfr_nan_p (mpc_realref (op)) || mpfr_nan_p (mpc_imagref (op))) in mpc_asin()
34 if (mpfr_inf_p (mpc_realref (op)) || mpfr_inf_p (mpc_imagref (op))) in mpc_asin()
36 mpfr_set_nan (mpc_realref (rop)); in mpc_asin()
39 else if (mpfr_zero_p (mpc_realref (op))) in mpc_asin()
41 mpfr_set (mpc_realref (rop), mpc_realref (op), GMP_RNDN); in mpc_asin()
46 mpfr_set_nan (mpc_realref (rop)); in mpc_asin()
53 if (mpfr_inf_p (mpc_realref (op)) || mpfr_inf_p (mpc_imagref (op))) in mpc_asin()
56 if (mpfr_inf_p (mpc_realref (op))) in mpc_asin()
60 inex_re = set_pi_over_2 (mpc_realref (rop), in mpc_asin()
61 (mpfr_signbit (mpc_realref (op)) ? -1 : 1), MPC_RND_RE (rnd)); in mpc_asin()
[all …]
H A Dtan.c37 if (mpfr_nan_p (mpc_realref (op))) in mpc_tan()
52 mpfr_set_nan (mpc_realref (rop)); in mpc_tan()
59 if (mpfr_cmp_ui (mpc_realref (op), 0) == 0) in mpc_tan()
69 mpfr_set_nan (mpc_realref (rop)); in mpc_tan()
74 else if (mpfr_inf_p (mpc_realref (op))) in mpc_tan()
82 const int sign_re = mpfr_signbit (mpc_realref (op)); in mpc_tan()
85 mpfr_set_ui (mpc_realref (rop), 0, MPC_RND_RE (rnd)); in mpc_tan()
86 mpfr_setsign (mpc_realref (rop), mpc_realref (rop), sign_re, GMP_RNDN); in mpc_tan()
98 mpfr_set_nan (mpc_realref (rop)); in mpc_tan()
114 mpfr_sin_cos (s, c, mpc_realref (op), GMP_RNDN); in mpc_tan()
[all …]
H A Dacos.c39 if (mpfr_nan_p (mpc_realref (op)) || mpfr_nan_p (mpc_imagref (op))) in mpc_acos()
41 if (mpfr_inf_p (mpc_realref (op)) || mpfr_inf_p (mpc_imagref (op))) in mpc_acos()
44 mpfr_set_nan (mpc_realref (rop)); in mpc_acos()
46 else if (mpfr_zero_p (mpc_realref (op))) in mpc_acos()
48 inex_re = set_pi_over_2 (mpc_realref (rop), +1, MPC_RND_RE (rnd)); in mpc_acos()
53 mpfr_set_nan (mpc_realref (rop)); in mpc_acos()
60 if (mpfr_inf_p (mpc_realref (op)) || mpfr_inf_p (mpc_imagref (op))) in mpc_acos()
62 if (mpfr_inf_p (mpc_realref (op))) in mpc_acos()
66 if (mpfr_sgn (mpc_realref (op)) > 0) in mpc_acos()
69 set_pi_over_2 (mpc_realref (rop), +1, MPC_RND_RE (rnd)); in mpc_acos()
[all …]
H A Ddiv.c28 int sign = MPFR_SIGNBIT (mpc_realref (w)); in mpc_div_zero()
33 mpfr_mul (mpc_realref (a), infty, mpc_realref (z), MPC_RND_RE (rnd)); in mpc_div_zero()
47 a = (mpfr_inf_p (mpc_realref (z)) ? MPFR_SIGNBIT (mpc_realref (z)) : 0); in mpc_div_inf_fin()
57 x = a * MPC_MPFR_SIGN (mpc_realref (w)) + b * MPC_MPFR_SIGN (mpc_imagref (w)); in mpc_div_inf_fin()
58 y = b * MPC_MPFR_SIGN (mpc_realref (w)) - a * MPC_MPFR_SIGN (mpc_imagref (w)); in mpc_div_inf_fin()
71 mpfr_add (sign, mpc_realref (w), mpc_imagref (w), GMP_RNDN); in mpc_div_inf_fin()
73 mpfr_sub (sign, mpc_realref (w), mpc_imagref (w), GMP_RNDN); in mpc_div_inf_fin()
77 mpfr_sub (sign, mpc_realref (w), mpc_imagref (w), GMP_RNDN); in mpc_div_inf_fin()
79 mpfr_add (sign, mpc_realref (w), mpc_imagref (w), GMP_RNDN); in mpc_div_inf_fin()
84 mpfr_sub (sign, mpc_imagref (w), mpc_realref (w), GMP_RNDN); in mpc_div_inf_fin()
[all …]
H A Dlog10.c36 prec = mpfr_get_prec ((flag == 0) ? mpc_realref (rop) : mpc_imagref (rop)); in mpc_log10_aux()
50 mpfr_atan2 (mpc_imagref (tmp), mpc_imagref (op), mpc_realref (op), in mpc_log10_aux()
61 mpfr_log (mpc_realref (tmp), mpc_realref (op), MPC_RND_RE (rnd)); in mpc_log10_aux()
62 mpfr_div (mpc_realref (tmp), mpc_realref (tmp), log10, GMP_RNDN); in mpc_log10_aux()
63 ok = mpfr_can_round (mpc_realref (tmp), prec - 2, GMP_RNDN, in mpc_log10_aux()
67 ret = mpfr_set (mpc_realref (rop), mpc_realref (tmp), in mpc_log10_aux()
103 if (mpfr_nan_p (mpc_realref (op))) in mpc_log10()
107 mpfr_set_inf (mpc_realref (rop), +1); in mpc_log10()
110 mpfr_set_nan (mpc_realref (rop)); in mpc_log10()
116 if (mpfr_inf_p (mpc_realref (op))) in mpc_log10()
[all …]
H A Dsin_cos.c40 if (mpfr_nan_p (mpc_realref (op_loc)) || mpfr_nan_p (mpc_imagref (op_loc))) { in mpc_sin_cos_nonfinite()
46 if (!mpfr_zero_p (mpc_realref (op_loc))) in mpc_sin_cos_nonfinite()
47 mpfr_set_nan (mpc_realref (rop_sin)); in mpc_sin_cos_nonfinite()
58 else if (mpfr_inf_p (mpc_realref (op_loc))) { in mpc_sin_cos_nonfinite()
59 mpfr_set_nan (mpc_realref (rop_sin)); in mpc_sin_cos_nonfinite()
71 else if (mpfr_zero_p (mpc_realref (op_loc))) { in mpc_sin_cos_nonfinite()
84 mpfr_sin_cos (s, c, mpc_realref (op_loc), GMP_RNDZ); in mpc_sin_cos_nonfinite()
85 mpfr_set_inf (mpc_realref (rop_sin), MPFR_SIGN (s)); in mpc_sin_cos_nonfinite()
93 if (mpfr_nan_p (mpc_realref (op_loc))) { in mpc_sin_cos_nonfinite()
101 mpfr_set_inf (mpc_realref (rop_cos), +1); in mpc_sin_cos_nonfinite()
[all …]
H A Dlog.c39 if (mpfr_nan_p (mpc_realref (op))) { in mpc_log()
41 mpfr_set_inf (mpc_realref (rop), +1); in mpc_log()
43 mpfr_set_nan (mpc_realref (rop)); in mpc_log()
48 if (mpfr_inf_p (mpc_realref (op))) in mpc_log()
49 mpfr_set_inf (mpc_realref (rop), +1); in mpc_log()
51 mpfr_set_nan (mpc_realref (rop)); in mpc_log()
56 inex_im = mpfr_atan2 (mpc_imagref (rop), mpc_imagref (op), mpc_realref (op), in mpc_log()
58 mpfr_set_inf (mpc_realref (rop), +1); in mpc_log()
64 re_cmp = mpfr_cmp_ui (mpc_realref (op), 0); in mpc_log()
68 inex_im = mpfr_atan2 (mpc_imagref (rop), mpc_imagref (op), mpc_realref (op), in mpc_log()
[all …]
H A Dpow.c132 if (mpfr_zero_p (mpc_realref(z))) in fix_sign()
140 mpfr_neg (mpc_realref(z), mpc_realref(z), GMP_RNDZ); in fix_sign()
180 int sign_rex = mpfr_signbit (mpc_realref(x)); in mpc_pow_exact()
182 int x_imag = mpfr_zero_p (mpc_realref(x)); in mpc_pow_exact()
186 if (mpc_realref (z) == y || mpc_imagref (z) == y) in mpc_pow_exact()
213 ec = mpfr_get_z_exp (c, mpc_realref(x)); in mpc_pow_exact()
405 ret = mpfr_set_z (mpc_realref(z), a, MPC_RND_RE(rnd)); in mpc_pow_exact()
407 mpfr_mul_2si (mpc_realref(z), mpc_realref(z), ed, MPC_RND_RE(rnd)); in mpc_pow_exact()
494 if (y_real && mpfr_zero_p (mpc_realref(y))) /* case y zero */ in mpc_pow()
496 if (x_real && mpfr_zero_p (mpc_realref(x))) in mpc_pow()
[all …]
H A Dmul.c36 int xrs = mpfr_signbit (mpc_realref (x)) ? -1 : 1; in mul_infinite()
38 int yrs = mpfr_signbit (mpc_realref (y)) ? -1 : 1; in mul_infinite()
47 if ( mpfr_nan_p (mpc_realref (x)) || mpfr_nan_p (mpc_imagref (x)) in mul_infinite()
48 || mpfr_nan_p (mpc_realref (y)) || mpfr_nan_p (mpc_imagref (y))) { in mul_infinite()
52 else if (mpfr_inf_p (mpc_realref (x))) { in mul_infinite()
54 u = ( mpfr_zero_p (mpc_realref (y)) in mul_infinite()
61 || (mpfr_inf_p (mpc_imagref (x)) && mpfr_zero_p (mpc_realref (y))) in mul_infinite()
62 || (mpfr_zero_p (mpc_imagref (x)) && mpfr_inf_p (mpc_realref (y))) in mul_infinite()
70 || (mpfr_zero_p (mpc_realref (x)) && mpfr_inf_p (mpc_realref (y))) in mul_infinite()
71 || (mpfr_inf_p (mpc_realref (y)) && xrs*yrs == xis*yis) in mul_infinite()
[all …]
H A Dfma.c49 mpfr_init2 (rea_reb, mpfr_get_prec (mpc_realref(a)) + mpfr_get_prec (mpc_realref(b))); in mpc_fma_naive()
50 mpfr_init2 (rea_imb, mpfr_get_prec (mpc_realref(a)) + mpfr_get_prec (mpc_imagref(b))); in mpc_fma_naive()
51 mpfr_init2 (ima_reb, mpfr_get_prec (mpc_imagref(a)) + mpfr_get_prec (mpc_realref(b))); in mpc_fma_naive()
54 mpfr_mul (rea_reb, mpc_realref(a), mpc_realref(b), GMP_RNDZ); /* exact */ in mpc_fma_naive()
55 mpfr_mul (rea_imb, mpc_realref(a), mpc_imagref(b), GMP_RNDZ); /* exact */ in mpc_fma_naive()
56 mpfr_mul (ima_reb, mpc_imagref(a), mpc_realref(b), GMP_RNDZ); /* exact */ in mpc_fma_naive()
63 pre13 = bound_prec_addsub (rea_reb, mpc_realref(c)); in mpc_fma_naive()
65 pre23 = bound_prec_addsub (ima_imb, mpc_realref(c)); in mpc_fma_naive()
71 inex_re = mpfr_add (mpc_realref(r), tmp, mpc_realref(c), MPC_RND_RE(rnd)); in mpc_fma_naive()
78 mpfr_add (tmp, rea_reb, mpc_realref(c), GMP_RNDZ); /* exact */ in mpc_fma_naive()
[all …]
H A Dexp.c33 if (mpfr_nan_p (mpc_realref (op)) || mpfr_nan_p (mpc_imagref (op))) in mpc_exp()
45 if (mpfr_inf_p (mpc_realref (op))) in mpc_exp()
47 if (mpfr_signbit (mpc_realref (op))) in mpc_exp()
51 mpfr_set_inf (mpc_realref (rop), +1); in mpc_exp()
56 mpfr_set_nan (mpc_realref (rop)); in mpc_exp()
67 inex_re = mpfr_exp (mpc_realref(rop), mpc_realref(op), MPC_RND_RE(rnd)); in mpc_exp()
72 if (mpfr_zero_p (mpc_realref (op))) in mpc_exp()
75 inex_re = mpfr_cos (mpc_realref (rop), mpc_imagref (op), MPC_RND_RE(rnd)); in mpc_exp()
81 if (mpfr_inf_p (mpc_realref (op))) in mpc_exp()
90 if (mpfr_signbit (mpc_realref (op))) in mpc_exp()
[all …]
H A Datan.c56 s_re = mpfr_signbit (mpc_realref (op)); in mpc_atan()
60 if (mpfr_nan_p (mpc_realref (op)) || mpfr_nan_p (mpc_imagref (op))) in mpc_atan()
62 if (mpfr_nan_p (mpc_realref (op))) in mpc_atan()
64 mpfr_set_nan (mpc_realref (rop)); in mpc_atan()
76 if (mpfr_inf_p (mpc_realref (op))) in mpc_atan()
78 inex_re = set_pi_over_2 (mpc_realref (rop), -s_re, MPC_RND_RE (rnd)); in mpc_atan()
83 mpfr_set_nan (mpc_realref (rop)); in mpc_atan()
90 if (mpfr_inf_p (mpc_realref (op)) || mpfr_inf_p (mpc_imagref (op))) in mpc_atan()
92 inex_re = set_pi_over_2 (mpc_realref (rop), -s_re, MPC_RND_RE (rnd)); in mpc_atan()
104 inex_re = mpfr_atan (mpc_realref (rop), mpc_realref (op), MPC_RND_RE (rnd)); in mpc_atan()
[all …]
H A Dmul_i.c35 mpfr_swap (mpc_realref (a), mpc_imagref (a)); in mpc_mul_i()
38 mpfr_set (mpc_realref (a), mpc_imagref (b), GMP_RNDN); in mpc_mul_i()
39 mpfr_set (mpc_imagref (a), mpc_realref (b), GMP_RNDN); in mpc_mul_i()
42 MPFR_CHANGE_SIGN (mpc_realref (a)); in mpc_mul_i()
56 inex_im = mpfr_set (mpc_imagref (a), mpc_realref (b), MPC_RND_IM (rnd)); in mpc_mul_i()
61 inex_im = mpfr_neg (mpc_imagref (a), mpc_realref (b), MPC_RND_IM (rnd)); in mpc_mul_i()
63 mpfr_clear (mpc_realref (a)); in mpc_mul_i()
64 mpc_realref (a)[0] = tmp [0]; in mpc_mul_i()
69 inex_re = mpfr_neg (mpc_realref (a), mpc_imagref (b), MPC_RND_RE (rnd)); in mpc_mul_i()
70 inex_im = mpfr_set (mpc_imagref (a), mpc_realref (b), MPC_RND_IM (rnd)); in mpc_mul_i()
[all …]
H A Dsqrt.c41 const int re_cmp = mpfr_cmp_ui (mpc_realref (b), 0), in mpc_sqrt()
58 mpfr_set_inf (mpc_realref (a), +1); in mpc_sqrt()
63 if (mpfr_inf_p (mpc_realref (b))) in mpc_sqrt()
65 if (mpfr_signbit (mpc_realref (b))) in mpc_sqrt()
71 mpfr_set_ui (mpc_realref (a), 0, GMP_RNDN); in mpc_sqrt()
78 mpfr_set_nan (mpc_realref (a)); in mpc_sqrt()
89 mpfr_set_inf (mpc_realref (a), +1); in mpc_sqrt()
107 if (mpfr_nan_p (mpc_realref (b)) || mpfr_nan_p (mpc_imagref (b))) in mpc_sqrt()
109 mpfr_set_nan (mpc_realref (a)); in mpc_sqrt()
127 inex_w = mpfr_sqrt (mpc_realref (a), mpc_realref (b), MPC_RND_RE (rnd)); in mpc_sqrt()
[all …]
H A Dsqr.c177 if (mpfr_nan_p (mpc_realref (op)) || mpfr_nan_p (mpc_imagref (op))) { in mpc_sqr()
178 mpfr_set_nan (mpc_realref (rop)); in mpc_sqr()
181 else if (mpfr_inf_p (mpc_realref (op))) { in mpc_sqr()
184 MPFR_SIGN (mpc_realref (op)) * MPFR_SIGN (mpc_imagref (op))); in mpc_sqr()
185 mpfr_set_nan (mpc_realref (rop)); in mpc_sqr()
192 MPFR_SIGN (mpc_realref (op)) * MPFR_SIGN (mpc_imagref (op))); in mpc_sqr()
193 mpfr_set_inf (mpc_realref (rop), +1); in mpc_sqr()
197 if (mpfr_zero_p (mpc_realref (op))) in mpc_sqr()
201 MPFR_SIGN (mpc_realref (op)) * MPFR_SIGN (mpc_imagref (op))); in mpc_sqr()
202 mpfr_set_inf (mpc_realref (rop), -1); in mpc_sqr()
[all …]
H A Dacosh.c36 if (mpfr_zero_p (mpc_realref (op)) && mpfr_nan_p (mpc_imagref (op))) in mpc_acosh()
38 mpfr_set_nan (mpc_realref (rop)); in mpc_acosh()
52 tmp[0] = mpc_realref (a)[0]; in mpc_acosh()
53 mpc_realref (a)[0] = mpc_imagref (a)[0]; in mpc_acosh()
64 tmp[0] = mpc_realref (a)[0]; in mpc_acosh()
65 mpc_realref (a)[0] = mpc_imagref (a)[0]; in mpc_acosh()
67 MPFR_CHANGE_SIGN (mpc_realref (a)); in mpc_acosh()
H A Dsinh.c34 mpc_realref (z)[0] = mpc_imagref (op)[0]; in mpc_sinh()
35 mpc_imagref (z)[0] = mpc_realref (op)[0]; in mpc_sinh()
36 mpc_realref (sin_z)[0] = mpc_imagref (rop)[0]; in mpc_sinh()
37 mpc_imagref (sin_z)[0] = mpc_realref (rop)[0]; in mpc_sinh()
42 mpc_realref (rop)[0] = mpc_imagref (sin_z)[0]; in mpc_sinh()
43 mpc_imagref (rop)[0] = mpc_realref (sin_z)[0]; in mpc_sinh()
H A Dtanh.c34 mpc_realref (z)[0] = mpc_imagref (op)[0]; in mpc_tanh()
35 mpc_imagref (z)[0] = mpc_realref (op)[0]; in mpc_tanh()
36 mpc_realref (tan_z)[0] = mpc_imagref (rop)[0]; in mpc_tanh()
37 mpc_imagref (tan_z)[0] = mpc_realref (rop)[0]; in mpc_tanh()
42 mpc_realref (rop)[0] = mpc_imagref (tan_z)[0]; in mpc_tanh()
43 mpc_imagref (rop)[0] = mpc_realref (tan_z)[0]; in mpc_tanh()
H A Dmul_fr.c29 if (c == mpc_realref (a)) in mpc_mul_fr()
33 real [0] = mpc_realref (a) [0]; in mpc_mul_fr()
35 inex_re = mpfr_mul (real, mpc_realref(b), c, MPC_RND_RE(rnd)); in mpc_mul_fr()
37 mpfr_set (mpc_realref (a), real, GMP_RNDN); /* exact */ in mpc_mul_fr()
39 if (c == mpc_realref (a)) in mpc_mul_fr()
H A Datanh.c31 mpc_realref (z)[0] = mpc_imagref (op)[0]; in mpc_atanh()
32 mpc_imagref (z)[0] = mpc_realref (op)[0]; in mpc_atanh()
33 MPFR_CHANGE_SIGN (mpc_realref (z)); in mpc_atanh()
42 tmp[0] = mpc_realref (a)[0]; in mpc_atanh()
43 mpc_realref (a)[0] = mpc_imagref (a)[0]; in mpc_atanh()
H A Dasinh.c32 mpc_realref (z)[0] = mpc_imagref (op)[0]; in mpc_asinh()
33 mpc_imagref (z)[0] = mpc_realref (op)[0]; in mpc_asinh()
34 MPFR_CHANGE_SIGN (mpc_realref (z)); in mpc_asinh()
45 tmp[0] = mpc_realref (a)[0]; in mpc_asinh()
46 mpc_realref (a)[0] = mpc_imagref (a)[0]; in mpc_asinh()
H A Dmpc-impl.h101 #define MPC_PREC_RE(x) (mpfr_get_prec(mpc_realref(x)))
108 #define mpc_inf_p(z) (mpfr_inf_p(mpc_realref(z))||mpfr_inf_p(mpc_imagref(z)))
111 #define mpc_zero_p(z) (mpfr_zero_p(mpc_realref(z))&&mpfr_zero_p(mpc_imagref(z)))
114 #define mpc_fin_p(z) (mpfr_number_p(mpc_realref(z))&&mpfr_number_p(mpc_imagref(z)))
116 …c_nan_p(z) ((mpfr_nan_p(mpc_realref(z)) && !mpfr_inf_p(mpc_imagref(z))) || (mpfr_nan_p(mpc_imagref…
H A Dcosh.c30 mpc_realref (z)[0] = mpc_imagref (op)[0]; in mpc_cosh()
31 mpc_imagref (z)[0] = mpc_realref (op)[0]; in mpc_cosh()
32 MPFR_CHANGE_SIGN (mpc_realref (z)); in mpc_cosh()
H A Dadd.c29 inex_re = mpfr_add (mpc_realref(a), mpc_realref(b), mpc_realref(c), MPC_RND_RE(rnd)); in mpc_add()
H A Dsub.c28 inex_re = mpfr_sub (mpc_realref(a), mpc_realref(b), mpc_realref(c), MPC_RND_RE(rnd)); in mpc_sub()

123