docs: private void doGenerateBrownianMotion() {
  		if(brownianIncrements != null) {
  			return;	// Nothing to do
  		}
  
  		// Create random number sequence generator
  		MersenneTwister mersenneTwister = new MersenneTwister(seed);
  
  		// Allocate memory
  		double[][][] brownianIncrementsArray = new double[timeDiscretization.getNumberOfTimeSteps()][numberOfFactors][numberOfPaths];
  
  		// Pre-calculate square roots of deltaT
  		double[] sqrtOfTimeStep = new double[timeDiscretization.getNumberOfTimeSteps()];
  		for(int timeIndex=0; timeIndex<sqrtOfTimeStep.length; timeIndex++) {
  			sqrtOfTimeStep[timeIndex] = Math.sqrt(timeDiscretization.getTimeStep(timeIndex));
  		}
  
  		/*
  		 * Generate normal distributed independent increments.
  		 *
  		 * The inner loop goes over time and factors.
  		 * MersenneTwister is known to generate "independent" increments in 623 dimensions.
  		 * Since we want to generate independent streams (paths), the loop over path is the outer loop.
  		 */
  		for(int path=0; path<numberOfPaths; path++) {
  			for(int timeIndex=0; timeIndex<timeDiscretization.getNumberOfTimeSteps(); timeIndex++) {
  				double sqrtDeltaT = sqrtOfTimeStep[timeIndex];
  				// Generate uncorrelated Brownian increment
  				for(int factor=0; factor<numberOfFactors; factor++) {
  					double uniformIncrement = mersenneTwister.nextDouble();
  					brownianIncrementsArray[timeIndex][factor][path] = net.finmath.functions.NormalDistribution.inverseCumulativeDistribution(uniformIncrement) * sqrtDeltaT;
  				}
  			}
  		}
  
  		// Allocate memory for RandomVariable wrapper objects.
  		brownianIncrements = new RandomVariableInterface[timeDiscretization.getNumberOfTimeSteps()][numberOfFactors];
  
  		// Wrap the values in RandomVariable objects
  		for(int timeIndex=0; timeIndex<timeDiscretization.getNumberOfTimeSteps(); timeIndex++) {
  			double time = timeDiscretization.getTime(timeIndex+1);
  			for(int factor=0; factor<numberOfFactors; factor++) {
  				brownianIncrements[timeIndex][factor] =
  						randomVariableFactory.createRandomVariable(time, brownianIncrementsArray[timeIndex][factor]);
  			}
  		}
  	}