如何准确复制mpfr_t？

Question

我只想知道如何使用mpfr_t 3.1.4制作一个精确的MPFR副本:相同的精度，完全相同的值，没有舍入。

void copy (mpfr_t x)
{
    mpfr_t y;
    mpfr_init2 (y, mpfr_get_prec (x));
    // Bummer: mpfr_set will round.
}

Answer 1

最好的方法是像您一样在设置精度之后使用mpfr_set。没有四舍五入，或者如果您愿意，舍入操作将不会修改值。如果值为NaN，则它的符号位将不会被复制，但您可能并不关心(如果您真的关心，仍然可以使用mpfr_copysign复制符号位)。

不要用mpfr_t或类似的东西复制memcpy。在此副本之后，您将得到相同的值，但两个数字都将共享这个值(因为其中一个字段是指向此值的指针)，而修改其中的一个数字将修改另一个数字，您将得到更不稳定的行为(可能的崩溃.)如果数字变化到一个特定的值(NaN，±Inf，±0)或其精度发生变化。