docs: 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;
    }