From 0b9e336ad7b2e9f2144edf00bcd7dbb35daa3d90 Mon Sep 17 00:00:00 2001
From: heister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Fri, 22 Nov 2013 14:42:37 +0000
Subject: [PATCH] rework r31759 to my liking and fix a few more places

git-svn-id: https://svn.dealii.org/trunk@31761 0785d39b-7218-0410-832d-ea1e28bc413d
---
 tests/aniso/anisotropic_crash.cc              |  2 +-
 tests/aniso/refine_and_coarsen_anisotropic.cc |  6 +++---
 tests/base/full_tensor_10.cc                  |  2 +-
 tests/base/geometry_info_1.cc                 |  2 +-
 tests/base/index_set_14.cc                    |  4 ++--
 tests/base/index_set_15.cc                    |  4 ++--
 tests/base/index_set_17.cc                    |  2 +-
 tests/base/index_set_18.cc                    |  2 +-
 tests/base/index_set_19.cc                    |  2 +-
 tests/base/index_set_20.cc                    |  4 ++--
 tests/base/index_set_22.cc                    |  2 +-
 tests/base/polynomial_lagrange_ops.cc         | 18 ++++++++---------
 tests/base/polynomial_lagrange_order.cc       |  2 +-
 tests/base/polynomial_lagrange_product.cc     |  2 +-
 tests/base/symmetric_tensor_19.cc             |  2 +-
 tests/base/symmetric_tensor_25.cc             |  2 +-
 tests/base/symmetric_tensor_26.cc             |  4 ++--
 tests/bits/chunk_sparse_matrix_15.cc          |  2 +-
 tests/bits/fe_q_constraints.cc                |  4 ++--
 tests/bits/find_cell_6.cc                     |  4 ++--
 tests/bits/full_matrix_56.cc                  |  2 +-
 tests/bits/joa_1.cc                           |  4 ++--
 tests/bits/point_value_hp_01.cc               |  2 +-
 tests/bits/point_value_hp_02.cc               |  2 +-
 tests/bits/rt_2.cc                            |  2 +-
 tests/bits/step-15.cc                         |  2 +-
 tests/deal.II/coarsening_02.cc                |  2 +-
 tests/deal.II/coarsening_02_1d.cc             |  2 +-
 tests/deal.II/constraint_graph.cc             |  2 +-
 tests/deal.II/constraints_local_to_global.cc  |  2 +-
 .../constraints_local_to_global_chunk.cc      |  2 +-
 tests/deal.II/dof_handler_number_cache.cc     |  2 +-
 tests/deal.II/dof_handler_number_cache_02.cc  |  2 +-
 tests/distributed_grids/2d_coarsening_02.cc   |  2 +-
 tests/distributed_grids/2d_coarsening_03.cc   |  2 +-
 tests/distributed_grids/2d_coarsening_05.cc   |  2 +-
 tests/distributed_grids/2d_refinement_02.cc   |  2 +-
 tests/distributed_grids/2d_refinement_03.cc   |  2 +-
 tests/distributed_grids/2d_refinement_05.cc   |  2 +-
 tests/distributed_grids/3d_coarsening_02.cc   |  2 +-
 tests/distributed_grids/3d_coarsening_03.cc   |  2 +-
 tests/distributed_grids/3d_coarsening_04.cc   |  2 +-
 tests/distributed_grids/3d_coarsening_05.cc   |  2 +-
 tests/distributed_grids/3d_refinement_02.cc   |  2 +-
 tests/distributed_grids/3d_refinement_03.cc   |  2 +-
 tests/distributed_grids/3d_refinement_05.cc   |  2 +-
 tests/distributed_grids/3d_refinement_08.cc   |  2 +-
 .../count_dofs_per_block_01.cc                |  2 +-
 .../dof_handler_number_cache.cc               |  2 +-
 .../distributed_grids/solution_transfer_03.cc |  4 ++--
 tests/fail/abf_approximation_01.cc            |  4 ++--
 tests/fail/hp-step-15.cc                      |  2 +-
 tests/fail/vectors_boundary_rhs_hp_02.cc      |  2 +-
 tests/fail/vectors_rhs_hp_02.cc               |  2 +-
 tests/fe/rt_8.cc                              |  2 +-
 tests/fe/rt_9.cc                              |  2 +-
 tests/fe/rt_approximation_01.cc               |  4 ++--
 tests/fe/rt_normal_02.cc                      |  2 +-
 tests/hp/count_dofs_per_block_01.cc           |  2 +-
 tests/hp/count_dofs_per_block_02.cc           |  2 +-
 tests/hp/crash_06.cc                          |  2 +-
 tests/hp/crash_10.cc                          |  4 ++--
 tests/hp/crash_16.cc                          |  2 +-
 tests/hp/create_laplace_matrix_01.cc          |  2 +-
 tests/hp/create_laplace_matrix_01b.cc         |  2 +-
 tests/hp/create_laplace_matrix_02.cc          |  2 +-
 tests/hp/create_laplace_matrix_02b.cc         |  2 +-
 tests/hp/create_laplace_matrix_03.cc          |  2 +-
 tests/hp/create_laplace_matrix_03b.cc         |  2 +-
 tests/hp/create_laplace_matrix_04.cc          |  2 +-
 tests/hp/create_laplace_matrix_04b.cc         |  2 +-
 .../create_laplace_matrix_constraints_02.cc   |  2 +-
 .../create_laplace_matrix_constraints_02b.cc  |  2 +-
 tests/hp/create_mass_matrix_01.cc             |  2 +-
 tests/hp/create_mass_matrix_01b.cc            |  2 +-
 tests/hp/create_mass_matrix_02.cc             |  2 +-
 tests/hp/create_mass_matrix_02b.cc            |  2 +-
 tests/hp/create_mass_matrix_03.cc             |  2 +-
 tests/hp/create_mass_matrix_03b.cc            |  2 +-
 tests/hp/create_mass_matrix_04.cc             |  2 +-
 tests/hp/create_mass_matrix_04b.cc            |  2 +-
 tests/hp/create_mass_matrix_05.cc             |  2 +-
 tests/hp/create_mass_matrix_constraints_02.cc |  2 +-
 .../hp/create_mass_matrix_constraints_02b.cc  |  2 +-
 tests/hp/dof_handler_number_cache.cc          |  2 +-
 tests/hp/get_active_fe_indices.cc             |  2 +-
 tests/hp/hp_constraints_common.h              |  8 ++++----
 tests/hp/hp_hanging_nodes_02.cc               |  2 +-
 tests/hp/integrate_difference.cc              |  2 +-
 tests/hp/integrate_difference_02.cc           |  2 +-
 tests/hp/matrices_hp.cc                       |  2 +-
 tests/hp/n_active_fe_indices.cc               |  2 +-
 tests/hp/random.cc                            |  2 +-
 tests/hp/solution_transfer.cc                 | 12 +++++------
 tests/hp/step-3b.cc                           |  2 +-
 tests/hp/step-3c.cc                           |  2 +-
 tests/hp/vectors_boundary_rhs_01.cc           |  2 +-
 tests/hp/vectors_boundary_rhs_02.cc           |  2 +-
 tests/hp/vectors_boundary_rhs_03.cc           |  2 +-
 tests/hp/vectors_boundary_rhs_hp_01.cc        |  2 +-
 tests/hp/vectors_boundary_rhs_hp_03.cc        |  2 +-
 tests/hp/vectors_point_source_01.cc           |  2 +-
 tests/hp/vectors_point_source_hp_01.cc        |  2 +-
 tests/hp/vectors_rhs_01.cc                    |  2 +-
 tests/hp/vectors_rhs_02.cc                    |  2 +-
 tests/hp/vectors_rhs_03.cc                    |  2 +-
 tests/hp/vectors_rhs_hp_01.cc                 |  2 +-
 tests/hp/vectors_rhs_hp_03.cc                 |  2 +-
 tests/lac/full_matrix.cc                      |  2 +-
 tests/lac/full_matrix_04.cc                   |  8 ++++----
 tests/lac/gmres_reorthogonalize_01.cc         |  2 +-
 tests/lac/gmres_reorthogonalize_03.cc         |  2 +-
 tests/lac/precondition_chebyshev_01.cc        |  2 +-
 tests/lac/sparse_ilu.cc                       |  2 +-
 tests/lac/sparse_ilu_t.cc                     |  2 +-
 tests/lac/sparse_matrix_Tmmult_01.cc          |  6 +++---
 tests/lac/sparse_matrix_Tmmult_02.cc          |  8 ++++----
 tests/lac/sparse_matrix_mmult_01.cc           |  6 +++---
 tests/lac/sparse_matrix_mmult_02.cc           |  8 ++++----
 tests/lac/sparse_mic.cc                       |  2 +-
 tests/lac/vector_accumulation.cc              |  4 ++--
 tests/matrix_free/copy.cc                     |  2 +-
 tests/matrix_free/evaluate_1d_shape.cc        |  6 +++---
 .../matrix_free/evaluate_1d_shape_evenodd.cc  |  6 +++---
 tests/matrix_free/get_functions_common.h      |  2 +-
 tests/matrix_free/get_functions_gl.cc         |  2 +-
 tests/matrix_free/get_functions_mappingq.cc   |  2 +-
 tests/matrix_free/get_functions_multife.cc    |  2 +-
 tests/matrix_free/get_functions_multife2.cc   |  2 +-
 .../get_functions_q_hierarchical.cc           |  2 +-
 tests/matrix_free/get_functions_rect.cc       |  2 +-
 tests/matrix_free/get_functions_variants.cc   |  2 +-
 tests/matrix_free/get_values_plain.cc         |  2 +-
 tests/matrix_free/integrate_functions.cc      |  4 ++--
 .../integrate_functions_multife.cc            |  8 ++++----
 .../integrate_functions_multife2.cc           |  8 ++++----
 tests/matrix_free/matrix_vector_10.cc         |  2 +-
 tests/matrix_free/matrix_vector_11.cc         |  2 +-
 tests/matrix_free/matrix_vector_12.cc         |  4 ++--
 tests/matrix_free/matrix_vector_13.cc         |  4 ++--
 tests/matrix_free/matrix_vector_common.h      |  2 +-
 tests/matrix_free/matrix_vector_curl.cc       |  2 +-
 tests/matrix_free/matrix_vector_div.cc        |  2 +-
 tests/matrix_free/matrix_vector_hp.cc         |  4 ++--
 tests/matrix_free/matrix_vector_mg.cc         |  4 ++--
 tests/matrix_free/matrix_vector_stokes.cc     |  2 +-
 .../matrix_vector_stokes_noflux.cc            |  4 ++--
 tests/matrix_free/no_index_initialize.cc      |  2 +-
 tests/matrix_free/thread_correctness.cc       |  2 +-
 tests/matrix_free/thread_correctness_hp.cc    |  6 +++---
 tests/matrix_free/update_mapping_only.cc      |  2 +-
 tests/mpi/constraint_matrix_trilinos_bug.cc   |  2 +-
 tests/mpi/count_dofs_per_block_02.cc          |  2 +-
 tests/mpi/crash_03.cc                         |  2 +-
 tests/mpi/dof_handler_number_cache.cc         |  2 +-
 tests/mpi/mg_05.cc                            |  2 +-
 tests/mpi/mg_06.cc                            |  2 +-
 tests/mpi/p4est_2d_constraintmatrix_03.cc     |  2 +-
 tests/mpi/p4est_2d_constraintmatrix_04.cc     |  6 +++---
 tests/mpi/p4est_2d_dofhandler_02.cc           |  2 +-
 tests/mpi/p4est_2d_dofhandler_03.cc           |  2 +-
 tests/mpi/p4est_2d_dofhandler_04.cc           |  2 +-
 tests/mpi/p4est_2d_refine_03.cc               |  2 +-
 tests/mpi/p4est_3d_constraintmatrix_01.cc     |  2 +-
 tests/mpi/p4est_3d_constraintmatrix_02.cc     |  2 +-
 tests/mpi/p4est_3d_constraintmatrix_03.cc     |  6 +++---
 tests/mpi/p4est_3d_refine_01.cc               |  2 +-
 tests/mpi/trilinos_matvec_01.cc               |  2 +-
 tests/mpi/trilinos_matvec_02.cc               |  2 +-
 tests/mpi/trilinos_matvec_03.cc               |  2 +-
 tests/petsc/66.cc                             |  6 +++---
 tests/petsc/67.cc                             |  6 +++---
 tests/petsc/68.cc                             |  8 ++++----
 tests/petsc/69.cc                             |  8 ++++----
 tests/petsc/slowness_04.cc                    | 20 +++++++++----------
 tests/tests.h                                 | 12 ++++++-----
 tests/trilinos/66.cc                          |  6 +++---
 tests/trilinos/67.cc                          |  6 +++---
 tests/trilinos/68.cc                          |  8 ++++----
 tests/trilinos/69.cc                          |  8 ++++----
 tests/trilinos/slowness_04.cc                 | 20 +++++++++----------
 181 files changed, 290 insertions(+), 288 deletions(-)

diff --git a/tests/aniso/anisotropic_crash.cc b/tests/aniso/anisotropic_crash.cc
index 9429c1d5b7..35f33a1635 100644
--- a/tests/aniso/anisotropic_crash.cc
+++ b/tests/aniso/anisotropic_crash.cc
@@ -55,7 +55,7 @@ int main()
   Triangulation<2>::active_cell_iterator cell = tri.begin_active(), end = tri.end();
   for ( ; cell != end; ++cell )
     {
-      switch (myrand()%4)
+      switch (Testing::rand()%4)
         {
           /// If a randomly drawn
           /// number is 0 or 1 we
diff --git a/tests/aniso/refine_and_coarsen_anisotropic.cc b/tests/aniso/refine_and_coarsen_anisotropic.cc
index be0f76b27a..85f1d215dd 100644
--- a/tests/aniso/refine_and_coarsen_anisotropic.cc
+++ b/tests/aniso/refine_and_coarsen_anisotropic.cc
@@ -81,10 +81,10 @@ void test_isotropic (int type, std::ostream *logfile)
     {
       cell = tria.begin_active();
       for (; cell!=endc; ++cell)
-        if (myrand()%5==0)
+        if (Testing::rand()%5==0)
           {
-            if (myrand()%2==0)
-              cell->set_refine_flag(ref_cases[myrand()%
+            if (Testing::rand()%2==0)
+              cell->set_refine_flag(ref_cases[Testing::rand()%
                                               RefinementCase<dim>::isotropic_refinement]);
           }
         else
diff --git a/tests/base/full_tensor_10.cc b/tests/base/full_tensor_10.cc
index 3f45efeeda..5d9d4598b6 100644
--- a/tests/base/full_tensor_10.cc
+++ b/tests/base/full_tensor_10.cc
@@ -33,7 +33,7 @@ void test ()
   // as in symmetric_tensor_10
   for (unsigned int i=0; i<dim; ++i)
     for (unsigned int j=0; j<dim; ++j)
-      t[i][j] = 1. * myrand() / RAND_MAX;
+      t[i][j] = 1. * Testing::rand() / RAND_MAX;
 
   for (unsigned int i=0; i<dim; ++i)
     for (unsigned int j=0; j<dim; ++j)
diff --git a/tests/base/geometry_info_1.cc b/tests/base/geometry_info_1.cc
index ca7e276a88..89938fb79c 100644
--- a/tests/base/geometry_info_1.cc
+++ b/tests/base/geometry_info_1.cc
@@ -27,7 +27,7 @@
 
 double rand_2 ()
 {
-  return 1.*myrand()/RAND_MAX*4-2.;
+  return 1.*Testing::rand()/RAND_MAX*4-2.;
 }
 
 
diff --git a/tests/base/index_set_14.cc b/tests/base/index_set_14.cc
index 3f794dcf82..a69ac9dfb6 100644
--- a/tests/base/index_set_14.cc
+++ b/tests/base/index_set_14.cc
@@ -37,8 +37,8 @@ void test ()
   // repetitions of previous ones
   for (unsigned int i=0; i<9*is1.size()/10; ++i)
     {
-      is1.add_index (myrand() % is1.size());
-      is2.add_index (myrand() % is2.size());
+      is1.add_index (Testing::rand() % is1.size());
+      is2.add_index (Testing::rand() % is2.size());
     }
 
   IndexSet is3 = is1 & is2;
diff --git a/tests/base/index_set_15.cc b/tests/base/index_set_15.cc
index 21c66c48fe..3c05906c65 100644
--- a/tests/base/index_set_15.cc
+++ b/tests/base/index_set_15.cc
@@ -37,8 +37,8 @@ void test ()
   // repetitions of previous ones
   for (unsigned int i=0; i<9*is1.size()/10; ++i)
     {
-      is1.add_index (myrand() % is1.size());
-      is2.add_index (myrand() % is2.size());
+      is1.add_index (Testing::rand() % is1.size());
+      is2.add_index (Testing::rand() % is2.size());
     }
 
   IndexSet is3 = is1;
diff --git a/tests/base/index_set_17.cc b/tests/base/index_set_17.cc
index cd5422a7b2..a495c13d7a 100644
--- a/tests/base/index_set_17.cc
+++ b/tests/base/index_set_17.cc
@@ -36,7 +36,7 @@ void test ()
   // repetitions of previous ones
   for (unsigned int i=0; i<9*is1.size()/10; ++i)
     {
-      is1.add_index (myrand() % is1.size());
+      is1.add_index (Testing::rand() % is1.size());
     }
 
   is1.compress();
diff --git a/tests/base/index_set_18.cc b/tests/base/index_set_18.cc
index 5c101775fb..1b837743ed 100644
--- a/tests/base/index_set_18.cc
+++ b/tests/base/index_set_18.cc
@@ -35,7 +35,7 @@ void test ()
   // set, some of which may be
   // repetitions of previous ones
   for (unsigned int i=0; i<9*is1.size()/10; ++i)
-    is1.add_index (myrand() % is1.size());
+    is1.add_index (Testing::rand() % is1.size());
 
   is1.compress();
 
diff --git a/tests/base/index_set_19.cc b/tests/base/index_set_19.cc
index 04476b4d71..8f5664acd9 100644
--- a/tests/base/index_set_19.cc
+++ b/tests/base/index_set_19.cc
@@ -35,7 +35,7 @@ void test ()
   // set, some of which may be
   // repetitions of previous ones
   for (unsigned int i=0; i<9*is1.size()/10; ++i)
-    is1.add_index (myrand() % is1.size());
+    is1.add_index (Testing::rand() % is1.size());
 
   std::vector<types::global_dof_index> indices;
   is1.fill_index_vector (indices);
diff --git a/tests/base/index_set_20.cc b/tests/base/index_set_20.cc
index 7e985b391d..0a3021a70b 100644
--- a/tests/base/index_set_20.cc
+++ b/tests/base/index_set_20.cc
@@ -86,8 +86,8 @@ void test()
       IndexSet b(size);
       for (unsigned int i=0; i<9*a.size()/10; ++i)
         {
-          a.add_index (myrand() % a.size());
-          b.add_index (myrand() % a.size());
+          a.add_index (Testing::rand() % a.size());
+          b.add_index (Testing::rand() % a.size());
         }
       testor(a, b, false);
     }
diff --git a/tests/base/index_set_22.cc b/tests/base/index_set_22.cc
index ea0ffe37b7..906d4eb1e1 100644
--- a/tests/base/index_set_22.cc
+++ b/tests/base/index_set_22.cc
@@ -35,7 +35,7 @@ void test ()
   // set, some of which may be
   // repetitions of previous ones
   for (unsigned int i=0; i<9*is1.size()/10; ++i)
-    is1.add_index (myrand() % is1.size());
+    is1.add_index (Testing::rand() % is1.size());
 
   std::vector<bool> zeros_and_ones(is1.size());
   is1.fill_binary_vector (zeros_and_ones);
diff --git a/tests/base/polynomial_lagrange_ops.cc b/tests/base/polynomial_lagrange_ops.cc
index 4a4d662f23..75898a59d3 100644
--- a/tests/base/polynomial_lagrange_ops.cc
+++ b/tests/base/polynomial_lagrange_ops.cc
@@ -37,8 +37,8 @@ void check_scale (const std::vector<Polynomial<double> > &p)
   for (unsigned int i=0; i<p.size(); ++i)
     {
       Polynomial<double> q = p[i];
-      double x = (double)myrand()/RAND_MAX;
-      double factor = 5.*(double)myrand()/RAND_MAX;
+      double x = (double)Testing::rand()/RAND_MAX;
+      double factor = 5.*(double)Testing::rand()/RAND_MAX;
       q.scale (factor);
       double value1 = p[i].value (factor * x);
       double value2 = q   .value (x);
@@ -62,8 +62,8 @@ void check_shift (const std::vector<Polynomial<double> > &p)
   for (unsigned int i=0; i<p.size(); ++i)
     {
       Polynomial<double> q = p[i];
-      double x = (double)myrand()/RAND_MAX;
-      double a = 10.*(-1.+2.*(double)myrand()/RAND_MAX);
+      double x = (double)Testing::rand()/RAND_MAX;
+      double a = 10.*(-1.+2.*(double)Testing::rand()/RAND_MAX);
       q.shift (a);
       double value1 = p[i].value (x+a);
       double value2 = q   .value (x);
@@ -83,8 +83,8 @@ void check_mult_scalar (const std::vector<Polynomial<double> > &p)
   for (unsigned int i=0; i<p.size(); ++i)
     {
       Polynomial<double> q = p[i];
-      double x = (double)myrand()/RAND_MAX;
-      double a = (double)myrand()/RAND_MAX;
+      double x = (double)Testing::rand()/RAND_MAX;
+      double a = (double)Testing::rand()/RAND_MAX;
       q *= a;
       double value1 = p[i].value (x) * a;
       double value2 = q   .value (x);
@@ -107,7 +107,7 @@ void check_mult (const std::vector<Polynomial<double> > &p)
         {
           Polynomial<double> q = p[i];
           q *= p[j];
-          double x = (double)myrand()/RAND_MAX;
+          double x = (double)Testing::rand()/RAND_MAX;
           double value1 = p[i].value (x) * p[j].value(x);
           double value2 = q   .value (x);
           if (std::fabs(value1-value2) > std::max(1e-13,1e-13*std::fabs(value1)))
@@ -139,7 +139,7 @@ void check_expand (const std::vector<Polynomial<double> > &p)
   for (unsigned int i=0; i<p.size(); ++i)
     {
       Polynomial<double> q = p[i];
-      double x = (double)myrand()/RAND_MAX;
+      double x = (double)Testing::rand()/RAND_MAX;
       q += zero;
       double value1 = p[i].value (x);
       double value2 = q   .value (x);
@@ -166,7 +166,7 @@ void check_mult_expand (const std::vector<Polynomial<double> > &p)
           Polynomial<double> q = p[i];
           q += zero;
           q *= p[j];
-          double x = (double)myrand()/RAND_MAX;
+          double x = (double)Testing::rand()/RAND_MAX;
           double value1 = p[i].value (x) * p[j].value(x);
           double value2 = q   .value (x);
           if (std::fabs(value1-value2) > std::max(1e-10,1e-10*std::fabs(value1)))
diff --git a/tests/base/polynomial_lagrange_order.cc b/tests/base/polynomial_lagrange_order.cc
index 9899f14182..0f55292bc5 100644
--- a/tests/base/polynomial_lagrange_order.cc
+++ b/tests/base/polynomial_lagrange_order.cc
@@ -68,7 +68,7 @@ void check_constant (const std::vector<Polynomial<double> > &p)
   deallog << "Representation of one at random points";
   for (unsigned int j=0; j<12; ++j)
     {
-      double x = (double)myrand()/RAND_MAX;
+      double x = (double)Testing::rand()/RAND_MAX;
       double value = 0;
       for (unsigned int i=0; i<p.size(); ++i)
         {
diff --git a/tests/base/polynomial_lagrange_product.cc b/tests/base/polynomial_lagrange_product.cc
index 47d3aced53..d11e59a832 100644
--- a/tests/base/polynomial_lagrange_product.cc
+++ b/tests/base/polynomial_lagrange_product.cc
@@ -46,7 +46,7 @@ void check_derivatives (const std::vector<Polynomial<double> > &p,
   Monomial<double> zero (0,0);
   for (unsigned int j=0; j<p.size(); ++j)
     {
-      double x = (double)myrand()/RAND_MAX;
+      double x = (double)Testing::rand()/RAND_MAX;
       p[j].value (x, values);
       Polynomial<double> q = p[j];
       q += zero;
diff --git a/tests/base/symmetric_tensor_19.cc b/tests/base/symmetric_tensor_19.cc
index 64d7b96a60..fe782feb38 100644
--- a/tests/base/symmetric_tensor_19.cc
+++ b/tests/base/symmetric_tensor_19.cc
@@ -112,7 +112,7 @@ void test ()
         for (unsigned int j=0; j<3; ++j)
           for (unsigned int k=0; k<3; ++k)
             for (unsigned int l=0; l<3; ++l)
-              A[i][j][k][l] = 1.*myrand()/RAND_MAX;
+              A[i][j][k][l] = 1.*Testing::rand()/RAND_MAX;
       break;
 
     default:
diff --git a/tests/base/symmetric_tensor_25.cc b/tests/base/symmetric_tensor_25.cc
index 46026f214f..e9bf54a9ce 100644
--- a/tests/base/symmetric_tensor_25.cc
+++ b/tests/base/symmetric_tensor_25.cc
@@ -33,7 +33,7 @@ void check ()
     {
       SymmetricTensor<2,dim> S;
       for (unsigned int i=0; i<S.n_independent_components; ++i)
-        S[S.unrolled_to_component_indices (i)] = myrand() % 10;
+        S[S.unrolled_to_component_indices (i)] = Testing::rand() % 10;
 
       deallog << "S = " << S << std::endl;
       deallog << "first invariant  = " << first_invariant(S) << std::endl;
diff --git a/tests/base/symmetric_tensor_26.cc b/tests/base/symmetric_tensor_26.cc
index 7a15ffa345..b7ed88f2a9 100644
--- a/tests/base/symmetric_tensor_26.cc
+++ b/tests/base/symmetric_tensor_26.cc
@@ -30,11 +30,11 @@ void check ()
 {
   SymmetricTensor<2,dim> S;
   for (unsigned int i=0; i<S.n_independent_components; ++i)
-    S[S.unrolled_to_component_indices (i)] = myrand() % 10;
+    S[S.unrolled_to_component_indices (i)] = Testing::rand() % 10;
 
   Tensor<1,dim> x;
   for (unsigned int i=0; i<dim; ++i)
-    x[i] = myrand() % 10;
+    x[i] = Testing::rand() % 10;
 
   deallog << "S = " << S << std::endl;
   deallog << "x = " << x << std::endl;
diff --git a/tests/bits/chunk_sparse_matrix_15.cc b/tests/bits/chunk_sparse_matrix_15.cc
index d938f327da..137861362b 100644
--- a/tests/bits/chunk_sparse_matrix_15.cc
+++ b/tests/bits/chunk_sparse_matrix_15.cc
@@ -52,7 +52,7 @@ void test (const unsigned int chunk_size)
       // the iterator and compare it with the exact value
       Vector<double> src(m.n()), dst(m.m()), dst_ref(m.m());
       for (unsigned int i=0; i<src.size(); ++i)
-        src(i) = (double)myrand()/RAND_MAX;
+        src(i) = (double)Testing::rand()/RAND_MAX;
       for (unsigned int i=0; i<m.m(); ++i)
         {
           double sum = 0;
diff --git a/tests/bits/fe_q_constraints.cc b/tests/bits/fe_q_constraints.cc
index 4e8262c5f2..cd0236a045 100644
--- a/tests/bits/fe_q_constraints.cc
+++ b/tests/bits/fe_q_constraints.cc
@@ -78,7 +78,7 @@ TestFunction<dim>::TestFunction (const unsigned int p_order) :
   for (unsigned int d = 0; d < dim; ++d)
     {
       for (unsigned int c = 0; c < p_order; ++c)
-        coeff[c] = (double) myrand() / (double) RAND_MAX;
+        coeff[c] = (double) Testing::rand() / (double) RAND_MAX;
       base.push_back (Polynomials::Polynomial<double> (coeff));
     }
 }
@@ -138,7 +138,7 @@ void TestFEQConstraints<dim>::refine_grid_random ()
   Vector<double> estimated_error_per_cell (n_cells);
 
   for (unsigned int i = 0; i < n_cells; ++i)
-    estimated_error_per_cell(i) = (double) myrand() / (double) RAND_MAX;
+    estimated_error_per_cell(i) = (double) Testing::rand() / (double) RAND_MAX;
 
   GridRefinement::refine_and_coarsen_fixed_number (triangulation,
                                                    estimated_error_per_cell,
diff --git a/tests/bits/find_cell_6.cc b/tests/bits/find_cell_6.cc
index 662dc7a394..95341c21a8 100644
--- a/tests/bits/find_cell_6.cc
+++ b/tests/bits/find_cell_6.cc
@@ -78,7 +78,7 @@ void check1 ()
     {
       for (int j=0; j<1000; ++j)
         {
-          Point<3> p ((myrand()%1000)/1000.0,(rand()%1000)/1000.0,(rand()%1000)/1000.0);
+          Point<3> p ((Testing::rand()%1000)/1000.0,(rand()%1000)/1000.0,(rand()%1000)/1000.0);
           if (!inside(tria, p))
             deallog << "NOT INSIDE" << std::endl;
           GridTools::find_active_cell_around_point(tria, p);
@@ -88,7 +88,7 @@ void check1 ()
            cell != tria.end(); ++cell)
         for (unsigned int f=0; f<GeometryInfo<3>::faces_per_cell; ++f)
           {
-            if (cell->face(f)->at_boundary() && (myrand()%5)==1)
+            if (cell->face(f)->at_boundary() && (Testing::rand()%5)==1)
               cell->set_refine_flag();
           }
 
diff --git a/tests/bits/full_matrix_56.cc b/tests/bits/full_matrix_56.cc
index b5c0546b7c..daf064bfaa 100644
--- a/tests/bits/full_matrix_56.cc
+++ b/tests/bits/full_matrix_56.cc
@@ -45,7 +45,7 @@ void random_matrix(FullMatrix<double> &A)
   for (unsigned int i=0; i<A.m(); ++i)
     for (unsigned int j=0; j<A.n(); ++j)
       {
-        double rnd = myrand();
+        double rnd = Testing::rand();
         rnd /= RAND_MAX;
         A(i,j) = (i==j) ? A.m()+rnd : rnd;
       }
diff --git a/tests/bits/joa_1.cc b/tests/bits/joa_1.cc
index d20bae63b0..8c5b45a5aa 100644
--- a/tests/bits/joa_1.cc
+++ b/tests/bits/joa_1.cc
@@ -1073,8 +1073,8 @@ void LaplaceProblem<dim>::run ()
   // points inside a radius of 0.9)
   for (int i=0; i<1000; i++)
     {
-      double r = sqrt((0.9*myrand()/RAND_MAX));
-      double phi = 2*3.14*(1.0*myrand()/RAND_MAX);
+      double r = sqrt((0.9*Testing::rand()/RAND_MAX));
+      double phi = 2*3.14*(1.0*Testing::rand()/RAND_MAX);
       double x = r*cos(phi);
       double y = r*sin(phi);
       Point<2> p(x,y);
diff --git a/tests/bits/point_value_hp_01.cc b/tests/bits/point_value_hp_01.cc
index 683840a2b3..b99d1564ed 100644
--- a/tests/bits/point_value_hp_01.cc
+++ b/tests/bits/point_value_hp_01.cc
@@ -123,7 +123,7 @@ check ()
   endc = dof_handler.end();
   for (; cell!=endc; ++cell)
     {
-      cell->set_active_fe_index (myrand() % fe.size());
+      cell->set_active_fe_index (Testing::rand() % fe.size());
     }
 
   dof_handler.distribute_dofs (fe);
diff --git a/tests/bits/point_value_hp_02.cc b/tests/bits/point_value_hp_02.cc
index 6d27e90f51..9b7d41d7cd 100644
--- a/tests/bits/point_value_hp_02.cc
+++ b/tests/bits/point_value_hp_02.cc
@@ -132,7 +132,7 @@ check ()
   endc = dof_handler.end();
   for (; cell!=endc; ++cell)
     {
-      cell->set_active_fe_index (myrand() % fe.size());
+      cell->set_active_fe_index (Testing::rand() % fe.size());
     }
 
   dof_handler.distribute_dofs (fe);
diff --git a/tests/bits/rt_2.cc b/tests/bits/rt_2.cc
index 2a41f81990..308b5d3b23 100644
--- a/tests/bits/rt_2.cc
+++ b/tests/bits/rt_2.cc
@@ -158,7 +158,7 @@ int main ()
   // values between 0 and 1.
   Vector<double> solution(dof_handler.n_dofs ());
   for (unsigned int i = 0; i < dof_handler.n_dofs (); ++i)
-    solution(i) = (double) myrand() / (double) RAND_MAX;
+    solution(i) = (double) Testing::rand() / (double) RAND_MAX;
 
   // Now check if the function is
   // continuous in normal direction
diff --git a/tests/bits/step-15.cc b/tests/bits/step-15.cc
index b74fac85e6..31b7279a8f 100644
--- a/tests/bits/step-15.cc
+++ b/tests/bits/step-15.cc
@@ -85,7 +85,7 @@ double InitializationValues::value (const Point<1> &p,
                                     const unsigned int) const
 {
   const double base = std::pow(p(0), 1./3.);
-  const double random = 2.*myrand()/RAND_MAX-1;
+  const double random = 2.*Testing::rand()/RAND_MAX-1;
   return std::max (base+.1*random, 0.);
 }
 
diff --git a/tests/deal.II/coarsening_02.cc b/tests/deal.II/coarsening_02.cc
index 7454364b96..df24077c2b 100644
--- a/tests/deal.II/coarsening_02.cc
+++ b/tests/deal.II/coarsening_02.cc
@@ -85,7 +85,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/deal.II/coarsening_02_1d.cc b/tests/deal.II/coarsening_02_1d.cc
index 7f2df7c629..35774e8473 100644
--- a/tests/deal.II/coarsening_02_1d.cc
+++ b/tests/deal.II/coarsening_02_1d.cc
@@ -85,7 +85,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/deal.II/constraint_graph.cc b/tests/deal.II/constraint_graph.cc
index 3370131cf4..ce14af6b0f 100644
--- a/tests/deal.II/constraint_graph.cc
+++ b/tests/deal.II/constraint_graph.cc
@@ -56,7 +56,7 @@ void test ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof_handler.begin_active();
        cell != dof_handler.end(); ++cell)
-    cell->set_active_fe_index (myrand() % fe_collection.size());
+    cell->set_active_fe_index (Testing::rand() % fe_collection.size());
 
   dof_handler.distribute_dofs(fe_collection);
 
diff --git a/tests/deal.II/constraints_local_to_global.cc b/tests/deal.II/constraints_local_to_global.cc
index 22d806f836..a171d265b8 100644
--- a/tests/deal.II/constraints_local_to_global.cc
+++ b/tests/deal.II/constraints_local_to_global.cc
@@ -90,7 +90,7 @@ void test ()
           if (counter % 42 == 0)
             local_mat(i,j) = 0;
           else
-            local_mat (i,j) = (double)myrand() / RAND_MAX;
+            local_mat (i,j) = (double)Testing::rand() / RAND_MAX;
       cell->get_dof_indices (local_dof_indices);
       constraints.distribute_local_to_global (local_mat, local_dof_indices,
                                               sparse);
diff --git a/tests/deal.II/constraints_local_to_global_chunk.cc b/tests/deal.II/constraints_local_to_global_chunk.cc
index 5d2bfa6b6d..b9ad731dac 100644
--- a/tests/deal.II/constraints_local_to_global_chunk.cc
+++ b/tests/deal.II/constraints_local_to_global_chunk.cc
@@ -94,7 +94,7 @@ void test (unsigned int chunk_size)
           if (counter % 42 == 0)
             local_mat(i,j) = 0;
           else
-            local_mat (i,j) = (double)myrand() / RAND_MAX;
+            local_mat (i,j) = (double)Testing::rand() / RAND_MAX;
       cell->get_dof_indices (local_dof_indices);
       constraints.distribute_local_to_global (local_mat, local_dof_indices,
                                               sparse);
diff --git a/tests/deal.II/dof_handler_number_cache.cc b/tests/deal.II/dof_handler_number_cache.cc
index 5e8ee38ae3..72744cd7ba 100644
--- a/tests/deal.II/dof_handler_number_cache.cc
+++ b/tests/deal.II/dof_handler_number_cache.cc
@@ -57,7 +57,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/deal.II/dof_handler_number_cache_02.cc b/tests/deal.II/dof_handler_number_cache_02.cc
index 463c8bbd85..16cf0da912 100644
--- a/tests/deal.II/dof_handler_number_cache_02.cc
+++ b/tests/deal.II/dof_handler_number_cache_02.cc
@@ -55,7 +55,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/distributed_grids/2d_coarsening_02.cc b/tests/distributed_grids/2d_coarsening_02.cc
index 78adf91af3..9c70b623a4 100644
--- a/tests/distributed_grids/2d_coarsening_02.cc
+++ b/tests/distributed_grids/2d_coarsening_02.cc
@@ -69,7 +69,7 @@ void test(std::ostream & /*out*/)
       // refine one-fifth of cells randomly
       std::vector<bool> flags (tr.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/distributed_grids/2d_coarsening_03.cc b/tests/distributed_grids/2d_coarsening_03.cc
index cccd24abaf..f41f35dccb 100644
--- a/tests/distributed_grids/2d_coarsening_03.cc
+++ b/tests/distributed_grids/2d_coarsening_03.cc
@@ -111,7 +111,7 @@ void test(std::ostream & /*out*/)
       // refine one-fifth of cells randomly
       std::vector<bool> flags (tr.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/distributed_grids/2d_coarsening_05.cc b/tests/distributed_grids/2d_coarsening_05.cc
index 4359953f18..ddc62a6f8e 100644
--- a/tests/distributed_grids/2d_coarsening_05.cc
+++ b/tests/distributed_grids/2d_coarsening_05.cc
@@ -64,7 +64,7 @@ void test(std::ostream & /*out*/)
       // refine one fifth of all cells each
       // time (but at least one)
       for (unsigned int i=0; i<tr.n_active_cells() / 5 + 1; ++i)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
 
       // refine tr and tr2
       unsigned int index=0;
diff --git a/tests/distributed_grids/2d_refinement_02.cc b/tests/distributed_grids/2d_refinement_02.cc
index 82a73c0394..6a59bfe155 100644
--- a/tests/distributed_grids/2d_refinement_02.cc
+++ b/tests/distributed_grids/2d_refinement_02.cc
@@ -71,7 +71,7 @@ void test(std::ostream & /*out*/)
     {
       std::vector<bool> flags (tr.n_active_cells(), false);
       {
-        const unsigned int x = myrand() % flags.size();
+        const unsigned int x = Testing::rand() % flags.size();
         deallog << "Refining cell " << x << std::endl;
         flags[x] = true;
       }
diff --git a/tests/distributed_grids/2d_refinement_03.cc b/tests/distributed_grids/2d_refinement_03.cc
index e9b78a1912..b4079a2621 100644
--- a/tests/distributed_grids/2d_refinement_03.cc
+++ b/tests/distributed_grids/2d_refinement_03.cc
@@ -83,7 +83,7 @@ void test(std::ostream & /*out*/)
     {
       std::vector<bool> flags (tr.n_active_cells());
       for (unsigned int j=0; j<flags.size(); ++j)
-        flags[j] = (myrand() < 0.2*RAND_MAX);
+        flags[j] = (Testing::rand() < 0.2*RAND_MAX);
 
       InterGridMap<Triangulation<dim> > intergrid_map;
       intergrid_map.make_mapping (tr, tr2);
diff --git a/tests/distributed_grids/2d_refinement_05.cc b/tests/distributed_grids/2d_refinement_05.cc
index 83a8ad1f1a..5c1ac4d15a 100644
--- a/tests/distributed_grids/2d_refinement_05.cc
+++ b/tests/distributed_grids/2d_refinement_05.cc
@@ -61,7 +61,7 @@ void test(std::ostream & /*out*/)
       // time (but at least one)
       for (unsigned int i=0; i<tr.n_active_cells() / 5 + 1; ++i)
         {
-          const unsigned int x = myrand() % flags.size();
+          const unsigned int x = Testing::rand() % flags.size();
           deallog << "Refining cell " << x << std::endl;
           flags[x] = true;
         }
diff --git a/tests/distributed_grids/3d_coarsening_02.cc b/tests/distributed_grids/3d_coarsening_02.cc
index a20ab8393a..d0995f5a34 100644
--- a/tests/distributed_grids/3d_coarsening_02.cc
+++ b/tests/distributed_grids/3d_coarsening_02.cc
@@ -69,7 +69,7 @@ void test(std::ostream & /*out*/)
       // refine one-fifth of cells randomly
       std::vector<bool> flags (tr.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/distributed_grids/3d_coarsening_03.cc b/tests/distributed_grids/3d_coarsening_03.cc
index 794e1f3542..c97427bb30 100644
--- a/tests/distributed_grids/3d_coarsening_03.cc
+++ b/tests/distributed_grids/3d_coarsening_03.cc
@@ -66,7 +66,7 @@ void test(std::ostream & /*out*/)
       // refine 10% of cells randomly
       std::vector<bool> flags (tr.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/10 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/distributed_grids/3d_coarsening_04.cc b/tests/distributed_grids/3d_coarsening_04.cc
index e366a3e87a..724f5975d3 100644
--- a/tests/distributed_grids/3d_coarsening_04.cc
+++ b/tests/distributed_grids/3d_coarsening_04.cc
@@ -70,7 +70,7 @@ void test(std::ostream & /*out*/)
       // refine one-tenth of cells randomly
       std::vector<bool> flags (tr.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/30 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/distributed_grids/3d_coarsening_05.cc b/tests/distributed_grids/3d_coarsening_05.cc
index 4b213982cb..bea8903030 100644
--- a/tests/distributed_grids/3d_coarsening_05.cc
+++ b/tests/distributed_grids/3d_coarsening_05.cc
@@ -64,7 +64,7 @@ void test(std::ostream & /*out*/)
       // refine one fifth of all cells each
       // time (but at least one)
       for (unsigned int i=0; i<tr.n_active_cells() / 5 + 1; ++i)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
 
       // refine tr and tr2
       unsigned int index=0;
diff --git a/tests/distributed_grids/3d_refinement_02.cc b/tests/distributed_grids/3d_refinement_02.cc
index 8a50d220c4..b9e101fc60 100644
--- a/tests/distributed_grids/3d_refinement_02.cc
+++ b/tests/distributed_grids/3d_refinement_02.cc
@@ -71,7 +71,7 @@ void test(std::ostream & /*out*/)
     {
       std::vector<bool> flags (tr.n_active_cells(), false);
       {
-        const unsigned int x = myrand() % flags.size();
+        const unsigned int x = Testing::rand() % flags.size();
         deallog << "Refining cell " << x << std::endl;
         flags[x] = true;
       }
diff --git a/tests/distributed_grids/3d_refinement_03.cc b/tests/distributed_grids/3d_refinement_03.cc
index 2b41275641..a38f8d9925 100644
--- a/tests/distributed_grids/3d_refinement_03.cc
+++ b/tests/distributed_grids/3d_refinement_03.cc
@@ -63,7 +63,7 @@ void test(std::ostream & /*out*/)
     {
       std::vector<bool> flags (tr.n_active_cells());
       for (unsigned int j=0; j<flags.size(); ++j)
-        flags[j] = (myrand() < 0.2*RAND_MAX);
+        flags[j] = (Testing::rand() < 0.2*RAND_MAX);
 
       InterGridMap<Triangulation<dim> > intergrid_map;
       intergrid_map.make_mapping (tr, tr2);
diff --git a/tests/distributed_grids/3d_refinement_05.cc b/tests/distributed_grids/3d_refinement_05.cc
index 03b7552abb..c8536fd96d 100644
--- a/tests/distributed_grids/3d_refinement_05.cc
+++ b/tests/distributed_grids/3d_refinement_05.cc
@@ -62,7 +62,7 @@ void test(std::ostream & /*out*/)
       // time (but at least one)
       for (unsigned int i=0; i<tr.n_active_cells() / 5 + 1; ++i)
         {
-          const unsigned int x = myrand() % flags.size();
+          const unsigned int x = Testing::rand() % flags.size();
           deallog << "Refining cell " << x << std::endl;
           flags[x] = true;
         }
diff --git a/tests/distributed_grids/3d_refinement_08.cc b/tests/distributed_grids/3d_refinement_08.cc
index 36ac26f8e2..86d6515fc3 100644
--- a/tests/distributed_grids/3d_refinement_08.cc
+++ b/tests/distributed_grids/3d_refinement_08.cc
@@ -63,7 +63,7 @@ void test(std::ostream & /*out*/)
       deallog << "Refining cells ... " << std::endl;
       for (unsigned int i=0; i<tr.n_active_cells() / 50 + 1; ++i)
         {
-          const unsigned int x = myrand() % flags.size();
+          const unsigned int x = Testing::rand() % flags.size();
           flags[x] = true;
         }
 
diff --git a/tests/distributed_grids/count_dofs_per_block_01.cc b/tests/distributed_grids/count_dofs_per_block_01.cc
index 22c373e0d3..99ae577644 100644
--- a/tests/distributed_grids/count_dofs_per_block_01.cc
+++ b/tests/distributed_grids/count_dofs_per_block_01.cc
@@ -60,7 +60,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/distributed_grids/dof_handler_number_cache.cc b/tests/distributed_grids/dof_handler_number_cache.cc
index 1aabc66112..846b097fea 100644
--- a/tests/distributed_grids/dof_handler_number_cache.cc
+++ b/tests/distributed_grids/dof_handler_number_cache.cc
@@ -61,7 +61,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/distributed_grids/solution_transfer_03.cc b/tests/distributed_grids/solution_transfer_03.cc
index a7f382835b..d40d269cc2 100644
--- a/tests/distributed_grids/solution_transfer_03.cc
+++ b/tests/distributed_grids/solution_transfer_03.cc
@@ -85,7 +85,7 @@ void test(std::ostream & /*out*/)
            cell = tr.begin();
            cell != tr.end(); ++cell)
         {
-          if (cell->has_children() && !(myrand() % 3))
+          if (cell->has_children() && !(Testing::rand() % 3))
             for (unsigned int c=0; c<GeometryInfo<dim>::max_children_per_cell; ++c)
               if (!cell->child(c)->has_children())
                 cell->child(c)->set_coarsen_flag();
@@ -94,7 +94,7 @@ void test(std::ostream & /*out*/)
            cell = tr.begin_active();
            cell != tr.end(); ++cell)
         {
-          if (!(myrand() % 3))
+          if (!(Testing::rand() % 3))
             {
               cell->clear_coarsen_flag();
               cell->set_refine_flag();
diff --git a/tests/fail/abf_approximation_01.cc b/tests/fail/abf_approximation_01.cc
index d98354cfcb..6419a1bc72 100644
--- a/tests/fail/abf_approximation_01.cc
+++ b/tests/fail/abf_approximation_01.cc
@@ -293,7 +293,7 @@ public:
       {
         for (unsigned int i=0; i < deg; ++i)
           {
-            double c = ((double) myrand()) / ((double) RAND_MAX + 1);
+            double c = ((double) Testing::rand()) / ((double) RAND_MAX + 1);
             coeff[i] = c;
           }
         polys.push_back (Polynomials::Polynomial<double> (coeff));
@@ -518,7 +518,7 @@ int main ()
     // Try arbitrary deformation ...
     for (unsigned int i=0; i < deformation.size (); ++i)
       {
-        double c = ((double) myrand()) / ((double) RAND_MAX + 1);
+        double c = ((double) Testing::rand()) / ((double) RAND_MAX + 1);
         deformation (i) = 0.35 * (c - 0.5);
       }
 
diff --git a/tests/fail/hp-step-15.cc b/tests/fail/hp-step-15.cc
index 6d050878a8..a76e365331 100644
--- a/tests/fail/hp-step-15.cc
+++ b/tests/fail/hp-step-15.cc
@@ -86,7 +86,7 @@ double InitializationValues::value (const Point<1> &p,
                                     const unsigned int) const
 {
   const double base = std::pow(p(0), 1./3.);
-  const double random = 2.*myrand()/RAND_MAX-1;
+  const double random = 2.*Testing::rand()/RAND_MAX-1;
   return std::max (base+.1*random, 0.);
 }
 
diff --git a/tests/fail/vectors_boundary_rhs_hp_02.cc b/tests/fail/vectors_boundary_rhs_hp_02.cc
index 590e3b7197..8dc62a8b24 100644
--- a/tests/fail/vectors_boundary_rhs_hp_02.cc
+++ b/tests/fail/vectors_boundary_rhs_hp_02.cc
@@ -91,7 +91,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/fail/vectors_rhs_hp_02.cc b/tests/fail/vectors_rhs_hp_02.cc
index 577b2870ff..71712cddc5 100644
--- a/tests/fail/vectors_rhs_hp_02.cc
+++ b/tests/fail/vectors_rhs_hp_02.cc
@@ -91,7 +91,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/fe/rt_8.cc b/tests/fe/rt_8.cc
index a4526be072..fa2c573579 100644
--- a/tests/fe/rt_8.cc
+++ b/tests/fe/rt_8.cc
@@ -87,7 +87,7 @@ test (const unsigned int degree)
 
   Vector<double> tmp1(dofs_per_cell), tmp2(dofs_per_cell);
   for (unsigned int i=0; i<dofs_per_cell; ++i)
-    tmp1(i) = 1.*myrand()/RAND_MAX;
+    tmp1(i) = 1.*Testing::rand()/RAND_MAX;
   cg.solve (mass_matrix, tmp2, tmp1, PreconditionIdentity());
 
   deallog << "Degree=" << degree
diff --git a/tests/fe/rt_9.cc b/tests/fe/rt_9.cc
index 5f2e30610d..28cc4c9ea6 100644
--- a/tests/fe/rt_9.cc
+++ b/tests/fe/rt_9.cc
@@ -78,7 +78,7 @@ test (const unsigned int degree)
 
   Vector<double> tmp1(dofs_per_cell), tmp2(dofs_per_cell);
   for (unsigned int i=0; i<dofs_per_cell; ++i)
-    tmp1(i) = 1.*myrand()/RAND_MAX;
+    tmp1(i) = 1.*Testing::rand()/RAND_MAX;
   cg.solve (mass_matrix, tmp2, tmp1, PreconditionIdentity());
 
   deallog << "Degree=" << degree
diff --git a/tests/fe/rt_approximation_01.cc b/tests/fe/rt_approximation_01.cc
index 9218074434..eb4a6af787 100644
--- a/tests/fe/rt_approximation_01.cc
+++ b/tests/fe/rt_approximation_01.cc
@@ -293,7 +293,7 @@ public:
       {
         for (unsigned int i=0; i < deg; ++i)
           {
-            double c = ((double) myrand()) / ((double) RAND_MAX + 1);
+            double c = ((double) Testing::rand()) / ((double) RAND_MAX + 1);
             coeff[i] = c;
           }
         polys.push_back (Polynomials::Polynomial<double> (coeff));
@@ -523,7 +523,7 @@ int main ()
     // Try arbitrary deformation ...
     for (unsigned int i=0; i < deformation.size (); ++i)
       {
-        double c = ((double) myrand()) / ((double) RAND_MAX + 1);
+        double c = ((double) Testing::rand()) / ((double) RAND_MAX + 1);
         deformation (i) = 0.35 * (c - 0.5);
       }
 
diff --git a/tests/fe/rt_normal_02.cc b/tests/fe/rt_normal_02.cc
index 1b6cc0c58a..5b18693c01 100644
--- a/tests/fe/rt_normal_02.cc
+++ b/tests/fe/rt_normal_02.cc
@@ -291,7 +291,7 @@ int main (int /*argc*/, char **/*argv*/)
 
   // Fill solution vector with random values between 0 and 1.
   for (unsigned int i = 0; i < dof_handler->n_dofs (); ++i)
-    solution(i) = (double) myrand() / (double) RAND_MAX;
+    solution(i) = (double) Testing::rand() / (double) RAND_MAX;
 
   // Now check if the function is continuous in normal
   // direction.
diff --git a/tests/hp/count_dofs_per_block_01.cc b/tests/hp/count_dofs_per_block_01.cc
index 125f6de7ff..e384875032 100644
--- a/tests/hp/count_dofs_per_block_01.cc
+++ b/tests/hp/count_dofs_per_block_01.cc
@@ -82,7 +82,7 @@ check ()
   cell = dof_handler.begin_active(),
   endc = dof_handler.end();
   for (; cell!=endc; ++cell)
-    cell->set_active_fe_index (myrand() % fe.size());
+    cell->set_active_fe_index (Testing::rand() % fe.size());
 
   dof_handler.distribute_dofs (fe);
 
diff --git a/tests/hp/count_dofs_per_block_02.cc b/tests/hp/count_dofs_per_block_02.cc
index a5adcfa998..1fbbdb4b28 100644
--- a/tests/hp/count_dofs_per_block_02.cc
+++ b/tests/hp/count_dofs_per_block_02.cc
@@ -82,7 +82,7 @@ check ()
   cell = dof_handler.begin_active(),
   endc = dof_handler.end();
   for (; cell!=endc; ++cell)
-    cell->set_active_fe_index (myrand() % fe.size());
+    cell->set_active_fe_index (Testing::rand() % fe.size());
 
   dof_handler.distribute_dofs (fe);
 
diff --git a/tests/hp/crash_06.cc b/tests/hp/crash_06.cc
index b688e5db7f..21f25fcd99 100644
--- a/tests/hp/crash_06.cc
+++ b/tests/hp/crash_06.cc
@@ -70,7 +70,7 @@ void do_check (const Triangulation<dim> &triangulation,
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof_handler.begin_active();
        cell != dof_handler.end(); ++cell)
-    cell->set_active_fe_index (myrand() % fe.size());
+    cell->set_active_fe_index (Testing::rand() % fe.size());
   dof_handler.distribute_dofs (fe);
 
   deallog << "n_dofs=" << dof_handler.n_dofs() << std::endl;
diff --git a/tests/hp/crash_10.cc b/tests/hp/crash_10.cc
index a868290e1c..be75be61a2 100644
--- a/tests/hp/crash_10.cc
+++ b/tests/hp/crash_10.cc
@@ -71,7 +71,7 @@ void test ()
   typename hp::DoFHandler<dim>::active_cell_iterator cell = dof_handler.begin_active (),
                                                      endc = dof_handler.end ();
   for (; cell != endc; ++cell)
-    cell->set_active_fe_index (myrand() % fe.size());
+    cell->set_active_fe_index (Testing::rand() % fe.size());
 
   dof_handler.distribute_dofs (fe);
 
@@ -113,7 +113,7 @@ int main ()
   // same value it has in
   // hp_hanging_nodes_02
   for (unsigned int i=0; i<64; ++i)
-    myrand();
+    Testing::rand();
   test<3> ();
 
   return 0;
diff --git a/tests/hp/crash_16.cc b/tests/hp/crash_16.cc
index 230627f4fe..e97612c11e 100644
--- a/tests/hp/crash_16.cc
+++ b/tests/hp/crash_16.cc
@@ -68,7 +68,7 @@ void test ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof_handler.begin_active();
        cell != dof_handler.end(); ++cell)
-    cell->set_active_fe_index (myrand() % fe.size());
+    cell->set_active_fe_index (Testing::rand() % fe.size());
   dof_handler.distribute_dofs (fe);
 
   // loop over all lines and make sure that
diff --git a/tests/hp/create_laplace_matrix_01.cc b/tests/hp/create_laplace_matrix_01.cc
index 9302e2e068..1bb40cd3b5 100644
--- a/tests/hp/create_laplace_matrix_01.cc
+++ b/tests/hp/create_laplace_matrix_01.cc
@@ -80,7 +80,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_laplace_matrix_01b.cc b/tests/hp/create_laplace_matrix_01b.cc
index 362d2fa730..d44f0d5ba0 100644
--- a/tests/hp/create_laplace_matrix_01b.cc
+++ b/tests/hp/create_laplace_matrix_01b.cc
@@ -81,7 +81,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_laplace_matrix_02.cc b/tests/hp/create_laplace_matrix_02.cc
index b4bc3c017b..58f755f8db 100644
--- a/tests/hp/create_laplace_matrix_02.cc
+++ b/tests/hp/create_laplace_matrix_02.cc
@@ -82,7 +82,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_laplace_matrix_02b.cc b/tests/hp/create_laplace_matrix_02b.cc
index dfcf4b1bc9..753c66f9ab 100644
--- a/tests/hp/create_laplace_matrix_02b.cc
+++ b/tests/hp/create_laplace_matrix_02b.cc
@@ -83,7 +83,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_laplace_matrix_03.cc b/tests/hp/create_laplace_matrix_03.cc
index 485f897043..d47c37e143 100644
--- a/tests/hp/create_laplace_matrix_03.cc
+++ b/tests/hp/create_laplace_matrix_03.cc
@@ -79,7 +79,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_laplace_matrix_03b.cc b/tests/hp/create_laplace_matrix_03b.cc
index 45701b331a..b337cf4b63 100644
--- a/tests/hp/create_laplace_matrix_03b.cc
+++ b/tests/hp/create_laplace_matrix_03b.cc
@@ -80,7 +80,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_laplace_matrix_04.cc b/tests/hp/create_laplace_matrix_04.cc
index bd28208994..c9c6ba4e00 100644
--- a/tests/hp/create_laplace_matrix_04.cc
+++ b/tests/hp/create_laplace_matrix_04.cc
@@ -81,7 +81,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_laplace_matrix_04b.cc b/tests/hp/create_laplace_matrix_04b.cc
index 7ffe72beaf..ff851a29ae 100644
--- a/tests/hp/create_laplace_matrix_04b.cc
+++ b/tests/hp/create_laplace_matrix_04b.cc
@@ -82,7 +82,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_laplace_matrix_constraints_02.cc b/tests/hp/create_laplace_matrix_constraints_02.cc
index 99b4c53635..a1922a1656 100644
--- a/tests/hp/create_laplace_matrix_constraints_02.cc
+++ b/tests/hp/create_laplace_matrix_constraints_02.cc
@@ -77,7 +77,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_laplace_matrix_constraints_02b.cc b/tests/hp/create_laplace_matrix_constraints_02b.cc
index 5c9a084566..156517288a 100644
--- a/tests/hp/create_laplace_matrix_constraints_02b.cc
+++ b/tests/hp/create_laplace_matrix_constraints_02b.cc
@@ -77,7 +77,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_01.cc b/tests/hp/create_mass_matrix_01.cc
index 0abe928230..4ce66bdcdd 100644
--- a/tests/hp/create_mass_matrix_01.cc
+++ b/tests/hp/create_mass_matrix_01.cc
@@ -80,7 +80,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_01b.cc b/tests/hp/create_mass_matrix_01b.cc
index 0a574b3001..5abac27879 100644
--- a/tests/hp/create_mass_matrix_01b.cc
+++ b/tests/hp/create_mass_matrix_01b.cc
@@ -81,7 +81,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_02.cc b/tests/hp/create_mass_matrix_02.cc
index 8ea98aa19d..e1f6550afd 100644
--- a/tests/hp/create_mass_matrix_02.cc
+++ b/tests/hp/create_mass_matrix_02.cc
@@ -82,7 +82,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_02b.cc b/tests/hp/create_mass_matrix_02b.cc
index 32102a9a76..c3372109f2 100644
--- a/tests/hp/create_mass_matrix_02b.cc
+++ b/tests/hp/create_mass_matrix_02b.cc
@@ -83,7 +83,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_03.cc b/tests/hp/create_mass_matrix_03.cc
index 1efcbfcb9e..8825edb228 100644
--- a/tests/hp/create_mass_matrix_03.cc
+++ b/tests/hp/create_mass_matrix_03.cc
@@ -79,7 +79,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_03b.cc b/tests/hp/create_mass_matrix_03b.cc
index c53a90d169..f7e86e8d39 100644
--- a/tests/hp/create_mass_matrix_03b.cc
+++ b/tests/hp/create_mass_matrix_03b.cc
@@ -82,7 +82,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_04.cc b/tests/hp/create_mass_matrix_04.cc
index 5b1016a6ac..4a018f3320 100644
--- a/tests/hp/create_mass_matrix_04.cc
+++ b/tests/hp/create_mass_matrix_04.cc
@@ -81,7 +81,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_04b.cc b/tests/hp/create_mass_matrix_04b.cc
index deb080e399..78f6d522fb 100644
--- a/tests/hp/create_mass_matrix_04b.cc
+++ b/tests/hp/create_mass_matrix_04b.cc
@@ -82,7 +82,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_05.cc b/tests/hp/create_mass_matrix_05.cc
index baea9e0be3..81e37d6de5 100644
--- a/tests/hp/create_mass_matrix_05.cc
+++ b/tests/hp/create_mass_matrix_05.cc
@@ -75,7 +75,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_constraints_02.cc b/tests/hp/create_mass_matrix_constraints_02.cc
index 8df741b047..6d3280ead3 100644
--- a/tests/hp/create_mass_matrix_constraints_02.cc
+++ b/tests/hp/create_mass_matrix_constraints_02.cc
@@ -77,7 +77,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/create_mass_matrix_constraints_02b.cc b/tests/hp/create_mass_matrix_constraints_02b.cc
index d073d58b4e..55e5fcb46f 100644
--- a/tests/hp/create_mass_matrix_constraints_02b.cc
+++ b/tests/hp/create_mass_matrix_constraints_02b.cc
@@ -77,7 +77,7 @@ check ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active();
        cell != dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/dof_handler_number_cache.cc b/tests/hp/dof_handler_number_cache.cc
index f1c2fe5df3..21ae525678 100644
--- a/tests/hp/dof_handler_number_cache.cc
+++ b/tests/hp/dof_handler_number_cache.cc
@@ -60,7 +60,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/hp/get_active_fe_indices.cc b/tests/hp/get_active_fe_indices.cc
index 262c077619..ae2090693e 100644
--- a/tests/hp/get_active_fe_indices.cc
+++ b/tests/hp/get_active_fe_indices.cc
@@ -57,7 +57,7 @@ void test ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof_handler.begin_active();
        cell != dof_handler.end(); ++cell)
-    cell->set_active_fe_index (myrand() % fe_collection.size());
+    cell->set_active_fe_index (Testing::rand() % fe_collection.size());
 
   dof_handler.distribute_dofs(fe_collection);
 
diff --git a/tests/hp/hp_constraints_common.h b/tests/hp/hp_constraints_common.h
index 6b580b9adb..818ca03e53 100644
--- a/tests/hp/hp_constraints_common.h
+++ b/tests/hp/hp_constraints_common.h
@@ -67,7 +67,7 @@ void do_check (const Triangulation<dim> &triangulation,
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof_handler.begin_active();
        cell != dof_handler.end(); ++cell)
-    cell->set_active_fe_index (rand() % fe.size());
+    cell->set_active_fe_index (Testing::rand() % fe.size());
   dof_handler.distribute_dofs (fe);
 
   deallog << "n_dofs=" << dof_handler.n_dofs() << std::endl;
@@ -125,7 +125,7 @@ void test_with_hanging_nodes_random (const hp::FECollection<dim> &fe)
            cell = triangulation.begin_active();
            cell != triangulation.end();
            ++cell)
-        if (rand() % 4 == 0)
+        if (Testing::rand() % 4 == 0)
           cell->set_refine_flag ();
       triangulation.execute_coarsening_and_refinement ();
     }
@@ -148,8 +148,8 @@ void test_with_hanging_nodes_random_aniso (const hp::FECollection<dim> &fe)
            cell = triangulation.begin_active();
            cell != triangulation.end();
            ++cell)
-        if (rand() % 4 == 0)
-          cell->set_refine_flag (RefinementCase<dim>(rand() % RefinementCase<dim>::isotropic_refinement + 1));
+        if (Testing::rand() % 4 == 0)
+          cell->set_refine_flag (RefinementCase<dim>(Testing::rand() % RefinementCase<dim>::isotropic_refinement + 1));
       triangulation.execute_coarsening_and_refinement ();
     }
 
diff --git a/tests/hp/hp_hanging_nodes_02.cc b/tests/hp/hp_hanging_nodes_02.cc
index 959e727aab..f94dd9d9f2 100644
--- a/tests/hp/hp_hanging_nodes_02.cc
+++ b/tests/hp/hp_hanging_nodes_02.cc
@@ -92,7 +92,7 @@ void run (bool random_p,
     {
       for (; cell != endc; ++cell)
         {
-          cell->set_active_fe_index ((int)(4.0 * (double) myrand() / (double) RAND_MAX));
+          cell->set_active_fe_index ((int)(4.0 * (double) Testing::rand() / (double) RAND_MAX));
         }
     }
   else
diff --git a/tests/hp/integrate_difference.cc b/tests/hp/integrate_difference.cc
index 91235de444..e0e8ae9120 100644
--- a/tests/hp/integrate_difference.cc
+++ b/tests/hp/integrate_difference.cc
@@ -67,7 +67,7 @@ void test ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof_handler.begin_active();
        cell != dof_handler.end(); ++cell)
-    cell->set_active_fe_index (myrand() % fe_collection.size());
+    cell->set_active_fe_index (Testing::rand() % fe_collection.size());
 
   dof_handler.distribute_dofs(fe_collection);
 
diff --git a/tests/hp/integrate_difference_02.cc b/tests/hp/integrate_difference_02.cc
index e0ebb1c5af..87770b6bd5 100644
--- a/tests/hp/integrate_difference_02.cc
+++ b/tests/hp/integrate_difference_02.cc
@@ -67,7 +67,7 @@ void test ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof_handler.begin_active();
        cell != dof_handler.end(); ++cell)
-    cell->set_active_fe_index (myrand() % fe_collection.size());
+    cell->set_active_fe_index (Testing::rand() % fe_collection.size());
 
   dof_handler.distribute_dofs(fe_collection);
 
diff --git a/tests/hp/matrices_hp.cc b/tests/hp/matrices_hp.cc
index 315dd67040..81eeae01e1 100644
--- a/tests/hp/matrices_hp.cc
+++ b/tests/hp/matrices_hp.cc
@@ -150,7 +150,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/n_active_fe_indices.cc b/tests/hp/n_active_fe_indices.cc
index 28f54d287c..b8d536fe8a 100644
--- a/tests/hp/n_active_fe_indices.cc
+++ b/tests/hp/n_active_fe_indices.cc
@@ -136,7 +136,7 @@ void test ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof_handler.begin_active();
        cell != dof_handler.end(); ++cell)
-    cell->set_active_fe_index (myrand() % fe_collection.size());
+    cell->set_active_fe_index (Testing::rand() % fe_collection.size());
 
   dof_handler.distribute_dofs(fe_collection);
 
diff --git a/tests/hp/random.cc b/tests/hp/random.cc
index e747cd5fbe..d3f3a3b2d8 100644
--- a/tests/hp/random.cc
+++ b/tests/hp/random.cc
@@ -56,7 +56,7 @@ void test ()
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof_handler.begin_active();
        cell != dof_handler.end(); ++cell)
-    cell->set_active_fe_index (myrand() % fe_collection.size());
+    cell->set_active_fe_index (Testing::rand() % fe_collection.size());
 
   dof_handler.distribute_dofs(fe_collection);
 
diff --git a/tests/hp/solution_transfer.cc b/tests/hp/solution_transfer.cc
index 64b94d10f3..ea0c2e8205 100644
--- a/tests/hp/solution_transfer.cc
+++ b/tests/hp/solution_transfer.cc
@@ -105,11 +105,11 @@ void transfer(std::ostream &out)
         if (counter < 15)
           cell->set_active_fe_index(1);
         else
-          cell->set_active_fe_index(myrand()%max_degree);
+          cell->set_active_fe_index(Testing::rand()%max_degree);
         if (counter < 15)
           celldg->set_active_fe_index(1);
         else
-          celldg->set_active_fe_index(myrand()%max_degree);
+          celldg->set_active_fe_index(Testing::rand()%max_degree);
       }
   }
 
@@ -140,9 +140,9 @@ void transfer(std::ostream &out)
     {
       if (counter > 120)
         cell->set_coarsen_flag();
-      else if (myrand() % 3 == 0)
+      else if (Testing::rand() % 3 == 0)
         cell->set_refine_flag();
-      else if (myrand() % 3 == 3)
+      else if (Testing::rand() % 3 == 3)
         cell->set_coarsen_flag();
     }
 
@@ -164,11 +164,11 @@ void transfer(std::ostream &out)
         if (counter > 20 && counter < 90)
           cell->set_active_fe_index(0);
         else
-          cell->set_active_fe_index(myrand()%max_degree);
+          cell->set_active_fe_index(Testing::rand()%max_degree);
         if (counter > 20 && counter < 90)
           celldg->set_active_fe_index(0);
         else
-          celldg->set_active_fe_index(myrand()%max_degree);
+          celldg->set_active_fe_index(Testing::rand()%max_degree);
       }
   }
 
diff --git a/tests/hp/step-3b.cc b/tests/hp/step-3b.cc
index 68a219c55b..9881bfac1b 100644
--- a/tests/hp/step-3b.cc
+++ b/tests/hp/step-3b.cc
@@ -110,7 +110,7 @@ void LaplaceProblem::make_grid_and_dofs ()
   unsigned int cell_no = 0;
   for (; cell != endc; ++cell)
     {
-      cell->set_active_fe_index (myrand() % fe.size());
+      cell->set_active_fe_index (Testing::rand() % fe.size());
       deallog << "Cell " << cell << " has fe_index=" << cell->active_fe_index()
               << std::endl;
       ++cell_no;
diff --git a/tests/hp/step-3c.cc b/tests/hp/step-3c.cc
index 7161040315..162034bc7c 100644
--- a/tests/hp/step-3c.cc
+++ b/tests/hp/step-3c.cc
@@ -109,7 +109,7 @@ void LaplaceProblem::make_grid_and_dofs ()
   unsigned int cell_no = 0;
   for (; cell != endc; ++cell)
     {
-      cell->set_active_fe_index (myrand() % fe.size());
+      cell->set_active_fe_index (Testing::rand() % fe.size());
       deallog << "Cell " << cell << " has fe_index=" << cell->active_fe_index()
               << std::endl;
       ++cell_no;
diff --git a/tests/hp/vectors_boundary_rhs_01.cc b/tests/hp/vectors_boundary_rhs_01.cc
index 8476ea4afd..117e29889b 100644
--- a/tests/hp/vectors_boundary_rhs_01.cc
+++ b/tests/hp/vectors_boundary_rhs_01.cc
@@ -91,7 +91,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_boundary_rhs_02.cc b/tests/hp/vectors_boundary_rhs_02.cc
index d7be80fb7d..bbfb624c29 100644
--- a/tests/hp/vectors_boundary_rhs_02.cc
+++ b/tests/hp/vectors_boundary_rhs_02.cc
@@ -90,7 +90,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_boundary_rhs_03.cc b/tests/hp/vectors_boundary_rhs_03.cc
index 79ccf747f1..d7fbfad646 100644
--- a/tests/hp/vectors_boundary_rhs_03.cc
+++ b/tests/hp/vectors_boundary_rhs_03.cc
@@ -89,7 +89,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_boundary_rhs_hp_01.cc b/tests/hp/vectors_boundary_rhs_hp_01.cc
index 131acc3168..a97062db21 100644
--- a/tests/hp/vectors_boundary_rhs_hp_01.cc
+++ b/tests/hp/vectors_boundary_rhs_hp_01.cc
@@ -92,7 +92,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_boundary_rhs_hp_03.cc b/tests/hp/vectors_boundary_rhs_hp_03.cc
index 336e53f4a5..ec379d42f4 100644
--- a/tests/hp/vectors_boundary_rhs_hp_03.cc
+++ b/tests/hp/vectors_boundary_rhs_hp_03.cc
@@ -90,7 +90,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_point_source_01.cc b/tests/hp/vectors_point_source_01.cc
index 85005a1a87..abfe4aeae0 100644
--- a/tests/hp/vectors_point_source_01.cc
+++ b/tests/hp/vectors_point_source_01.cc
@@ -68,7 +68,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_point_source_hp_01.cc b/tests/hp/vectors_point_source_hp_01.cc
index 29720607b1..1b45bc34dc 100644
--- a/tests/hp/vectors_point_source_hp_01.cc
+++ b/tests/hp/vectors_point_source_hp_01.cc
@@ -68,7 +68,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_rhs_01.cc b/tests/hp/vectors_rhs_01.cc
index 398549e9fe..7cb5704ecb 100644
--- a/tests/hp/vectors_rhs_01.cc
+++ b/tests/hp/vectors_rhs_01.cc
@@ -91,7 +91,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_rhs_02.cc b/tests/hp/vectors_rhs_02.cc
index e703fac95a..2f5d15b73e 100644
--- a/tests/hp/vectors_rhs_02.cc
+++ b/tests/hp/vectors_rhs_02.cc
@@ -90,7 +90,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_rhs_03.cc b/tests/hp/vectors_rhs_03.cc
index a27d07c444..52b21c16db 100644
--- a/tests/hp/vectors_rhs_03.cc
+++ b/tests/hp/vectors_rhs_03.cc
@@ -89,7 +89,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_rhs_hp_01.cc b/tests/hp/vectors_rhs_hp_01.cc
index 3a4f215c83..2dee753110 100644
--- a/tests/hp/vectors_rhs_hp_01.cc
+++ b/tests/hp/vectors_rhs_hp_01.cc
@@ -92,7 +92,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/hp/vectors_rhs_hp_03.cc b/tests/hp/vectors_rhs_hp_03.cc
index 40d54f393b..31dab39f12 100644
--- a/tests/hp/vectors_rhs_hp_03.cc
+++ b/tests/hp/vectors_rhs_hp_03.cc
@@ -90,7 +90,7 @@ check ()
   hp::DoFHandler<dim> dof(tr);
   for (typename hp::DoFHandler<dim>::active_cell_iterator
        cell = dof.begin_active(); cell!=dof.end(); ++cell)
-    cell->set_active_fe_index (myrand() % element.size());
+    cell->set_active_fe_index (Testing::rand() % element.size());
 
   dof.distribute_dofs(element);
 
diff --git a/tests/lac/full_matrix.cc b/tests/lac/full_matrix.cc
index 0876f633d7..507d7deb65 100644
--- a/tests/lac/full_matrix.cc
+++ b/tests/lac/full_matrix.cc
@@ -37,7 +37,7 @@ void random_matrix(FullMatrix<double> &A)
   for (unsigned int i=0; i<A.m(); ++i)
     for (unsigned int j=0; j<A.n(); ++j)
       {
-        double rnd = myrand();
+        double rnd = Testing::rand();
         rnd /= RAND_MAX;
         A(i,j) = (i==j) ? A.m()+rnd : rnd;
       }
diff --git a/tests/lac/full_matrix_04.cc b/tests/lac/full_matrix_04.cc
index 9a7e1d269b..9de8c6159b 100644
--- a/tests/lac/full_matrix_04.cc
+++ b/tests/lac/full_matrix_04.cc
@@ -52,16 +52,16 @@ void test (const unsigned int n, const unsigned int m)
   for (unsigned int i=0; i<n; ++i)
     {
       for (unsigned int j=0; j<n; ++j)
-        A(i,j) = myrand();
+        A(i,j) = Testing::rand();
       for (unsigned int j=0; j<m; ++j)
         {
-          B(i,j) = Bt(j,i) = myrand();
-          D(j,i) = Dt(i,j) = myrand();
+          B(i,j) = Bt(j,i) = Testing::rand();
+          D(j,i) = Dt(i,j) = Testing::rand();
         }
     }
   for (unsigned int i=0; i<m; ++i)
     for (unsigned int j=0; j<m; ++j)
-      C(i,j) = myrand();
+      C(i,j) = Testing::rand();
 
   // Compare first Schur complement
   // with mmult.
diff --git a/tests/lac/gmres_reorthogonalize_01.cc b/tests/lac/gmres_reorthogonalize_01.cc
index 49fe200b8a..6968b5d5fa 100644
--- a/tests/lac/gmres_reorthogonalize_01.cc
+++ b/tests/lac/gmres_reorthogonalize_01.cc
@@ -36,7 +36,7 @@ void test (unsigned int variant)
   FullMatrix<number> matrix(n, n);
   for (unsigned int i=0; i<n; ++i)
     for (unsigned int j=0; j<n; ++j)
-      matrix(i,j) = -0.1 + 0.2 * myrand()/RAND_MAX;
+      matrix(i,j) = -0.1 + 0.2 * Testing::rand()/RAND_MAX;
 
   // put diagonal entries of different strengths. these are very challenging
   // for GMRES and will usually take a lot of iterations until the Krylov
diff --git a/tests/lac/gmres_reorthogonalize_03.cc b/tests/lac/gmres_reorthogonalize_03.cc
index b6830f91a4..26dd1b61e5 100644
--- a/tests/lac/gmres_reorthogonalize_03.cc
+++ b/tests/lac/gmres_reorthogonalize_03.cc
@@ -35,7 +35,7 @@ void test (unsigned int variant)
   FullMatrix<number> matrix(n, n);
   for (unsigned int i=0; i<n; ++i)
     for (unsigned int j=0; j<n; ++j)
-      matrix(i,j) = -0.1 + 0.2 * myrand()/RAND_MAX;
+      matrix(i,j) = -0.1 + 0.2 * Testing::rand()/RAND_MAX;
 
   // put diagonal entries of different strengths. these are very challenging
   // for GMRES and will usually take a lot of iterations until the Krylov
diff --git a/tests/lac/precondition_chebyshev_01.cc b/tests/lac/precondition_chebyshev_01.cc
index 23982d6cd8..bb4b86a5b9 100644
--- a/tests/lac/precondition_chebyshev_01.cc
+++ b/tests/lac/precondition_chebyshev_01.cc
@@ -55,7 +55,7 @@ check()
 
   Vector<double> in (size), out(size);
   for (unsigned int i=0; i<size; ++i)
-    in(i) = (double)myrand()/RAND_MAX;
+    in(i) = (double)Testing::rand()/RAND_MAX;
 
   PreconditionChebyshev<FullMatrixModified,Vector<double> > prec;
   PreconditionChebyshev<FullMatrixModified,Vector<double> >::AdditionalData
diff --git a/tests/lac/sparse_ilu.cc b/tests/lac/sparse_ilu.cc
index f91a5dddbb..126c8b5a4a 100644
--- a/tests/lac/sparse_ilu.cc
+++ b/tests/lac/sparse_ilu.cc
@@ -103,7 +103,7 @@ int main()
           for (unsigned int i=0; i<3; ++i)
             {
               for (unsigned int j=0; j<dim; ++j)
-                v(j) = 1. * myrand()/RAND_MAX;
+                v(j) = 1. * Testing::rand()/RAND_MAX;
 
               A.vmult (tmp1, v);
               ilu.vmult (tmp2, tmp1);
diff --git a/tests/lac/sparse_ilu_t.cc b/tests/lac/sparse_ilu_t.cc
index 8044a18a96..453f71744c 100644
--- a/tests/lac/sparse_ilu_t.cc
+++ b/tests/lac/sparse_ilu_t.cc
@@ -104,7 +104,7 @@ int main()
           for (unsigned int i=0; i<3; ++i)
             {
               for (unsigned int j=0; j<dim; ++j)
-                v(j) = 1. * myrand()/RAND_MAX;
+                v(j) = 1. * Testing::rand()/RAND_MAX;
 
               A.Tvmult (tmp1, v);
               ilu.Tvmult (tmp2, tmp1);
diff --git a/tests/lac/sparse_matrix_Tmmult_01.cc b/tests/lac/sparse_matrix_Tmmult_01.cc
index c9db6880d4..823e1cc309 100644
--- a/tests/lac/sparse_matrix_Tmmult_01.cc
+++ b/tests/lac/sparse_matrix_Tmmult_01.cc
@@ -56,9 +56,9 @@ void test (const unsigned int n)
   for (unsigned int i=0; i<n; ++i)
     {
       for (unsigned int j=0; j<n; ++j)
-        A.set(i,j,myrand());
+        A.set(i,j,Testing::rand());
       for (unsigned int j=0; j<n; ++j)
-        B.set(i,j,myrand());
+        B.set(i,j,Testing::rand());
     }
 
   // now form the matrix-matrix product and
@@ -67,7 +67,7 @@ void test (const unsigned int n)
 
   Vector<double> x(n), y(n), z(n), tmp(n);
   for (unsigned int j=0; j<n; ++j)
-    x(j) = myrand();
+    x(j) = Testing::rand();
 
   // then test for correctness
   C.vmult (y, x);
diff --git a/tests/lac/sparse_matrix_Tmmult_02.cc b/tests/lac/sparse_matrix_Tmmult_02.cc
index 83aff6f4e8..c639bff45f 100644
--- a/tests/lac/sparse_matrix_Tmmult_02.cc
+++ b/tests/lac/sparse_matrix_Tmmult_02.cc
@@ -56,14 +56,14 @@ void test (const unsigned int n)
   for (unsigned int i=0; i<n; ++i)
     {
       for (unsigned int j=0; j<n; ++j)
-        A.set(i,j,myrand());
+        A.set(i,j,Testing::rand());
       for (unsigned int j=0; j<n; ++j)
-        B.set(i,j,myrand());
+        B.set(i,j,Testing::rand());
     }
 
   Vector<double> v(n);
   for (unsigned int j=0; j<n; ++j)
-    v(j) = myrand();
+    v(j) = Testing::rand();
 
   // now form the matrix-matrix product and
   // initialize a test rhs
@@ -71,7 +71,7 @@ void test (const unsigned int n)
 
   Vector<double> x(n), y(n), z(n), tmp(n);
   for (unsigned int j=0; j<n; ++j)
-    x(j) = myrand();
+    x(j) = Testing::rand();
 
   // then test for correctness
   C.vmult (y, x);
diff --git a/tests/lac/sparse_matrix_mmult_01.cc b/tests/lac/sparse_matrix_mmult_01.cc
index de2070609e..20ac941be9 100644
--- a/tests/lac/sparse_matrix_mmult_01.cc
+++ b/tests/lac/sparse_matrix_mmult_01.cc
@@ -56,9 +56,9 @@ void test (const unsigned int n)
   for (unsigned int i=0; i<n; ++i)
     {
       for (unsigned int j=0; j<n; ++j)
-        A.set(i,j,myrand());
+        A.set(i,j,Testing::rand());
       for (unsigned int j=0; j<n; ++j)
-        B.set(i,j,myrand());
+        B.set(i,j,Testing::rand());
     }
 
   // now form the matrix-matrix product and
@@ -67,7 +67,7 @@ void test (const unsigned int n)
 
   Vector<double> x(n), y(n), z(n), tmp(n);
   for (unsigned int j=0; j<n; ++j)
-    x(j) = myrand();
+    x(j) = Testing::rand();
 
   // then test for correctness
   C.vmult (y, x);
diff --git a/tests/lac/sparse_matrix_mmult_02.cc b/tests/lac/sparse_matrix_mmult_02.cc
index fc24fdff97..f94e783fc3 100644
--- a/tests/lac/sparse_matrix_mmult_02.cc
+++ b/tests/lac/sparse_matrix_mmult_02.cc
@@ -56,14 +56,14 @@ void test (const unsigned int n)
   for (unsigned int i=0; i<n; ++i)
     {
       for (unsigned int j=0; j<n; ++j)
-        A.set(i,j,myrand());
+        A.set(i,j,Testing::rand());
       for (unsigned int j=0; j<n; ++j)
-        B.set(i,j,myrand());
+        B.set(i,j,Testing::rand());
     }
 
   Vector<double> v(n);
   for (unsigned int j=0; j<n; ++j)
-    v(j) = myrand();
+    v(j) = Testing::rand();
 
   // now form the matrix-matrix product and
   // initialize a test rhs
@@ -71,7 +71,7 @@ void test (const unsigned int n)
 
   Vector<double> x(n), y(n), z(n), tmp(n);
   for (unsigned int j=0; j<n; ++j)
-    x(j) = myrand();
+    x(j) = Testing::rand();
 
   // then test for correctness
   C.vmult (y, x);
diff --git a/tests/lac/sparse_mic.cc b/tests/lac/sparse_mic.cc
index f1a03a00da..9dacf4b804 100644
--- a/tests/lac/sparse_mic.cc
+++ b/tests/lac/sparse_mic.cc
@@ -105,7 +105,7 @@ int main()
           for (unsigned int i=0; i<3; ++i)
             {
               for (unsigned int j=0; j<dim; ++j)
-                v(j) = 1. * myrand()/RAND_MAX;
+                v(j) = 1. * Testing::rand()/RAND_MAX;
 
               A.vmult (tmp1, v);
               mic.vmult (tmp2, tmp1);
diff --git a/tests/lac/vector_accumulation.cc b/tests/lac/vector_accumulation.cc
index 0c5cd1fafe..f2a65dd581 100644
--- a/tests/lac/vector_accumulation.cc
+++ b/tests/lac/vector_accumulation.cc
@@ -34,10 +34,10 @@ void check_norms ()
 {
   for (unsigned int test=0; test<20; ++test)
     {
-      const unsigned int size = myrand() % 100000;
+      const unsigned int size = Testing::rand() % 100000;
       Vector<number> vec (size);
       for (unsigned int i=0; i<size; ++i)
-        vec(i) = static_cast<number>(myrand())/static_cast<number>(RAND_MAX);
+        vec(i) = static_cast<number>(Testing::rand())/static_cast<number>(RAND_MAX);
       const typename Vector<number>::real_type norm = vec.l2_norm();
       for (unsigned int i=0; i<30; ++i)
         Assert (vec.l2_norm() == norm, ExcInternalError());
diff --git a/tests/matrix_free/copy.cc b/tests/matrix_free/copy.cc
index 6eebe79bd1..1ebee2cf6e 100644
--- a/tests/matrix_free/copy.cc
+++ b/tests/matrix_free/copy.cc
@@ -96,7 +96,7 @@ void sub_test()
         {
           if (constraints.is_constrained(i))
             continue;
-          const double entry = myrand()/(double)RAND_MAX;
+          const double entry = Testing::rand()/(double)RAND_MAX;
           in_dist(i) = entry;
         }
 
diff --git a/tests/matrix_free/evaluate_1d_shape.cc b/tests/matrix_free/evaluate_1d_shape.cc
index d908646d28..44cc266fe5 100644
--- a/tests/matrix_free/evaluate_1d_shape.cc
+++ b/tests/matrix_free/evaluate_1d_shape.cc
@@ -36,7 +36,7 @@ void test()
   for (unsigned int i=0; i<(M+1)/2; ++i)
     for (unsigned int j=0; j<N; ++j)
       {
-        shape[i][j] = -1. + 2. * (double)myrand()/RAND_MAX;
+        shape[i][j] = -1. + 2. * (double)Testing::rand()/RAND_MAX;
         if (type == 1)
           shape[M-1-i][N-1-j] = -shape[i][j];
         else
@@ -53,7 +53,7 @@ void test()
 
   double x[N], x_ref[N], y[M], y_ref[M];
   for (unsigned int i=0; i<N; ++i)
-    x[i] = (double)myrand()/RAND_MAX;
+    x[i] = (double)Testing::rand()/RAND_MAX;
 
   // compute reference
   for (unsigned int i=0; i<M; ++i)
@@ -82,7 +82,7 @@ void test()
 
 
   for (unsigned int i=0; i<M; ++i)
-    y[i] = (double)myrand()/RAND_MAX;
+    y[i] = (double)Testing::rand()/RAND_MAX;
 
   // compute reference
   for (unsigned int i=0; i<N; ++i)
diff --git a/tests/matrix_free/evaluate_1d_shape_evenodd.cc b/tests/matrix_free/evaluate_1d_shape_evenodd.cc
index a4b78a84fb..4db1704541 100644
--- a/tests/matrix_free/evaluate_1d_shape_evenodd.cc
+++ b/tests/matrix_free/evaluate_1d_shape_evenodd.cc
@@ -37,7 +37,7 @@ void test()
   for (unsigned int i=0; i<(M+1)/2; ++i)
     for (unsigned int j=0; j<N; ++j)
       {
-        shape[i][j] = -1. + 2. * (double)myrand()/RAND_MAX;
+        shape[i][j] = -1. + 2. * (double)Testing::rand()/RAND_MAX;
         if (type == 1)
           shape[M-1-i][N-1-j] = -shape[i][j];
         else
@@ -67,7 +67,7 @@ void test()
 
   double x[N], x_ref[N], y[M], y_ref[M];
   for (unsigned int i=0; i<N; ++i)
-    x[i] = (double)myrand()/RAND_MAX;
+    x[i] = (double)Testing::rand()/RAND_MAX;
 
   // compute reference
   for (unsigned int i=0; i<M; ++i)
@@ -88,7 +88,7 @@ void test()
 
 
   for (unsigned int i=0; i<M; ++i)
-    y[i] = (double)myrand()/RAND_MAX;
+    y[i] = (double)Testing::rand()/RAND_MAX;
 
   // compute reference
   for (unsigned int i=0; i<N; ++i)
diff --git a/tests/matrix_free/get_functions_common.h b/tests/matrix_free/get_functions_common.h
index 6e23e6e2ba..aba500ec68 100644
--- a/tests/matrix_free/get_functions_common.h
+++ b/tests/matrix_free/get_functions_common.h
@@ -231,7 +231,7 @@ void do_test (const DoFHandler<dim> &dof,
     {
       if (constraints.is_constrained(i))
         continue;
-      const double entry = rand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       solution(i) = entry;
     }
 
diff --git a/tests/matrix_free/get_functions_gl.cc b/tests/matrix_free/get_functions_gl.cc
index af8131c811..c585c8e235 100644
--- a/tests/matrix_free/get_functions_gl.cc
+++ b/tests/matrix_free/get_functions_gl.cc
@@ -120,7 +120,7 @@ void test ()
     {
       if (constraints.is_constrained(i))
         continue;
-      const double entry = myrand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       solution(i) = entry;
     }
   constraints.distribute(solution);
diff --git a/tests/matrix_free/get_functions_mappingq.cc b/tests/matrix_free/get_functions_mappingq.cc
index 20e8195b6a..27e1179631 100644
--- a/tests/matrix_free/get_functions_mappingq.cc
+++ b/tests/matrix_free/get_functions_mappingq.cc
@@ -67,7 +67,7 @@ void test ()
     {
       if (constraints.is_constrained(i))
         continue;
-      const double entry = myrand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       solution(i) = entry;
     }
 
diff --git a/tests/matrix_free/get_functions_multife.cc b/tests/matrix_free/get_functions_multife.cc
index 8b0c421dd2..a6d8743e5c 100644
--- a/tests/matrix_free/get_functions_multife.cc
+++ b/tests/matrix_free/get_functions_multife.cc
@@ -249,7 +249,7 @@ void test ()
       {
         if (constraints[no]->is_constrained(i))
           continue;
-        const double entry = myrand()/(double)RAND_MAX;
+        const double entry = Testing::rand()/(double)RAND_MAX;
         src[no](i) = entry;
       }
 
diff --git a/tests/matrix_free/get_functions_multife2.cc b/tests/matrix_free/get_functions_multife2.cc
index 289f54dc6a..b4e8cd81f5 100644
--- a/tests/matrix_free/get_functions_multife2.cc
+++ b/tests/matrix_free/get_functions_multife2.cc
@@ -300,7 +300,7 @@ void test ()
       {
         if (constraints[no]->is_constrained(i))
           continue;
-        const double entry = myrand()/(double)RAND_MAX;
+        const double entry = Testing::rand()/(double)RAND_MAX;
         src[no](i) = entry;
       }
 
diff --git a/tests/matrix_free/get_functions_q_hierarchical.cc b/tests/matrix_free/get_functions_q_hierarchical.cc
index e360aba21d..ce3c796e20 100644
--- a/tests/matrix_free/get_functions_q_hierarchical.cc
+++ b/tests/matrix_free/get_functions_q_hierarchical.cc
@@ -123,7 +123,7 @@ void test ()
     {
       if (constraints.is_constrained(i))
         continue;
-      const double entry = myrand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       solution(i) = entry;
     }
   constraints.distribute(solution);
diff --git a/tests/matrix_free/get_functions_rect.cc b/tests/matrix_free/get_functions_rect.cc
index d99ffd9842..e31f4a5a79 100644
--- a/tests/matrix_free/get_functions_rect.cc
+++ b/tests/matrix_free/get_functions_rect.cc
@@ -89,7 +89,7 @@ void test ()
     {
       if (constraints.is_constrained(i))
         continue;
-      const double entry = myrand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       solution(i) = entry;
     }
   constraints.distribute(solution);
diff --git a/tests/matrix_free/get_functions_variants.cc b/tests/matrix_free/get_functions_variants.cc
index 27665233d6..5a95affffb 100644
--- a/tests/matrix_free/get_functions_variants.cc
+++ b/tests/matrix_free/get_functions_variants.cc
@@ -169,7 +169,7 @@ void test ()
   // create vector with random entries
   for (unsigned int i=0; i<dof.n_dofs(); ++i)
     {
-      const double entry = myrand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       solution_dist(i) = entry;
     }
 
diff --git a/tests/matrix_free/get_values_plain.cc b/tests/matrix_free/get_values_plain.cc
index 4da355bb68..3cb6681f43 100644
--- a/tests/matrix_free/get_values_plain.cc
+++ b/tests/matrix_free/get_values_plain.cc
@@ -120,7 +120,7 @@ void do_test (const DoFHandler<dim> &dof,
     {
       if (constraints.is_constrained(i))
         continue;
-      const double entry = myrand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       solution(i) = entry;
     }
 
diff --git a/tests/matrix_free/integrate_functions.cc b/tests/matrix_free/integrate_functions.cc
index 096b293c6d..f42acc2e60 100644
--- a/tests/matrix_free/integrate_functions.cc
+++ b/tests/matrix_free/integrate_functions.cc
@@ -107,9 +107,9 @@ operator () (const MatrixFree<dim,Number> &data,
           // and their gradients
           for (unsigned int q=0; q<n_q_points; ++q)
             {
-              values[q][j] = myrand()/(double)RAND_MAX;
+              values[q][j] = Testing::rand()/(double)RAND_MAX;
               for (unsigned int d=0; d<dim; ++d)
-                gradients[q*dim+d][j] = -1. + 2. * (myrand()/(double)RAND_MAX);
+                gradients[q*dim+d][j] = -1. + 2. * (Testing::rand()/(double)RAND_MAX);
             }
           fe_val.reinit (data.get_cell_iterator(cell,j));
           data.get_cell_iterator(cell,j)->get_dof_indices(dof_indices);
diff --git a/tests/matrix_free/integrate_functions_multife.cc b/tests/matrix_free/integrate_functions_multife.cc
index 904e9bdbf8..bef0320048 100644
--- a/tests/matrix_free/integrate_functions_multife.cc
+++ b/tests/matrix_free/integrate_functions_multife.cc
@@ -123,9 +123,9 @@ operator () (const MatrixFree<dim,Number> &data,
           // and their gradients
           for (unsigned int q=0; q<n_q_points0; ++q)
             {
-              values0[q][j] = myrand()/(double)RAND_MAX;
+              values0[q][j] = Testing::rand()/(double)RAND_MAX;
               for (unsigned int d=0; d<dim; ++d)
-                gradients0[q*dim+d][j] = -1. + 2. * (myrand()/(double)RAND_MAX);
+                gradients0[q*dim+d][j] = -1. + 2. * (Testing::rand()/(double)RAND_MAX);
             }
           fe_val0.reinit (data.get_cell_iterator(cell,j,0));
           data.get_cell_iterator(cell,j,0)->get_dof_indices(dof_indices0);
@@ -149,9 +149,9 @@ operator () (const MatrixFree<dim,Number> &data,
 
           for (unsigned int q=0; q<n_q_points1; ++q)
             {
-              values1[q][j] = myrand()/(double)RAND_MAX;
+              values1[q][j] = Testing::rand()/(double)RAND_MAX;
               for (unsigned int d=0; d<dim; ++d)
-                gradients1[q*dim+d][j] = -1. + 2. * (myrand()/(double)RAND_MAX);
+                gradients1[q*dim+d][j] = -1. + 2. * (Testing::rand()/(double)RAND_MAX);
             }
           for (unsigned int i=0; i<dofs_per_cell1; ++i)
             {
diff --git a/tests/matrix_free/integrate_functions_multife2.cc b/tests/matrix_free/integrate_functions_multife2.cc
index 7a182c0808..804f875952 100644
--- a/tests/matrix_free/integrate_functions_multife2.cc
+++ b/tests/matrix_free/integrate_functions_multife2.cc
@@ -126,9 +126,9 @@ operator () (const MatrixFree<dim,Number> &data,
           // and their gradients
           for (unsigned int q=0; q<n_q_points0; ++q)
             {
-              values0[q][j] = myrand()/(double)RAND_MAX;
+              values0[q][j] = Testing::rand()/(double)RAND_MAX;
               for (unsigned int d=0; d<dim; ++d)
-                gradients0[q*dim+d][j] = -1. + 2. * (myrand()/(double)RAND_MAX);
+                gradients0[q*dim+d][j] = -1. + 2. * (Testing::rand()/(double)RAND_MAX);
             }
           fe_val0.reinit (data.get_cell_iterator(cell,j,0));
           data.get_cell_iterator(cell,j,0)->get_dof_indices(dof_indices0);
@@ -152,9 +152,9 @@ operator () (const MatrixFree<dim,Number> &data,
 
           for (unsigned int q=0; q<n_q_points1; ++q)
             {
-              values1[q][j] = myrand()/(double)RAND_MAX;
+              values1[q][j] = Testing::rand()/(double)RAND_MAX;
               for (unsigned int d=0; d<dim; ++d)
-                gradients1[q*dim+d][j] = -1. + 2. * (myrand()/(double)RAND_MAX);
+                gradients1[q*dim+d][j] = -1. + 2. * (Testing::rand()/(double)RAND_MAX);
             }
           for (unsigned int i=0; i<dofs_per_cell1; ++i)
             {
diff --git a/tests/matrix_free/matrix_vector_10.cc b/tests/matrix_free/matrix_vector_10.cc
index 6be8ed7321..b460be662d 100644
--- a/tests/matrix_free/matrix_vector_10.cc
+++ b/tests/matrix_free/matrix_vector_10.cc
@@ -129,7 +129,7 @@ void test ()
         owned_set.nth_index_in_set (i);
       if (constraints.is_constrained(glob_index))
         continue;
-      in.local_element(i) = (double)myrand()/RAND_MAX;
+      in.local_element(i) = (double)Testing::rand()/RAND_MAX;
     }
 
   mf.vmult (out, in);
diff --git a/tests/matrix_free/matrix_vector_11.cc b/tests/matrix_free/matrix_vector_11.cc
index daedda3923..1d111f2a4b 100644
--- a/tests/matrix_free/matrix_vector_11.cc
+++ b/tests/matrix_free/matrix_vector_11.cc
@@ -122,7 +122,7 @@ void test ()
         owned_set.nth_index_in_set (i);
       if (constraints.is_constrained(glob_index))
         continue;
-      in.local_element(i) = (double)myrand()/RAND_MAX;
+      in.local_element(i) = (double)Testing::rand()/RAND_MAX;
     }
 
   mf.vmult (ref, in);
diff --git a/tests/matrix_free/matrix_vector_12.cc b/tests/matrix_free/matrix_vector_12.cc
index 60b18c3a94..0c497fe15b 100644
--- a/tests/matrix_free/matrix_vector_12.cc
+++ b/tests/matrix_free/matrix_vector_12.cc
@@ -188,8 +188,8 @@ void test ()
         owned_set.nth_index_in_set (i);
       if (constraints.is_constrained(glob_index))
         continue;
-      in[0].local_element(i) = (double)myrand()/RAND_MAX;
-      in[1].local_element(i) = (double)myrand()/RAND_MAX;
+      in[0].local_element(i) = (double)Testing::rand()/RAND_MAX;
+      in[1].local_element(i) = (double)Testing::rand()/RAND_MAX;
     }
 
   mf.vmult (out, in);
diff --git a/tests/matrix_free/matrix_vector_13.cc b/tests/matrix_free/matrix_vector_13.cc
index dcbf078bd5..c777c8c79c 100644
--- a/tests/matrix_free/matrix_vector_13.cc
+++ b/tests/matrix_free/matrix_vector_13.cc
@@ -191,8 +191,8 @@ void test ()
         owned_set.nth_index_in_set (i);
       if (constraints.is_constrained(glob_index))
         continue;
-      in.block(0).local_element(i) = (double)myrand()/RAND_MAX;
-      in.block(1).local_element(i) = (double)myrand()/RAND_MAX;
+      in.block(0).local_element(i) = (double)Testing::rand()/RAND_MAX;
+      in.block(1).local_element(i) = (double)Testing::rand()/RAND_MAX;
     }
 
   mf.vmult (out, in);
diff --git a/tests/matrix_free/matrix_vector_common.h b/tests/matrix_free/matrix_vector_common.h
index f1c5b3b796..40e715a5ac 100644
--- a/tests/matrix_free/matrix_vector_common.h
+++ b/tests/matrix_free/matrix_vector_common.h
@@ -82,7 +82,7 @@ void do_test (const DoFHandler<dim> &dof,
     {
       if (constraints.is_constrained(i))
         continue;
-      const double entry = rand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       in(i) = entry;
       in_dist(i) = entry;
     }
diff --git a/tests/matrix_free/matrix_vector_curl.cc b/tests/matrix_free/matrix_vector_curl.cc
index 36fcc82a45..773ae0fc43 100644
--- a/tests/matrix_free/matrix_vector_curl.cc
+++ b/tests/matrix_free/matrix_vector_curl.cc
@@ -239,7 +239,7 @@ void test ()
   for (unsigned int i=0; i<dim; ++i)
     for (unsigned int j=0; j<system_rhs.block(i).size(); ++j)
       {
-        const double val = -1. + 2.*(double)myrand()/double(RAND_MAX);
+        const double val = -1. + 2.*(double)Testing::rand()/double(RAND_MAX);
         system_rhs.block(i)(j) = val;
       }
   constraints.condense(system_rhs);
diff --git a/tests/matrix_free/matrix_vector_div.cc b/tests/matrix_free/matrix_vector_div.cc
index c4af939d56..736e9a8c81 100644
--- a/tests/matrix_free/matrix_vector_div.cc
+++ b/tests/matrix_free/matrix_vector_div.cc
@@ -243,7 +243,7 @@ void test ()
   for (unsigned int i=0; i<dim; ++i)
     for (unsigned int j=0; j<system_rhs.block(i).size(); ++j)
       {
-        const double val = -1. + 2.*(double)myrand()/double(RAND_MAX);
+        const double val = -1. + 2.*(double)Testing::rand()/double(RAND_MAX);
         system_rhs.block(i)(j) = val;
       }
   constraints.condense(system_rhs);
diff --git a/tests/matrix_free/matrix_vector_hp.cc b/tests/matrix_free/matrix_vector_hp.cc
index 19b31a5fe1..b359397e68 100644
--- a/tests/matrix_free/matrix_vector_hp.cc
+++ b/tests/matrix_free/matrix_vector_hp.cc
@@ -137,7 +137,7 @@ void test ()
     endc = dof.end();
     for (; cell!=endc; ++cell)
       {
-        const unsigned int fe_index = myrand() % max_degree;
+        const unsigned int fe_index = Testing::rand() % max_degree;
         cell->set_active_fe_index (fe_index);
       }
   }
@@ -222,7 +222,7 @@ void test ()
   for (unsigned int i=0; i<dof.n_dofs(); ++i)
     {
       if (constraints.is_constrained(i) == false)
-        src(i) = (double)myrand()/RAND_MAX;
+        src(i) = (double)Testing::rand()/RAND_MAX;
     }
 
   // now perform matrix-vector product and check
diff --git a/tests/matrix_free/matrix_vector_mg.cc b/tests/matrix_free/matrix_vector_mg.cc
index 1bbc885654..4eb2670a22 100644
--- a/tests/matrix_free/matrix_vector_mg.cc
+++ b/tests/matrix_free/matrix_vector_mg.cc
@@ -177,7 +177,7 @@ void test ()
   for (unsigned int i=0; i<dof.n_dofs(); ++i)
     {
       if (constraints.is_constrained(i) == false)
-        src(i) = (double)myrand()/RAND_MAX;
+        src(i) = (double)Testing::rand()/RAND_MAX;
     }
 
   // now perform matrix-vector product and check
@@ -198,7 +198,7 @@ void test ()
       for (unsigned int i=0; i<dof.n_dofs(level); ++i)
         {
           if (mg_constraints[level].is_constrained(i) == false)
-            src(i) = (double)myrand()/RAND_MAX;
+            src(i) = (double)Testing::rand()/RAND_MAX;
         }
 
       // now perform matrix-vector product and check
diff --git a/tests/matrix_free/matrix_vector_stokes.cc b/tests/matrix_free/matrix_vector_stokes.cc
index fd8d549a3d..46e1463c3a 100644
--- a/tests/matrix_free/matrix_vector_stokes.cc
+++ b/tests/matrix_free/matrix_vector_stokes.cc
@@ -268,7 +268,7 @@ void test ()
   for (unsigned int i=0; i<dim+1; ++i)
     for (unsigned int j=0; j<system_rhs.block(i).size(); ++j)
       {
-        const double val = -1. + 2.*(double)myrand()/double(RAND_MAX);
+        const double val = -1. + 2.*(double)Testing::rand()/double(RAND_MAX);
         system_rhs.block(i)(j) = val;
         vec1[i](j) = val;
       }
diff --git a/tests/matrix_free/matrix_vector_stokes_noflux.cc b/tests/matrix_free/matrix_vector_stokes_noflux.cc
index 3d01a123e8..1860f6797e 100644
--- a/tests/matrix_free/matrix_vector_stokes_noflux.cc
+++ b/tests/matrix_free/matrix_vector_stokes_noflux.cc
@@ -283,13 +283,13 @@ void test ()
   for (unsigned int j=0; j<system_rhs.block(0).size(); ++j)
     if (constraints_u.is_constrained(j) == false)
       {
-        const double val = -1 + 2.*(double)myrand()/double(RAND_MAX);
+        const double val = -1 + 2.*(double)Testing::rand()/double(RAND_MAX);
         system_rhs.block(0)(j) = val;
       }
   for (unsigned int j=0; j<system_rhs.block(1).size(); ++j)
     if (constraints_p.is_constrained(j) == false)
       {
-        const double val = -1 + 2.*(double)myrand()/double(RAND_MAX);
+        const double val = -1 + 2.*(double)Testing::rand()/double(RAND_MAX);
         system_rhs.block(1)(j) = val;
       }
 
diff --git a/tests/matrix_free/no_index_initialize.cc b/tests/matrix_free/no_index_initialize.cc
index 4b4846e7c7..2c29799abb 100644
--- a/tests/matrix_free/no_index_initialize.cc
+++ b/tests/matrix_free/no_index_initialize.cc
@@ -126,7 +126,7 @@ void do_test (const DoFHandler<dim> &dof,
     {
       if (constraints.is_constrained(i))
         continue;
-      const double entry = myrand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       solution(i) = entry;
     }
 
diff --git a/tests/matrix_free/thread_correctness.cc b/tests/matrix_free/thread_correctness.cc
index a108ee3758..edbd65c2ac 100644
--- a/tests/matrix_free/thread_correctness.cc
+++ b/tests/matrix_free/thread_correctness.cc
@@ -110,7 +110,7 @@ void sub_test()
         {
           if (constraints.is_constrained(i))
             continue;
-          const double entry = myrand()/(double)RAND_MAX;
+          const double entry = Testing::rand()/(double)RAND_MAX;
           in_dist(i) = entry;
         }
 
diff --git a/tests/matrix_free/thread_correctness_hp.cc b/tests/matrix_free/thread_correctness_hp.cc
index 6c0130e573..1cb0a8cb88 100644
--- a/tests/matrix_free/thread_correctness_hp.cc
+++ b/tests/matrix_free/thread_correctness_hp.cc
@@ -107,7 +107,7 @@ void do_test (const unsigned int parallel_option)
       cell = tria.begin_active (),
       endc = tria.end();
       for (; cell!=endc; ++cell)
-        if (myrand() % (7-i) == 0)
+        if (Testing::rand() % (7-i) == 0)
           cell->set_refine_flag();
       tria.execute_coarsening_and_refinement();
     }
@@ -131,7 +131,7 @@ void do_test (const unsigned int parallel_option)
     endc = dof.end();
     for (; cell!=endc; ++cell)
       {
-        const unsigned int fe_index = myrand() % max_degree;
+        const unsigned int fe_index = Testing::rand() % max_degree;
         cell->set_active_fe_index (fe_index);
       }
   }
@@ -189,7 +189,7 @@ void do_test (const unsigned int parallel_option)
       for (unsigned int i=0; i<dof.n_dofs(); ++i)
         {
           if (constraints.is_constrained(i) == false)
-            src(i) = (double)myrand()/RAND_MAX;
+            src(i) = (double)Testing::rand()/RAND_MAX;
         }
 
       // now perform 30 matrix-vector products in
diff --git a/tests/matrix_free/update_mapping_only.cc b/tests/matrix_free/update_mapping_only.cc
index b24169b7ca..d24f62df7d 100644
--- a/tests/matrix_free/update_mapping_only.cc
+++ b/tests/matrix_free/update_mapping_only.cc
@@ -58,7 +58,7 @@ void test ()
     {
       if (constraints.is_constrained(i))
         continue;
-      const double entry = myrand()/(double)RAND_MAX;
+      const double entry = Testing::rand()/(double)RAND_MAX;
       solution(i) = entry;
     }
 
diff --git a/tests/mpi/constraint_matrix_trilinos_bug.cc b/tests/mpi/constraint_matrix_trilinos_bug.cc
index 4a11e07d39..86b61ee225 100644
--- a/tests/mpi/constraint_matrix_trilinos_bug.cc
+++ b/tests/mpi/constraint_matrix_trilinos_bug.cc
@@ -68,7 +68,7 @@ void test()
       endc = tr.end();
 
       for (; cell!=endc; ++cell)
-        if (myrand()%42==1)
+        if (Testing::rand()%42==1)
           cell->set_refine_flag ();
 
       tr.execute_coarsening_and_refinement ();
diff --git a/tests/mpi/count_dofs_per_block_02.cc b/tests/mpi/count_dofs_per_block_02.cc
index dbc6c47502..958357ff26 100644
--- a/tests/mpi/count_dofs_per_block_02.cc
+++ b/tests/mpi/count_dofs_per_block_02.cc
@@ -60,7 +60,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/mpi/crash_03.cc b/tests/mpi/crash_03.cc
index cea95adca7..7cac982c0b 100644
--- a/tests/mpi/crash_03.cc
+++ b/tests/mpi/crash_03.cc
@@ -54,7 +54,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/mpi/dof_handler_number_cache.cc b/tests/mpi/dof_handler_number_cache.cc
index 9a07498f78..58d7f29536 100644
--- a/tests/mpi/dof_handler_number_cache.cc
+++ b/tests/mpi/dof_handler_number_cache.cc
@@ -62,7 +62,7 @@ void test()
       // refine one-fifth of cells randomly
       std::vector<bool> flags (triangulation.n_active_cells(), false);
       for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
-        flags[myrand() % flags.size()] = true;
+        flags[Testing::rand() % flags.size()] = true;
       // make sure there's at least one that
       // will be refined
       flags[0] = true;
diff --git a/tests/mpi/mg_05.cc b/tests/mpi/mg_05.cc
index 7cb6335de5..fc90013800 100644
--- a/tests/mpi/mg_05.cc
+++ b/tests/mpi/mg_05.cc
@@ -101,7 +101,7 @@ void test()
         endc = tr.end();
       
       for (; cell!=endc; ++cell)
-	if (myrand()%42==1)
+	if (Testing::rand()%42==1)
 	  cell->set_refine_flag ();
       
       tr.execute_coarsening_and_refinement ();
diff --git a/tests/mpi/mg_06.cc b/tests/mpi/mg_06.cc
index 7e342327ea..2413116931 100644
--- a/tests/mpi/mg_06.cc
+++ b/tests/mpi/mg_06.cc
@@ -103,7 +103,7 @@ void test()
         endc = tr.end();
       
       for (; cell!=endc; ++cell)
-	if (myrand()%42==1)
+	if (Testing::rand()%42==1)
 	  cell->set_refine_flag ();
       
       tr.execute_coarsening_and_refinement ();
diff --git a/tests/mpi/p4est_2d_constraintmatrix_03.cc b/tests/mpi/p4est_2d_constraintmatrix_03.cc
index 77a2878c97..420c802021 100644
--- a/tests/mpi/p4est_2d_constraintmatrix_03.cc
+++ b/tests/mpi/p4est_2d_constraintmatrix_03.cc
@@ -123,7 +123,7 @@ void test()
       endc = tr.end();
 
       for (; cell!=endc; ++cell)
-        if (myrand()%42==1)
+        if (Testing::rand()%42==1)
           cell->set_refine_flag ();
 
       tr.execute_coarsening_and_refinement ();
diff --git a/tests/mpi/p4est_2d_constraintmatrix_04.cc b/tests/mpi/p4est_2d_constraintmatrix_04.cc
index eca7083716..ccbd947518 100644
--- a/tests/mpi/p4est_2d_constraintmatrix_04.cc
+++ b/tests/mpi/p4est_2d_constraintmatrix_04.cc
@@ -160,7 +160,7 @@ void test()
         endc = tr.end();
 
         for (; cell!=endc; ++cell)
-          if (myrand()%42==1)
+          if (Testing::rand()%42==1)
             cell->set_refine_flag ();
 
         tr.execute_coarsening_and_refinement ();
@@ -200,9 +200,9 @@ void test()
         for (; cell!=endc; ++cell)
           if (!cell->is_artificial() && !cell->is_ghost())
             {
-              if (myrand()%12==1)
+              if (Testing::rand()%12==1)
                 cell->set_refine_flag ();
-              else if (myrand()%7==1)
+              else if (Testing::rand()%7==1)
                 cell->set_coarsen_flag ();
             }
       }
diff --git a/tests/mpi/p4est_2d_dofhandler_02.cc b/tests/mpi/p4est_2d_dofhandler_02.cc
index 4de08b4f54..d2eab70e0d 100644
--- a/tests/mpi/p4est_2d_dofhandler_02.cc
+++ b/tests/mpi/p4est_2d_dofhandler_02.cc
@@ -59,7 +59,7 @@ void test()
           std::vector<bool> flags (tr.n_locally_owned_active_cells(), false);
           for (unsigned int i=0; i<tr.n_locally_owned_active_cells() / 5 + 1; ++i)
             {
-              const unsigned int x = myrand() % flags.size();
+              const unsigned int x = Testing::rand() % flags.size();
               flags[x] = true;
             }
 
diff --git a/tests/mpi/p4est_2d_dofhandler_03.cc b/tests/mpi/p4est_2d_dofhandler_03.cc
index 4e89fe123a..d60c5dc1dd 100644
--- a/tests/mpi/p4est_2d_dofhandler_03.cc
+++ b/tests/mpi/p4est_2d_dofhandler_03.cc
@@ -60,7 +60,7 @@ void test()
           std::vector<bool> flags (tr.n_locally_owned_active_cells(), false);
           for (unsigned int i=0; i<tr.n_locally_owned_active_cells() / 5 + 1; ++i)
             {
-              const unsigned int x = myrand() % flags.size();
+              const unsigned int x = Testing::rand() % flags.size();
               flags[x] = true;
             }
 
diff --git a/tests/mpi/p4est_2d_dofhandler_04.cc b/tests/mpi/p4est_2d_dofhandler_04.cc
index 93ebb77409..e4d50ce573 100644
--- a/tests/mpi/p4est_2d_dofhandler_04.cc
+++ b/tests/mpi/p4est_2d_dofhandler_04.cc
@@ -60,7 +60,7 @@ void test()
           std::vector<bool> flags (tr.n_locally_owned_active_cells(), false);
           for (unsigned int i=0; i<tr.n_locally_owned_active_cells() / 5 + 1; ++i)
             {
-              const unsigned int x = myrand() % flags.size();
+              const unsigned int x = Testing::rand() % flags.size();
               flags[x] = true;
             }
 
diff --git a/tests/mpi/p4est_2d_refine_03.cc b/tests/mpi/p4est_2d_refine_03.cc
index c1b1acf6cf..672b8d3d6d 100644
--- a/tests/mpi/p4est_2d_refine_03.cc
+++ b/tests/mpi/p4est_2d_refine_03.cc
@@ -60,7 +60,7 @@ void test()
           // will be refined.
           for (unsigned int i=0; i<tr.n_active_cells() / 5 + 1; ++i)
             {
-              const unsigned int x = myrand() % flags.size();
+              const unsigned int x = Testing::rand() % flags.size();
               flags[x] = true;
             }
 
diff --git a/tests/mpi/p4est_3d_constraintmatrix_01.cc b/tests/mpi/p4est_3d_constraintmatrix_01.cc
index 23261eb4f6..3cf00b8c8e 100644
--- a/tests/mpi/p4est_3d_constraintmatrix_01.cc
+++ b/tests/mpi/p4est_3d_constraintmatrix_01.cc
@@ -60,7 +60,7 @@ void test()
       endc = tr.end();
 
       for (; cell!=endc; ++cell)
-        if (myrand()%42==1)
+        if (Testing::rand()%42==1)
           cell->set_refine_flag ();
 
       tr.execute_coarsening_and_refinement ();
diff --git a/tests/mpi/p4est_3d_constraintmatrix_02.cc b/tests/mpi/p4est_3d_constraintmatrix_02.cc
index 7745179056..6299420a9a 100644
--- a/tests/mpi/p4est_3d_constraintmatrix_02.cc
+++ b/tests/mpi/p4est_3d_constraintmatrix_02.cc
@@ -60,7 +60,7 @@ void test()
       endc = tr.end();
 
       for (; cell!=endc; ++cell)
-        if (myrand()%42==1)
+        if (Testing::rand()%42==1)
           cell->set_refine_flag ();
 
       tr.execute_coarsening_and_refinement ();
diff --git a/tests/mpi/p4est_3d_constraintmatrix_03.cc b/tests/mpi/p4est_3d_constraintmatrix_03.cc
index 1b47188f1b..c0b1b41b55 100644
--- a/tests/mpi/p4est_3d_constraintmatrix_03.cc
+++ b/tests/mpi/p4est_3d_constraintmatrix_03.cc
@@ -159,7 +159,7 @@ void test()
         endc = tr.end();
 
         for (; cell!=endc; ++cell)
-          if (myrand()%42==1)
+          if (Testing::rand()%42==1)
             cell->set_refine_flag ();
 
         tr.execute_coarsening_and_refinement ();
@@ -199,9 +199,9 @@ void test()
         for (; cell!=endc; ++cell)
           if (!cell->is_artificial() && !cell->is_ghost())
             {
-              if (myrand()%12==1)
+              if (Testing::rand()%12==1)
                 cell->set_refine_flag ();
-              else if (myrand()%7==1)
+              else if (Testing::rand()%7==1)
                 cell->set_coarsen_flag ();
             }
       }
diff --git a/tests/mpi/p4est_3d_refine_01.cc b/tests/mpi/p4est_3d_refine_01.cc
index 84eb720cd6..05f499e711 100644
--- a/tests/mpi/p4est_3d_refine_01.cc
+++ b/tests/mpi/p4est_3d_refine_01.cc
@@ -60,7 +60,7 @@ void test()
           // will be refined.
           for (unsigned int i=0; i<tr.n_active_cells() / 5 + 1; ++i)
             {
-              const unsigned int x = myrand() % flags.size();
+              const unsigned int x = Testing::rand() % flags.size();
               flags[x] = true;
             }
 
diff --git a/tests/mpi/trilinos_matvec_01.cc b/tests/mpi/trilinos_matvec_01.cc
index 94ea3e4a44..4ded88c338 100644
--- a/tests/mpi/trilinos_matvec_01.cc
+++ b/tests/mpi/trilinos_matvec_01.cc
@@ -100,7 +100,7 @@ void test ()
   for (unsigned int i=0; i<col_partitioning.n_elements(); ++i)
     {
       const unsigned int global_index = col_partitioning.nth_index_in_set(i);
-      dx(global_index) = (double)myrand()/RAND_MAX;
+      dx(global_index) = (double)Testing::rand()/RAND_MAX;
       x(global_index)  = dx(global_index);
     }
   dy = 1.;
diff --git a/tests/mpi/trilinos_matvec_02.cc b/tests/mpi/trilinos_matvec_02.cc
index 6c601200e6..ef3e17b155 100644
--- a/tests/mpi/trilinos_matvec_02.cc
+++ b/tests/mpi/trilinos_matvec_02.cc
@@ -100,7 +100,7 @@ void test ()
   for (unsigned int i=0; i<row_partitioning.n_elements(); ++i)
     {
       const unsigned int global_index = row_partitioning.nth_index_in_set(i);
-      dy(global_index) = (double)myrand()/RAND_MAX;
+      dy(global_index) = (double)Testing::rand()/RAND_MAX;
       y(global_index)  = dy(global_index);
     }
   dx = 1.;
diff --git a/tests/mpi/trilinos_matvec_03.cc b/tests/mpi/trilinos_matvec_03.cc
index 7649ecb009..bd92dc2c9c 100644
--- a/tests/mpi/trilinos_matvec_03.cc
+++ b/tests/mpi/trilinos_matvec_03.cc
@@ -94,7 +94,7 @@ void test ()
   for (unsigned int i=0; i<col_partitioning.n_elements(); ++i)
     {
       const unsigned int global_index = col_partitioning.nth_index_in_set(i);
-      dx(global_index) = (double)myrand()/RAND_MAX;
+      dx(global_index) = (double)Testing::rand()/RAND_MAX;
       x(global_index)  = dx(global_index);
     }
   dy = 1.;
diff --git a/tests/petsc/66.cc b/tests/petsc/66.cc
index 6826926027..0c950d7cd1 100644
--- a/tests/petsc/66.cc
+++ b/tests/petsc/66.cc
@@ -40,7 +40,7 @@ void test (PETScWrappers::MatrixBase &m)
     {
       if (i>=5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.add (i,i-5, s);
           norm_sqr += s*s;
           ++nnz;
@@ -48,13 +48,13 @@ void test (PETScWrappers::MatrixBase &m)
 
       if (i<N-5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.add (i,i+5, s);
           norm_sqr += s*s;
           ++nnz;
         }
 
-      const double s = myrand();
+      const double s = Testing::rand();
       m.add (i,i,s);
       norm_sqr += s*s;
       ++nnz;
diff --git a/tests/petsc/67.cc b/tests/petsc/67.cc
index 34febf2126..3f5ce10a71 100644
--- a/tests/petsc/67.cc
+++ b/tests/petsc/67.cc
@@ -42,7 +42,7 @@ void test (PETScWrappers::MatrixBase &m)
     {
       if (i>=5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.add (i,i-5, s);
           norm_sqr += s*s;
           ++nnz;
@@ -50,13 +50,13 @@ void test (PETScWrappers::MatrixBase &m)
 
       if (i<N-5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.add (i,i+5, s);
           norm_sqr += s*s;
           ++nnz;
         }
 
-      const double s = myrand();
+      const double s = Testing::rand();
       m.add (i,i,s);
       norm_sqr += s*s;
       ++nnz;
diff --git a/tests/petsc/68.cc b/tests/petsc/68.cc
index 802565b7c7..d33074117c 100644
--- a/tests/petsc/68.cc
+++ b/tests/petsc/68.cc
@@ -40,7 +40,7 @@ void test (PETScWrappers::MatrixBase &m)
     {
       if (i>=5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.add (i,i-5, s);
           norm_sqr += s*s;
           ++nnz;
@@ -48,13 +48,13 @@ void test (PETScWrappers::MatrixBase &m)
 
       if (i<N-5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.add (i,i+5, s);
           norm_sqr += s*s;
           ++nnz;
         }
 
-      const double s = myrand();
+      const double s = Testing::rand();
       m.add (i,i,s);
       norm_sqr += s*s;
       ++nnz;
@@ -72,7 +72,7 @@ void test (PETScWrappers::MatrixBase &m)
 
   // now remove the entries of row N/2. set
   // diagonal entries to rnd
-  const double rnd = myrand();
+  const double rnd = Testing::rand();
   for (unsigned int i=0; i<N; ++i)
     {
       const double s = m.el(N/2,i);
diff --git a/tests/petsc/69.cc b/tests/petsc/69.cc
index 1d4e91eca8..5ec674841c 100644
--- a/tests/petsc/69.cc
+++ b/tests/petsc/69.cc
@@ -42,7 +42,7 @@ void test (PETScWrappers::MatrixBase &m)
     {
       if (i>=5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.add (i,i-5, s);
           norm_sqr += s*s;
           ++nnz;
@@ -50,13 +50,13 @@ void test (PETScWrappers::MatrixBase &m)
 
       if (i<N-5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.add (i,i+5, s);
           norm_sqr += s*s;
           ++nnz;
         }
 
-      const double s = myrand();
+      const double s = Testing::rand();
       m.add (i,i,s);
       norm_sqr += s*s;
       ++nnz;
@@ -74,7 +74,7 @@ void test (PETScWrappers::MatrixBase &m)
 
   // now remove the entries of rows N/2 and
   // N/3. set diagonal entries to rnd
-  const double rnd = myrand();
+  const double rnd = Testing::rand();
   for (size_type i=0; i<N; ++i)
     {
       const double s = m.el(N/2,i);
diff --git a/tests/petsc/slowness_04.cc b/tests/petsc/slowness_04.cc
index 0b75bfa048..0a616f8d4d 100644
--- a/tests/petsc/slowness_04.cc
+++ b/tests/petsc/slowness_04.cc
@@ -53,7 +53,7 @@ void test ()
       {
         // pick a random element among the
         // unused indices
-        const unsigned int k = myrand() % (N-i);
+        const unsigned int k = Testing::rand() % (N-i);
         permutation[i] = unused_indices[k];
 
         // then swap this used element to the
@@ -75,26 +75,26 @@ void test ()
         const unsigned int j=permutation[j_];
 
         const unsigned int global = i*N+j;
-        matrix.add(global, global, myrand());
+        matrix.add(global, global, Testing::rand());
         if (j>0)
           {
-            matrix.add(global-1, global, myrand());
-            matrix.add(global, global-1, myrand());
+            matrix.add(global-1, global, Testing::rand());
+            matrix.add(global, global-1, Testing::rand());
           }
         if (j<N-1)
           {
-            matrix.add(global+1, global, myrand());
-            matrix.add(global, global+1, myrand());
+            matrix.add(global+1, global, Testing::rand());
+            matrix.add(global, global+1, Testing::rand());
           }
         if (i>0)
           {
-            matrix.add(global-N, global, myrand());
-            matrix.add(global, global-N, myrand());
+            matrix.add(global-N, global, Testing::rand());
+            matrix.add(global, global-N, Testing::rand());
           }
         if (i<N-1)
           {
-            matrix.add(global+N, global, myrand());
-            matrix.add(global, global+N, myrand());
+            matrix.add(global+N, global, Testing::rand());
+            matrix.add(global, global+N, Testing::rand());
           }
       }
   matrix.compress (VectorOperation::add);
diff --git a/tests/tests.h b/tests/tests.h
index 471f68f559..6e10fe988b 100644
--- a/tests/tests.h
+++ b/tests/tests.h
@@ -36,7 +36,10 @@
 // Cygwin has a different implementation for rand() which causes many tests to fail.
 // This here is a reimplementation that gives the same sequence of numbers as a program
 // that uses rand() on a typical linux machine.
-int myrand(bool reseed=false, int seed=1) throw()
+// we put this into a namespace to not conflict with stdlib
+namespace Testing
+{
+int rand(bool reseed=false, int seed=1) throw()
 {
   static int r[32];
   static int k;
@@ -76,12 +79,11 @@ int myrand(bool reseed=false, int seed=1) throw()
 }
 
 // reseed our random number generator
-void mysrand(int seed) throw()
+void srand(int seed) throw()
 {
-  myrand(true, seed);
+  rand(true, seed);
+}
 }
-
-
 
 // given the name of a file, copy it to deallog
 // and then delete it
diff --git a/tests/trilinos/66.cc b/tests/trilinos/66.cc
index 764c49bd50..562adb36d0 100644
--- a/tests/trilinos/66.cc
+++ b/tests/trilinos/66.cc
@@ -41,7 +41,7 @@ void test (TrilinosWrappers::SparseMatrix &m)
     {
       if (i>=5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.set (i,i-5, s);
           norm_sqr += s*s;
           ++nnz;
@@ -49,13 +49,13 @@ void test (TrilinosWrappers::SparseMatrix &m)
 
       if (i<N-5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.set (i,i+5, s);
           norm_sqr += s*s;
           ++nnz;
         }
 
-      const double s = myrand();
+      const double s = Testing::rand();
       m.set (i,i,s);
       norm_sqr += s*s;
       ++nnz;
diff --git a/tests/trilinos/67.cc b/tests/trilinos/67.cc
index 23adcbb7bd..831139c5bf 100644
--- a/tests/trilinos/67.cc
+++ b/tests/trilinos/67.cc
@@ -41,7 +41,7 @@ void test (TrilinosWrappers::SparseMatrix &m)
     {
       if (i>=5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.set (i,i-5, s);
           norm_sqr += s*s;
           ++nnz;
@@ -49,13 +49,13 @@ void test (TrilinosWrappers::SparseMatrix &m)
 
       if (i<N-5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.set (i,i+5, s);
           norm_sqr += s*s;
           ++nnz;
         }
 
-      const double s = myrand();
+      const double s = Testing::rand();
       m.set (i,i,s);
       norm_sqr += s*s;
       ++nnz;
diff --git a/tests/trilinos/68.cc b/tests/trilinos/68.cc
index c0cb3860b4..54662507fd 100644
--- a/tests/trilinos/68.cc
+++ b/tests/trilinos/68.cc
@@ -41,7 +41,7 @@ void test (TrilinosWrappers::SparseMatrix &m)
     {
       if (i>=5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.set (i,i-5, s);
           norm_sqr += s*s;
           ++nnz;
@@ -49,13 +49,13 @@ void test (TrilinosWrappers::SparseMatrix &m)
 
       if (i<N-5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.set (i,i+5, s);
           norm_sqr += s*s;
           ++nnz;
         }
 
-      const double s = myrand();
+      const double s = Testing::rand();
       m.set (i,i,s);
       norm_sqr += s*s;
       ++nnz;
@@ -73,7 +73,7 @@ void test (TrilinosWrappers::SparseMatrix &m)
 
   // now remove the entries of row N/2. set
   // diagonal entries to rnd
-  const double rnd = myrand();
+  const double rnd = Testing::rand();
   for (unsigned int i=0; i<N; ++i)
     {
       const double s = m.el(N/2,i);
diff --git a/tests/trilinos/69.cc b/tests/trilinos/69.cc
index 6acb1f7c1b..d536357df7 100644
--- a/tests/trilinos/69.cc
+++ b/tests/trilinos/69.cc
@@ -41,7 +41,7 @@ void test (TrilinosWrappers::SparseMatrix &m)
     {
       if (i>=5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.set (i,i-5, s);
           norm_sqr += s*s;
           ++nnz;
@@ -49,13 +49,13 @@ void test (TrilinosWrappers::SparseMatrix &m)
 
       if (i<N-5)
         {
-          const double s = myrand();
+          const double s = Testing::rand();
           m.set (i,i+5, s);
           norm_sqr += s*s;
           ++nnz;
         }
 
-      const double s = myrand();
+      const double s = Testing::rand();
       m.set (i,i,s);
       norm_sqr += s*s;
       ++nnz;
@@ -73,7 +73,7 @@ void test (TrilinosWrappers::SparseMatrix &m)
 
   // now remove the entries of rows N/2 and
   // N/3. set diagonal entries to rnd
-  const double rnd = myrand();
+  const double rnd = Testing::rand();
   for (unsigned int i=0; i<N; ++i)
     {
       const double s = m.el(N/2,i);
diff --git a/tests/trilinos/slowness_04.cc b/tests/trilinos/slowness_04.cc
index e462b5154d..8e1760ca96 100644
--- a/tests/trilinos/slowness_04.cc
+++ b/tests/trilinos/slowness_04.cc
@@ -54,7 +54,7 @@ void test ()
       {
         // pick a random element among the
         // unused indices
-        const unsigned int k = myrand() % (N-i);
+        const unsigned int k = Testing::rand() % (N-i);
         permutation[i] = unused_indices[k];
 
         // then swap this used element to the
@@ -75,26 +75,26 @@ void test ()
         const unsigned int j=permutation[j_];
 
         const unsigned int global = i*N+j;
-        matrix.set(global, global, myrand());
+        matrix.set(global, global, Testing::rand());
         if (j>0)
           {
-            matrix.set(global-1, global, myrand());
-            matrix.set(global, global-1, myrand());
+            matrix.set(global-1, global, Testing::rand());
+            matrix.set(global, global-1, Testing::rand());
           }
         if (j<N-1)
           {
-            matrix.set(global+1, global, myrand());
-            matrix.set(global, global+1, myrand());
+            matrix.set(global+1, global, Testing::rand());
+            matrix.set(global, global+1, Testing::rand());
           }
         if (i>0)
           {
-            matrix.set(global-N, global, myrand());
-            matrix.set(global, global-N, myrand());
+            matrix.set(global-N, global, Testing::rand());
+            matrix.set(global, global-N, Testing::rand());
           }
         if (i<N-1)
           {
-            matrix.set(global+N, global, myrand());
-            matrix.set(global, global+N, myrand());
+            matrix.set(global+N, global, Testing::rand());
+            matrix.set(global, global+N, Testing::rand());
           }
       }
   matrix.compress (VectorOperation::insert);
-- 
2.39.5