blackvariancecurve.cpp File Reference

#include "toplevelfixture.hpp"
#include <boost/make_shared.hpp>
#include <boost/test/unit_test.hpp>
#include <ql/quotes/simplequote.hpp>
#include <ql/termstructures/volatility/equityfx/blackvariancecurve.hpp>
#include <ql/time/calendars/target.hpp>
#include <ql/time/daycounters/actualactual.hpp>
#include <qle/termstructures/blackvariancecurve3.hpp>

Go to the source code of this file.

Functions
	BOOST_AUTO_TEST_CASE (testBlackVarianceCurve)

Function Documentation

BOOST_AUTO_TEST_CASE()

BOOST_AUTO_TEST_CASE

(

testBlackVarianceCurve

)

Definition at line 37 of file blackvariancecurve.cpp.

                                             {
 
    BOOST_TEST_MESSAGE("Testing QuantExt::BlackVarianceCurve3...");
 
    SavedSettings backup;
    Settings::instance().evaluationDate() = Date(1, Dec, 2015);
    Date today = Settings::instance().evaluationDate();
 
    Natural settlementDays = 0;
    Calendar cal = TARGET();
    BusinessDayConvention bdc = Following;
    DayCounter dc = ActualActual(ActualActual::ISDA);
 
    vector<Time> times;
    vector<Date> dates;
    vector<Volatility> vols;
    vector<QuantLib::ext::shared_ptr<SimpleQuote> > simpleQuotes;
    vector<Handle<Quote> > quotes;
 
    Size numYears = 10;
    for (Size i = 1; i < numYears; i++) {
 
        Volatility vol = 0.1 + (0.01 * i); // 11% at 1Y, 12% at 2Y
        vols.push_back(vol);
 
        simpleQuotes.push_back(QuantLib::ext::make_shared<SimpleQuote>(vol));
        quotes.push_back(Handle<Quote>(simpleQuotes.back()));
 
        dates.push_back(Date(1, Dec, today.year() + i));
        times.push_back(dc.yearFraction(today, dates.back()));
    }
 
    // Build a QuantLib::BlackVarianceCurve
    BlackVarianceCurve bvcBase(today, dates, vols, dc);
    bvcBase.enableExtrapolation();
 
    // Build a QuantExt::BlackVarianceCurve3
    BlackVarianceCurve3 bvcTest(settlementDays, cal, bdc, dc, times, quotes);
    bvcTest.enableExtrapolation();
 
    Real strike = 1.0; // this is all ATM so we don't care
 
    // Check that bvcTest returns the expected values
    for (Size i = 0; i < times.size(); ++i) {
        BOOST_CHECK_CLOSE(bvcTest.blackVol(times[i], strike), vols[i], 1e-12);
        BOOST_CHECK_CLOSE(bvcTest.blackVol(dates[i], strike), vols[i], 1e-12);
    }
 
    // Now check that they give the same vols (including extrapolation)
    for (Time t = 0.1; t < numYears + 10.0; t += 0.1) {
        BOOST_CHECK_CLOSE(bvcBase.blackVol(t, strike), bvcTest.blackVol(t, strike), 1e-12);
    }
 
    // Now double the quotes
    for (Size i = 0; i < simpleQuotes.size(); ++i) {
        simpleQuotes[i]->setValue(simpleQuotes[i]->value() * 2.0);
    }
    // and check again
    for (Time t = 0.1; t < numYears + 10.0; t += 0.1) {
        BOOST_CHECK_CLOSE(bvcBase.blackVol(t, strike), 0.5 * bvcTest.blackVol(t, strike), 1e-12);
    }
}