QuantLib: a free/open-source library for quantitative finance
fully annotated source code - version 1.34
Loading...
Searching...
No Matches
ninepointlinearop.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) 2008 Andreas Gaida
5 Copyright (C) 2008 Ralph Schreyer
6 Copyright (C) 2008 Klaus Spanderen
7
8 This file is part of QuantLib, a free-software/open-source library
9 for financial quantitative analysts and developers - http://quantlib.org/
10
11 QuantLib is free software: you can redistribute it and/or modify it
12 under the terms of the QuantLib license. You should have received a
13 copy of the license along with this program; if not, please email
14 <quantlib-dev@lists.sf.net>. The license is also available online at
15 <http://quantlib.org/license.shtml>.
16
17 This program is distributed in the hope that it will be useful, but WITHOUT
18 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
19 FOR A PARTICULAR PURPOSE. See the license for more details.
20*/
21
22/*! \file ninepointlinearop.hpp
23 \brief nine point linear operator
24*/
25
26#ifndef quantlib_nine_point_linear_op_hpp
27#define quantlib_nine_point_linear_op_hpp
28
31#include <memory>
32
33namespace QuantLib {
34
35 class FdmMesher;
36
38 public:
40 const ext::shared_ptr<FdmMesher>& mesher);
45 ~NinePointLinearOp() override = default;
46
47 Array apply(const Array& r) const override;
48 NinePointLinearOp mult(const Array& u) const;
49
50 void swap(NinePointLinearOp& m) noexcept;
51
52 SparseMatrix toMatrix() const override;
53
54 protected:
55 NinePointLinearOp() = default;
56
58 std::unique_ptr<Size[]> i00_, i10_, i20_;
59 std::unique_ptr<Size[]> i01_, i21_;
60 std::unique_ptr<Size[]> i02_, i12_, i22_;
61 std::unique_ptr<Real[]> a00_, a10_, a20_;
62 std::unique_ptr<Real[]> a01_, a11_, a21_;
63 std::unique_ptr<Real[]> a02_, a12_, a22_;
64
65 ext::shared_ptr<FdmMesher> mesher_;
66 };
67
68
70 swap(m);
71 }
72
74 NinePointLinearOp temp(m);
75 swap(temp);
76 return *this;
77 }
78
80 swap(m);
81 return *this;
82 }
83
84}
85
86#endif
1-D array used in linear algebra.
Definition: array.hpp:52
std::unique_ptr< Size[]> i12_
SparseMatrix toMatrix() const override
std::unique_ptr< Size[]> i20_
std::unique_ptr< Real[]> a20_
~NinePointLinearOp() override=default
std::unique_ptr< Size[]> i02_
std::unique_ptr< Size[]> i21_
std::unique_ptr< Real[]> a11_
NinePointLinearOp & operator=(const NinePointLinearOp &m)
std::unique_ptr< Size[]> i01_
std::unique_ptr< Real[]> a00_
std::unique_ptr< Real[]> a02_
std::unique_ptr< Size[]> i10_
Array apply(const Array &r) const override
NinePointLinearOp mult(const Array &u) const
ext::shared_ptr< FdmMesher > mesher_
std::unique_ptr< Real[]> a22_
std::unique_ptr< Real[]> a10_
void swap(NinePointLinearOp &m) noexcept
std::unique_ptr< Real[]> a12_
std::unique_ptr< Real[]> a21_
std::unique_ptr< Size[]> i00_
std::unique_ptr< Size[]> i22_
std::unique_ptr< Real[]> a01_
linear operator to model a multi dimensinal pde system
std::size_t Size
size of a container
Definition: types.hpp:58
Definition: any.hpp:35
void swap(Array &v, Array &w) noexcept
Definition: array.hpp:903
boost::numeric::ublas::compressed_matrix< Real > SparseMatrix
ext::shared_ptr< YieldTermStructure > r
typedef for boost sparse matrix class