Valuing American options and Flexible Forwards contracts in time-dependent models

Battle Count: 72

The paper is primarily relevant to FX derivatives desks and structured products teams rather than algorithmic trading. It provides a fast and accurate pricing engine for Flexible Forward contracts (a common FX hedging instrument) and American options under stochastic volatility. The 1-2 second pricing speed makes it suitable for real-time risk management and trade pricing. However, it does not address trading strategy development, signal generation, or execution optimization. The early exercise boundary computation is critical for hedging and risk management of American-style positions. The volatility skew modeling is directly relevant to FX options trading desks.

The implementation requires: (1) recursive Riccati equation solvers with matrix propagators for the joint CF, (2) double cosine series expansions with careful truncation interval selection, (3) DCT-based expectation computation, (4) nonlinear Fredholm-Volterra integral equation solving via fixed-point iteration, (5) DSINC basis construction with conditional CF decomposition, (6) tilted CIR transform evaluation, (7) penalty-iteration ADI-MCS finite-difference benchmark, (8) handling of Feller condition violations, Gibbs oscillations, ITM options, and time-dependent parameters. The paper provides detailed algorithms but the interplay between multiple numerical methods, special functions, and boundary conditions makes this a highly complex implementation. The 432-case validation suite and multiple benchmark comparisons add further engineering overhead.

The paper provides detailed algorithmic pseudocode (Algorithms 1-3), complete parameter tables (Tables 1-6), explicit formulas for all components, and mentions that supporting Python code is available on GitHub. The model specification, Riccati recursion, COS/DSINC implementations, and FD benchmark are all described in sufficient detail for reproduction. However, the GitHub URL is not explicitly provided in the extract, and some implementation details (e.g., specific solver tolerances, grid construction parameters) require careful reading. The 432-case validation suite parameters are fully specified.