public void testMonteCarloLookback()
{
double tolerance = 0.1;
DayCounter dc = new Actual360();
Date today = Date.Today;
double strike = 90;
double t = 1;
double t1 = 0.25;
Date exDate = today + Convert.ToInt32(t * 360 + 0.5);
Exercise exercise = new EuropeanExercise(exDate);
SimpleQuote spot = new SimpleQuote(0.0);
SimpleQuote qRate = new SimpleQuote(0.0);
YieldTermStructure qTS = Utilities.flatRate(today, qRate, dc);
SimpleQuote rRate = new SimpleQuote(0.0);
YieldTermStructure rTS = Utilities.flatRate(today, rRate, dc);
SimpleQuote vol = new SimpleQuote(0.0);
BlackVolTermStructure volTS = Utilities.flatVol(today, vol, dc);
spot.setValue(100);
qRate.setValue(0);
rRate.setValue(0.06);
vol.setValue(0.1);
BlackScholesMertonProcess stochProcess =
new BlackScholesMertonProcess(
new Handle <Quote>(spot),
new Handle <YieldTermStructure>(qTS),
new Handle <YieldTermStructure>(rTS),
new Handle <BlackVolTermStructure>(volTS));
Option.Type[] types = new Option.Type[] { Option.Type.Call, Option.Type.Put };
for (int i = 0; i < types.Length; i++)
{
Option.Type type = types[i];
StrikedTypePayoff payoff = new PlainVanillaPayoff(type, strike);
/**
* Partial Fixed
* **/
Date lookbackStart = today + Convert.ToInt32(t1 * 360 + 0.5);
ContinuousPartialFixedLookbackOption partialFixedLookback
= new ContinuousPartialFixedLookbackOption(lookbackStart,
payoff,
exercise);
IPricingEngine engine =
new AnalyticContinuousPartialFixedLookbackEngine(stochProcess);
partialFixedLookback.setPricingEngine(engine);
double analytical = partialFixedLookback.NPV();
IPricingEngine mcpartialfixedengine =
new MakeMCLookbackEngine <ContinuousPartialFixedLookbackOption.Arguments, ContinuousPartialFixedLookbackOption.Results,
PseudoRandom, Statistics>(stochProcess)
.withSteps(2000)
.withAntitheticVariate()
.withSeed(1)
.withAbsoluteTolerance(tolerance)
.value();
partialFixedLookback.setPricingEngine(mcpartialfixedengine);
double monteCarlo = partialFixedLookback.NPV();
double diff = Math.Abs(analytical - monteCarlo);
if (diff > tolerance)
{
REPORT_FAILURE_MC("Partial Fixed", type, analytical, monteCarlo, tolerance);
}
/**
* Fixed
* **/
double minMax = 100;
ContinuousFixedLookbackOption fixedLookback
= new ContinuousFixedLookbackOption(minMax,
payoff,
exercise);
IPricingEngine analyticalfixedengine =
new AnalyticContinuousFixedLookbackEngine(stochProcess);
fixedLookback.setPricingEngine(analyticalfixedengine);
analytical = fixedLookback.NPV();
IPricingEngine mcfixedengine =
new MakeMCLookbackEngine <ContinuousFixedLookbackOption.Arguments, ContinuousFixedLookbackOption.Results,
PseudoRandom, Statistics>(stochProcess)
.withSteps(2000)
.withAntitheticVariate()
.withSeed(1)
.withAbsoluteTolerance(tolerance)
.value();
fixedLookback.setPricingEngine(mcfixedengine);
monteCarlo = fixedLookback.NPV();
diff = Math.Abs(analytical - monteCarlo);
if (diff > tolerance)
{
REPORT_FAILURE_MC("Fixed", type, analytical, monteCarlo, tolerance);
}
/**
* Partial Floating
* **/
double lambda = 1;
Date lookbackEnd = today + Convert.ToInt32(t1 * 360 + 0.5);
FloatingTypePayoff floatingPayoff = new FloatingTypePayoff(type);
ContinuousPartialFloatingLookbackOption partialFloating
= new ContinuousPartialFloatingLookbackOption(minMax,
lambda,
lookbackEnd,
floatingPayoff,
exercise);
IPricingEngine analyticalpartialFloatingengine =
new AnalyticContinuousPartialFloatingLookbackEngine(stochProcess);
partialFloating.setPricingEngine(analyticalpartialFloatingengine);
analytical = partialFloating.NPV();
IPricingEngine mcpartialfloatingengine =
new MakeMCLookbackEngine <ContinuousPartialFloatingLookbackOption.Arguments, ContinuousPartialFloatingLookbackOption.Results,
PseudoRandom, Statistics>
(stochProcess)
.withSteps(2000)
.withAntitheticVariate()
.withSeed(1)
.withAbsoluteTolerance(tolerance)
.value();
partialFloating.setPricingEngine(mcpartialfloatingengine);
monteCarlo = partialFloating.NPV();
diff = Math.Abs(analytical - monteCarlo);
if (diff > tolerance)
{
REPORT_FAILURE_MC("Partial Floating", type, analytical, monteCarlo, tolerance);
}
/**
* Floating
* **/
ContinuousFloatingLookbackOption floating =
new ContinuousFloatingLookbackOption(minMax,
floatingPayoff,
exercise);
IPricingEngine analyticalFloatingengine =
new AnalyticContinuousFloatingLookbackEngine(stochProcess);
floating.setPricingEngine(analyticalFloatingengine);
analytical = floating.NPV();
IPricingEngine mcfloatingengine =
new MakeMCLookbackEngine <ContinuousFloatingLookbackOption.Arguments, ContinuousFloatingLookbackOption.Results,
PseudoRandom, Statistics>
(stochProcess)
.withSteps(2000)
.withAntitheticVariate()
.withSeed(1)
.withAbsoluteTolerance(tolerance)
.value();
floating.setPricingEngine(mcfloatingengine);
monteCarlo = floating.NPV();
diff = Math.Abs(analytical - monteCarlo);
if (diff > tolerance)
{
REPORT_FAILURE_MC("Floating", type, analytical, monteCarlo, tolerance);
}
}
}
public void testAnalyticContinuousPartialFixedLookback()
{
// Testing analytic continuous fixed-strike lookback options
LookbackOptionData[] values =
{
// data from "Option Pricing Formulas, Second Edition", Haug, 2006, pg.148
//type, strike, minmax, s, q, r, t, v, l, t1, result, tol
new LookbackOptionData(Option.Type.Call, 90, 0, 100, 0, 0.06, 1, 0.1, 0, 0.25, 20.2845, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 90, 0, 100, 0, 0.06, 1, 0.1, 0, 0.5, 19.6239, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 90, 0, 100, 0, 0.06, 1, 0.1, 0, 0.75, 18.6244, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 110, 0, 100, 0, 0.06, 1, 0.1, 0, 0.25, 4.0432, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 110, 0, 100, 0, 0.06, 1, 0.1, 0, 0.5, 3.958, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 110, 0, 100, 0, 0.06, 1, 0.1, 0, 0.75, 3.7015, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 90, 0, 100, 0, 0.06, 1, 0.2, 0, 0.25, 27.5385, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 90, 0, 100, 0, 0.06, 1, 0.2, 0, 0.5, 25.8126, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 90, 0, 100, 0, 0.06, 1, 0.2, 0, 0.75, 23.4957, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 110, 0, 100, 0, 0.06, 1, 0.2, 0, 0.25, 11.4895, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 110, 0, 100, 0, 0.06, 1, 0.2, 0, 0.5, 10.8995, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 110, 0, 100, 0, 0.06, 1, 0.2, 0, 0.75, 9.8244, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 90, 0, 100, 0, 0.06, 1, 0.3, 0, 0.25, 35.4578, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 90, 0, 100, 0, 0.06, 1, 0.3, 0, 0.5, 32.7172, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 90, 0, 100, 0, 0.06, 1, 0.3, 0, 0.75, 29.1473, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 110, 0, 100, 0, 0.06, 1, 0.3, 0, 0.25, 19.725, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 110, 0, 100, 0, 0.06, 1, 0.3, 0, 0.5, 18.4025, 1.0e-4),
new LookbackOptionData(Option.Type.Call, 110, 0, 100, 0, 0.06, 1, 0.3, 0, 0.75, 16.2976, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 90, 0, 100, 0, 0.06, 1, 0.1, 0, 0.25, 0.4973, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 90, 0, 100, 0, 0.06, 1, 0.1, 0, 0.5, 0.4632, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 90, 0, 100, 0, 0.06, 1, 0.1, 0, 0.75, 0.3863, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 110, 0, 100, 0, 0.06, 1, 0.1, 0, 0.25, 12.6978, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 110, 0, 100, 0, 0.06, 1, 0.1, 0, 0.5, 10.9492, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 110, 0, 100, 0, 0.06, 1, 0.1, 0, 0.75, 9.1555, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 90, 0, 100, 0, 0.06, 1, 0.2, 0, 0.25, 4.5863, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 90, 0, 100, 0, 0.06, 1, 0.2, 0, 0.5, 4.1925, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 90, 0, 100, 0, 0.06, 1, 0.2, 0, 0.75, 3.5831, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 110, 0, 100, 0, 0.06, 1, 0.2, 0, 0.25, 19.0255, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 110, 0, 100, 0, 0.06, 1, 0.2, 0, 0.5, 16.9433, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 110, 0, 100, 0, 0.06, 1, 0.2, 0, 0.75, 14.6505, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 90, 0, 100, 0, 0.06, 1, 0.3, 0, 0.25, 9.9348, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 90, 0, 100, 0, 0.06, 1, 0.3, 0, 0.5, 9.1111, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 90, 0, 100, 0, 0.06, 1, 0.3, 0, 0.75, 7.9267, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 110, 0, 100, 0, 0.06, 1, 0.3, 0, 0.25, 25.2112, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 110, 0, 100, 0, 0.06, 1, 0.3, 0, 0.5, 22.8217, 1.0e-4),
new LookbackOptionData(Option.Type.Put, 110, 0, 100, 0, 0.06, 1, 0.3, 0, 0.75, 20.0566, 1.0e-4)
};
DayCounter dc = new Actual360();
Date today = Date.Today;
SimpleQuote spot = new SimpleQuote(0.0);
SimpleQuote qRate = new SimpleQuote(0.0);
YieldTermStructure qTS = Utilities.flatRate(today, qRate, dc);
SimpleQuote rRate = new SimpleQuote(0.0);
YieldTermStructure rTS = Utilities.flatRate(today, rRate, dc);
SimpleQuote vol = new SimpleQuote(0.0);
BlackVolTermStructure volTS = Utilities.flatVol(today, vol, dc);
for (int i = 0; i < values.Length; i++)
{
Date exDate = today + Convert.ToInt32(values[i].t * 360 + 0.5);
Exercise exercise = new EuropeanExercise(exDate);
spot.setValue(values[i].s);
qRate.setValue(values[i].q);
rRate.setValue(values[i].r);
vol.setValue(values[i].v);
StrikedTypePayoff payoff = new PlainVanillaPayoff(values[i].type, values[i].strike);
BlackScholesMertonProcess stochProcess = new BlackScholesMertonProcess(
new Handle <Quote>(spot),
new Handle <YieldTermStructure>(qTS),
new Handle <YieldTermStructure>(rTS),
new Handle <BlackVolTermStructure>(volTS));
IPricingEngine engine = new AnalyticContinuousPartialFixedLookbackEngine(stochProcess);
Date lookbackStart = today + Convert.ToInt32(values[i].t1 * 360 + 0.5);
ContinuousPartialFixedLookbackOption option = new ContinuousPartialFixedLookbackOption(lookbackStart,
payoff, exercise);
option.setPricingEngine(engine);
double calculated = option.NPV();
double expected = values[i].result;
double error = Math.Abs(calculated - expected);
if (error > values[i].tol)
{
REPORT_FAILURE_FIXED("value", values[i].minmax, payoff, exercise,
values[i].s, values[i].q, values[i].r, today,
values[i].v, expected, calculated, error,
values[i].tol);
}
}
}
