7得票1回答
使用scipy的quad例程对具有奇点的函数进行积分集成

我正在使用scipy.integrate v0.19.1中的quad函数来积分具有类似于根号的奇异性的函数,这些函数在积分区间的每个端点处都存在,例如: In [1]: quad(lambda x: 1/sqrt(1-x**2), -1, 1) 我使用来自numpy v1.12.0的sq...

7得票3回答
NumPy向量化与积分集成

我有一个向量 ,希望创建另一个同样长度的向量,其第k个分量为 问题是:我们如何进行向量化以提高速度? NumPy的vectorize()实际上是一个循环,所以不算。 Veedrac指出:" 没有办法将一个纯Python函数应用于NumPy数组的每个元素而不调用它那么多次 "。由于我使用...