/* The computation of z with magnitude of x and sign of y: z = (-1)^signbit(y) * abs(x), i.e. with the same sign bit as y, even if z is a NaN. Note: This function implements copysign from the IEEE-754 standard when no rounding occurs (e.g. if PREC(z) >= PREC(x)).  */
int mpfr_copysign(mpfr_ptr z, mpfr_srcptr x, mpfr_srcptr y, mpfr_rnd_t rnd_mode)
/* The computation of z with magnitude of x and sign of y: z = (-1)^signbit(y) * abs(x), i.e. with the same sign bit as y, even if z is a NaN. Note: This function implements copysign from the IEEE-754 standard when no rounding occurs (e.g. if PREC(z) >= PREC(x)).  */
int mpfr_copysign(mpfr_ptr z, mpfr_srcptr x, mpfr_srcptr y, mpfr_rnd_t rnd_mode)
{
  return mpfr_set4 (z, x, rnd_mode, MPFR_SIGN (y));
}