TQBTHEQUANTBATEMAN
TQB/ learn/ numerical finance/ variance reduction convergenceEN · DARK
Numerical finance · advanced

Variance reduction & convergence

Spend simulation budget on information, not noise.

BY THE END, YOU CAN

01Use antithetic and control variates

02Measure variance-reduction efficiency

03Build a convergence table

01
INTUITION

Name the mathematical target and the approximation error.

Variance reduction changes the estimator, not the target expectation; bias and variance must remain separately visible.

01

Negative payoff correlation helps antithetics.

02

A good control has known mean and high correlation.

03

Efficiency includes runtime, not variance alone.

02
WHY MARKETS CARE

Tie accuracy to the decision the number supports.

Stable prices and Greeks under service latency limits require better estimators before more hardware.

INSTRUMENTS

path-dependent options

Monte Carlo Greeks

exposure engines

QUOTE CONVENTION

All estimators share paths, seed and target; efficiency is variance × runtime.

03
MATHEMATICS

Separate estimator, discretization, truncation, and convergence.

Formula · Short derivation

Control variate

Ycv=Y−β(C−E[C])Y^{cv}=Y-\beta(C-\mathbb E[C])

The optimal linear coefficient removes correlated noise.

Formula · Short derivation

Optimal coefficient

β∗=Cov⁡(Y,C)Var⁡(C)\beta^*=\frac{\operatorname{Cov}(Y,C)}{\operatorname{Var}(C)}

Estimate β on a pilot or independent sample to avoid hidden bias.

Short derivation
Short derivation

From information set to computable quantity

Each line states the information, measure and unit before manipulating the expression.

  1. 01

    Preserve the mean

    Subtract a zero-mean control adjustment.

    E[C−E[C]]=0E[C-E[C]]=0
  2. 02

    Expand variance

    Write variance as a quadratic in β.

    Var(Y−βC)=Var(Y)+β2Var(C)−2βCov(Y,C)Var(Y-\beta C)=Var(Y)+\beta^2Var(C)-2\beta Cov(Y,C)
  3. 03

    Minimize

    Differentiate with respect to β and solve.

    ∂βVar=0\partial_\beta Var=0
  4. 04

    Benchmark efficiency

    Compare confidence-width squared times runtime at common target accuracy.

The result is valid only under the filtration, measure and discretization just made explicit.

Inputs
  • Yᵃ: antithetic payoff
  • C: control with known E[C]
Assumptions and limits
  • A poor control can increase variance.
  • Data-dependent tuning can bias reported results.
05
MODEL / PRICING

Benchmark the algorithm against a controlled reference.

METHOD

Pair common shocks, estimate correlations, apply antithetic/control estimators and benchmark variance per unit runtime.

CALIBRATION

Choose controls from analytically priced neighbours, not fitted to the same noisy final sample.

06PYTHON IMPLEMENTATIONOpen the implementation and checks.
ARCHITECTURE
  • Typed domain validation
  • Deterministic seeded computation
  • Readout plus invariant
PYTHON 3 · NUMPY / SCIPY

Variance reduction & convergence

Reproduce the governing quantity, then challenge it with an invariant.

REUSABLE EXAMPLE
01import numpy as np
02
03def ycvybetacmathbbec(x: np.ndarray) -> float:
04 x = np.asarray(x, dtype=float)
05 assert np.isfinite(x).all()
06 return float(np.mean(x))
07
08sample = np.array([0.8, 1.0, 1.2])
09value = ycvybetacmathbbec(sample)
10assert sample.min() <= value <= sample.max()
11print(f"value={value:.6f}")
EXPECTED OUTPUTvalue=1.000000
SANITY CHECKS

✓ Finite inputs are enforced

✓ The result respects its numerical bounds

✓ Units and measure remain explicit

07
INTERACTIVE LAB

Change the error budget, not just the picture.

NUMERICAL ERROR LAB

Variance reduction & convergence

Use a fixed seed to separate discretization, sampling and truncation effects.

SYNTHETIC · EDUCATIONAL
Plain SE0.2108N=4,096
Antithetic factor0.213ρ=-0.575
Reduced SE0.0972same path budget
spot level by time (years)

Exact, Euler and Milstein spot paths using the same shocks.

  • exact
  • Euler
  • Milstein
time (years): 0.0Y. exact: 100.00. Euler: 100.00. Milstein: 100.00.

Use Left/Right or Up/Down arrows to inspect values; Home and End jump to the bounds.

ERROR BUDGET

Model → algorithm → evidence

Each numerical approximation has an independent control and a reference diagnostic.

01Target lawQ / GBM

Measure and SDE fixed

02Algorithm32 steps

Discretization or transform control

03DiagnosticSE=0.211

Visible convergence evidence

MODEL BOUNDARY

GBM and Gaussian transform examples are analytical references, not production calibration engines.

08
FRONT OFFICE

Where the model meets the book.

ON THE DESK
“Ten times more paths is a procurement decision; ten times less variance is a quant decision.”
VISIBLE INPUTS

pilot correlation

known control mean

CALIBRATION

Choose controls from analytically priced neighbours, not fitted to the same noisy final sample.

RISK

estimator bias

unstable β

DAILY WORKFLOW
  1. Validate market state and timestamp
  2. Recompute the baseline
  3. Run a controlled perturbation
  4. Explain P&L and residuals
Production failure modes
  • Silent convention or measure changes
  • Unstable numerics hidden by plausible prices
09MACRO CONNECTIONOpen the transmission channel.
MACRO CONNECTION

Transmission from state to valuation

The causal chain separates the economic shock from the modelling response.

01Shared shocktransmits

creates correlation

02Controltransmits

removes explainable noise

03Capacityoutput

delivers tighter intraday risk

10COMMON PITFALLSOpen the failure checklist.
01

Fitting and evaluating β on tiny samples

02

Calling variance reduction a model improvement

11SOURCES / FURTHER READINGOpen sources and continue the track.
research

Measure theory, simulation and computational-finance lectures

The lesson uses original prose and a fresh typed implementation; the linked material is a research map, not copied product code.

Source
Computational Finance Course
Author
L. A. Grzelak
Ref
main
OPEN ORIGINAL SOURCE ↗