问题标签 [quad]

我需要使用 ctypes 函数来减少 python 中 quad 的运行时间。这是我最初的问题original question，但现在我知道我需要遵循什么路径。我需要按照与此处类似问题链接中相同的步骤进行操作。

但是，在我的情况下，将在数值积分中处理的函数是调用另一个 python 函数。像这样：

其中 cut_off 只是我在代码中决定的一个标量，而 kde_bad 是使用 KernelDensity 创建的内核对象。

所以我的问题是我需要如何在 C 中指定函数？相当于：

任何输入表示赞赏！

因此，我正在开发一种表观温度计算工具，以集成到具有许多不同功能的数据分析程序中。我的“测试代码”在我的主程序之外运行良好，并给出了正确的结果。有效的代码：

但是，当我将此代码集成到我的主 tkinter GUI 循环中时，我定义的被积函数仅返回 0.0 的值。不起作用的代码：

我认为这可能是一个四舍五入的问题，但使用 bigfloat 没有效果。我真的很困惑为什么它会在我的主代码之外工作，然后在集成时以这种方式失败。我一直在搜索我的 tkinter gui 窗口代码的其余部分，以寻找任何可能影响它的东西，但我很困惑。任何帮助，甚至猜测，都非常感谢。如果需要，我可以发布更多我的主要代码，但老实说，这是一团糟，因为我一直在使用 pyinstaller 将其编译为 .exe，并且它需要在同一个文件中定义所有非模块函数。为了回答@jsbueno，我对 rad_t_conv 的调用来自我主窗口中的一个按钮。

类定义如下：

然后称为：

合并函数为：

我现在很困惑我写的一个小函数。我想将以下函数嵌入到我的程序中：

其中是泊松分布的累积分布函数的卷积。

如果我做对了，我会像这样计算卷积，其中 k 是间隔中的事件数，Mu 只是平均值，例如 5，j 是卷积数：

函数 D 如下所示：

这是在 python 中完成的，如下所示：

现在我遇到的棘手部分是，如何让积分与 cdf 一起工作。目前我正在尝试这样做，但我不知道如何继续：

如果有人知道如何做到这一点，我将不胜感激:)

编辑： 为了检查 cdf 的卷积是否正确，我这样做了：

我知道这里有一些类似的问题，但似乎都没有真正解决我的问题。

我的代码如下所示：

这些功能确实有效，我的问题：如何绘制 d_L 与 z ？就像我在 d_L 的定义中有这个积分函数显然是一个问题，它取决于 z 和一些 args=(omega_m, )。

我有一个向量，并希望制作另一个相同长度的向量，其第 k 个分量是

问题是：我们如何将其矢量化以提高速度？NumPy vectorize() 实际上是一个for循环，所以不算。

Veedrac 指出“如果不多次调用它，就无法将纯 Python 函数应用于 NumPy 数组的每个元素”。由于我使用的是 NumPy 函数而不是“纯 Python”函数，我想可以进行矢量化，但我不知道如何。

附带说明一下，Cython for 循环是否总是比 NumPy 向量化更快？

我不断收到错误消息“提供的函数没有返回有效的浮点数。” ，我认为问题在于我尝试将数组传递给 quad 函数。我考虑过用类似的方法评估数组每个条目的积分：

它不起作用，但它是正确的方法吗？还有其他/更好的建议吗？

我想将这段 Mathematica 代码重写为 Python，但当时我很困惑，请帮助我！非常感谢。也许集成的结果与数组不同，但是，我不知道！

图片：我想用python重写的mathematica代码

那时我对python代码感到困扰......

我想把积分结果画成如下图，结果好像是一张白纸，是什么原因？请帮我！

输出，但情节是空白的！

I am receiving a ValueError when integrating using scipy.integrate.quad. Here is my simplified code:

The last line gives the exception:

Can someone please propose how to do this integration properly? Thanks

我需要执行与 python 的集成，但其中一个限制是变量，而不是数字（从 0 到 z）。

我尝试了以下方法：

我遇到的问题是dl(z,y,a)给了我一个数组，所以每当我想绘制或评估它时，我都会获得以下信息：

ValueError：具有多个元素的数组的真值不明确。

我不知道是否有任何解决方案。在此先感谢