docs: private double getIntegratedDriftAdjustment(int timeIndex) {
  		double integratedDriftAdjustment = 0;
  		for(int i=1; i<=timeIndex; i++) {
  			double t = getProcess().getTime(i-1);
  			double t2 = getProcess().getTime(i);
  
  			int timeIndexVolatilityModel = volatilityModel.getTimeDiscretization().getTimeIndex(t);
  			if(timeIndexVolatilityModel < 0)
  			{
  				timeIndexVolatilityModel = -timeIndexVolatilityModel-2;	// Get timeIndex corresponding to previous point
  			}
  			double meanReversion = volatilityModel.getMeanReversion(timeIndexVolatilityModel).doubleValue();
  
  			integratedDriftAdjustment += getShortRateConditionalVariance(0, t) * getB(t,t2)/(t2-t) * (t2-t) - integratedDriftAdjustment * meanReversion * (t2-t) * getB(t,t2)/(t2-t);
  		}
  		return integratedDriftAdjustment;
  	}