SegmentBasisFunctionGenerator.java

  1. package org.drip.spline.bspline;

  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>SegmentBasisFunctionGenerator</i> generates B Spline Functions of different order. It provides the
  82.  * following functionality:
  83.  *  <ul>
  84.  *      <li>
  85.  *          Create a Tension Monic B Spline Basis Function.
  86.  *      </li>
  87.  *      <li>
  88.  *          Construct a Sequence of Monic Basis Functions.
  89.  *      </li>
  90.  *      <li>
  91.  *          Create a sequence of B Splines of the specified order from the given inputs.
  92.  *      </li>
  93.  *  </ul>
  94.  *
  95.  * <br><br>
  96.  *  <ul>
  97.  *      <li><b>Module </b> = <a href = "https://github.com/lakshmiDRIP/DROP/tree/master/ComputationalCore.md">Computational Core Module</a></li>
  98.  *      <li><b>Library</b> = <a href = "https://github.com/lakshmiDRIP/DROP/tree/master/SplineBuilderLibrary.md">Spline Builder Library</a></li>
  99.  *      <li><b>Project</b> = <a href = "https://github.com/lakshmiDRIP/DROP/tree/master/src/main/java/org/drip/spline/README.md">Basis Splines and Linear Compounders across a Broad Family of Spline Basis Functions</a></li>
  100.  *      <li><b>Package</b> = <a href = "https://github.com/lakshmiDRIP/DROP/tree/master/src/main/java/org/drip/spline/bspline/README.md">de Boor Rational/Exponential/Tension B-Splines</a></li>
  101.  *  </ul>
  102.  * <br><br>
  103.  *
  104.  * @author Lakshmi Krishnamurthy
  105.  */

  106. public class SegmentBasisFunctionGenerator {

  107.     /**
  108.      * Create a Tension Monic B Spline Basis Function
  109.      *
  110.      * @param strHatType The Primitive Hat Type
  111.      * @param strShapeControlType Type of the Shape Controller to be used - NONE, LINEAR/QUADRATIC Rational
  112.      * @param adblPredictorOrdinate Array of Predictor Ordinates
  113.      * @param iDerivOrder The Derivative Order
  114.      * @param dblTension Tension
  115.      *
  116.      * @return The Tension Monic B Spline Basis Function Instance
  117.      */

  118.     public static final org.drip.spline.bspline.SegmentBasisFunction Monic (
  119.         final java.lang.String strHatType,
  120.         final java.lang.String strShapeControlType,
  121.         final double[] adblPredictorOrdinate,
  122.         final int iDerivOrder,
  123.         final double dblTension)
  124.     {
  125.         org.drip.spline.bspline.TensionBasisHat[] aTBH =
  126.             org.drip.spline.bspline.BasisHatPairGenerator.GenerateHatPair (strHatType, strShapeControlType,
  127.                 adblPredictorOrdinate[0], adblPredictorOrdinate[1], adblPredictorOrdinate[2], iDerivOrder,
  128.                     dblTension);

  129.         if (null == aTBH || 2 != aTBH.length) return null;

  130.         try {
  131.             return new org.drip.spline.bspline.SegmentMonicBasisFunction (aTBH[0], aTBH[1]);
  132.         } catch (java.lang.Exception e) {
  133.             e.printStackTrace();
  134.         }

  135.         return null;
  136.     }

  137.     /**
  138.      * Construct a Sequence of Monic Basis Functions
  139.      *
  140.      * @param strHatType The Primitive Hat Type
  141.      * @param strShapeControlType Type of the Shape Controller to be used - NONE, LINEAR/QUADRATIC Rational
  142.      * @param adblPredictorOrdinate Array of Predictor Ordinates
  143.      * @param iDerivOrder The Derivative Order
  144.      * @param dblTension Tension
  145.      *
  146.      * @return Sequence of Tension Monic B Spline Basis Functions
  147.      */

  148.     public static final org.drip.spline.bspline.SegmentBasisFunction[] MonicSequence (
  149.         final java.lang.String strHatType,
  150.         final java.lang.String strShapeControlType,
  151.         final double[] adblPredictorOrdinate,
  152.         final int iDerivOrder,
  153.         final double dblTension)
  154.     {
  155.         if (null == adblPredictorOrdinate) return null;

  156.         int iNumMonic = adblPredictorOrdinate.length - 2;
  157.         org.drip.spline.bspline.SegmentBasisFunction[] aSBFMonic = new
  158.             org.drip.spline.bspline.SegmentBasisFunction[iNumMonic];

  159.         if (0 >= iNumMonic) return null;

  160.         for (int i = 0; i < iNumMonic; ++i) {
  161.             TensionBasisHat[] aTBH = BasisHatPairGenerator.GenerateHatPair (strHatType, strShapeControlType,
  162.                 adblPredictorOrdinate[i], adblPredictorOrdinate[i + 1], adblPredictorOrdinate[i + 2],
  163.                     iDerivOrder, dblTension);

  164.             if (null == aTBH || 2 != aTBH.length) return null;

  165.             try {
  166.                 aSBFMonic[i] = new org.drip.spline.bspline.SegmentMonicBasisFunction (aTBH[0], aTBH[1]);
  167.             } catch (java.lang.Exception e) {
  168.                 e.printStackTrace();

  169.                 return null;
  170.             }
  171.         }

  172.         return aSBFMonic;
  173.     }

  174.     /**
  175.      * Create a sequence of B Splines of the specified order from the given inputs.
  176.      *
  177.      * @param iTargetBSplineOrder Desired B Spline Order
  178.      * @param aSBFPrev Array of Segment Basis Functions
  179.      *
  180.      * @return The sequence of B Splines of the specified order.
  181.      */

  182.     public static final org.drip.spline.bspline.SegmentBasisFunction[] MulticSequence (
  183.         final int iTargetBSplineOrder,
  184.         final org.drip.spline.bspline.SegmentBasisFunction[] aSBFPrev)
  185.     {
  186.         if (2 >= iTargetBSplineOrder || null == aSBFPrev) return null;

  187.         int iNumSBF = aSBFPrev.length - 1;
  188.         org.drip.spline.bspline.SegmentBasisFunction[] aSBF = new
  189.             org.drip.spline.bspline.SegmentBasisFunction[iNumSBF];

  190.         if (1 > iNumSBF) return null;

  191.         for (int i = 0; i < iNumSBF; ++i) {
  192.             try {
  193.                 aSBF[i] = new org.drip.spline.bspline.SegmentMulticBasisFunction (aSBFPrev[i],
  194.                     aSBFPrev[i + 1]);
  195.             } catch (java.lang.Exception e) {
  196.                 e.printStackTrace();

  197.                 return null;
  198.             }
  199.         }

  200.         return iTargetBSplineOrder == aSBF[0].bSplineOrder() ? aSBF : MulticSequence (iTargetBSplineOrder,
  201.             aSBF);
  202.     }
  203. }