/**
* Copyright (C) 2014 - present by OpenGamma Inc. and the OpenGamma group of companies
*
* Please see distribution for license.
*/
package com.opengamma.analytics.financial.interestrate.payments.provider;
import java.util.ArrayList;
import java.util.HashMap;
import java.util.List;
import java.util.Map;
import com.opengamma.analytics.financial.interestrate.payments.derivative.CouponIborAverageFixingDatesCompounding;
import com.opengamma.analytics.financial.provider.description.interestrate.MulticurveProviderInterface;
import com.opengamma.analytics.financial.provider.sensitivity.multicurve.ForwardSensitivity;
import com.opengamma.analytics.financial.provider.sensitivity.multicurve.MulticurveSensitivity;
import com.opengamma.analytics.financial.provider.sensitivity.multicurve.MultipleCurrencyMulticurveSensitivity;
import com.opengamma.analytics.financial.provider.sensitivity.multicurve.SimplyCompoundedForwardSensitivity;
import com.opengamma.util.ArgumentChecker;
import com.opengamma.util.money.MultipleCurrencyAmount;
import com.opengamma.util.tuple.DoublesPair;
/**
*
*/
public final class CouponIborAverageFixingDatesCompoundingDiscountingMethod {
/**
* The method unique instance.
*/
private static final CouponIborAverageFixingDatesCompoundingDiscountingMethod INSTANCE = new CouponIborAverageFixingDatesCompoundingDiscountingMethod();
/**
* Return the unique instance of the class.
* @return The instance.
*/
public static CouponIborAverageFixingDatesCompoundingDiscountingMethod getInstance() {
return INSTANCE;
}
/**
* Private constructor.
*/
private CouponIborAverageFixingDatesCompoundingDiscountingMethod() {
}
/**
* Compute the present value of a Ibor average coupon by discounting.
* @param coupon The coupon.
* @param multicurves The multi-curve provider.
* @return The present value.
*/
public MultipleCurrencyAmount presentValue(final CouponIborAverageFixingDatesCompounding coupon, final MulticurveProviderInterface multicurves) {
ArgumentChecker.notNull(coupon, "Coupon");
ArgumentChecker.notNull(multicurves, "Multi-curves provider");
final int nPeriods = coupon.getFixingTime().length;
double payOff = coupon.getInvestmentFactor();
for (int i = 0; i < nPeriods; ++i) {
double forward = ((i == 0) ? coupon.getAmountAccrued() : 0.0);
final int nDates = coupon.getFixingTime()[i].length;
for (int j = 0; j < nDates; ++j) {
final double forward1 = multicurves.getSimplyCompoundForwardRate(coupon.getIndex(), coupon.getFixingPeriodStartTime()[i][j], coupon.getFixingPeriodEndTime()[i][j],
coupon.getFixingPeriodAccrualFactor()[i][j]);
forward += coupon.getWeight()[i][j] * forward1;
}
payOff *= (1.0 + forward * coupon.getPaymentAccrualFactors()[i]);
}
final double df = multicurves.getDiscountFactor(coupon.getCurrency(), coupon.getPaymentTime());
final double pv = coupon.getNotional() * (payOff - 1.0) * df;
return MultipleCurrencyAmount.of(coupon.getCurrency(), pv);
}
/**
* Compute the present value sensitivity to yield for discounting curve and forward rate (in index convention) for forward curve.
* @param coupon The coupon.
* @param multicurve The multi-curve provider.
* @return The present value curve sensitivity.
*/
public MultipleCurrencyMulticurveSensitivity presentValueCurveSensitivity(final CouponIborAverageFixingDatesCompounding coupon, final MulticurveProviderInterface multicurve) {
ArgumentChecker.notNull(coupon, "Coupon");
ArgumentChecker.notNull(multicurve, "Curves");
final int nPeriods = coupon.getFixingTime().length;
final int[] nDates = new int[nPeriods];
final double[] fctr = new double[nPeriods];
double payOff = coupon.getInvestmentFactor();
for (int i = 0; i < nPeriods; ++i) {
double forward = ((i == 0) ? coupon.getAmountAccrued() : 0.0);
nDates[i] = coupon.getFixingTime()[i].length;
for (int j = 0; j < nDates[i]; ++j) {
final double forward1 = multicurve.getSimplyCompoundForwardRate(coupon.getIndex(), coupon.getFixingPeriodStartTime()[i][j], coupon.getFixingPeriodEndTime()[i][j],
coupon.getFixingPeriodAccrualFactor()[i][j]);
forward += coupon.getWeight()[i][j] * forward1;
}
fctr[i] = (1.0 + forward * coupon.getPaymentAccrualFactors()[i]);
payOff *= fctr[i];
}
final double df = multicurve.getDiscountFactor(coupon.getCurrency(), coupon.getPaymentTime());
// Backward sweep.
final double payOffBar = coupon.getNotional() * df;
final double dfBar = coupon.getNotional() * (payOff - 1.0);
final double[][] forwardBars = new double[nPeriods][];
for (int i = 0; i < nPeriods; ++i) {
forwardBars[i] = new double[nDates[i]];
for (int j = 0; j < nDates[i]; ++j) {
forwardBars[i][j] = coupon.getWeight()[i][j] * coupon.getPaymentAccrualFactors()[i] * payOff / fctr[i] * payOffBar;
}
}
final Map<String, List<DoublesPair>> mapDsc = new HashMap<>();
final List<DoublesPair> listDiscounting = new ArrayList<>();
listDiscounting.add(DoublesPair.of(coupon.getPaymentTime(), -coupon.getPaymentTime() * df * dfBar));
mapDsc.put(multicurve.getName(coupon.getCurrency()), listDiscounting);
final Map<String, List<ForwardSensitivity>> mapFwd = new HashMap<>();
final List<ForwardSensitivity> listForward = new ArrayList<>();
for (int i = 0; i < nPeriods; ++i) {
for (int j = 0; j < nDates[i]; ++j) {
listForward.add(new SimplyCompoundedForwardSensitivity(coupon.getFixingPeriodStartTime()[i][j],
coupon.getFixingPeriodEndTime()[i][j], coupon.getFixingPeriodAccrualFactor()[i][j], forwardBars[i][j]));
}
}
mapFwd.put(multicurve.getName(coupon.getIndex()), listForward);
return MultipleCurrencyMulticurveSensitivity.of(coupon.getCurrency(), MulticurveSensitivity.of(mapDsc, mapFwd));
}
}