mpmath

我试图找到一个表达式的拉普拉斯逆，在声明时除一个变量外，所有变量都已定义： from numpy import * import mpmath as mp p0 = 1 E = 2 c= 3 L = 4 x = 2.5 t = linspace(1,5,10) ulaplace = [] def U(s): return(c*p0*(-exp(L*s/c) + exp(s*(L + 2... »

最近我想编写一些代码来进行拉普拉斯逆变换。 我使用 mpmath.invertlaplace 函数。 它工作得很好，但我想删除 for 循环以获得这个结果。 这是一个例子。 import mpmath as mp def func = lambda s: 1 / (s + k) Time = mp.linspace(1,10,100) result = [] for i in Ti... »

我试图在 python 中描述Riemann Xi function 以进行可视化。这是一个整个函数，没有孔或极。在将数学转换为一些简单的 Python 代码时，我注意到有一个无穷大*零取消，很明显，Python 无法处理。 我当前的代码是这样的，使用 mpmath 库： def xi(s): return 1/2*s*(s-1)*cmath.pi**(-s/2)*mpmath.gamma(... »

我想使用 Python 对数据数组执行数值拉普拉斯逆变换。 我在mpmath 中找到了一个名为 invertlaplace 的算法，但它只接受 lambda 函数。... »

我需要对 g(u)jn(u) 类型进行积分，其中 g(u) 是一个没有零的平滑函数，而贝塞尔函数中的 jn(u) 是具有无穷大零的，但出现以下错误： TypeError: Cannot cast array data from dtype('O') to dtype('float64') according to the rule 'safe' 首先我需要将变量 x 更改为变量 u 并在新变... »

我想解决 Python 中多精度的广义特征值问题（ A.C = (lam).B.C 其中 A 和 B 是 3000x3000 矩阵，C 是 3000x1 向量。(lam) 是特征值。）。所以我已经安装了 MPMATH。但我找不到用于广义特征值问题的 Python 程序。如何在 Python 中解决多精度的广义特征值问题？... »

我是 sympy 的新手，但对 python 仍然很幼稚……我想解一个三角方程，找到它的零点。 （一旦有了语法，我就会使用更复杂的函数。） 我还找不到正确的语法。以下是我在 Spyder (Python 2.7) 的 iPython 控制台上尝试的内容： from sympy.solvers import solve from sympy import Symbol x = Symbol('x... »

我发现自己想使用作为 mpmath 包的一部分提供的近似值，但对它们应该做什么感到困惑： http://docs.sympy.org/dev/modules/mpmath/calculus/approximation.html sympy 表达式和 sympy.mpmath 表达式之间到底有什么区别？ 如果我想要一个符号表达式的泰勒近似值而不了解 mpmath 包在做什么，我可以执行以下操作... »

我想使用 mpmath.quad 来计算一个通过 sympy 简化的函数的积分值。 然后我用下面的代码 from sympy import * import mpmath as mp r = symbols("r") f = 4**(-r)*r*(0.115391565000863*2**r + 693.147180559945)*exp(-0.000332949677174323*2**r)... »

我希望这个问题不是太初级。我已经广泛搜索了解决方案，但还没有找到。 我最近开始使用带有 Sympy 的 Jupyter Notebook 来记笔记并在我的微积分 II 课上做作业（这真是一个巨大的好处！）。 但是，我唯一的问题是我无法弄清楚如何配置绘图的大小（即像素尺寸）。 直接使用matplotlib 很容易（特别是matplotlib.pyplot.figure()），但我使用Sympy... »

这是我第一次尝试将 JIT 用于 python，这是我想要加快速度的用例。我读了一些关于 numba 的内容，它看起来很简单，但是下面的代码没有提供任何加速。请原谅我可能犯的任何明显错误。 我也尝试按照 cython 的基本教程的建议进行操作，但在时间上也没有区别。 http://docs.cython.org/src/tutorial/cython_tutorial.html 我猜我必须做一... »

我想使用sympy中的jacobDN功能，所以我下载了它并python setup.py install它，成功了。 当我想像documentation 那样使用它时： >>> from sympy.mpmath import * Traceback (most recent call last): File "<stdin>", line 1, in &lt... »

我有以下代码，我需要在其中求解表达式以找到根。需要求解omega的表达式。 import numpy as np from sympy import Symbol,lambdify import scipy from mpmath import findroot, exp eta = 1.5 tau = 5 /1000 omega = Symbol("omega") Tf = exp(1j * ... »

我可以使用 scipy quad 和 nquad 进行涉及可变积分限制的四重积分。问题是，当无法达到所要求的容差时，使用的默认精度会引发错误。使用 mpmath 积分器，我可以通过设置 mp.dps = 任意来定义任意精度，但是我看不到限制是否以及如何像使用 nquad 一样变得可变。 Mpmath 还在 quadgl 中使用 Gauss-Legendre 方法提供了非常快速的执行，这是非常可取的... »

我是 SymPy 的新手，希望有人能指出优化我的代码的方法。 我需要对一个包含非常高小数位 (150–300) 的表达式进行数值计算，并且每个参数集需要 30 秒或更长时间——考虑到要计算的参数空间，这非常长。 我在积分处理中使用了lambdify 和mpmath 后端和meijerg=True，它显着降低了运行时间。还有其他可以使用的方法吗？理想情况下，将评估时间推到 1 秒以下会很棒。我的代码... »

我的 Windows 10 机器上安装了 2 个版本的 python。 3.7 和 3.4。我需要为我的旧 python 发行版 3.4 安装一些包。对于某些软件包，有一个可用的可执行文件，我可以安装 matplotlib 和 numpy 之类的软件包。对于Scipy 和mpmath，我找不到相同的。也在尝试： C:\Python34> ./scripts/pip.exe install s... »

我有 4 个输入变量（浮点数）： Xmax Xmin 百分比 模式 我想为 s 求解以下（相当长的）方程： > (1/2+1/2*erf((log(Xmax)-(log(mode)+s**2))/(sqrt(2)*s))-(1/2+1/2*erf((log(Xmin)-(log(mode)+s**2))/(sqrt(2)*s))) - percentage == 0 我想使用 m... »

我正在尝试使用 mpmath.findroot 的多维牛顿法以高精度数值求解方程组。这是一个示例系统： def f(x_0, x_1, x_2, x_3, x_4, x_5, y_0, y_1, y_2, y_3, y_4, y_5, l_0, l_1,l_2, l_3, l_4, l_5, l_6, l_7, l_8, l_9, l_10, l_11, l_12, l_13, l_14): ... »

我有一些计算涉及快速爆炸的阶乘，因此我决定使用任意精度库mpmath。 我的代码如下所示： import numpy as np import mpmath as mp import time a = np.linspace( 0, 100e-2, 100 ) b = np.linspace( 0, np.pi ) c = np.arange( 30 ) t = t... »

我有： import numpy as np from mpmath import * mpf(np.array(range(0,600))) 但它不会让我这样做： TypeError: cannot create mpf from array 那我该怎么办？ 基本上我将使用这个数组，并根据情况（例如1.35626567e1084 或6.2345252e-2732）将元素乘以令人难... »