具有mpfr系数的对称多项式？

Question

我想使用Sympy的多项式，但我也想使用更高的系数。

只是在做这似乎给了我多项式与sympy.core.numbers.float系数。

import sympy
from sympy import Poly
from sympy.abc import x
from gmpy2 import mpfr, get_context

get_context().precision = 150

#float64 can't tell this from 1.0
one_and_change = mpfr('1.0000000000000000000000000000000000001')
#mpfr('1.0000000000000000000000000000000000001000000005',150)

p = [one_and_change]
px = Poly(p, x)

print(px)
# Poly(1.0, x, domain='RR')
print(px.is_one)
# True
print(type(px.all_coeffs()[0]))
# <class 'sympy.core.numbers.Float'>

我也尝试了sympy.mpmath.mpf ，结果相同。

这也不起作用：[ 1 ]

domain = sympy.polys.domains.realfield.RealField(150)
px = Poly(p, x, domain=domain)
print(type(px.all_coeffs()[0]))
# <class 'sympy.core.numbers.Float'>

Answer 1

有一些障碍：

• gmpy.mpfr没有._sympy_方法，因此它将在中间步骤中转换为float。
• sympy.Poly默认情况下，将sympy.polys.domain.RR用于浮点系数，并将用于转换。
• RR加载时使用53作为其精度，并忽略mpmath.mp.precision 。
• 转换时， RR使用RR.precision并忽略参数的精度。

解：

1. 系数必须是Sympy类型（例如sympy.Float ，具有扩展的精度）。
2. 对于域，请使用以下之一：
   • 设置sympy.polys.domain.RR._context.prec = my_precision 。
   • 通过domain='EX' ，不进行转换。
   • 传递自定义域（例如sympy.polys.domains.realfield.RealField(150) ）。