HwCG Class Reference

#include <ored/scripting/models/hwcg.hpp>

Collaboration diagram for HwCG:

Public Member Functions
	HwCG (const std::string &qualifier, QuantExt::ComputationGraph &g, const std::function< QuantLib::ext::shared_ptr< IrLgm1fParametrization >()> &p, std::vector< std::pair< std::size_t, std::function< double(void)> > > &modelParameters, const bool sloppySimDates=false, const std::set< Date > &effSimDates={})
QuantLib::ext::shared_ptr< IrLgm1fParametrization >	parametrization () const
std::size_t	numeraire (const Date &d, const std::size_t x, const Handle< YieldTermStructure > &discountCurve=Handle< YieldTermStructure >(), const std::string &discountCurveId="default") const
std::size_t	discountBond (const Date &d, Date e, const std::size_t x, const Handle< YieldTermStructure > &discountCurve=Handle< YieldTermStructure >(), const std::string &discountCurveId="default") const
std::size_t	reducedDiscountBond (const Date &d, const Date &e, const std::size_t x, const Handle< YieldTermStructure > &discountCurve=Handle< YieldTermStructure >(), const std::string &discountCurveId="default") const
std::size_t	fixing (const QuantLib::ext::shared_ptr< InterestRateIndex > &index, const Date &fixingDate, const Date &t, const std::size_t x) const

Private Attributes
std::string	qualifier_
QuantExt::ComputationGraph &	g_
std::function< QuantLib::ext::shared_ptr< IrLgm1fParametrization >()>	p_
std::vector< std::pair< std::size_t, std::function< double(void)> > > &	modelParameters_
bool	sloppySimDates_
std::set< Date >	effSimDates_

Detailed Description

Definition at line 36 of file hwcg.hpp.

Constructor & Destructor Documentation

HwCG()

HwCG	(	const std::string &	qualifier,
		QuantExt::ComputationGraph &	g,
		const std::function< QuantLib::ext::shared_ptr< IrLgm1fParametrization >()> &	p,
		std::vector< std::pair< std::size_t, std::function< double(void)> > > &	modelParameters,
		const bool	sloppySimDates = false,
		const std::set< Date > &	effSimDates = {}
	)

Definition at line 38 of file hwcg.hpp.

                                                                              {})
        : qualifier_(qualifier), g_(g), p_(p), modelParameters_(modelParameters), sloppySimDates_(sloppySimDates),
          effSimDates_(effSimDates) {}

Member Function Documentation

parametrization()

QuantLib::ext::shared_ptr< IrLgm1fParametrization > parametrization

(

)

const

Definition at line 45 of file hwcg.hpp.

{ return p_(); }

numeraire()

std::size_t numeraire	(	const Date &	d,
		const std::size_t	x,
		const Handle< YieldTermStructure > &	discountCurve = Handle<YieldTermStructure>(),
		const std::string &	discountCurveId = "default"
	)		const

Definition at line 34 of file hwcg.cpp.

                                                                    {
    std::string id = "__hw_" + qualifier_ + "_N_" + ore::data::to_string(d) + "_" + discountCurveId;
    std::size_t n;
    if (n = cg_var(g_, id, ComputationGraph::VarDoesntExist::Nan); n == ComputationGraph::nan) {
        auto p(p_);
        Date ds = getSloppyDate(d, sloppySimDates_, effSimDates_);
        Real t = p()->termStructure()->timeFromReference(d);
        // Real ts = p()->termStructure()->timeFromReference(ds);
        std::size_t I = cg_var(g_, "__hw_" + qualifier_ + "_I_" + ore::data::to_string(ds));
 
        std::string id_P0t = "__dsc_" + ore::data::to_string(ds) + "_" + discountCurveId;
        std::size_t P0t = addModelParameter(g_, modelParameters_, id_P0t, [p, discountCurve, t] {
            return (discountCurve.empty() ? p()->termStructure() : discountCurve)->discount(t);
        });
        n = cg_div(g_, cg_exp(g_, I), P0t);
        g_.setVariable(id, n);
    }
    return n;
}

Here is the call graph for this function:

Here is the caller graph for this function:

discountBond()

std::size_t discountBond	(	const Date &	d,
		Date	e,
		const std::size_t	x,
		const Handle< YieldTermStructure > &	discountCurve = Handle<YieldTermStructure>(),
		const std::string &	discountCurveId = "default"
	)		const

Definition at line 55 of file hwcg.cpp.

                                                                       {
    e = std::max(d, e);
    std::string id =
        "__hw_" + qualifier_ + "_P_" + ore::data::to_string(d) + "_" + ore::data::to_string(e) + "_" + discountCurveId;
    std::size_t n;
    if (n = cg_var(g_, id, ComputationGraph::VarDoesntExist::Nan), n == ComputationGraph::nan) {
        auto p = p_;
        Date ds = getSloppyDate(d, sloppySimDates_, effSimDates_);
        Date es = getSloppyDate(e, sloppySimDates_, effSimDates_);
        Real t = p()->termStructure()->timeFromReference(d);
        Real T = p()->termStructure()->timeFromReference(e);
        Real ts = p()->termStructure()->timeFromReference(ds);
        Real Ts = p()->termStructure()->timeFromReference(es);
        std::string id_P0t = "__dsc_" + ore::data::to_string(ds) + "_" + discountCurveId;
        std::string id_P0T = "__dsc_" + ore::data::to_string(es) + "_" + discountCurveId;
 
        std::size_t y = cg_var(g_, "__hw_" + qualifier_ + "_y_" + ore::data::to_string(ds));
        std::size_t g = cg_const(g_, Ts - ts);
        std::size_t exponent = cg_negative(
            g_, cg_add(g_, cg_mult(g_, g, x), cg_mult(g_, cg_const(g_, 0.5), cg_mult(g_, y, cg_mult(g_, g, g)))));
 
        std::size_t P0t = addModelParameter(g_, modelParameters_, id_P0t, [p, discountCurve, t] {
            return (discountCurve.empty() ? p()->termStructure() : discountCurve)->discount(t);
        });
        std::size_t P0T = addModelParameter(g_, modelParameters_, id_P0T, [p, discountCurve, T] {
            return (discountCurve.empty() ? p()->termStructure() : discountCurve)->discount(T);
        });
 
        n = cg_mult(g_, cg_div(g_, P0T, P0t), cg_exp(g_, exponent));
        g_.setVariable(id, n);
    }
    return n;
}

Here is the call graph for this function:

Here is the caller graph for this function:

reducedDiscountBond()

std::size_t reducedDiscountBond	(	const Date &	d,
		const Date &	e,
		const std::size_t	x,
		const Handle< YieldTermStructure > &	discountCurve = Handle<YieldTermStructure>(),
		const std::string &	discountCurveId = "default"
	)		const

Definition at line 91 of file hwcg.cpp.

                                                                              {
    std::string id =
        "__hw_" + qualifier_ + "_Pr_" + ore::data::to_string(d) + "_" + ore::data::to_string(e) + "_" + discountCurveId;
    std::size_t n;
    if (n = cg_var(g_, id, ComputationGraph::VarDoesntExist::Nan), n == ComputationGraph::nan) {
        n = cg_div(g_, discountBond(d, e, x, discountCurve, discountCurveId),
                   numeraire(d, x, discountCurve, discountCurveId));
        g_.setVariable(id, n);
    }
    return n;
}

Here is the call graph for this function:

fixing()

std::size_t fixing	(	const QuantLib::ext::shared_ptr< InterestRateIndex > &	index,
		const Date &	fixingDate,
		const Date &	t,
		const std::size_t	x
	)		const

Definition at line 106 of file hwcg.cpp.

                                                  {
 
    std::string id =
        "__irFix_" + index->name() + "_" + ore::data::to_string(fixingDate) + "_" + ore::data::to_string(t);
    std::size_t n;
    if (n = cg_var(g_, id, ComputationGraph::VarDoesntExist::Nan); n == ComputationGraph::nan) {
 
        // handle case where fixing is deterministic
 
        Date today = Settings::instance().evaluationDate();
        if (fixingDate <= today) {
 
            // handle historical fixing
 
            n = addModelParameter(g_, modelParameters_, id,
                                  [index, fixingDate]() { return index->fixing(fixingDate); });
 
        } else if (auto ibor = QuantLib::ext::dynamic_pointer_cast<IborIndex>(index)) {
 
            // Ibor Index
 
            Date d1 = std::max(t, ibor->valueDate(fixingDate));
            Date d2 = ibor->maturityDate(d1);
 
            Time dt = ibor->dayCounter().yearFraction(d1, d2);
 
            std::size_t disc1 = discountBond(t, d1, x, ibor->forwardingTermStructure(), "fwd_" + index->name());
            std::size_t disc2 = discountBond(t, d2, x, ibor->forwardingTermStructure(), "fwd_" + index->name());
 
            n = cg_div(g_, cg_subtract(g_, cg_div(g_, disc1, disc2), cg_const(g_, 1.0)), cg_const(g_, dt));
 
        } else {
            QL_FAIL("HwCG::fixing(): only ibor indices handled so far, index = " << index->name());
        }
        g_.setVariable(id, n);
    }
    return n;
}

Here is the call graph for this function:

Member Data Documentation

qualifier_

std::string qualifier_

private

Definition at line 64 of file hwcg.hpp.

g_

QuantExt::ComputationGraph& g_

private

Definition at line 65 of file hwcg.hpp.

p_

std::function<QuantLib::ext::shared_ptr<IrLgm1fParametrization>()> p_

private

Definition at line 66 of file hwcg.hpp.

modelParameters_

std::vector<std::pair<std::size_t, std::function<double(void)> > >& modelParameters_

private

Definition at line 67 of file hwcg.hpp.

sloppySimDates_

bool sloppySimDates_

private

Definition at line 68 of file hwcg.hpp.

effSimDates_

std::set<Date> effSimDates_

private

Definition at line 69 of file hwcg.hpp.