QuantLib: a free/open-source library for quantitative finance
fully annotated source code - version 1.34
Loading...
Searching...
No Matches
bsmrndcalculator.cpp
Go to the documentation of this file.
1/* -*- mode: c++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
2
3/*
4 Copyright (C) 2015 Johannes Göttker-Schnetmann
5 Copyright (C) 2015 Klaus Spanderen
6
7 This file is part of QuantLib, a free-software/open-source library
8 for financial quantitative analysts and developers - http://quantlib.org/
9
10 QuantLib is free software: you can redistribute it and/or modify it
11 under the terms of the QuantLib license. You should have received a
12 copy of the license along with this program; if not, please email
13 <quantlib-dev@lists.sf.net>. The license is also available online at
14 <http://quantlib.org/license.shtml>.
15
16 This program is distributed in the hope that it will be useful, but WITHOUT
17 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
18 FOR A PARTICULAR PURPOSE. See the license for more details.
19*/
20
21/*! \file bsmrndcalculator.hpp
22 \brief risk neutral terminal density calculator for the
23 Black-Scholes-Merton model with constant volatility
24*/
25
29#include <cmath>
30#include <utility>
31
32namespace QuantLib {
33
34 BSMRNDCalculator::BSMRNDCalculator(ext::shared_ptr<GeneralizedBlackScholesProcess> process)
35 : process_(std::move(process)) {}
36
37 std::pair<Real, Volatility>
39 const Volatility stdDev =
40 process_->blackVolatility()->blackVol(t, std::exp(x))*std::sqrt(t);
41 const Real mean = std::log(process_->x0()) - 0.5*stdDev*stdDev
42 + std::log( process_->dividendYield()->discount(t)
43 / process_->riskFreeRate()->discount(t));
44
45 return std::make_pair(mean, stdDev);
46 }
47
49 std::pair<Real, Volatility> p = distributionParams(x, t);
50 return NormalDistribution(p.first, p.second)(x);
51 }
52
54 std::pair<Real, Volatility> p = distributionParams(x, t);
55 return CumulativeNormalDistribution(p.first, p.second)(x);
56 }
57
59 std::pair<Real, Volatility> p = distributionParams(x, t);
60 return InvCumulativeNormalDistribution(p.first, p.second)(x);
61 }
62}
Black-Scholes processes.
risk neutral terminal density calculator for the Black-Scholes-Merton model with constant volatility
BSMRNDCalculator(ext::shared_ptr< GeneralizedBlackScholesProcess > process)
Real invcdf(Real q, Time t) const override
const ext::shared_ptr< GeneralizedBlackScholesProcess > process_
Real pdf(Real x, Time t) const override
Real cdf(Real x, Time t) const override
std::pair< Real, Volatility > distributionParams(Real x, Time t) const
Cumulative normal distribution function.
Normal distribution function.
const DefaultType & t
Real Time
continuous quantity with 1-year units
Definition: types.hpp:62
QL_REAL Real
real number
Definition: types.hpp:50
Real Volatility
volatility
Definition: types.hpp:78
Definition: any.hpp:35
InverseCumulativeNormal InvCumulativeNormalDistribution
STL namespace.
normal, cumulative and inverse cumulative distributions