Function randMatrix

src/utilities/data/Matrix.cpp:256–286

generates a Matrix of M*N points randomly drawn between and including a and b.

/// generates a Matrix of M*N points randomly drawn between and including a and b.
Matrix randMatrix(double a, double b, unsigned M, unsigned N) {
  Matrix result(M, N);

  // handle degenerate case
  if (equal(a, b)) {
    for (unsigned i = 0; i < M; ++i) {
      for (unsigned j = 0; j < N; ++j) {
        result(i, j) = a;
      }
    }
    return result;
  }

  // ETH@20100120. What library does this come from? The user should be able to seed the
  // generator independently of this function.
  // seed random number generator
  thread_local std::minstd_rand generator(42u);

  // define distribution
  boost::uniform_real<> dist(a, b);

  // create a generator
  boost::variate_generator<std::minstd_rand&, boost::uniform_real<>> uniformGenerator(generator, dist);

  for (unsigned i = 0; i < M; ++i) {
    for (unsigned j = 0; j < N; ++j) {
      result(i, j) = uniformGenerator();
    }
  }
  return result;
}

/// sum
double sum(const Matrix& matrix) {