DiscountCurveJacobianRegressorSet.java

  1. package org.drip.regression.curvejacobian;

  2. /*
  3.  * -*- mode: java; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
  4.  */

  5. /*!
  6.  * Copyright (C) 2020 Lakshmi Krishnamurthy
  7.  * Copyright (C) 2019 Lakshmi Krishnamurthy
  8.  * Copyright (C) 2018 Lakshmi Krishnamurthy
  9.  * Copyright (C) 2017 Lakshmi Krishnamurthy
  10.  * Copyright (C) 2016 Lakshmi Krishnamurthy
  11.  * Copyright (C) 2015 Lakshmi Krishnamurthy
  12.  * Copyright (C) 2014 Lakshmi Krishnamurthy
  13.  * Copyright (C) 2013 Lakshmi Krishnamurthy
  14.  *
  15.  *  This file is part of DROP, an open-source library targeting analytics/risk, transaction cost analytics,
  16.  *      asset liability management analytics, capital, exposure, and margin analytics, valuation adjustment
  17.  *      analytics, and portfolio construction analytics within and across fixed income, credit, commodity,
  18.  *      equity, FX, and structured products. It also includes auxiliary libraries for algorithm support,
  19.  *      numerical analysis, numerical optimization, spline builder, model validation, statistical learning,
  20.  *      and computational support.
  21.  *  
  22.  *      https://lakshmidrip.github.io/DROP/
  23.  *  
  24.  *  DROP is composed of three modules:
  25.  *  
  26.  *  - DROP Product Core - https://lakshmidrip.github.io/DROP-Product-Core/
  27.  *  - DROP Portfolio Core - https://lakshmidrip.github.io/DROP-Portfolio-Core/
  28.  *  - DROP Computational Core - https://lakshmidrip.github.io/DROP-Computational-Core/
  29.  *
  30.  *  DROP Product Core implements libraries for the following:
  31.  *  - Fixed Income Analytics
  32.  *  - Loan Analytics
  33.  *  - Transaction Cost Analytics
  34.  *
  35.  *  DROP Portfolio Core implements libraries for the following:
  36.  *  - Asset Allocation Analytics
  37.  *  - Asset Liability Management Analytics
  38.  *  - Capital Estimation Analytics
  39.  *  - Exposure Analytics
  40.  *  - Margin Analytics
  41.  *  - XVA Analytics
  42.  *
  43.  *  DROP Computational Core implements libraries for the following:
  44.  *  - Algorithm Support
  45.  *  - Computation Support
  46.  *  - Function Analysis
  47.  *  - Model Validation
  48.  *  - Numerical Analysis
  49.  *  - Numerical Optimizer
  50.  *  - Spline Builder
  51.  *  - Statistical Learning
  52.  *
  53.  *  Documentation for DROP is Spread Over:
  54.  *
  55.  *  - Main                     => https://lakshmidrip.github.io/DROP/
  56.  *  - Wiki                     => https://github.com/lakshmiDRIP/DROP/wiki
  57.  *  - GitHub                   => https://github.com/lakshmiDRIP/DROP
  58.  *  - Repo Layout Taxonomy     => https://github.com/lakshmiDRIP/DROP/blob/master/Taxonomy.md
  59.  *  - Javadoc                  => https://lakshmidrip.github.io/DROP/Javadoc/index.html
  60.  *  - Technical Specifications => https://github.com/lakshmiDRIP/DROP/tree/master/Docs/Internal
  61.  *  - Release Versions         => https://lakshmidrip.github.io/DROP/version.html
  62.  *  - Community Credits        => https://lakshmidrip.github.io/DROP/credits.html
  63.  *  - Issues Catalog           => https://github.com/lakshmiDRIP/DROP/issues
  64.  *  - JUnit                    => https://lakshmidrip.github.io/DROP/junit/index.html
  65.  *  - Jacoco                   => https://lakshmidrip.github.io/DROP/jacoco/index.html
  66.  *
  67.  *  Licensed under the Apache License, Version 2.0 (the "License");
  68.  *      you may not use this file except in compliance with the License.
  69.  *  
  70.  *  You may obtain a copy of the License at
  71.  *      http://www.apache.org/licenses/LICENSE-2.0
  72.  *  
  73.  *  Unless required by applicable law or agreed to in writing, software
  74.  *      distributed under the License is distributed on an "AS IS" BASIS,
  75.  *      WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  76.  *  
  77.  *  See the License for the specific language governing permissions and
  78.  *      limitations under the License.
  79.  */

  80. /**
  81.  * <i>DiscountCurveJacobianRegressorSet</i> implements the regression analysis for the full discount curve
  82.  * (built from cash/future/swap) Sensitivity Jacobians. Specifically, it computes the PVDF micro-Jack.
  83.  *
  84.  * <br><br>
  85.  *  <ul>
  86.  *      <li><b>Module </b> = <a href = "https://github.com/lakshmiDRIP/DROP/tree/master/ComputationalCore.md">Computational Core Module</a></li>
  87.  *      <li><b>Library</b> = <a href = "https://github.com/lakshmiDRIP/DROP/tree/master/ComputationSupportLibrary.md">Computation Support</a></li>
  88.  *      <li><b>Project</b> = <a href = "https://github.com/lakshmiDRIP/DROP/tree/master/src/main/java/org/drip/regression/README.md">Regression Engine Core and the Unit Regressors</a></li>
  89.  *      <li><b>Package</b> = <a href = "https://github.com/lakshmiDRIP/DROP/tree/master/src/main/java/org/drip/regression/curvejacobian/README.md">Curve Jacobian Reconciliation Regression Engine</a></li>
  90.  *  </ul>
  91.  * <br><br>
  92.  *
  93.  * @author Lakshmi Krishnamurthy
  94.  */

  95. public class DiscountCurveJacobianRegressorSet implements org.drip.regression.core.RegressorSet {
  96.     private java.lang.String _strRegressionScenario =
  97.         "org.drip.analytics.definition.IRSDiscountCurve.CompPVDFJacobian";

  98.     private java.util.List<org.drip.regression.core.UnitRegressor> _setRegressors = new
  99.         java.util.ArrayList<org.drip.regression.core.UnitRegressor>();

  100.     @Override public java.util.List<org.drip.regression.core.UnitRegressor> getRegressorSet()
  101.     {
  102.         return _setRegressors;
  103.     }

  104.     @Override public boolean setupRegressors()
  105.     {
  106.         try {
  107.             _setRegressors.add (new org.drip.regression.core.UnitRegressionExecutor ("DiscountCurveJacobian",
  108.                 _strRegressionScenario)
  109.             {
  110.                 org.drip.analytics.date.JulianDate dtStart = null;
  111.                 org.drip.state.discount.MergedDiscountForwardCurve dcIRS = null;
  112.                 org.drip.numerical.differentiation.WengertJacobian wjPVDF = null;
  113.                 org.drip.numerical.differentiation.WengertJacobian aWJComp[] = null;
  114.                 org.drip.product.definition.CalibratableComponent aCompCalib[] = null;

  115.                 org.drip.param.market.LatentStateFixingsContainer lsfc = new
  116.                     org.drip.param.market.LatentStateFixingsContainer();

  117.                 @Override public boolean preRegression()
  118.                 {
  119.                     int NUM_DC_INSTR = 15;
  120.                     int aiDate[] = new int[NUM_DC_INSTR];
  121.                     double adblRate[] = new double[NUM_DC_INSTR];
  122.                     double adblCompCalibValue[] = new double[NUM_DC_INSTR];
  123.                     aWJComp = new org.drip.numerical.differentiation.WengertJacobian[NUM_DC_INSTR];
  124.                     java.lang.String astrCalibMeasure[] = new java.lang.String[NUM_DC_INSTR];
  125.                     aCompCalib = new org.drip.product.definition.CalibratableComponent[NUM_DC_INSTR];

  126.                     if (null == (dtStart = org.drip.analytics.date.DateUtil.CreateFromYMD (2011, 4, 6)))
  127.                         return false;

  128.                     aiDate[0] = dtStart.addDays ((int)(365.25 * 1 + 2)).julian(); // 4Y

  129.                     aiDate[1] = dtStart.addDays ((int)(365.25 * 2 + 2)).julian(); // 5Y

  130.                     aiDate[2] = dtStart.addDays ((int)(365.25 * 3 + 2)).julian(); // 6Y

  131.                     aiDate[3] = dtStart.addDays ((int)(365.25 * 7 + 2)).julian(); // 7Y

  132.                     aiDate[4] = dtStart.addDays ((int)(365.25 * 8 + 2)).julian(); // 8Y

  133.                     aiDate[5] = dtStart.addDays ((int)(365.25 * 9 + 2)).julian(); // 9Y

  134.                     aiDate[6] = dtStart.addDays ((int)(365.25 * 10 + 2)).julian(); // 10Y

  135.                     aiDate[7] = dtStart.addDays ((int)(365.25 * 11 + 2)).julian(); // 11Y

  136.                     aiDate[8] = dtStart.addDays ((int)(365.25 * 12 + 2)).julian(); // 12Y

  137.                     aiDate[9] = dtStart.addDays ((int)(365.25 * 15 + 2)).julian(); // 15Y

  138.                     aiDate[10] = dtStart.addDays ((int)(365.25 * 20 + 2)).julian(); // 20Y

  139.                     aiDate[11] = dtStart.addDays ((int)(365.25 * 25 + 2)).julian(); // 25Y

  140.                     aiDate[12] = dtStart.addDays ((int)(365.25 * 30 + 2)).julian(); // 30Y

  141.                     aiDate[13] = dtStart.addDays ((int)(365.25 * 40 + 2)).julian(); // 40Y

  142.                     aiDate[14] = dtStart.addDays ((int)(365.25 * 50 + 2)).julian(); // 50Y

  143.                     adblCompCalibValue[0] = .0166;
  144.                     adblCompCalibValue[1] = .0206;
  145.                     adblCompCalibValue[2] = .0241;
  146.                     adblCompCalibValue[3] = .0269;
  147.                     adblCompCalibValue[4] = .0292;
  148.                     adblCompCalibValue[5] = .0311;
  149.                     adblCompCalibValue[6] = .0326;
  150.                     adblCompCalibValue[7] = .0340;
  151.                     adblCompCalibValue[8] = .0351;
  152.                     adblCompCalibValue[9] = .0375;
  153.                     adblCompCalibValue[10] = .0393;
  154.                     adblCompCalibValue[11] = .0402;
  155.                     adblCompCalibValue[12] = .0407;
  156.                     adblCompCalibValue[13] = .0409;
  157.                     adblCompCalibValue[14] = .0409;
  158.                     org.drip.param.period.CompositePeriodSetting cpsFixed = null;
  159.                     org.drip.param.period.CompositePeriodSetting cpsFloating = null;
  160.                     org.drip.param.period.UnitCouponAccrualSetting ucasFixed = null;
  161.                     org.drip.param.period.ComposableFixedUnitSetting cfusFixed = null;
  162.                     org.drip.param.period.ComposableFloatingUnitSetting cfusFloating = null;

  163.                     try {
  164.                         ucasFixed = new org.drip.param.period.UnitCouponAccrualSetting (2, "Act/360", false,
  165.                             "Act/360", false, "USD", true,
  166.                                 org.drip.analytics.support.CompositePeriodBuilder.ACCRUAL_COMPOUNDING_RULE_GEOMETRIC);

  167.                         cfusFloating = new org.drip.param.period.ComposableFloatingUnitSetting ("3M",
  168.                             org.drip.analytics.support.CompositePeriodBuilder.EDGE_DATE_SEQUENCE_SINGLE,
  169.                                 null, org.drip.state.identifier.ForwardLabel.Standard ("USD-3M"),
  170.                                     org.drip.analytics.support.CompositePeriodBuilder.REFERENCE_PERIOD_IN_ADVANCE,
  171.                             0.);

  172.                         cfusFixed = new org.drip.param.period.ComposableFixedUnitSetting ("6M",
  173.                             org.drip.analytics.support.CompositePeriodBuilder.EDGE_DATE_SEQUENCE_REGULAR,
  174.                                 null, 0., 0., "USD");

  175.                         cpsFloating = new org.drip.param.period.CompositePeriodSetting (4, "3M", "USD", null,
  176.                             -1., null, null, null, null);

  177.                         cpsFixed = new org.drip.param.period.CompositePeriodSetting (2, "6M", "USD", null,
  178.                             1., null, null, null, null);
  179.                     } catch (java.lang.Exception e) {
  180.                         e.printStackTrace();

  181.                         return false;
  182.                     }

  183.                     for (int i = 0; i < NUM_DC_INSTR; ++i) {
  184.                         adblRate[i] = 0.01;
  185.                         astrCalibMeasure[i] = "Rate";

  186.                         try {
  187.                             java.util.List<java.lang.Integer> lsFixedStreamEdgeDate =
  188.                                 org.drip.analytics.support.CompositePeriodBuilder.BackwardEdgeDates (dtStart,
  189.                                     new org.drip.analytics.date.JulianDate (aiDate[i]), "6M", null,
  190.                                         org.drip.analytics.support.CompositePeriodBuilder.SHORT_STUB);

  191.                             java.util.List<java.lang.Integer> lsFloatingStreamEdgeDate =
  192.                                 org.drip.analytics.support.CompositePeriodBuilder.BackwardEdgeDates (dtStart,
  193.                                     new org.drip.analytics.date.JulianDate (aiDate[i]), "3M", null,
  194.                                         org.drip.analytics.support.CompositePeriodBuilder.SHORT_STUB);

  195.                             org.drip.product.rates.Stream floatingStream = new org.drip.product.rates.Stream
  196.                                 (org.drip.analytics.support.CompositePeriodBuilder.FloatingCompositeUnit
  197.                                     (lsFloatingStreamEdgeDate, cpsFloating, cfusFloating));

  198.                             org.drip.product.rates.Stream fixedStream = new org.drip.product.rates.Stream
  199.                                 (org.drip.analytics.support.CompositePeriodBuilder.FixedCompositeUnit
  200.                                     (lsFixedStreamEdgeDate, cpsFixed, ucasFixed, cfusFixed));


  201.                             aCompCalib[i] = new org.drip.product.rates.FixFloatComponent (fixedStream,
  202.                                 floatingStream, null);
  203.                         } catch (java.lang.Exception e) {
  204.                             e.printStackTrace();

  205.                             return false;
  206.                         }
  207.                     }

  208.                     lsfc.add (dtStart.addDays (2), org.drip.state.identifier.ForwardLabel.Standard
  209.                         ("USD-6M"), 0.0042);

  210.                     return null != (dcIRS =
  211.                         org.drip.state.creator.ScenarioDiscountCurveBuilder.NonlinearBuild (dtStart, "USD",
  212.                             aCompCalib, adblCompCalibValue, astrCalibMeasure, lsfc));
  213.                 }

  214.                 @Override public boolean execRegression()
  215.                 {
  216.                     for (int i = 0; i < aCompCalib.length; ++i) {
  217.                         try {
  218.                             if (null == (aWJComp[i] = aCompCalib[i].jackDDirtyPVDManifestMeasure (new
  219.                                 org.drip.param.valuation.ValuationParams (dtStart, dtStart, "USD"), null,
  220.                                     org.drip.param.creator.MarketParamsBuilder.Create (dcIRS, null, null,
  221.                                         null, null, null, lsfc), null)))
  222.                                 return false;
  223.                         } catch (java.lang.Exception e) {
  224.                             e.printStackTrace();

  225.                             return false;
  226.                         }
  227.                     }

  228.                     return null != (wjPVDF = dcIRS.compJackDPVDManifestMeasure (dtStart));
  229.                 }

  230.                 @Override public boolean postRegression (
  231.                     final org.drip.regression.core.RegressionRunDetail rnvd)
  232.                 {
  233.                     for (int i = 0; i < aCompCalib.length; ++i) {
  234.                         if (!rnvd.set ("PVDFMicroJack_" + aCompCalib[i].name(), aWJComp[i].displayString()))
  235.                             return false;
  236.                     }

  237.                     return rnvd.set ("CompPVDFJacobian", wjPVDF.displayString());
  238.                 }
  239.             });
  240.         } catch (java.lang.Exception e) {
  241.             e.printStackTrace();
  242.         }

  243.         return true;
  244.     }

  245.     @Override public java.lang.String getSetName()
  246.     {
  247.         return _strRegressionScenario;
  248.     }
  249. }