QuantLib: a free/open-source library for quantitative finance
fully annotated source code - version 1.34
Loading...
Searching...
No Matches
brownianbridge.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) 2006 StatPro Italia srl
5
6 This file is part of QuantLib, a free-software/open-source library
7 for financial quantitative analysts and developers - http://quantlib.org/
8
9 QuantLib is free software: you can redistribute it and/or modify it
10 under the terms of the QuantLib license. You should have received a
11 copy of the license along with this program; if not, please email
12 <quantlib-dev@lists.sf.net>. The license is also available online at
13 <http://quantlib.org/license.shtml>.
14
15 This program is distributed in the hope that it will be useful, but WITHOUT
16 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
17 FOR A PARTICULAR PURPOSE. See the license for more details.
18*/
19
20// ===========================================================================
21// NOTE: The following copyright notice applies to the original code,
22//
23// Copyright (C) 2002 Peter Jäckel "Monte Carlo Methods in Finance".
24// All rights reserved.
25//
26// Permission to use, copy, modify, and distribute this software is freely
27// granted, provided that this notice is preserved.
28// ===========================================================================
29
31
32namespace QuantLib {
33
35 : size_(steps), t_(size_), sqrtdt_(size_),
36 bridgeIndex_(size_), leftIndex_(size_), rightIndex_(size_),
37 leftWeight_(size_), rightWeight_(size_), stdDev_(size_) {
38 for (Size i=0; i<size_; ++i)
39 t_[i] = static_cast<Time>(i+1);
40 initialize();
41 }
42
43 BrownianBridge::BrownianBridge(const std::vector<Time>& times)
44 : size_(times.size()), t_(times), sqrtdt_(size_),
45 bridgeIndex_(size_), leftIndex_(size_), rightIndex_(size_),
46 leftWeight_(size_), rightWeight_(size_), stdDev_(size_) {
47 initialize();
48 }
49
51 : size_(timeGrid.size()-1), t_(size_), sqrtdt_(size_),
52 bridgeIndex_(size_), leftIndex_(size_), rightIndex_(size_),
53 leftWeight_(size_), rightWeight_(size_), stdDev_(size_) {
54 for (Size i=0; i<size_; ++i)
55 t_[i] = timeGrid[i+1];
56 initialize();
57 }
58
59
61
62 sqrtdt_[0] = std::sqrt(t_[0]);
63 for (Size i=1; i<size_; ++i)
64 sqrtdt_[i] = std::sqrt(t_[i]-t_[i-1]);
65
66 // map is used to indicate which points are already constructed.
67 // If map[i] is zero, path point i is yet unconstructed.
68 // map[i]-1 is the index of the variate that constructs
69 // the path point # i.
70 std::vector<Size> map(size_, 0);
71
72 // The first point in the construction is the global step.
73 map[size_-1] = 1;
74 // The global step is constructed from the first variate.
75 bridgeIndex_[0] = size_-1;
76 // The variance of the global step
77 stdDev_[0] = std::sqrt(t_[size_-1]);
78 // The global step to the last point in time is special.
79 leftWeight_[0] = rightWeight_[0] = 0.0;
80 for (Size j=0, i=1; i<size_; ++i) {
81 // Find the next unpopulated entry in the map.
82 while (map[j] != 0U)
83 ++j;
84 Size k = j;
85 // Find the next populated entry in the map from there.
86 while (map[k] == 0U)
87 ++k;
88 // l-1 is now the index of the point to be constructed next.
89 Size l = j + ((k-1-j)>>1);
90 map[l] = i;
91 // The i-th Gaussian variate will be used to set point l-1.
92 bridgeIndex_[i] = l;
93 leftIndex_[i] = j;
94 rightIndex_[i] = k;
95 if (j != 0) {
96 leftWeight_[i]= (t_[k]-t_[l])/(t_[k]-t_[j-1]);
97 rightWeight_[i] = (t_[l]-t_[j-1])/(t_[k]-t_[j-1]);
98 stdDev_[i] =
99 std::sqrt(((t_[l]-t_[j-1])*(t_[k]-t_[l]))
100 /(t_[k]-t_[j-1]));
101 } else {
102 leftWeight_[i] = (t_[k]-t_[l])/t_[k];
103 rightWeight_[i] = t_[l]/t_[k];
104 stdDev_[i] = std::sqrt(t_[l]*(t_[k]-t_[l])/t_[k]);
105 }
106 j=k+1;
107 if (j>=size_)
108 j=0; // wrap around
109 }
110 }
111
112}
113
Browian bridge.
std::vector< Size > leftIndex_
std::vector< Real > sqrtdt_
std::vector< Size > rightIndex_
std::vector< Time > t_
std::vector< Real > stdDev_
std::vector< Real > rightWeight_
std::vector< Size > bridgeIndex_
std::vector< Real > leftWeight_
time grid class
Definition: timegrid.hpp:43
Real Time
continuous quantity with 1-year units
Definition: types.hpp:62
std::size_t Size
size of a container
Definition: types.hpp:58
Definition: any.hpp:35
Size size_
Definition: pseudosqrt.cpp:76