All Downloads are FREE. Search and download functionalities are using the official Maven repository.

net.finmath.montecarlo.interestrate.models.covariance.AbstractShortRateVolatilityModelParametric Maven / Gradle / Ivy

/*
 * (c) Copyright Christian P. Fries, Germany. Contact: [email protected].
 *
 * Created on 20.05.2006
 */
package net.finmath.montecarlo.interestrate.models.covariance;

import java.text.DecimalFormat;
import java.text.Format;
import java.util.ArrayList;
import java.util.Arrays;
import java.util.HashMap;
import java.util.Map;
import java.util.concurrent.Callable;
import java.util.concurrent.ExecutionException;
import java.util.concurrent.ExecutorService;
import java.util.concurrent.Future;
import java.util.concurrent.FutureTask;
import java.util.logging.Level;
import java.util.logging.Logger;

import net.finmath.exception.CalculationException;
import net.finmath.montecarlo.BrownianMotion;
import net.finmath.montecarlo.BrownianMotionFromMersenneRandomNumbers;
import net.finmath.montecarlo.interestrate.CalibrationProduct;
import net.finmath.montecarlo.interestrate.LIBORMonteCarloSimulationFromLIBORModel;
import net.finmath.montecarlo.interestrate.ShortRateModel;
import net.finmath.montecarlo.interestrate.models.HullWhiteModel;
import net.finmath.montecarlo.process.EulerSchemeFromProcessModel;
import net.finmath.optimizer.Optimizer;
import net.finmath.optimizer.Optimizer.ObjectiveFunction;
import net.finmath.optimizer.OptimizerFactory;
import net.finmath.optimizer.OptimizerFactoryLevenbergMarquardt;
import net.finmath.optimizer.SolverException;
import net.finmath.stochastic.RandomVariable;
import net.finmath.time.TimeDiscretization;

/**
 * Base class for parametric volatility models, see also {@link AbstractShortRateVolatilityModel}.
 *
 * Parametric models feature a parameter vector which can be inspected
 * and modified for calibration purposes.
 *
 * The parameter vector may have zero length, which indicated that the model
 * is not calibrateable.
 *
 * This class includes the implementation of a generic calibration algorithm.
 * If you provide an arbitrary list of calibration products, the class can return
 * a new instance where the parameters are chosen such that the (weighted) root-mean-square
 * error of the difference of the value of the calibration products and given target
 * values is minimized.
 *
 * @author Ruben Duarte
 * @version 1.1
 */
public abstract class AbstractShortRateVolatilityModelParametric extends AbstractShortRateVolatilityModel implements ShortRateVolatilityModelParametric, ShortRateVolatilityModelCalibrateable {

	private static final long serialVersionUID = 7015719361182945464L;

	private static final Logger logger = Logger.getLogger("net.finmath");

	/**
	 * Constructor consuming time discretization.
	 *
	 * @param timeDiscretization The vector of simulation time discretization points.
	 */
	public AbstractShortRateVolatilityModelParametric(final TimeDiscretization timeDiscretization) {
		super(timeDiscretization);
	}

	/**
	 * Get the parameters of determining this parametric
	 * volatility model. The parameters are usually free parameters
	 * which may be used in calibration.
	 *
	 * @return Parameter vector.
	 */
	@Override
	public abstract RandomVariable[]	getParameter();


	@Override
	public abstract Object clone();

	@Override
	public double[] getParameterAsDouble() {
		final RandomVariable[] parameters = getParameter();
		final double[] parametersAsDouble = new double[parameters.length];
		for(int i=0; iStrings):
	 * 
    *
  • brownianMotion: Under this key an object implementing {@link net.finmath.montecarlo.BrownianMotion} may be provided. If so, this Brownian motion is used to build the valuation model.
  • *
  • maxIterations: Under this key an object of type Integer may be provided specifying the maximum number of iterations.
  • *
  • accuracy: Under this key an object of type Double may be provided specifying the desired accuracy. Note that this is understood in the sense that the solver will stop if the iteration does not improve by more than this number.
  • *
* * @param calibrationModel The LIBOR market model to be used for calibrations (specifies forward curve and tenor discretization). * @param calibrationProducts The array of calibration products. * @param calibrationParameters A map of type Map<String, Object> specifying some (optional) calibration parameters. * @return A new parametric model of the same type than this one, but with calibrated parameters. * @throws CalculationException Thrown if calibration has failed. */ /* (non-Javadoc) * @see net.finmath.montecarlo.interestrate.models.covariance.ShortRateVolatilityModelCalibrateable#getCloneCalibrated(net.finmath.montecarlo.interestrate.ShortRateModel, net.finmath.montecarlo.interestrate.CalibrationProduct[], java.util.Map) */ @Override public AbstractShortRateVolatilityModelParametric getCloneCalibrated(final ShortRateModel calibrationModel, final CalibrationProduct[] calibrationProducts, final Map calibrationParameters) throws CalculationException { return getCloneCalibratedLegazy(calibrationModel, calibrationProducts, calibrationParameters); // if(calibrationParameters == null) { // calibrationParameters = new HashMap<>(); // } // int numberOfPaths = (Integer)calibrationParameters.getOrDefault("numberOfPaths", 2000); // int seed = (Integer)calibrationParameters.getOrDefault("seed", 31415); // int numberOfFactors = 2; // Hull-White model implementation uses two factors for exact discretization. // int maxIterations = (Integer)calibrationParameters.getOrDefault("maxIterations", 400); // double accuracy = (Double)calibrationParameters.getOrDefault("accuracy", 1E-7); // final BrownianMotion brownianMotion = (BrownianMotion)calibrationParameters.getOrDefault("brownianMotion", new BrownianMotionFromMersenneRandomNumbers(getTimeDiscretization(), numberOfFactors, numberOfPaths, seed)); // // RandomVariable[] initialParameters = this.getParameter(); // RandomVariable[] lowerBound = new RandomVariable[initialParameters.length]; // RandomVariable[] upperBound = new RandomVariable[initialParameters.length]; // RandomVariable[] parameterStep = new RandomVariable[initialParameters.length]; // Arrays.fill(lowerBound, new RandomVariableFromDoubleArray(Double.NEGATIVE_INFINITY)); // Arrays.fill(upperBound, new RandomVariableFromDoubleArray(Double.POSITIVE_INFINITY)); // Double parameterStepParameter = (Double)calibrationParameters.get("parameterStep"); // Arrays.fill(parameterStep, new RandomVariableFromDoubleArray(parameterStepParameter != null ? parameterStepParameter.doubleValue() : 1E-4)); // // /* // * We allow for 2 simultaneous calibration models. // * Note: In the case of a Monte-Carlo calibration, the memory requirement is that of // * one model with 2 times the number of paths. In the case of an analytic calibration // * memory requirement is not the limiting factor. // */ // int numberOfThreads = 2; // Object optimizerFactory = calibrationParameters.getOrDefault("optimizerFactory", new OptimizerFactoryLevenbergMarquardt(maxIterations, accuracy, numberOfThreads)); // OptimizerFactory optimizerFactoryParameter = (OptimizerFactory)calibrationParameters.get("optimizerFactory"); // // final ExecutorService executor = null; // // int numberOfThreadsForProductValuation = 2 * Math.max(2, Runtime.getRuntime().availableProcessors()); // // Executors.newFixedThreadPool(numberOfThreadsForProductValuation); // // ObjectiveFunction calibrationError = new ObjectiveFunction() { // // Calculate model values for given parameters // @Override // public void setValues(double[] parameters, double[] values) throws SolverException { // // AbstractShortRateVolatilityModelParametric calibrationVolatilityModel = AbstractShortRateVolatilityModelParametric.this.getCloneWithModifiedParameters(parameters); // // // Create a HullWhiteModel with the new volatility structure. // // TODO the case has be removed after the interface has been refactored: // HullWhiteModel model = (HullWhiteModel)calibrationModel.getCloneWithModifiedVolatilityModel(calibrationVolatilityModel); // EulerSchemeFromProcessModel process = new EulerSchemeFromProcessModel(brownianMotion); // final LIBORMonteCarloSimulationFromLIBORModel modelMonteCarloSimulation = new LIBORMonteCarloSimulationFromLIBORModel(model, process); // // ArrayList> valueFutures = new ArrayList<>(calibrationProducts.length); // for(int calibrationProductIndex=0; calibrationProductIndex worker = new Callable() { // @Override // public RandomVariable call() { // try { // return calibrationProducts[workerCalibrationProductIndex].getProduct().getValue(0.0, modelMonteCarloSimulation).sub(calibrationProducts[workerCalibrationProductIndex].getTargetValue()).mult(calibrationProducts[workerCalibrationProductIndex].getWeight()); // } catch (CalculationException e) { // // We do not signal exceptions to keep the solver working and automatically exclude non-working calibration products. // return null; // } catch (Exception e) { // // We do not signal exceptions to keep the solver working and automatically exclude non-working calibration products. // return null; // } // } // }; // if(executor != null) { // Future valueFuture = executor.submit(worker); // valueFutures.add(calibrationProductIndex, valueFuture); // } // else { // FutureTask valueFutureTask = new FutureTask<>(worker); // valueFutureTask.run(); // valueFutures.add(calibrationProductIndex, valueFutureTask); // } // } // for(int calibrationProductIndex=0; calibrationProductIndexStrings): *
    *
  • brownianMotion: Under this key an object implementing {@link net.finmath.montecarlo.BrownianMotion} may be provided. If so, this Brownian motion is used to build the valuation model.
  • *
  • maxIterations: Under this key an object of type Integer may be provided specifying the maximum number of iterations.
  • *
  • accuracy: Under this key an object of type Double may be provided specifying the desired accuracy. Note that this is understood in the sense that the solver will stop if the iteration does not improve by more than this number.
  • *
* * @param calibrationModel The LIBOR market model to be used for calibrations (specifies forward curve and tenor discretization). * @param calibrationProducts The array of calibration products. * @param calibrationParameters A map of type Map<String, Object> specifying some (optional) calibration parameters. * @return A new parametric model of the same type than this one, but with calibrated parameters. * @throws CalculationException Thrown if calibration has failed. */ public AbstractShortRateVolatilityModelParametric getCloneCalibratedLegazy(final ShortRateModel calibrationModel, final CalibrationProduct[] calibrationProducts, Map calibrationParameters) throws CalculationException { if(calibrationParameters == null) { calibrationParameters = new HashMap<>(); } final Integer numberOfPathsParameter = (Integer)calibrationParameters.get("numberOfPaths"); final Integer seedParameter = (Integer)calibrationParameters.get("seed"); final Integer maxIterationsParameter = (Integer)calibrationParameters.get("maxIterations"); final Double parameterStepParameter = (Double)calibrationParameters.get("parameterStep"); final Double accuracyParameter = (Double)calibrationParameters.get("accuracy"); final BrownianMotion brownianMotionParameter = (BrownianMotion)calibrationParameters.get("brownianMotion"); final double[] initialParameters = this.getParameterAsDouble(); final double[] lowerBound = new double[initialParameters.length]; final double[] upperBound = new double[initialParameters.length]; final double[] parameterStep = new double[initialParameters.length]; final double[] zero = new double[calibrationProducts.length]; Arrays.fill(lowerBound, Double.NEGATIVE_INFINITY); Arrays.fill(upperBound, Double.POSITIVE_INFINITY); Arrays.fill(parameterStep, parameterStepParameter != null ? parameterStepParameter.doubleValue() : 1E-4); Arrays.fill(zero, 0); /* * We allow for 2 simultaneous calibration models. * Note: In the case of a Monte-Carlo calibration, the memory requirement is that of * one model with 2 times the number of paths. In the case of an analytic calibration * memory requirement is not the limiting factor. */ final int numberOfThreads = 2; final OptimizerFactory optimizerFactoryParameter = (OptimizerFactory)calibrationParameters.get("optimizerFactory"); final int numberOfPaths = numberOfPathsParameter != null ? numberOfPathsParameter.intValue() : 2000; final int seed = seedParameter != null ? seedParameter.intValue() : 31415; final int numberOfFactors = 2; // Hull-White model implementation uses two factors for exact discretization. final int maxIterations = maxIterationsParameter != null ? maxIterationsParameter.intValue() : 400; final double accuracy = accuracyParameter != null ? accuracyParameter.doubleValue() : 1E-7; final BrownianMotion brownianMotion = brownianMotionParameter != null ? brownianMotionParameter : new BrownianMotionFromMersenneRandomNumbers(getTimeDiscretization(), numberOfFactors, numberOfPaths, seed); final OptimizerFactory optimizerFactory = optimizerFactoryParameter != null ? optimizerFactoryParameter : new OptimizerFactoryLevenbergMarquardt(maxIterations, accuracy, numberOfThreads); final ExecutorService executor = null; // int numberOfThreadsForProductValuation = 2 * Math.max(2, Runtime.getRuntime().availableProcessors()); // Executors.newFixedThreadPool(numberOfThreadsForProductValuation); final ObjectiveFunction calibrationError = new ObjectiveFunction() { // Calculate model values for given parameters @Override public void setValues(final double[] parameters, final double[] values) throws SolverException { final ShortRateVolatilityModelParametric calibrationVolatilityModel = AbstractShortRateVolatilityModelParametric.this.getCloneWithModifiedParameters(parameters); // Create a HullWhiteModel with the new volatility structure. // TODO the case has be removed after the interface has been refactored: final HullWhiteModel model = (HullWhiteModel)calibrationModel.getCloneWithModifiedVolatilityModel(calibrationVolatilityModel); final EulerSchemeFromProcessModel process = new EulerSchemeFromProcessModel(model, brownianMotion); final LIBORMonteCarloSimulationFromLIBORModel modelMonteCarloSimulation = new LIBORMonteCarloSimulationFromLIBORModel(model, process); final ArrayList> valueFutures = new ArrayList<>(calibrationProducts.length); for(int calibrationProductIndex=0; calibrationProductIndex worker = new Callable() { @Override public RandomVariable call() { try { return calibrationProducts[workerCalibrationProductIndex].getProduct().getValue(0.0, modelMonteCarloSimulation).sub(calibrationProducts[workerCalibrationProductIndex].getTargetValue()).mult(calibrationProducts[workerCalibrationProductIndex].getWeight()); } catch (final Exception e) { // We do not signal exceptions to keep the solver working and automatically exclude non-working calibration products. return null; } } }; if(executor != null) { final Future valueFuture = executor.submit(worker); valueFutures.add(calibrationProductIndex, valueFuture); } else { final FutureTask valueFutureTask = new FutureTask<>(worker); valueFutureTask.run(); valueFutures.add(calibrationProductIndex, valueFutureTask); } } for(int calibrationProductIndex=0; calibrationProductIndex




© 2015 - 2024 Weber Informatics LLC | Privacy Policy