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
|
/* -*- mode: c++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
/*
Copyright (C) 2003 Ferdinando Ametrano
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.
*/
/*! \file drifttermstructure.hpp
\brief Drift term structure
*/
#ifndef quantlib_drift_term_structure_hpp
#define quantlib_drift_term_structure_hpp
#include <ql/termstructures/yield/zeroyieldstructure.hpp>
#include <ql/termstructures/volatility/equityfx/blackvoltermstructure.hpp>
namespace QuantLib {
//! Drift term structure
/*! Drift term structure for modelling the common drift term:
riskFreeRate - dividendYield - 0.5*vol*vol
\note This term structure will remain linked to the original
structures, i.e., any changes in the latters will be
reflected in this structure as well.
*/
class DriftTermStructure : public ZeroYieldStructure {
public:
DriftTermStructure(const Handle<YieldTermStructure>& riskFreeTS,
const Handle<YieldTermStructure>& dividendTS,
const Handle<BlackVolTermStructure>& blackVolTS);
//! \name YieldTermStructure interface
//@{
DayCounter dayCounter() const { return riskFreeTS_->dayCounter(); }
Calendar calendar() const;
const Date& referenceDate() const;
Date maxDate() const { return maxDate_; }
//@}
protected:
//! returns the discount factor as seen from the evaluation date
Rate zeroYieldImpl(Time) const;
private:
Handle<YieldTermStructure> riskFreeTS_, dividendTS_;
Handle<BlackVolTermStructure> blackVolTS_;
Real underlyingLevel_;
Date maxDate_;
};
// inline definitions
inline DriftTermStructure::DriftTermStructure(
const Handle<YieldTermStructure>& riskFreeTS,
const Handle<YieldTermStructure>& dividendTS,
const Handle<BlackVolTermStructure>& blackVolTS)
: ZeroYieldStructure(riskFreeTS->dayCounter()),
riskFreeTS_(riskFreeTS),
dividendTS_(dividendTS),
blackVolTS_(blackVolTS) {
registerWith(riskFreeTS_);
registerWith(dividendTS_);
registerWith(blackVolTS_);
maxDate_ = std::min(dividendTS_->maxDate(),
riskFreeTS_->maxDate());
maxDate_ = std::min(maxDate_, blackVolTS_->maxDate());
}
inline Calendar DriftTermStructure::calendar() const {
return riskFreeTS_->calendar();
}
inline const Date& DriftTermStructure::referenceDate() const {
// warning: here it is assumed that all TS have the same referenceDate
// It should be QL_REQUIREd
return riskFreeTS_->referenceDate();
}
inline Rate DriftTermStructure::zeroYieldImpl(Time t) const {
// warning: here it is assumed that
// a) all TS have the same daycount.
// b) all TS have the same referenceDate
// It should be QL_REQUIREd
return riskFreeTS_->zeroRate(t, Continuous, NoFrequency, true)
- dividendTS_->zeroRate(t, Continuous, NoFrequency, true)
- 0.5 * blackVolTS_->blackVol(t, underlyingLevel_, true)
* blackVolTS_->blackVol(t, underlyingLevel_, true);
}
}
#endif
|
