Dear all,
It nothing big, but I thing adding method fill_random(const double& a=-1., const double &b=1.) to vector classes that will fill a vector with random values (from range [a,b]) would be useful. I'm using it in most of my programs, it's great for testing solvers, especially multigrid smoothers.
Best,
Michal