final double[][] res = new double[n][];
for (int i = 0; i < n; i++) {
final int m = strikes[i].length;
res[i] = new double[m];
for (int j = 0; j < m; j++) {
final BlackVolatilitySurfaceMoneyness bumpedSurface = _surfaceInterpolator.getBumpedVolatilitySurface(marketData, i, j, SHIFT);
final LocalVolatilitySurfaceMoneyness bumpedLV = _dupireCalculator.getLocalVolatility(bumpedSurface);
final PDETerminalResults1D pdeResBumped = _pdeCalculator.runPDESolver(bumpedLV, option);
for (int k = 0; k < 4; k++) {
vols[k] = BlackFormulaRepository.impliedVolatility(pdeResBumped.getFunctionValue(index + k), 1.0, moneyness[k],
expiry, option.isCall());