public void testFlatTermVolatilityStripping2()
{
// Testing forward/forward vol stripping from flat term vol
// surface using OptionletStripper2 class...");
CommonVars vars = new CommonVars();
Settings.setEvaluationDate(Date.Today);
vars.setFlatTermVolCurve();
vars.setFlatTermVolSurface();
IborIndex iborIndex = new Euribor6M(vars.yieldTermStructure);
// optionletstripper1
OptionletStripper1 optionletStripper1 = new OptionletStripper1(vars.flatTermVolSurface,
iborIndex, null, vars.accuracy);
StrippedOptionletAdapter strippedOptionletAdapter1 = new StrippedOptionletAdapter(optionletStripper1);
Handle <OptionletVolatilityStructure> vol1 = new Handle <OptionletVolatilityStructure>(strippedOptionletAdapter1);
vol1.link.enableExtrapolation();
// optionletstripper2
OptionletStripper optionletStripper2 = new OptionletStripper2(optionletStripper1, vars.flatTermVolCurve);
StrippedOptionletAdapter strippedOptionletAdapter2 = new StrippedOptionletAdapter(optionletStripper2);
Handle <OptionletVolatilityStructure> vol2 = new Handle <OptionletVolatilityStructure>(strippedOptionletAdapter2);
vol2.link.enableExtrapolation();
// consistency check: diff(stripped vol1-stripped vol2)
for (int strikeIndex = 0; strikeIndex < vars.strikes.Count; ++strikeIndex)
{
for (int tenorIndex = 0; tenorIndex < vars.optionTenors.Count; ++tenorIndex)
{
double strippedVol1 = vol1.link.volatility(vars.optionTenors[tenorIndex], vars.strikes[strikeIndex], true);
double strippedVol2 = vol2.link.volatility(vars.optionTenors[tenorIndex], vars.strikes[strikeIndex], true);
// vol from flat vol surface (for comparison only)
double flatVol = vars.flatTermVolSurface.volatility(vars.optionTenors[tenorIndex], vars.strikes[strikeIndex], true);
double error = Math.Abs(strippedVol1 - strippedVol2);
if (error > vars.tolerance)
{
QAssert.Fail("\noption tenor: " + vars.optionTenors[tenorIndex] +
"\nstrike: " + vars.strikes[strikeIndex] +
"\nstripped vol1: " + strippedVol1 +
"\nstripped vol2: " + strippedVol2 +
"\nflat vol: " + flatVol +
"\nerror: " + error +
"\ntolerance: " + vars.tolerance);
}
}
}
}