void OptionletStripperTest::testFlatTermVolatilityStripping2() {
BOOST_TEST_MESSAGE(
"Testing forward/forward vol stripping from flat term vol "
"surface using OptionletStripper2 class...");
CommonVars vars;
Settings::instance().evaluationDate() = Date::todaysDate();
vars.setFlatTermVolCurve();
vars.setFlatTermVolSurface();
shared_ptr<IborIndex> iborIndex(new Euribor6M(vars.yieldTermStructure));
// optionletstripper1
shared_ptr<OptionletStripper1> optionletStripper1(new
OptionletStripper1(vars.flatTermVolSurface,
iborIndex,
Null<Rate>(),
vars.accuracy));
boost::shared_ptr<StrippedOptionletAdapter> strippedOptionletAdapter1(new
StrippedOptionletAdapter(optionletStripper1));
Handle<OptionletVolatilityStructure> vol1(strippedOptionletAdapter1);
vol1->enableExtrapolation();
// optionletstripper2
shared_ptr<OptionletStripper> optionletStripper2(new
OptionletStripper2(optionletStripper1, vars.flatTermVolCurve));
shared_ptr<StrippedOptionletAdapter> strippedOptionletAdapter2(new
StrippedOptionletAdapter(optionletStripper2));
Handle<OptionletVolatilityStructure> vol2(strippedOptionletAdapter2);
vol2->enableExtrapolation();
// consistency check: diff(stripped vol1-stripped vol2)
for (Size strikeIndex=0; strikeIndex<vars.strikes.size(); ++strikeIndex) {
for (Size tenorIndex=0; tenorIndex<vars.optionTenors.size(); ++tenorIndex) {
Volatility strippedVol1 = vol1->volatility(vars.optionTenors[tenorIndex],
vars.strikes[strikeIndex], true);
Volatility strippedVol2 = vol2->volatility(vars.optionTenors[tenorIndex],
vars.strikes[strikeIndex], true);
// vol from flat vol surface (for comparison only)
Volatility flatVol = vars.flatTermVolSurface->volatility(vars.optionTenors[tenorIndex],
vars.strikes[strikeIndex], true);
Real error = std::fabs(strippedVol1-strippedVol2);
if (error>vars.tolerance)
BOOST_FAIL("\noption tenor: " << vars.optionTenors[tenorIndex] <<
"\nstrike: " << io::rate(vars.strikes[strikeIndex]) <<
"\nstripped vol1: " << io::rate(strippedVol1) <<
"\nstripped vol2: " << io::rate(strippedVol2) <<
"\nflat vol: " << io::rate(flatVol) <<
"\nerror: " << io::rate(error) <<
"\ntolerance: " << io::rate(vars.tolerance));
}
}
}
