问题标签 [mpmath]

I would like to convert an mpmath function to a function that can work on numpy arrays. Let's say I have for example the following

However, now the values in the array A are of the mpmath type mpf. So what is the fastest/best way to take a function in mpmath, and convert it to a function that can act on numpy arrays, but that returns standard float and not mpmath floats? Or just convert an array of mpf to numpy floats? the function float() works to convert numbers, but not arrays.

我正在使用 mpmath python 库在某些计算期间获得精度，但我需要将结果转换为 numpy 本机类型。

更准确地说，我需要在 numpy.ndarray （包含浮点类型）中转换一个 mpmath 矩阵（包含 mpf 对象类型）。

我用原始方法解决了这个问题：

我的问题是：

有更好/更优雅/更容易/更聪明的方法吗？

更新：

反复阅读 mpmath 文档，我发现这个非常有用的方法：tolist()，它可以按如下方式使用：

它看起来更好更优雅（不需要 for 循环）

有更好的方法吗？我的第二个解决方案是舍入还是砍掉多余的数字？

我正在寻找用于 python 的 Jacobi theta 函数的导数的实现。我找到了这个

http://mpmath.googlecode.com/svn/trunk/doc/build/functions/elliptic.html#jtheta

但是，导数是针对错误参数计算的。我想要关于 q 的导数，而不是 z。有谁知道获得这个的最佳方法？仅以有限差分对导数进行数值计算并不能很好地工作（我使用的 ODE 求解器在执行此操作时会崩溃）。

（实际上我只希望它也适用于特殊情况 z=0，如果可能的话。对于 E4 Eisenstein 系列的相同问题也足够了）

谢谢！

您好 StackExchange 的人们，

我正在将库 mpmath 用于一个项目，我需要一些有关数组的帮助，即通过函数传递它们。

我用 mpf 元素创建数组没有问题。这段代码摘录让我可以做到这一点：

但是，如果我在下一行调用：

它给了我一个错误。有没有办法通过函数传递带有 mpf 元素的数组？

编辑：我尝试在每个 r1 和 r2 函数中使用矢量化，如下所示：

但它一直不起作用。

Numpy 是一个用于高效数值数组的库。

mpmath 在 gmpy 的支持下，是一个用于高效多精度数的库。

如何有效地将它们组合在一起？还是仅使用带有 mpmath 数字的 Numpy 数组已经很有效？

要求“与本机浮点数一样高效”是没有意义的，但您可以要求它接近等效 C 代码（或者，如果不这样做，Java/C# 代码）的效率。特别是，一个高效的多精度数字数组意味着您可以进行矢量化操作，而不必__add__在 Global Interpreter 中查找一百万次。

编辑：致亲密的选民：我的问题是关于将它们组合在一起的有效方法。可能重复中的答案特别指出，幼稚的方法效率不高。

拥有一个 dtype=object 的 numpy 数组可能会产生误导，因为强大的 numpy 机制使得使用标准 dtypes 的操作超级快，现在由默认对象的 python 运算符处理，这意味着速度不会存在不再

我是 sympy 的新手，但对 python 仍然很天真……我想求解一个三角方程，找到它的零点。（一旦有了语法，我就会使用更复杂的函数。）

我还找不到正确的语法。这是我在 Spyder (Python 2.7) 的 iPython 控制台上尝试的内容：

我收到了这个错误：

回溯（最近一次通话最后）：

文件“”，第 1 行，在求解（sin（x），x）中

NameError：名称'sin'未定义

好的，所以我需要正确引用正弦函数。

根据 sympy 文档，我认为这是 in mpath，但这不起作用：

ImportError：没有名为 mpmath 的模块

如何加载/访问 mpmath 或其他方式来获取正弦函数？

我目前正在查看 Meijerg 的 mpmath文档。它说

mpmath.meijerg（A_S，B_S，Z，r = 1，** Kwargs）评估Meijer G功能，定义为 适当的轮廓L（请参阅参考文献）。

• 有 p 个元素 a_j。参数 a_s 应该是一对列表，第一个包含 n 个元素 a_1, ..., a_n，第二个包含 pn 个元素 a_{n+1}, ... a_p。
• 有 q 个元素 b_j。参数 b_s 应该是一对列表，第一个包含 m 个元素 b_1, ..., b_m，第二个包含 qm 个元素 b_{m+1}, ... b_q。

给定的第二个示例有系列 a、b 的四个参数：

例如，对于第一行，我不清楚这是一对每个都包含一个数字的集合，[[1],[1]]还是一个包含两个元素和一个空集合的集合[[1, 1], []}

令人困惑的是，没有提及，他们做了混合。他们将第二个示例中的函数称为

这意味着前两个数字属于第一个集合，但后两个数字分别属于一个集合。有人可以向我解释为什么吗？这背后的逻辑是什么？

我正在尝试使用mpmath.polyroots整数系数来查找简单多项式的根，该多项式x*(x-4)**3在展开时具有 的系数向量[1, -12, 48, 64, 0]。以下代码失败：

出现错误：

增加步数没有帮助。预期的答案显然是[0,4,4,4]。

如果存在多重性，mpmath 是否无法找到多项式的根？我该如何解决这个问题？

mpmath声称支持“任意精度浮点运算”。

但是 。. .

我错过了什么吗？或者这是一个基本的限制mpmath？

@jonrsharpe 建议问题是我已经提交float了erf. 但是，下面的代码表明这不是问题：

我希望这个问题不是太初级。我已经广泛搜索了一个解决方案，但还没有找到一个。

我最近开始使用带有 Sympy 的 Jupyter Notebook 来记笔记并在我的微积分 II 课上做作业（这是一个巨大的好处！）。

但是，我唯一的问题是我无法弄清楚如何配置绘图的大小（即像素尺寸）。

matplotlib直接使用（特别是）很容易matplotlib.pyplot.figure()，但我正在使用该Sympy.mpmath.plot模块，因为Sympy对于我们在本课程中所做的符号操作来说效果更好。我知道Sympy有它自己的绘图模块，但mpmath到目前为止似乎更容易使用（当然，除了这个问题）。

但是，我查看了mpmath文档并反复搜索问题，但没有解决方案。

如何更改使用mpmathAPI 绘制函数产生的图像大小？