Questions tagged «mathematica»

3

1
替换C ++中的Mathematica的QuasiMonteCarlo集成
我有一个Mathematica程序,使用该QuasiMonteCarlo方法在3或4维中执行一些积分。问题是,要花费一些令人讨厌的时间,以致于某些计算无法在我们的HPC群集上提供的最大工作时间内完成。因此,我正在考虑用C ++重写该程序,我怀疑它将大大提高它的速度。 我查看了GSL文档,虽然有关于准随机序列和常规MC集成的部分,但没有任何东西可以将它们组合在一起。同样,一两个Google搜索并没有发现任何看起来像是广受信任的实现。在C ++中,经过良好测试的QMC集成实现方式有哪些选择? 为了保持一致性,如果可以选择,我宁愿使用Mathematica实现的类似于Halton-Hammersley-Wozniakowski方法的方法。

1
从实部数值恢复解析连续的虚部
我的情况。 我有一个通过复杂积分定义的复杂变量函数。我感兴趣的是此函数在虚轴上的值。我可以在以下功能区上对该函数进行数值访问:。形式上,积分表达式在该域之外是发散的,因此我需要解析继续。总结一下我的情况,ž = (X ,Ý )∈ (- ∞ ,∞ )× [ - 1 ,1 ]F(z)f(z)f(z)ž= (x ,y)∈ (- ∞ ,∞ )× [ - 1 ,1 ]z=(x,y)∈(−∞,∞)×[−1,1]z=(x,y)\in (-\infty,\infty)\times[-1,1] 这是我从数字上了解的功能区上的信息:F(z)f(z)f(z) 它同时围绕虚轴和实轴对称。 在处衰减为零。ř È (ž)→ ∞Re(z)→∞Re(z)\rightarrow\infty 它在附近爆炸。我不知道这可能是极点或分支点。我怀疑这种奇异性的性质(以及解析延续的所有其他孤立奇异性)取决于此函数的特定参数化(有关详细信息,请参见下面的积分)ξž= ± 我z=±iz=\pm iξξ\xi 实际上,绘制时它看起来与或非常相似。这是真实部分的图:1 /(1 + z 2 )2 n塞奇2(z)sech2(z)\text{sech}^2(z)1 /(1 + z2)2 n1/(1+z2)2n1/(1+z^2)^{2n} 我的问题是,鉴于我对该函数有大量的信息(在该功能区上对其进行总数值访问),我是否可以通过某种方式沿虚轴数值计算对该函数的近似值?我正在使用Mathematica。 我对虚轴上的值感兴趣的原因是,我需要评估此函数的以下傅立叶变换: …
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.