pith. sign in

arxiv: 2605.22427 · v3 · pith:UCV6JWMWnew · submitted 2026-05-21 · 💱 q-fin.CP · q-fin.PR

Faster Monotone Implied Volatility Solver

Fabien Le Floc'h This is my paper

Pith reviewed 2026-06-30 16:28 UTC · model grok-4.3

classification 💱 q-fin.CP q-fin.PR
keywords implied volatilityBlack-Scholesnumerical solvermonotone convergenceEuler-Chebyshev iterationHalley methodoption pricing
0
0 comments X

The pith

ThiopheneIV inverts Black-Scholes implied volatility with a proven monotone core using a lower-bound seed plus six targeted iterations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents ThiopheneIV as a Black-Scholes implied-volatility solver built around a monotone iteration scheme. It begins with the Choi-Huh-Su L3 lower-bound seed, then applies three Euler-Chebyshev steps on the lower branch and three Halley steps on the upper branch. The authors prove that in exact arithmetic the seed always lies below the root and both iteration maps increase monotonically without overshooting. Practical double-precision code must add parity normalisation, Bachelier-limit handling, saturated-price treatment, and an optional Jäckel-Newton polish to keep the observed accuracy. The work shows close agreement with multiprecision references across standard and adversarial test cases while remaining faster than several recent alternatives.

Core claim

ThiopheneIV starts from the Choi-Huh-Su L3 lower-bound seed and applies three Euler-Chebyshev steps followed by three Halley steps; the paper proves that in exact arithmetic the seed lies below the root and both maps increase monotonically without overshooting, while the double-precision implementation adds explicit guards for parity, Bachelier limits, saturation, and optional polishing to match multiprecision Black reference prices at low latency.

What carries the argument

The monotone core formed by the Choi-Huh-Su L3 lower-bound seed together with the Euler-Chebyshev lower-branch map and the Halley upper-branch map, whose monotonicity without overshoot is proved in exact arithmetic.

If this is right

  • The solver matches multiprecision Black reference prices across standard grids, market-like data, high-volatility regimes, and adversarial corners.
  • It runs at low latency while providing explicit production guards for boundary cases.
  • Direct comparisons place its accuracy and speed against Jäckel's Let's Be Rational and other recent solvers.
  • A convergence proof supplies a clean iteration core, yet robust production use still requires careful boundary handling and choice of pricing objective.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The separation between the exact-arithmetic proof and the double-precision guards implies that similar monotone schemes for other pricing models would also need dedicated boundary logic before deployment.
  • The optional Jäckel-Newton polish suggests that the core six-step method may be tuned further for speed-accuracy trade-offs in specific market regimes.
  • Emphasis on matching a chosen pricing objective rather than a universal reference points to potential differences in solver behaviour when pricing the same option under different conventions.

Load-bearing premise

The double-precision implementation with its parity normalisation, Bachelier-limit handling, saturated-price treatment, and optional polish continues to respect the exact-arithmetic monotonicity and accuracy properties.

What would settle it

A single test point where the solver returns an implied volatility that, when substituted back into the Black-Scholes formula, produces a price differing from the input price by more than a few units in the last place of double precision.

Figures

Figures reproduced from arXiv: 2605.22427 by Fabien Le Floc'h.

Figure 1
Figure 1. Figure 1: One-ulp price monotonicity diagnostic for 𝐹∗∕𝐾∗ = 0.955 and 𝑣 = 0.05. Each curve shows the local OTM-price step produced by advancing the total volatility by one double-precision value. Negative points mark local downward steps despite the monotonicity of the real-valued Black price. 3 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 1
Figure 1. Figure 1: visualises two of these one-ulp sweeps in total volatility. In these practical cases, the beta￾normalised erf cx/log formula improves on the textbook CDF path, while the branched reference remains cleaner still. 0 100 200 300 400 500 One-ulp step index 100 101 102 103 104 Absolute price error (ulps) CDF tail cancellation ∗∕∗ = 0.925, = 2, = 15% 0 100 200 300 400 500 One-ulp step index Commons erfcx improve… view at source ↗
Figure 2
Figure 2. Figure 2: Fixed-𝑥 comparison of ThiopheneIV and ThiopheneIV+ against Jäckel’s implementation. Each panel varies the reference total volatility at fixed log-moneyness and reports relative total-volatility error. by rounding, objective evaluation mismatch, and ulp quantisation; vega maps the remaining few-volatility-ulp bias into sub-ulp price changes, so the repriced residual can look mean-zero. 5.3. Operation Costs … view at source ↗
Figure 2
Figure 2. Figure 2: One-ulp price-step diagnostic for 𝐹∗∕𝐾∗ = 0.955 and 𝑣 = 0.05. Each curve shows the local OTM-price step produced by advancing total volatility by one double-precision value. Negative points mark local downward steps despite the monotonicity of the real-valued Black price [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Fixed-𝑥 convergence slices for the Choi L3 seed and the Euler–Chebyshev chain. Each panel fixes log-moneyness and varies the reference total volatility. Curves show relative total-volatility error after the seed and after one, two, and three full-precision Euler–Chebyshev steps; values are floored at 10−18 only for log-scale plotting. 0 100 200 300 400 500 One-ulp step index −150 −125 −100 −75 −50 −25 0 25… view at source ↗
Figure 3
Figure 3. Figure 3: summarises the input-domain branch map. The guards detailed in Appendix B handle adversarial inputs rather than only smooth benchmark grids. 10−14 10−12 10−10 10−8 10−6 10−4 10−2 1 Distance from ATM = − (log scale; left edge approximates = 0) 10−12 10−9 10−6 10−3 1∕2 0.999 1 − 10−6 OTM-normalised price (logit scale) lower-tail objective Choi L3 seed ln() 3 Euler–Chebyshev steps rst guard tiny near-ATM Bach… view at source ↗
Figure 4
Figure 4. Figure 4: Local round-trip diagnostics for the optional Jäckel–Newton polish. The left panel measures implied-volatility error when the input price is generated by the expanded Jäckel reference price; the right panel measures the relative price residual after inversion and repricing with the same reference formula. 9 [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 4
Figure 4. Figure 4: Fixed-𝑥 comparison of ThiopheneIV and ThiopheneIV+ against Jäckel’s implementation. Each panel varies the reference total volatility at fixed log-moneyness and reports relative total-volatility error. 5.2. Jäckel–Newton Polish Diagnostic The cutoff 𝑐 ≤ 1∕2 used by ThiopheneIV+ is not the boundary of Jäckel’s expanded price-evaluation regions. Let 𝑐J denote the Jäckel expansion threshold derived in Appendix… view at source ↗
Figure 5
Figure 5. Figure 5: Fixed-𝑥 convergence slices for the Choi L3 seed and the ThiopheneIV refinement chain. Each panel fixes log-moneyness and varies the reference total volatility. Curves show relative total-volatility error after the seed and after one, two, and three full-precision refinement steps; values are floored at 10−18 only for log-scale plotting. Thus the polish above 𝑐J is mainly a choice to align with Jäckel’s dou… view at source ↗
Figure 6
Figure 6. Figure 6: Local round-trip diagnostics for the optional Jäckel–Newton polish. The left panel measures implied-volatility error when the input price is generated by the expanded Jäckel reference price; the right panel measures the relative price residual after inversion and repricing with the same reference formula [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
read the original abstract

We present ThiopheneIV, a Black-Scholes implied-volatility solver with a monotone core and explicit production guards. The solver starts from the simple Choi-Huh-Su L3 lower-bound seed and applies three Euler-Chebyshev steps on a lower branch and three Halley steps on the remaining upper branch. We prove that, in exact arithmetic, the seed lies below the root and both maps increase monotonically without overshooting. We also detail the practical challenges encountered for a double-precision implementation: parity normalisation, microscopic Bachelier-limit handling, saturated price treatment, and an optional J\"ackel-Newton polish. Across standard grids, market-like data, high-volatility cases, and adversarial corners, ThiopheneIV agrees closely with multiprecision Black reference prices at low latency. We provide detailed comparisons with recent solvers, including J\"ackel's Let's Be Rational. The broader lesson is that a convergence proof gives a clean core, but robust production inversion still depends on boundary handling and on the pricing objective one chooses to match.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The paper presents ThiopheneIV, a Black-Scholes implied-volatility solver whose core begins with the Choi-Huh-Su L3 lower-bound seed and applies three Euler-Chebyshev steps on the lower branch followed by three Halley steps on the upper branch. It proves that, in exact arithmetic, the seed lies below the root and both iteration maps are strictly increasing and never overshoot. The manuscript also describes the additional machinery required for a double-precision implementation (parity normalisation, microscopic Bachelier-limit handling, saturated-price treatment, and an optional Jäckel-Newton polish) and reports close agreement with multiprecision reference prices across standard grids, market data, high-volatility regimes, and adversarial corners, together with comparisons against Jäckel’s Let’s Be Rational solver.

Significance. If the exact-arithmetic proof is correct and the floating-point implementation can be shown to inherit the monotonicity and no-overshoot guarantees, the work supplies a theoretically grounded, low-latency solver that is directly usable in production pricing engines. The explicit discussion of the boundary-handling interventions needed to move from the proved regime to double precision is itself a useful contribution to the literature on robust numerical inversion.

major comments (1)
  1. [double-precision implementation section] The section describing the double-precision implementation and production guards: the monotonicity and no-overshoot properties are proved only in exact arithmetic, yet the production solver introduces parity normalisation, Bachelier-limit handling, saturated-price treatment, and an optional Jäckel-Newton polish. No separate argument or exhaustive floating-point analysis is supplied showing that these interventions preserve the proved properties once rounding occurs. Because the headline accuracy and performance claims rest on the implemented solver behaving as the exact proof predicts, this transfer constitutes a load-bearing gap.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for highlighting this important point regarding the distinction between the exact-arithmetic proof and the floating-point implementation. We provide a point-by-point response below.

read point-by-point responses

  1. Referee: [double-precision implementation section] The section describing the double-precision implementation and production guards: the monotonicity and no-overshoot properties are proved only in exact arithmetic, yet the production solver introduces parity normalisation, Bachelier-limit handling, saturated-price treatment, and an optional Jäckel-Newton polish. No separate argument or exhaustive floating-point analysis is supplied showing that these interventions preserve the proved properties once rounding occurs. Because the headline accuracy and performance claims rest on the implemented solver behaving as the exact proof predicts, this transfer constitutes a load-bearing gap.

    Authors: We acknowledge that the monotonicity and no-overshoot properties are proved exclusively in exact arithmetic. The double-precision implementation includes several production guards—parity normalisation, Bachelier-limit handling, saturated-price treatment, and an optional Jäckel-Newton polish—to address numerical issues that arise in IEEE 754 floating-point arithmetic, such as underflow, overflow, and loss of precision near boundaries. These guards are necessary to ensure the solver remains robust and accurate in practice. We do not supply a formal floating-point analysis demonstrating that these interventions preserve the exact-arithmetic properties under all rounding conditions. Instead, the accuracy of the implemented solver is validated through extensive empirical testing against multiprecision reference Black-Scholes prices across a variety of test cases, including standard grids, market data, high-volatility regimes, and adversarial corners. In these tests, the solver exhibits the expected behavior without observed instances of non-monotonicity or overshooting. We agree that explicitly noting this distinction would improve the manuscript. We will revise the double-precision implementation section to clarify that the theoretical guarantees apply to the exact-arithmetic core, while the floating-point version relies on empirical validation for its practical performance and accuracy claims. The performance claims are supported by direct timing measurements, independent of the proof. revision: partial

Circularity Check

0 steps flagged

No significant circularity; mathematical proof stands independently of implementation details

full rationale

The paper's core contribution is an explicit proof in exact arithmetic that the Choi-Huh-Su L3 seed lies below the root and that the three Euler-Chebyshev steps (lower branch) plus three Halley steps (upper branch) are strictly increasing and non-overshooting. This proof is presented as a self-contained mathematical argument and does not rely on any fitted parameters renamed as predictions, self-citations that bear the central load, or definitions that smuggle the result into the inputs. The double-precision guards (parity normalisation, Bachelier-limit handling, saturated-price treatment, optional Jäckel-Newton polish) are described separately as practical additions outside the proved regime; they do not retroactively alter or define the exact-arithmetic claim. The manuscript further anchors its accuracy claims to external multiprecision Black reference prices and comparisons with other solvers, keeping the derivation independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review based solely on abstract; no explicit free parameters, axioms, or invented entities are stated beyond the implicit domain assumption that the Black-Scholes price is strictly monotone in volatility.

axioms (1)
  • domain assumption Black-Scholes call/put price is strictly increasing in volatility for fixed strike and maturity
    Required for the root to be unique and for the monotonicity proof to apply; implicit in any implied-volatility solver.

pith-pipeline@v0.9.1-grok · 5699 in / 1279 out tokens · 39379 ms · 2026-06-30T16:28:57.342503+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. A Practical Guide to Strip Caplet Volatilities

    q-fin.CP 2026-05 unverdicted novelty 6.0

    Interpolation scheme and node placement drive instability in caplet stripping; continuous flat-linear kernels, midpoint nodes, and exponential reparametrization or Hyman splines deliver stable positive curves with neg...

  2. A Practical Guide to Strip Caplet Volatilities

    q-fin.CP 2026-05 unverdicted novelty 5.0

    Constructive algorithm for arbitrage-free caplet volatility stripping via time-value interpolation, data outlier correction, and new compact-kernel bootstrap interpolants.

  3. Implying Volatility: How Fast Can We Go?

    q-fin.CP 2026-05 unverdicted novelty 4.0

    FlashIV is a new Black-Scholes implied volatility solver using input normalization, erfcx residual, and fixed Householder refinement that runs faster than Jäckel's Let's Be Rational while staying close to its reference price.

Reference graph

Works this paper leans on

12 extracted references · 6 canonical work pages · cited by 2 Pith papers

  1. [1]

    The Pricing of Commodity Contracts.J

    Black, F. The Pricing of Commodity Contracts.J. Financ. Econ.1976,3, 167–179. https://doi.org/10.1016/0304- 405X(76)90024-6. 25

    work page doi:10.1016/0304- 1976

  2. [2]

    https://doi.org/10.1002/fut.70024

    Cui,Z.;Liu,Y.;Yao,Y.TighterBoundsforImpliedVolatilityWiththeDiracDeltaFamilyMethod.J.FuturesMarkets 2025,45,e70024. https://doi.org/10.1002/fut.70024

    work page doi:10.1002/fut.70024 2025

  3. [3]

    Let's Be Rational

    Jäckel,P.Let’sBeRational.Wilmott2017,2015,40–53. https://doi.org/10.1002/wilm.10395

    work page doi:10.1002/wilm.10395 2015

  4. [4]

    Availableonline: https://chasethedevil

    LeFloc’h,F.MonotonicityoftheBlack–ScholesOptionPricesinPractice. Availableonline: https://chasethedevil. github.io/post/vol_monotonicity_in_practice/(accessedon19May2026)

  5. [5]

    Tighter ‘Uniform Bounds for Black–Scholes Implied Volatility’ and the Applications to Root-Finding.Oper

    Choi, J.; Huh, J.; Su, N. Tighter ‘Uniform Bounds for Black–Scholes Implied Volatility’ and the Applications to Root-Finding.Oper. Res. Lett.2024,57,107189. https://doi.org/10.1016/j.orl.2024.107189

    work page doi:10.1016/j.orl.2024.107189 2024

  6. [6]

    Birnbaum,Z.W.AnInequalityforMill’sRatio.Ann. Math. Statist.1942,13,245–246. https://doi.org/10.1214/aoms/ 1177731611

    work page doi:10.1214/aoms/ 1942

  7. [7]

    On Information and Sufficiency,

    Sampford, M.R. Some Inequalities on Mill’s Ratio and Related Functions.Ann. Math. Statist.1953,24, 130–132. https://doi.org/10.1214/aoms/1177729093

    work page doi:10.1214/aoms/1177729093 1953

  8. [8]

    Abramowitz,M.;Stegun,I.A.HandbookofMathematicalFunctions;NationalBureauofStandards,AppliedMathe- maticsSeries55: Washington,DC,USA,1964

    1964

  9. [9]

    Traub,J.F.IterativeMethodsfortheSolutionofEquations;ChelseaPublishingCompany: NewYork,NY,USA,1982

    1982

  10. [10]

    Householder,A.S.TheNumericalTreatmentofaSingleNonlinearEquation;McGraw-Hill: NewYork,NY,USA,1970

    1970

  11. [11]

    Healy,J.AppliedQuantitativeFinanceforEquityDerivatives,4thed.;self-published,2024

    2024

  12. [12]

    Availableonline: https://ssrn.com/abstract=2420757 (accessedon16May2026)

    LeFloc’h,F.FastandAccurateAnalyticBasisPointVolatility. Availableonline: https://ssrn.com/abstract=2420757 (accessedon16May2026). 26