/usr/include/ql/math/matrixutilities/pseudosqrt.hpp is in libquantlib0-dev 1.7.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | /* -*- mode: c++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
/*
Copyright (C) 2003, 2004 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 pseudosqrt.hpp
\brief pseudo square root of a real symmetric matrix
*/
#ifndef quantlib_pseudo_sqrt_hpp
#define quantlib_pseudo_sqrt_hpp
#include <ql/math/matrix.hpp>
namespace QuantLib {
//! algorithm used for matricial pseudo square root
struct SalvagingAlgorithm {
enum Type { None, Spectral, Hypersphere, LowerDiagonal, Higham };
};
//! Returns the pseudo square root of a real symmetric matrix
/*! Given a matrix \f$ M \f$, the result \f$ S \f$ is defined
as the matrix such that \f$ S S^T = M. \f$
If the matrix is not positive semi definite, it can
return an approximation of the pseudo square root
using a (user selected) salvaging algorithm.
For more information see: R. Rebonato and P. Jäckel, The most
general methodology to create a valid correlation matrix for
risk management and option pricing purposes, The Journal of
Risk, 2(2), Winter 1999/2000.
http://www.rebonato.com/correlationmatrix.pdf
Revised and extended in "Monte Carlo Methods in Finance",
by Peter Jäckel, Chapter 6.
\pre the given matrix must be symmetric.
\relates Matrix
\warning Higham algorithm only works for correlation matrices.
\test
- the correctness of the results is tested by reproducing
known good data.
- the correctness of the results is tested by checking
returned values against numerical calculations.
*/
const Disposable<Matrix> pseudoSqrt(
const Matrix&,
SalvagingAlgorithm::Type = SalvagingAlgorithm::None);
//! Returns the rank-reduced pseudo square root of a real symmetric matrix
/*! The result matrix has rank<=maxRank. If maxRank>=size, then the
specified percentage of eigenvalues out of the eigenvalues' sum is
retained.
If the input matrix is not positive semi definite, it can return an
approximation of the pseudo square root using a (user selected)
salvaging algorithm.
\pre the given matrix must be symmetric.
\relates Matrix
*/
const Disposable<Matrix> rankReducedSqrt(const Matrix&,
Size maxRank,
Real componentRetainedPercentage,
SalvagingAlgorithm::Type);
}
#endif
|