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
|
/* -*- mode: c++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
/*
Copyright (C) 2001, 2002, 2003 Sadruddin Rejeb
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/models/calibrationhelper.hpp>
#include <ql/math/solvers1d/brent.hpp>
namespace QuantLib {
class CalibrationHelper::ImpliedVolatilityHelper {
public:
ImpliedVolatilityHelper(const CalibrationHelper& helper,
Real value)
: helper_(helper), value_(value) {}
Real operator()(Volatility x) const {
return value_ - helper_.blackPrice(x);
}
private:
const CalibrationHelper& helper_;
Real value_;
};
Volatility CalibrationHelper::impliedVolatility(Real targetValue,
Real accuracy,
Size maxEvaluations,
Volatility minVol,
Volatility maxVol) const {
ImpliedVolatilityHelper f(*this,targetValue);
Brent solver;
solver.setMaxEvaluations(maxEvaluations);
return solver.solve(f,accuracy,volatility_->value(),minVol,maxVol);
}
Real CalibrationHelper::calibrationError() {
if (calibrateVolatility_) {
const Real lowerPrice = blackPrice(0.001);
const Real upperPrice = blackPrice(10);
const Real modelPrice = modelValue();
Volatility implied;
if (modelPrice <= lowerPrice)
implied = 0.001;
else
if (modelPrice >= upperPrice)
implied = 10.0;
else
implied = this->impliedVolatility(
modelPrice, 1e-12, 5000, 0.001, 10);
return implied - volatility_->value();
}
else {
return std::fabs(marketValue() - modelValue())/marketValue();
}
}
}
|
