1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
/* -*- mode: c++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
/*
Copyright (C) 2018 Roy Zywina
Copyright (C) 2019 Eisuke Tani
This file is part of QuantLib, a free-software/open-source library
for financial quantitative analysts and developers - http://quantlib.org/
QuantLib is free software: you can redistribute it and/or modify it
under the terms of the QuantLib license. You should have received a
copy of the license along with this program; if not, please email
<quantlib-dev@lists.sf.net>. The license is also available online at
<http://quantlib.org/license.shtml>.
This program is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
FOR A PARTICULAR PURPOSE. See the license for more details.
*/
#include <ql/instruments/overnightindexfuture.hpp>
#include <ql/indexes/indexmanager.hpp>
#include <ql/event.hpp>
#include <utility>
namespace QuantLib {
OvernightIndexFuture::OvernightIndexFuture(ext::shared_ptr<OvernightIndex> overnightIndex,
const Date& valueDate,
const Date& maturityDate,
Handle<Quote> convexityAdjustment,
RateAveraging::Type averagingMethod)
: overnightIndex_(std::move(overnightIndex)), valueDate_(valueDate),
maturityDate_(maturityDate), convexityAdjustment_(std::move(convexityAdjustment)),
averagingMethod_(averagingMethod) {
QL_REQUIRE(overnightIndex_, "null overnight index");
registerWith(overnightIndex_);
}
Real OvernightIndexFuture::averagedRate() const {
Date today = Settings::instance().evaluationDate();
Calendar calendar = overnightIndex_->fixingCalendar();
DayCounter dayCounter = overnightIndex_->dayCounter();
Handle<YieldTermStructure> forwardCurve = overnightIndex_->forwardingTermStructure();
Real avg = 0;
Date d1 = valueDate_;
const TimeSeries<Real>& history = IndexManager::instance()
.getHistory(overnightIndex_->name());
Real fwd;
while (d1 < maturityDate_) {
Date d2 = calendar.advance(d1, 1, Days);
if (d1 < today) {
fwd = history[d1];
QL_REQUIRE(fwd != Null<Real>(), "missing rate on " <<
d1 << " for index " << overnightIndex_->name());
} else {
fwd = forwardCurve->forwardRate(d1, d2, dayCounter, Simple).rate();
}
avg += fwd * dayCounter.yearFraction(d1, d2);
d1 = d2;
}
return avg / dayCounter.yearFraction(valueDate_, maturityDate_);
}
Real OvernightIndexFuture::compoundedRate() const {
Date today = Settings::instance().evaluationDate();
Calendar calendar = overnightIndex_->fixingCalendar();
DayCounter dayCounter = overnightIndex_->dayCounter();
Handle<YieldTermStructure> forwardCurve = overnightIndex_->forwardingTermStructure();
Real prod = 1;
if (today > valueDate_) {
// can't value on a weekend inside reference period because we
// won't know the reset rate until start of next business day.
// user can supply an estimate if they really want to do this
today = calendar.adjust(today);
// for valuations inside the reference period, index quotes
// must have been populated in the history
const TimeSeries<Real>& history = IndexManager::instance()
.getHistory(overnightIndex_->name());
Date d1 = valueDate_;
while (d1 < today) {
Real r = history[d1];
QL_REQUIRE(r != Null<Real>(), "missing rate on " <<
d1 << " for index " << overnightIndex_->name());
Date d2 = calendar.advance(d1, 1, Days);
prod *= 1 + r * dayCounter.yearFraction(d1, d2);
d1 = d2;
}
}
DiscountFactor forwardDiscount = forwardCurve->discount(maturityDate_);
if (valueDate_ > today) {
forwardDiscount /= forwardCurve->discount(valueDate_);
}
prod /= forwardDiscount;
return (prod - 1) / dayCounter.yearFraction(valueDate_, maturityDate_);
}
Real OvernightIndexFuture::rate() const {
switch (averagingMethod_) {
case RateAveraging::Simple:
return averagedRate();
break;
case RateAveraging::Compound:
return compoundedRate();
break;
default:
QL_FAIL("unknown compounding convention (" << Integer(averagingMethod_) << ")");
}
}
bool OvernightIndexFuture::isExpired() const {
return detail::simple_event(maturityDate_).hasOccurred();
}
Real OvernightIndexFuture::convexityAdjustment() const {
return convexityAdjustment_.empty() ? 0.0 : convexityAdjustment_->value();
}
void OvernightIndexFuture::performCalculations() const {
Rate R = convexityAdjustment() + rate();
NPV_ = 100.0 * (1.0 - R);
}
}
|
