Quantlib python Heston 模型：生成路径，得到 "Boost assertion failed: px != 0"

Question

我正在尝试使用 GaussianPathGenerator 和 Quantlib 中的 HestonProcess 生成底层路径 python。但它给了我“RuntimeError：Boost assertion failed：px！= 0”。

如果我用 BS 或 HullWhite 进程替换 Heston 进程，路径会生成正常。有谁知道为什么 Heston 那个不工作？

非常感谢您的帮助！

将 QuantLib 导入为 ql

today = ql.Date(21, 10, 2020)
daycount = ql.Actual360()
calendar = ql.TARGET()
r_ts = ql.YieldTermStructureHandle(ql.FlatForward(today, 0.02, daycount))
d_ts = ql.YieldTermStructureHandle(ql.FlatForward(today, 0.01, daycount))
v0, kappa, theta, sigma, rho = 0.03, 10, 0.03, 0.4, 0.3
hs_process = ql.HestonProcess(r_ts, d_ts, ql.QuoteHandle(ql.SimpleQuote(100)), v0, kappa, theta, sigma, rho)

bs_process = ql.BlackScholesProcess(ql.QuoteHandle(ql.SimpleQuote(100)), r_ts, 
                                    ql.BlackVolTermStructureHandle(ql.BlackConstantVol(0, calendar, 0.20, daycount)))

timestep = 50
length = 1.0
rng = ql.GaussianRandomSequenceGenerator(ql.UniformRandomSequenceGenerator(timestep, ql.UniformRandomGenerator()))
seq = ql.GaussianPathGenerator(hs_process, length, timestep, rng, False)

path_1 = seq.next()

RuntimeError                              Traceback (most recent call last)
<ipython-input-23-e8c522a62dd1> in <module>
----> 1 path1 = seq.next()

~\Anaconda3\lib\site-packages\QuantLib\QuantLib.py in next(self)
  22328 
  22329     def next(self):
> 22330         return _QuantLib.GaussianPathGenerator_next(self)
  22331 
  22332     def antithetic(self):

RuntimeError: Boost assertion failed: px != 0

Answer 1

Heston 过程是二维的（它演化了底层和的波动性），因此需要将其传递给 GaussianMultiPathGenerator。不幸的是，包装器中的 SWIG 机制没有捕捉到类型不匹配；它只是试图将过程转换为一维过程，当转换失败时会导致空指针。

出于同样的原因，您必须使用 2*timesteps 初始化 GaussianRandomSequenceGenerator，因为它必须为过程中的两个变量提供足够的随机数。

总而言之，您的代码应为：

rng = ql.GaussianRandomSequenceGenerator(ql.UniformRandomSequenceGenerator(2*timestep, ql.UniformRandomGenerator()))
times = list(ql.TimeGrid(length, timestep))
seq = ql.GaussianMultiPathGenerator(hs_process, times, rng)

调用 path = seq.next() 将 return 一个 Multipath 实例； path.value()[0] 将是底层证券的路径，path.value()[1] 是其波动率的路径。

Quantlib python Heston 模型：生成路径，得到 "Boost assertion failed: px != 0"

Quantlib python Heston model: generate path, get "Boost assertion failed: px != 0"

python

montecarlo

quantlib