/usr/include/ql/methods/finitedifferences/stepcondition.hpp is in libquantlib0-dev 1.4-2.
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 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 | /* -*- mode: c++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
/*
Copyright (C) 2000, 2001, 2002, 2003 RiskMap srl
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 stepcondition.hpp
\brief conditions to be applied at every time step
*/
#ifndef quantlib_step_condition_hpp
#define quantlib_step_condition_hpp
#include <ql/math/array.hpp>
#include <ql/instruments/payoffs.hpp>
namespace QuantLib {
//! condition to be applied at every time step
/*! \ingroup findiff */
template <class array_type>
class StepCondition {
public:
virtual ~StepCondition() {}
virtual void applyTo(array_type& a, Time t) const = 0;
};
/* Abstract base class which allows step conditions to use both
payoff and array functions */
template <class array_type>
class CurveDependentStepCondition :
public StepCondition<array_type> {
public:
void applyTo(Array &a, Time) const {
#pragma omp parallel for
for (Size i = 0; i < a.size(); i++) {
a[i] =
applyToValue(a[i], getValue(a,i));
}
}
protected:
CurveDependentStepCondition(Option::Type type, Real strike)
: curveItem_(new PayoffWrapper(type, strike)) {};
CurveDependentStepCondition(const Payoff *p)
: curveItem_(new PayoffWrapper(p)) {};
CurveDependentStepCondition(const array_type & a)
: curveItem_(new ArrayWrapper(a)) {};
class CurveWrapper;
boost::shared_ptr<CurveWrapper> curveItem_;
Real getValue(const array_type &a, Size index) const {
return curveItem_->getValue(a, index);
}
virtual Real applyToValue(Real, Real) const {
QL_FAIL("not yet implemented");
}
class CurveWrapper {
public:
virtual ~CurveWrapper() {}
virtual Real getValue(const array_type &a,
int i) = 0;
};
class ArrayWrapper : public CurveWrapper {
private:
array_type value_;
public:
ArrayWrapper (const array_type &a)
: value_(a) {}
Real getValue(const array_type&, int i) {
return value_[i];
}
};
class PayoffWrapper : public CurveWrapper {
private:
boost::shared_ptr<Payoff> payoff_;
public:
PayoffWrapper (const Payoff * p)
: payoff_(p) {};
PayoffWrapper (Option::Type type, Real strike)
: payoff_(new PlainVanillaPayoff(type, strike)) {};
Real getValue(const array_type &a,
int i) {
return (*payoff_)(a[i]);
}
};
};
//! %null step condition
/*! \ingroup findiff */
template <class array_type>
class NullCondition : public StepCondition<array_type> {
public:
void applyTo(array_type&, Time) const {}
};
}
#endif
|