QuantLib
: a free/open-source library for quantitative finance
fully annotated source code - version 1.34
Loading...
Searching...
No Matches
ql
math
matrixutilities
tqreigendecomposition.hpp
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) 2005 Klaus Spanderen
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
/*! \file tqreigendecomposition.hpp
21
\brief tridiag. QR eigen decomposition with explicite shift aka Wilkinson
22
*/
23
24
#ifndef quantlib_tqr_eigen_decomposition_hpp
25
#define quantlib_tqr_eigen_decomposition_hpp
26
27
#include <
ql/math/array.hpp
>
28
#include <
ql/math/matrix.hpp
>
29
30
namespace
QuantLib
{
31
32
//! tridiag. QR eigen decomposition with explicite shift aka Wilkinson
33
/*! References:
34
35
Wilkinson, J.H. and Reinsch, C. 1971, Linear Algebra, vol. II of
36
Handbook for Automatic Computation (New York: Springer-Verlag)
37
38
"Numerical Recipes in C", 2nd edition,
39
Press, Teukolsky, Vetterling, Flannery,
40
41
\test the correctness of the result is tested by checking it
42
against known good values.
43
*/
44
class
TqrEigenDecomposition
{
45
public
:
46
enum
EigenVectorCalculation
{
WithEigenVector
,
47
WithoutEigenVector
,
48
OnlyFirstRowEigenVector
};
49
50
enum
ShiftStrategy
{
NoShift
,
51
Overrelaxation
,
52
CloseEigenValue
};
53
54
TqrEigenDecomposition
(
const
Array
& diag,
55
const
Array
& sub,
56
EigenVectorCalculation
calc =
WithEigenVector
,
57
ShiftStrategy
strategy =
CloseEigenValue
);
58
59
const
Array
&
eigenvalues
()
const
{
return
d_
; }
60
const
Matrix
&
eigenvectors
()
const
{
return
ev_
; }
61
62
Size
iterations
()
const
{
return
iter_
; }
63
64
private
:
65
bool
offDiagIsZero
(
Size
k,
Array
& e);
66
67
Size
iter_
= 0;
68
Array
d_
;
69
Matrix
ev_
;
70
};
71
72
}
73
74
75
#endif
array.hpp
1-D array used in linear algebra.
QuantLib::Array
1-D array used in linear algebra.
Definition:
array.hpp:52
QuantLib::Matrix
Matrix used in linear algebra.
Definition:
matrix.hpp:41
QuantLib::TqrEigenDecomposition
tridiag. QR eigen decomposition with explicite shift aka Wilkinson
Definition:
tqreigendecomposition.hpp:44
QuantLib::TqrEigenDecomposition::ev_
Matrix ev_
Definition:
tqreigendecomposition.hpp:69
QuantLib::TqrEigenDecomposition::iter_
Size iter_
Definition:
tqreigendecomposition.hpp:67
QuantLib::TqrEigenDecomposition::d_
Array d_
Definition:
tqreigendecomposition.hpp:68
QuantLib::TqrEigenDecomposition::ShiftStrategy
ShiftStrategy
Definition:
tqreigendecomposition.hpp:50
QuantLib::TqrEigenDecomposition::Overrelaxation
@ Overrelaxation
Definition:
tqreigendecomposition.hpp:51
QuantLib::TqrEigenDecomposition::NoShift
@ NoShift
Definition:
tqreigendecomposition.hpp:50
QuantLib::TqrEigenDecomposition::CloseEigenValue
@ CloseEigenValue
Definition:
tqreigendecomposition.hpp:52
QuantLib::TqrEigenDecomposition::EigenVectorCalculation
EigenVectorCalculation
Definition:
tqreigendecomposition.hpp:46
QuantLib::TqrEigenDecomposition::WithoutEigenVector
@ WithoutEigenVector
Definition:
tqreigendecomposition.hpp:47
QuantLib::TqrEigenDecomposition::OnlyFirstRowEigenVector
@ OnlyFirstRowEigenVector
Definition:
tqreigendecomposition.hpp:48
QuantLib::TqrEigenDecomposition::WithEigenVector
@ WithEigenVector
Definition:
tqreigendecomposition.hpp:46
QuantLib::TqrEigenDecomposition::offDiagIsZero
bool offDiagIsZero(Size k, Array &e)
Definition:
tqreigendecomposition.cpp:144
QuantLib::TqrEigenDecomposition::eigenvectors
const Matrix & eigenvectors() const
Definition:
tqreigendecomposition.hpp:60
QuantLib::TqrEigenDecomposition::eigenvalues
const Array & eigenvalues() const
Definition:
tqreigendecomposition.hpp:59
QuantLib::TqrEigenDecomposition::iterations
Size iterations() const
Definition:
tqreigendecomposition.hpp:62
QuantLib::Size
std::size_t Size
size of a container
Definition:
types.hpp:58
matrix.hpp
matrix used in linear algebra.
QuantLib
Definition:
any.hpp:35
Generated by
Doxygen
1.9.5