public void LessComplexCurve()
{
var startDate = new DateTime(2016, 05, 20);
var depoTenors = new Frequency[] { 3.Months() };
var OISdepoTenors = new Frequency[] { 1.Bd() };
double[] depoPricesZAR = { 0.06 };
string[] FRATenors = { "3x6", "6x9", "9x12", "12x15", "15x18", "18x21", "21x24" };
double[] FRAPricesZAR = { 0.065, 0.07, 0.075, 0.077, 0.08, 0.081, 0.082 };
var ZARpillarDatesDepo = depoTenors.Select(x => startDate.AddPeriod(RollType.MF, _jhb, x)).ToArray();
var ZARpillarDatesFRA = FRATenors.Select(x => startDate.AddPeriod(RollType.MF, _jhb, new Frequency(x.Split('x')[1] + "M"))).ToArray();
var ZARpillarDates3m = ZARpillarDatesDepo.Union(ZARpillarDatesFRA).Distinct().OrderBy(x => x).ToArray();
var ZARdepos = new IrSwap[depoTenors.Length];
var ZARFRAs = new ForwardRateAgreement[FRATenors.Length];
var FIC = new FundingInstrumentCollection();
for (var i = 0; i < FRATenors.Length; i++)
{
ZARFRAs[i] = new ForwardRateAgreement(startDate, FRATenors[i], FRAPricesZAR[i], _zar3m, SwapPayReceiveType.Payer, FraDiscountingType.Isda, "ZAR.JIBAR.3M", "ZAR.JIBAR.3M")
{
SolveCurve = "ZAR.JIBAR.3M"
};
FIC.Add(ZARFRAs[i]);
}
for (var i = 0; i < depoTenors.Length; i++)
{
ZARdepos[i] = new IrSwap(startDate, depoTenors[i], _zar3m, depoPricesZAR[i], SwapPayReceiveType.Payer, "ZAR.JIBAR.3M", "ZAR.JIBAR.3M")
{
SolveCurve = "ZAR.JIBAR.3M"
};
FIC.Add(ZARdepos[i]);
}
var ZARcurve3m = new IrCurve(ZARpillarDates3m, new double[ZARpillarDates3m.Length], startDate, "ZAR.JIBAR.3M", Interpolator1DType.LinearFlatExtrap)
{
SolveStage = 0
};
var engine = new FundingModel(startDate, new IrCurve[] { ZARcurve3m });
var S = new NewtonRaphsonMultiCurveSolverStagedWithAnalyticJacobian();
//var S = new NewtonRaphsonMultiCurveSolverStaged();
S.Solve(engine, FIC);
foreach (var ins in FIC)
{
var pv = ins.Pv(engine, false);
Assert.Equal(0.0, pv, 7);
}
}
public void ComplexCurve()
{
var startDate = new DateTime(2016, 05, 20);
var depoTenors = new Frequency[] { 3.Months() };
var OISdepoTenors = new Frequency[] { 1.Bd() };
double[] depoPricesZAR = { 0.06 };
double[] depoPricesUSD = { 0.01 };
double[] OISdepoPricesZAR = { 0.055 };
double[] OISdepoPricesUSD = { 0.009 };
string[] FRATenors = { "3x6", "6x9", "9x12", "12x15", "15x18", "18x21", "21x24" };
double[] FRAPricesZAR = { 0.065, 0.07, 0.075, 0.077, 0.08, 0.081, 0.082 };
double[] FRAPricesUSD = { 0.012, 0.013, 0.014, 0.015, 0.016, 0.017, 0.018 };
Frequency[] swapTenors = { 3.Years(), 4.Years(), 5.Years(), 6.Years(), 7.Years(), 8.Years(), 9.Years(), 10.Years(), 12.Years(), 15.Years(), 20.Years(), 25.Years(), 30.Years() };
double[] swapPricesZAR = { 0.08, 0.083, 0.085, 0.087, 0.089, 0.091, 0.092, 0.093, 0.094, 0.097, 0.099, 0.099, 0.099 };
double[] swapPricesUSD = { 0.017, 0.018, 0.019, 0.020, 0.021, 0.022, 0.023, 0.024, 0.025, 0.026, 0.027, 0.028, 0.03 };
Frequency[] oisTenors = { 3.Months(), 6.Months(), 1.Years(), 18.Months(), 2.Years(), 3.Years(), 4.Years(), 5.Years(), 6.Years(), 7.Years(), 8.Years(), 9.Years(), 10.Years(), 12.Years(), 15.Years(), 20.Years(), 25.Years(), 30.Years() };
double[] oisPricesZAR = { 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004 };
double[] oisPricesUSD = { 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002, 0.002 };
var ZARpillarDatesDepo = depoTenors.Select(x => startDate.AddPeriod(RollType.MF, JHB, x)).ToArray();
var ZARpillarDatesFRA = FRATenors.Select(x => startDate.AddPeriod(RollType.MF, JHB, new Frequency(x.Split('x')[1] + "M"))).ToArray();
var ZARpillarDatesSwap = swapTenors.Select(x => startDate.AddPeriod(RollType.MF, JHB, x)).ToArray();
var ZARpillarDates3m = ZARpillarDatesDepo.Union(ZARpillarDatesSwap).Union(ZARpillarDatesFRA).Distinct().OrderBy(x => x).ToArray();
var ZARpillarDatesDepoOIS = OISdepoTenors.Select(x => startDate.AddPeriod(RollType.MF, JHB, x)).ToArray();
var ZARpillarDatesOISSwap = oisTenors.Select(x => startDate.AddPeriod(RollType.MF, JHB, x)).ToArray();
var ZARpillarDatesOIS = ZARpillarDatesDepoOIS.Union(ZARpillarDatesOISSwap).Distinct().OrderBy(x => x).ToArray();
var USDpillarDatesDepo = depoTenors.Select(x => startDate.AddPeriod(RollType.MF, _usd, x)).ToArray();
var USDpillarDatesFRA = FRATenors.Select(x => startDate.AddPeriod(RollType.MF, _usd, new Frequency(x.Split('x')[1] + "M"))).ToArray();
var USDpillarDatesSwap = swapTenors.Select(x => startDate.AddPeriod(RollType.MF, _usd, x)).ToArray();
var USDpillarDates3m = USDpillarDatesDepo.Union(USDpillarDatesSwap).Union(USDpillarDatesFRA).Distinct().OrderBy(x => x).ToArray();
var USDpillarDatesDepoOIS = OISdepoTenors.Select(x => startDate.AddPeriod(RollType.MF, _usd, x)).ToArray();
var USDpillarDatesOISSwap = oisTenors.Select(x => startDate.AddPeriod(RollType.MF, _usd, x)).ToArray();
var USDpillarDatesOIS = USDpillarDatesDepoOIS.Union(USDpillarDatesOISSwap).Distinct().OrderBy(x => x).ToArray();
var ZARswaps = new IrSwap[swapTenors.Length];
var ZARdepos = new IrSwap[depoTenors.Length];
var ZARdeposOIS = new IrSwap[OISdepoTenors.Length];
var ZARoisSwaps = new IrBasisSwap[oisTenors.Length];
var ZARFRAs = new ForwardRateAgreement[FRATenors.Length];
var USDswaps = new IrSwap[swapTenors.Length];
var USDdepos = new IrSwap[depoTenors.Length];
var USDdeposOIS = new IrSwap[OISdepoTenors.Length];
var USDoisSwaps = new IrBasisSwap[oisTenors.Length];
var USDFRAs = new ForwardRateAgreement[FRATenors.Length];
var FIC = new FundingInstrumentCollection(TestProviderHelper.CurrencyProvider);
for (var i = 0; i < FRATenors.Length; i++)
{
ZARFRAs[i] = new ForwardRateAgreement(startDate, FRATenors[i], FRAPricesZAR[i], _zar3m, SwapPayReceiveType.Payer, FraDiscountingType.Isda, "ZAR.JIBAR.3M", "ZAR.DISC.CSA_ZAR")
{
SolveCurve = "ZAR.JIBAR.3M"
};
FIC.Add(ZARFRAs[i]);
USDFRAs[i] = new ForwardRateAgreement(startDate, FRATenors[i], FRAPricesUSD[i], usd3m, SwapPayReceiveType.Payer, FraDiscountingType.Isda, "USD.LIBOR.3M", "USD.DISC.CSA_USD")
{
SolveCurve = "USD.LIBOR.3M"
};
FIC.Add(USDFRAs[i]);
}
for (var i = 0; i < oisTenors.Length; i++)
{
ZARoisSwaps[i] = new IrBasisSwap(startDate, oisTenors[i], oisPricesZAR[i], true, zaron, _zar3m, "ZAR.JIBAR.3M", "ZAR.DISC.CSA_ZAR", "ZAR.DISC.CSA_ZAR")
{
SolveCurve = "ZAR.DISC.CSA_ZAR"
};
FIC.Add(ZARoisSwaps[i]);
USDoisSwaps[i] = new IrBasisSwap(startDate, oisTenors[i], oisPricesUSD[i], true, usdon, usd3m, "USD.LIBOR.3M", "USD.DISC.CSA_USD", "USD.DISC.CSA_USD")
{
SolveCurve = "USD.DISC.CSA_USD"
};
FIC.Add(USDoisSwaps[i]);
}
for (var i = 0; i < swapTenors.Length; i++)
{
ZARswaps[i] = new IrSwap(startDate, swapTenors[i], _zar3m, swapPricesZAR[i], SwapPayReceiveType.Payer, "ZAR.JIBAR.3M", "ZAR.DISC.CSA_ZAR")
{
SolveCurve = "ZAR.JIBAR.3M"
};
FIC.Add(ZARswaps[i]);
USDswaps[i] = new IrSwap(startDate, swapTenors[i], usd3m, swapPricesUSD[i], SwapPayReceiveType.Payer, "USD.LIBOR.3M", "USD.DISC.CSA_USD")
{
SolveCurve = "USD.LIBOR.3M"
};
FIC.Add(USDswaps[i]);
}
for (var i = 0; i < depoTenors.Length; i++)
{
ZARdepos[i] = new IrSwap(startDate, depoTenors[i], _zar3m, depoPricesZAR[i], SwapPayReceiveType.Payer, "ZAR.JIBAR.3M", "ZAR.DISC.CSA_ZAR")
{
SolveCurve = "ZAR.JIBAR.3M"
};
FIC.Add(ZARdepos[i]);
USDdepos[i] = new IrSwap(startDate, depoTenors[i], usd3m, depoPricesUSD[i], SwapPayReceiveType.Payer, "USD.LIBOR.3M", "USD.DISC.CSA_USD")
{
SolveCurve = "USD.LIBOR.3M"
};
FIC.Add(USDdepos[i]);
}
for (var i = 0; i < OISdepoTenors.Length; i++)
{
ZARdeposOIS[i] = new IrSwap(startDate, OISdepoTenors[i], zaron, OISdepoPricesZAR[i], SwapPayReceiveType.Payer, "ZAR.DISC.CSA_ZAR", "ZAR.DISC.CSA_ZAR")
{
SolveCurve = "ZAR.DISC.CSA_ZAR"
};
FIC.Add(ZARdeposOIS[i]);
USDdeposOIS[i] = new IrSwap(startDate, OISdepoTenors[i], usdon, OISdepoPricesUSD[i], SwapPayReceiveType.Payer, "USD.DISC.CSA_USD", "USD.DISC.CSA_USD")
{
SolveCurve = "USD.DISC.CSA_USD"
};
FIC.Add(USDdeposOIS[i]);
}
var ZARcurve3m = new IrCurve(ZARpillarDates3m, new double[ZARpillarDates3m.Length], startDate, "ZAR.JIBAR.3M", Interpolator1DType.LinearFlatExtrap, ccyZar)
{
SolveStage = 0
};
var ZARcurveOIS = new IrCurve(ZARpillarDatesOIS, new double[ZARpillarDatesOIS.Length], startDate, "ZAR.DISC.CSA_ZAR", Interpolator1DType.LinearFlatExtrap, ccyZar)
{
SolveStage = 0
};
var USDcurve3m = new IrCurve(USDpillarDates3m, new double[USDpillarDates3m.Length], startDate, "USD.LIBOR.3M", Interpolator1DType.LinearFlatExtrap, ccyUsd)
{
SolveStage = 1
};
var USDcurveOIS = new IrCurve(USDpillarDatesOIS, new double[USDpillarDatesOIS.Length], startDate, "USD.DISC.CSA_USD", Interpolator1DType.LinearFlatExtrap, ccyUsd)
{
SolveStage = 1
};
var engine = new FundingModel(startDate, new IrCurve[] { ZARcurve3m, ZARcurveOIS, USDcurve3m, USDcurveOIS }, TestProviderHelper.CurrencyProvider, TestProviderHelper.CalendarProvider);
var ZARcurve3m0 = new IrCurve(ZARpillarDates3m, new double[ZARpillarDates3m.Length], startDate, "ZAR.JIBAR.3M", Interpolator1DType.LinearFlatExtrap, ccyZar)
{
SolveStage = 0
};
var ZARcurveOIS0 = new IrCurve(ZARpillarDatesOIS, new double[ZARpillarDatesOIS.Length], startDate, "ZAR.DISC.CSA_ZAR", Interpolator1DType.LinearFlatExtrap, ccyZar)
{
SolveStage = 0
};
var USDcurve3m0 = new IrCurve(USDpillarDates3m, new double[USDpillarDates3m.Length], startDate, "USD.LIBOR.3M", Interpolator1DType.LinearFlatExtrap, ccyUsd)
{
SolveStage = 1
};
var USDcurveOIS0 = new IrCurve(USDpillarDatesOIS, new double[USDpillarDatesOIS.Length], startDate, "USD.DISC.CSA_USD", Interpolator1DType.LinearFlatExtrap, ccyUsd)
{
SolveStage = 1
};
var engine0 = new FundingModel(startDate, new IrCurve[] { ZARcurve3m0, ZARcurveOIS0, USDcurve3m0, USDcurveOIS0 }, TestProviderHelper.CurrencyProvider, TestProviderHelper.CalendarProvider);
var S = new NewtonRaphsonMultiCurveSolverStagedWithAnalyticJacobian()
{
Tollerance = IsCoverageOnly ? 1 : 0.00000001,
MaxItterations = IsCoverageOnly ? 1 : 100,
};
var S0 = new NewtonRaphsonMultiCurveSolverStaged()
{
Tollerance = IsCoverageOnly ? 1 : 0.00000001,
MaxItterations = IsCoverageOnly ? 1 : 100,
};
S.Solve(engine, FIC);
S0.Solve(engine0, FIC);
if (!IsCoverageOnly)
{
foreach (var ins in FIC)
{
var pv = ins.Pv(engine, false);
Assert.Equal(0.0, pv, 7);
}
foreach (var curve in engine.Curves)
{
var otherCurve = engine0.Curves[curve.Key];
Assert.Equal(curve.Value.NumberOfPillars, otherCurve.NumberOfPillars);
var otherRates = otherCurve.GetRates();
var rates = curve.Value.GetRates();
for (var i = 0; i < otherRates.Length; i++)
{
Assert.Equal(otherRates[i], rates[i], 10);
}
}
}
}