code: protected static double[][] randomInitialSolution(final int size, final int dim, Random random) {
      double[][] sol = new double[size][dim];
      for(int i = 0; i < size; i++) {
        for(int j = 0; j < dim; j++) {
          sol[i][j] = random.nextGaussian() * INITIAL_SOLUTION_SCALE;
        }
      }
      return sol;
    }
docs: Generate a random initial solution.
  
  @param size Data set size
  @param dim Output dimensionality
  @param random Random generator
  @return Initial solution matrix