/**
* Copyright (C) 2012 - present by OpenGamma Inc. and the OpenGamma group of companies
*
* Please see distribution for license.
*/
package com.opengamma.analytics.financial.provider.curve.issuer;
import java.util.Set;
import com.opengamma.analytics.financial.interestrate.InstrumentDerivative;
import com.opengamma.analytics.financial.provider.description.interestrate.IssuerProviderDiscount;
import com.opengamma.analytics.financial.provider.sensitivity.issuer.AbstractParameterSensitivityIssuerMatrixCalculator;
import com.opengamma.analytics.math.function.Function1D;
import com.opengamma.analytics.math.matrix.DoubleMatrix1D;
import com.opengamma.analytics.math.matrix.DoubleMatrix2D;
/**
* Function computing the Jacobian of the error of valuation produce by a array representing the curve parameters.
*/
public class IssuerDiscountFinderJacobian extends Function1D<DoubleMatrix1D, DoubleMatrix2D> {
/**
* The instrument parameter sensitivity calculator.
*/
private final AbstractParameterSensitivityIssuerMatrixCalculator _parameterSensitivityCalculator;
/**
* The data required for curve building.
*/
private final IssuerDiscountBuildingData _data;
/**
* Constructor.
* @param parameterSensitivityCalculator The instrument parameter sensitivity calculator.
* @param data The data required for curve building.
*/
public IssuerDiscountFinderJacobian(final AbstractParameterSensitivityIssuerMatrixCalculator parameterSensitivityCalculator, final IssuerDiscountBuildingData data) {
_parameterSensitivityCalculator = parameterSensitivityCalculator;
_data = data;
}
@Override
public DoubleMatrix2D evaluate(DoubleMatrix1D x) {
final IssuerProviderDiscount provider = _data.getKnownData().copy();
IssuerProviderDiscount newCurves = _data.getGeneratorMarket().evaluate(x);
provider.setAll(newCurves);
Set<String> curvesSet = _data.getGeneratorMarket().getCurvesList();
final int nbParameters = _data.getNumberOfInstruments();
final double[][] res = new double[nbParameters][nbParameters];
for (int loopinstrument = 0; loopinstrument < _data.getNumberOfInstruments(); loopinstrument++) {
InstrumentDerivative deriv = _data.getInstrument(loopinstrument);
res[loopinstrument] = _parameterSensitivityCalculator.calculateSensitivity(deriv, provider, curvesSet).getData();
}
return new DoubleMatrix2D(res);
}
}