23#include <ql/exercise.hpp>
24#include <ql/methods/finitedifferences/utilities/escroweddividendadjustment.hpp>
25#include <ql/methods/finitedifferences/meshers/fdmblackscholesmesher.hpp>
26#include <ql/methods/finitedifferences/meshers/fdmmeshercomposite.hpp>
27#include <ql/methods/finitedifferences/solvers/fdmblackscholessolver.hpp>
28#include <ql/methods/finitedifferences/stepconditions/fdmstepconditioncomposite.hpp>
29#include <ql/methods/finitedifferences/utilities/fdmshoutloginnervaluecalculator.hpp>
30#include <ql/pricingengines/vanilla/fdblackscholesshoutengine.hpp>
31#include <ql/processes/blackscholesprocess.hpp>
36 ext::shared_ptr<GeneralizedBlackScholesProcess> process,
41 : process_(
std::move(process)), explicitDividends_(false),
42 tGrid_(tGrid), xGrid_(xGrid), dampingSteps_(dampingSteps),
43 schemeDesc_(schemeDesc) {
48 ext::shared_ptr<GeneralizedBlackScholesProcess> process,
54 : process_(
std::move(process)), dividends_(
std::move(dividends)), explicitDividends_(true),
55 tGrid_(tGrid), xGrid_(xGrid), dampingSteps_(dampingSteps),
56 schemeDesc_(schemeDesc) {
63 QL_DEPRECATED_DISABLE_WARNING
65 QL_DEPRECATED_ENABLE_WARNING
69 const Date settlementDate =
process_->riskFreeRate()->referenceDate();
71 const auto escrowedDividendAdj =
72 ext::make_shared<EscrowedDividendAdjustment>(
76 [&](
Date d){ return process_->time(d); },
79 const Real divAdj = escrowedDividendAdj
80 ->dividendAdjustment(
process_->time(settlementDate));
82 QL_REQUIRE(
process_->x0() + divAdj > 0.0,
83 "spot minus dividends becomes negative");
86 ext::dynamic_pointer_cast<PlainVanillaPayoff>(
arguments_.payoff);
88 QL_REQUIRE(payoff,
"non plain vanilla payoff given");
92 const auto mesher = ext::make_shared<FdmMesherComposite>(
93 ext::make_shared<FdmBlackScholesMesher>(
96 std::pair<Real, Real>(payoff->strike(), 0.1),
97 emptyDividendSchedule,
98 ext::shared_ptr<FdmQuantoHelper>(),
101 const auto innerValuecalculator =
102 ext::make_shared<FdmShoutLogInnerValueCalculator>(
104 escrowedDividendAdj, maturity, payoff, mesher, 0);
107 for (
const auto& cf: passedDividends)
108 zeroDividendSchedule.push_back(
109 ext::make_shared<FixedDividend>(0.0, cf->date()));
111 const auto conditions =
113 zeroDividendSchedule,
115 innerValuecalculator,
116 process_->riskFreeRate()->referenceDate(),
117 process_->riskFreeRate()->dayCounter());
121 conditions, innerValuecalculator,
125 ext::make_shared<FdmBlackScholesSolver>(
132 results_.delta = solver->deltaAt(spot);
133 results_.gamma = solver->gammaAt(spot);
134 results_.theta = solver->thetaAt(spot);
DividendSchedule cashFlow
void calculate() const override
const ext::shared_ptr< GeneralizedBlackScholesProcess > process_
DividendSchedule dividends_
FdBlackScholesShoutEngine(ext::shared_ptr< GeneralizedBlackScholesProcess >, Size tGrid=100, Size xGrid=100, Size dampingSteps=0, const FdmSchemeDesc &schemeDesc=FdmSchemeDesc::Douglas())
const FdmSchemeDesc schemeDesc_
static ext::shared_ptr< FdmStepConditionComposite > vanillaComposite(const DividendSchedule &schedule, const ext::shared_ptr< Exercise > &exercise, const ext::shared_ptr< FdmMesher > &mesher, const ext::shared_ptr< FdmInnerValueCalculator > &calculator, const Date &refDate, const DayCounter &dayCounter)
DividendVanillaOption::results results_
DividendVanillaOption::arguments arguments_
Shared handle to an observable.
template class providing a null value for a given type.
std::pair< iterator, bool > registerWith(const ext::shared_ptr< Observable > &)
Real Time
continuous quantity with 1-year units
std::size_t Size
size of a container
std::vector< ext::shared_ptr< Dividend > > DividendSchedule
OperatorTraits< FdmLinearOp >::bc_set FdmBoundaryConditionSet