]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
rework r31759 to my liking and fix a few more places
authorheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 22 Nov 2013 14:42:37 +0000 (14:42 +0000)
committerheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 22 Nov 2013 14:42:37 +0000 (14:42 +0000)
git-svn-id: https://svn.dealii.org/trunk@31761 0785d39b-7218-0410-832d-ea1e28bc413d

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

index 9429c1d5b7a55727dac7fcfff9adb200dc6a4c70..35f33a1635261dd48c2c8bfc9a0f25f78ed3454e 100644 (file)
@@ -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
index be0f76b27a2d057acac90f6b94c7116608635c31..85f1d215ddda335a9fe9126228e83fc14a560202 100644 (file)
@@ -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
index 3f45efeedadf547e47ce77a1744cdda24d08d44f..5d9d4598b6384df1d869452d57a236e36b921f51 100644 (file)
@@ -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)
index ca7e276a883ac9c79e54641a7c82da953481fbb6..89938fb79cd19b05efd916e916dcb0699d60c20a 100644 (file)
@@ -27,7 +27,7 @@
 
 double rand_2 ()
 {
-  return 1.*myrand()/RAND_MAX*4-2.;
+  return 1.*Testing::rand()/RAND_MAX*4-2.;
 }
 
 
index 3f794dcf827078b01b221da30317577c597fee1f..a69ac9dfb676a327efe20eae4f47eec2276dbcb9 100644 (file)
@@ -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;
index 21c66c48fe36e5ecddc7143834c061b83db786cd..3c05906c657f9c2da2ee74526c86b65e74455163 100644 (file)
@@ -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;
index cd5422a7b24141dd3807c4525010552a877edec2..a495c13d7a12140b9f1a13a690c966ae40c2a442 100644 (file)
@@ -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();
index 5c101775fb947fa0b732e4ed36342d02a21207ba..1b837743edc625ea2d8d0b257be9d0e331537e23 100644 (file)
@@ -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();
 
index 04476b4d71a992818e184b5d51a63b68b14319a5..8f5664acd98b8725b33091b1bfd480b81ec259d3 100644 (file)
@@ -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);
index 7e985b391d5cbb12e137cf6bb29ee3831891cdb5..0a3021a70b599bcf86bca9413f996196ad6eb737 100644 (file)
@@ -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);
     }
index ea0ffe37b7cdf4bbc9d0fd7f5392574a28481721..906d4eb1e1ce1e5900e173c0bc8866e0d8446f9f 100644 (file)
@@ -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);
index 4a4d662f23477e8a2d18cd5529ec904132010d6e..75898a59d3dce0d751cec65c499d35881ef109a4 100644 (file)
@@ -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)))
index 9899f14182aeada391661249339fe82eec14e353..0f55292bc5aacee5c0df7c6217c038224cbabca5 100644 (file)
@@ -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)
         {
index 47d3aced5310b4665eace5a1b46c7077adefadcd..d11e59a83288110c9f3b9acc7efa10056228dfe9 100644 (file)
@@ -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;
index 64d7b96a60ef663b0100cce21a2cd4398ec6643b..fe782feb38434340bc49a892109e26858cccab15 100644 (file)
@@ -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:
index 46026f214f88c310fe8b3673d7d8f8ccee51c659..e9bf54a9ce9fa96d3d78487fd07b2d557ab664a1 100644 (file)
@@ -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;
index 7a15ffa345e63513f504646cf31bee0ea42c002f..b7ed88f2a95729a90f36ace329cd410370ff8345 100644 (file)
@@ -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;
index d938f327da054069e6867bf9288b7099f8093f59..137861362b6d6c53baea2a3f57f579e2094deee1 100644 (file)
@@ -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;
index 4e8262c5f24bae92ca0f83cc4aef4b2d919ff4d0..cd0236a045b3d36ff3fa84a6f44fb0167cdd6e1b 100644 (file)
@@ -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,
index 662dc7a3941f44e9d41239ab8692cc3a4a9ba5f7..95341c21a838338b3cfe2f68c7aca5a563ca036b 100644 (file)
@@ -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();
           }
 
index b5c0546b7c2b3388ff65529c84af9e5f52f24ec9..daf064bfaa463d9c86185004d6101b8a69759104 100644 (file)
@@ -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;
       }
index d20bae63b0ca081b0abc40d8b188e4f14d2d0fa3..8c5b45a5aa5396ba496b99c565f3957c54da990c 100644 (file)
@@ -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);
index 683840a2b396d26ba00b20390519bc455cc2f325..b99d1564ede0572deb8be5f46be6c2a040b6199c 100644 (file)
@@ -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);
index 6d27e90f5114dde4b2910820c54cfe4451a971ec..9b7d41d7cd22b74a68a51310bb91343ef3d77334 100644 (file)
@@ -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);
index 2a41f819903f6458f8dde77e5f46dee7f0461988..308b5d3b23fa55ece686c943ac66e1a0558baab0 100644 (file)
@@ -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
index b74fac85e64475899c63d10f939bdf8cdd85bf81..31b7279a8ff4fee0b1fac526e20550e795f43f29 100644 (file)
@@ -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.);
 }
 
index 7454364b966110c17fa1a5c3dbe57e81133364d6..df24077c2bd9d8573bffcd69f4b33ed8201a18ea 100644 (file)
@@ -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;
index 7f2df7c629576f235b96523b9e2c0a4a3f35a1ed..35774e8473565e3c321dbe116e98a41e1a172037 100644 (file)
@@ -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;
index 3370131cf43650d246f67eb14a9ade7f050ddc1c..ce14af6b0f72231f9eda72b1511861658230a081 100644 (file)
@@ -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);
 
index 22d806f83682a8fea896c9c22809ab82756d3158..a171d265b82038093a4af1b6e17ee4bcce018914 100644 (file)
@@ -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);
index 5d2bfa6b6d21270d4b2a7246099a40c57d491af1..b9ad731dac9f399db15c95d78c8a5d36f504b663 100644 (file)
@@ -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);
index 5e8ee38ae3264100aeefea108730fa215af0c807..72744cd7ba752343a72399af0b98ea624c037084 100644 (file)
@@ -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;
index 463c8bbd85c69949ac41f45714847ced4c7821c6..16cf0da912627818b683b85fcac1ca8b3c5eb0c8 100644 (file)
@@ -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;
index 78adf91af32a668de08f6d695cc61fe78907e975..9c70b623a43d80f4d6801400607d85830b843f46 100644 (file)
@@ -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;
index cccd24abaf4fb57c4fa7a56ccba4eea9def4faca..f41f35dccbe79eb0e86677c0639e1254b6fb3cd3 100644 (file)
@@ -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;
index 4359953f1849f82c860811a617c66ad9a4e09313..ddc62a6f8e17545500435f5a238ffd6676298cbe 100644 (file)
@@ -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;
index 82a73c0394ff1bb18bbde3b4bb23e693fdedd324..6a59bfe15516f3412b64fd12c774d248ed3e64c6 100644 (file)
@@ -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;
       }
index e9b78a19121d5209355ee4d9200f4c2ea26b68b0..b4079a26215a158bc38b36fde8eb0cf6b9e206bb 100644 (file)
@@ -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);
index 83a8ad1f1a37c3b2fd6bf9d228ea483351180ca0..5c1ac4d15a1e0123bd0e7ee2b5105cf9bc5ab7f8 100644 (file)
@@ -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;
         }
index a20ab8393ae3d2c6811c6da92db21657a447b13f..d0995f5a345caa6c7658d4e7b1a4a5fee3b0e974 100644 (file)
@@ -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;
index 794e1f35422dcc923bea80995c8624979112d2ae..c97427bb3075df71ab052e391d4973615045ac65 100644 (file)
@@ -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;
index e366a3e87aaa0bf1fcc4f41412456682566af64e..724f5975d3e9b710f5a8bc2130aae9797c4609bd 100644 (file)
@@ -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;
index 4b213982cbf9a89c5722ba69a861662dcd3c2bed..bea89030308f336133c557984efb27e856057335 100644 (file)
@@ -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;
index 8a50d220c4c78ac51c85f24a6c1b3734d7ce0c2d..b9e101fc60d1243074216628c0aaae4c18c7afc4 100644 (file)
@@ -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;
       }
index 2b41275641402d61e361a210f9be8081d42c0285..a38f8d9925910ef74c4243741239a901816ad3f4 100644 (file)
@@ -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);
index 03b7552abb890d6752ed3e827cf36c043f717732..c8536fd96d076d196fbcf509f7e606c079965a73 100644 (file)
@@ -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;
         }
index 36ac26f8e212664d745784400659f56dbc246160..86d6515fc3b37eac5cd715a201eded0d4c47226c 100644 (file)
@@ -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;
         }
 
index 22c373e0d39b4c71e338b4614ce1852ddb6b038e..99ae577644ff1322a6466ab5fbc3f1d0a0144413 100644 (file)
@@ -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;
index 1aabc661121b62a31de2204d7fffde449a5a2ebf..846b097fea2dcf120d97bbe8849832badab0cde3 100644 (file)
@@ -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;
index a7f382835bb6267ba7c3eecd7e6a9aad694440f5..d40d269cc219ff4d25aae459b5f4d4dbb0cc73b2 100644 (file)
@@ -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();
index d98354cfcb1c9ae098cd2e95c6c41074e6bbb733..6419a1bc72cc86cd46bc21a967a5718c0e7db6d7 100644 (file)
@@ -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);
       }
 
index 6d050878a8862a2cd60d757f18e4c29c7c9b6419..a76e36533125f8b686f74bfa990831df856efc45 100644 (file)
@@ -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.);
 }
 
index 590e3b71973de07827ea27b38f1542dab7eca607..8dc62a8b24e4a41eb551b7e97c7bbae03e6f5865 100644 (file)
@@ -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);
 
index 577b2870ffc7bca13d9b2268916104870368d927..71712cddc5a9bccaaac277ae68f0705ea9ff3a55 100644 (file)
@@ -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);
 
index a4526be072085bd6d18834bfa813af0d5860704e..fa2c573579eaf0728877f5a3308214fa8587cce3 100644 (file)
@@ -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
index 5f2e30610d1747c5afe03da2b3a992d4297e1f56..28cc4c9ea66776eb178c69b44266901553fc95c6 100644 (file)
@@ -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
index 9218074434b60c3493b00794aa6994783c83aa6d..eb4a6af7876da9b0197ee3f0b7e1c452cecc1785 100644 (file)
@@ -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);
       }
 
index 1b6cc0c58a324f7c1692b39f3dd770dc0bc26370..5b18693c01050c40a084a8cce63363dfb859564c 100644 (file)
@@ -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.
index 125f6de7ff78ecde521bdd2e3f759cf155ceab4b..e384875032261d00f40542ba5298188b2c1c1c7e 100644 (file)
@@ -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);
 
index a5adcfa9982300db2db8fb97bafed83def5f6b3f..1fbbdb4b28d3371a78a71664d71cc808ebf53fb6 100644 (file)
@@ -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);
 
index b688e5db7f5391a318737d7abb3f6b2ca7ada6b6..21f25fcd9904a584a59a2d8ea10ae82f0ef5bcce 100644 (file)
@@ -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;
index a868290e1c7eedbf192690eb45750946cc585ff3..be75be61a2f4f3bcc7a2b9f9b445bfdf50684fdc 100644 (file)
@@ -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;
index 230627f4fe86b0778c85d0b4a9dc2ccaf6edb838..e97612c11e6bc57ab445a36b8bfb9ade46f6de36 100644 (file)
@@ -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
index 9302e2e068be7e4e686640f7ca5cf69aa21c1625..1bb40cd3b53a1a36a5f473b0997d1555d9486a72 100644 (file)
@@ -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);
 
index 362d2fa7301ba9f57b73139226ca1b15ba592b23..d44f0d5ba01f0e84de07ce4a8bd481567a6c72d1 100644 (file)
@@ -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);
 
index b4bc3c017bf48f6145999ce89f5afe229cca72fd..58f755f8db367b7b6d8414ef0767c12f1a968afc 100644 (file)
@@ -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);
 
index dfcf4b1bc9524ecc8af419f9ce6818811eb14dad..753c66f9ab85a23dcbaa6fb9c6c04646b2c11918 100644 (file)
@@ -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);
 
index 485f897043551aa070aaaa13e7c92b5cfaaa7ffd..d47c37e1434dc9fc3bdc1521d1eb2c4b36a62caf 100644 (file)
@@ -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);
 
index 45701b331ae78b923aedb7bad6ccb68a5a6efbbd..b337cf4b63b00291dc626fc3f203933b36b2fb6c 100644 (file)
@@ -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);
 
index bd28208994238c1cfc1178c570db6585190e7549..c9c6ba4e00dcdc9431493b08e4c12d890b52e788 100644 (file)
@@ -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);
 
index 7ffe72beafc2add560da0f7b3aeca120cea01bb4..ff851a29ae8e02e6bb592d3cb64813ae5fbf4f75 100644 (file)
@@ -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);
 
index 99b4c53635a8a7d6ab9f9a1414783e46d76a92d4..a1922a1656e5a1aab40a5d47dd36a2a1714d8551 100644 (file)
@@ -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);
 
index 5c9a084566cfb4173cccf0f0391c2363d7703170..156517288ac17d5fbc45ab07dcb3e46845469562 100644 (file)
@@ -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);
 
index 0abe9282307b95bfc9a0642bc7f7e938f331c75e..4ce66bdcddd93d31bf87a5b18b5b93c7a0c883ac 100644 (file)
@@ -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);
 
index 0a574b30011cde1e5eee27a8fb4f90e59537109c..5abac27879fd97f1863eff3e426b8742d1ed7075 100644 (file)
@@ -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);
 
index 8ea98aa19ddbc9a389928906ebf419cb70f33fb5..e1f6550afd1c5e7be27638452e1c3e8e1ae684c8 100644 (file)
@@ -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);
 
index 32102a9a76568d0346ffec77764b6defd1abc6b7..c3372109f213f9933a77b1503c1ed61271ace37b 100644 (file)
@@ -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);
 
index 1efcbfcb9ef7e8d3274cb1543cca3e158aeae0ec..8825edb228b225d67ceb3bacbe70f673d044625e 100644 (file)
@@ -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);
 
index c53a90d16955c0c2b547d3ddc21f6124fb61491e..f7e86e8d39b4d885d79b76d75c28414424ccaf31 100644 (file)
@@ -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);
 
index 5b1016a6ac6b8e5ebba25e14f35915940281ec41..4a018f3320bd217062b5b9da6502681d6b3fde75 100644 (file)
@@ -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);
 
index deb080e399720584bea8444a593e0a38be8d48d4..78f6d522fb78c5912f05637b1bd71c048263b944 100644 (file)
@@ -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);
 
index baea9e0be39bdb174ecf8c9298e539abcaaf9fa2..81e37d6de5d5e01d297a58f74d246b1d80799c3e 100644 (file)
@@ -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);
 
index 8df741b04753f887f33ba32950f97eb153adefcc..6d3280ead3ce9ca93c532613e02fbb39e527f4c1 100644 (file)
@@ -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);
 
index d073d58b4e0b933f22bddf1a8c89159a6c58c380..55e5fcb46f543ce3cd2736dffadc343cd44418f9 100644 (file)
@@ -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);
 
index f1c2fe5df382d8c4ad01f83eacf0ee3da5d6d36a..21ae5256785309df8c45dfac106f2b07b8181158 100644 (file)
@@ -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;
index 262c077619b3b4b171d888b4695a6e34552082c5..ae2090693e080b5afdae98247a06a934c4d682ec 100644 (file)
@@ -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);
 
index 6b580b9adbd31d09f30ed5e8cc7ff5051281e21d..818ca03e5391d152c9c2d8c6662faf5ba76dac0d 100644 (file)
@@ -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 ();
     }
 
index 959e727aabaffaa43cc115018a13616bddd2b3bd..f94dd9d9f2ffd3a47b4ae7b816d23929e1b15653 100644 (file)
@@ -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
index 91235de444098eeaf1c6c6c18882bbc4dba80581..e0e8ae91206e1df42059b8daa13191e0c9353d13 100644 (file)
@@ -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);
 
index e0ebb1c5af7d6684721303a21e08c39bfbea78f3..87770b6bd54d7f6063191d685960d7868de939c5 100644 (file)
@@ -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);
 
index 315dd670402521e7507293cff0effe74891eeba0..81eeae01e1780a0e14d0ea6a21cf5f603e239912 100644 (file)
@@ -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);
 
index 28f54d287caf49b5f9f84b4fd3aa757cdf47d95d..b8d536fe8a95a7d215349b15cc709e6403013297 100644 (file)
@@ -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);
 
index e747cd5fbe545858e08cb484d2cd663de6a5a28a..d3f3a3b2d8ffa4148ba8ce50a2e1e11e51bd9d42 100644 (file)
@@ -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);
 
index 64b94d10f345b47abc8d9d655b08de800307efa0..ea0c2e820564461df4b2941f758fbd9305a7dd27 100644 (file)
@@ -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);
       }
   }
 
index 68a219c55b82b120dc46fe75e1c829ca98b8e518..9881bfac1ba64750df6d0eb2d3b4b508e31917ad 100644 (file)
@@ -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;
index 716104031500695af73b9ebb6c816c0a9965232e..162034bc7c0512aa528e07e0b6b4182186e103d0 100644 (file)
@@ -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;
index 8476ea4afd550723b47821c65d964cff912f77c4..117e29889b9427e443d4d22d42bc7eaf24f596e8 100644 (file)
@@ -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);
 
index d7be80fb7d5aa1ca4da1c29cca76146766f3953d..bbfb624c29b93b8f380a51f4470e1f7710ca5052 100644 (file)
@@ -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);
 
index 79ccf747f1596e46b4ed64d701518b9707557166..d7fbfad6464662a73f612f0226168f13c6228c57 100644 (file)
@@ -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);
 
index 131acc3168b448d1c9cd2058490063d38f7f796a..a97062db2111a472ef95e6bb261ab655e630db14 100644 (file)
@@ -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);
 
index 336e53f4a50e46c3523ffa55e9c5353324c23d0f..ec379d42f4112255a952b055ef103a545c9a64e0 100644 (file)
@@ -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);
 
index 85005a1a87920ae2a0f5f5480b0c36a6d2e7c6fa..abfe4aeae0d7e1126cd6466bf9560fc967cddde7 100644 (file)
@@ -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);
 
index 29720607b1552523331376d485f1fc9d624fa843..1b45bc34dc76740fa1016866909f5cd82f3616f2 100644 (file)
@@ -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);
 
index 398549e9fe21384b9738b67b473614e8bedd5377..7cb5704ecb118f4af4be4cfc75bafb36fd61a038 100644 (file)
@@ -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);
 
index e703fac95af0e6a2f5695be6b1430a4765a143f6..2f5d15b73e763e1c465290f86ced068fe83ade7f 100644 (file)
@@ -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);
 
index a27d07c4447e3c0df611040e89b0e8b60885a610..52b21c16db2ae090369755b38dca2309cb7bb164 100644 (file)
@@ -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);
 
index 3a4f215c83cbbe2862a8520e90c05188438aac1c..2dee753110feda9df3e351042cab1573c6d416b1 100644 (file)
@@ -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);
 
index 40d54f393b9b8bbd0b6e7377cd78492bda27d3a2..31dab39f128e80599ca967404a41f6bc45ee32c5 100644 (file)
@@ -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);
 
index 0876f633d74fa13858f7d40663e9866f2d79723c..507d7deb65af9e830d116fe126fd21cad4af0ced 100644 (file)
@@ -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;
       }
index 9a7e1d269b874c27759214f2475e12ca91a71b6c..9de8c6159b8695756ae1cbf5df97901309062445 100644 (file)
@@ -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.
index 49fe200b8a5b96ac9d4a7282f7b2645a83db2e64..6968b5d5fa55f9eee21f4137a2c58d157eb66158 100644 (file)
@@ -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
index b6830f91a428a5e2e1f832d340d6c04d9db0e172..26dd1b61e5280836b5528b97ed7153bd083313bf 100644 (file)
@@ -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
index 23982d6cd85c720694d6d34bce6e44257be9be2a..bb4b86a5b91cf0d2c5a18ee0fc56da29fee6c757 100644 (file)
@@ -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
index f91a5dddbbafd1b956bd0e586cab0d53476ac668..126c8b5a4aa3c95ec9cdb6095efae4a456244974 100644 (file)
@@ -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);
index 8044a18a964528d7a86d5b146dbc8a5dfa3df988..453f71744c0aa0f84fe79ad46206f2e3a0214012 100644 (file)
@@ -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);
index c9db6880d4572633ca4eaec52b4b790efdc5ec21..823e1cc3094e0bfa23bc4d98ec2c18f282e1b876 100644 (file)
@@ -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);
index 83aff6f4e809d5374a4dfb3fbf2abef479f1b9c7..c639bff45f067edbb675a448b98b4436387b5a91 100644 (file)
@@ -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);
index de2070609e4be34aeefb95842718551baa6b5416..20ac941be9c1e3eb1276f54c4f35bc7cd920ca57 100644 (file)
@@ -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);
index fc24fdff979615ce1c141d95c3b918ea57120482..f94e783fc327780cc3deb33b441db31b523c05ae 100644 (file)
@@ -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);
index f1a03a00da317a488b67166aa62ff4356d810f5f..9dacf4b804460f9cba4145870ac274b8f4cee363 100644 (file)
@@ -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);
index 0c5cd1fafe44006486333fc1550db037ada19c11..f2a65dd581e57de1c23067c97ebfa2d08d73a37e 100644 (file)
@@ -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());
index 6eebe79bd1ba5d3a240a0edb43ef457c73e0d4a1..1ebee2cf6e9178b867a5b0e238535cd2cebe6b27 100644 (file)
@@ -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;
         }
 
index d908646d282c485a1243b54b341d20343f37c922..44cc266fe552fdd28b0e7a914b27aeb07e12f0ab 100644 (file)
@@ -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)
index a4b78a84fbc60e7fa18096322564b69b8c92b697..4db1704541c6507004f1d37d50d804b15b6185ff 100644 (file)
@@ -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)
index 6e23e6e2ba142c6162fd13dfcc6275f8db77f8f2..aba500ec682de61828979c2a395f2829057a4ce6 100644 (file)
@@ -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;
     }
 
index af8131c811d505bd3b6b9ae9fb4ed482e8cdb9fb..c585c8e235ab7509df5e1a8478acc2674ad08f48 100644 (file)
@@ -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);
index 20e8195b6a8edcca508fa16b3638f27149f4d59b..27e1179631cf448e52a75f1f219518e07fbbd939 100644 (file)
@@ -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;
     }
 
index 8b0c421dd279a051a572d82a53109da0c45a1398..a6d8743e5c1451ed67350dafab9b046cd8105b73 100644 (file)
@@ -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;
       }
 
index 289f54dc6aaf2599f4deba3fb32676106d3c8001..b4e8cd81f58cf98b8dc2eeb080385b0f2358ea31 100644 (file)
@@ -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;
       }
 
index e360aba21d6bbd8ffd3e535a2ad7fea8ffa0a2a5..ce3c796e2093548d8cf5e3693170c95915dff11f 100644 (file)
@@ -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);
index d99ffd984231b3add6116f190635a7e8433fae8f..e31f4a5a798b96b1d977d55e9183edb90acc6eaa 100644 (file)
@@ -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);
index 27665233d6bc29ef396a13161bc1e8fcabdfe714..5a95affffb0fefedd0aca8122f69ccabe5c1bd47 100644 (file)
@@ -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;
     }
 
index 4da355bb68f771f71b201dc365a5023394161dc2..3cb6681f43e9e495aa03edac9c2d249a3ed71e37 100644 (file)
@@ -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;
     }
 
index 096b293c6d0d4428c718655c3e95a5db94efe655..f42acc2e60f77cd33bcafa8fba2de1091394b603 100644 (file)
@@ -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);
index 904e9bdbf86695f8eceea29a1cf9b122812bacba..bef032004836834ddf091f7ed8e1c0ce61f47f69 100644 (file)
@@ -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)
             {
index 7a182c08086f32350084117734a3960e8086a548..804f875952e8f020c56f2dc9455a0f539f349624 100644 (file)
@@ -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)
             {
index 6be8ed73219e12d78afd9a0f488681151f459ffd..b460be662d6fac395340dca46bd981534d78521c 100644 (file)
@@ -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);
index daedda39233fbb8e9adbb9dee1b85c8ae0cd52af..1d111f2a4b9ad193f0b94610d330d54a898df2cf 100644 (file)
@@ -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);
index 60b18c3a94242389d27f22d7461812ee236a3278..0c497fe15bbfc39bb6781e47ba8b975c8aebe42a 100644 (file)
@@ -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);
index dcbf078bd56751aa7fe8c9a67800d0bbbde594b1..c777c8c79cce67a6cb28ca4cbf1692cc4917961c 100644 (file)
@@ -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);
index f1c5b3b7960debdc355924591cce8ff01f361194..40e715a5acadc8028246353eff6522bce7be58bc 100644 (file)
@@ -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;
     }
index 36fcc82a4589615154c779cdf549305d89284775..773ae0fc43035512ca76245060949cf5b95a2fe2 100644 (file)
@@ -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);
index c4af939d566221fe8c7b1c26df8dd07bbcde2b73..736e9a8c81167467ebf34d696ed9fc0bb2d0d691 100644 (file)
@@ -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);
index 19b31a5fe181846802466d689a0b80d02422a094..b359397e68def8533b1c68de832e533dbbf7c6db 100644 (file)
@@ -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
index 1bbc885654718d708aad90d219749d29a8b784cd..4eb2670a22d668a39165ce05dff01b7c28d2e631 100644 (file)
@@ -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
index fd8d549a3dfe4d98132d3c413cb58e64ab381972..46e1463c3ac94b21210ccf540552e15f35b0e0d8 100644 (file)
@@ -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;
       }
index 3d01a123e87dc63f6dc374fd95c17c65061d0f58..1860f6797ef216771bfd45ce0a457acaaceb8ee4 100644 (file)
@@ -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;
       }
 
index 4b4846e7c744de7ef0676314b58af6404ecbd837..2c29799abb151c39da1f614aca79cf4c0634c539 100644 (file)
@@ -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;
     }
 
index a108ee375827b15b6295760f10dda1e16e34d6cd..edbd65c2accf32e78875ab29f3aa3cc8346217e2 100644 (file)
@@ -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;
         }
 
index 6c0130e5738c139e14eb6db0fb18bab6d32cab8b..1cb0a8cb88398490a35ba767541fd2cdadb034ff 100644 (file)
@@ -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
index b24169b7ca30989c1c34483a4ae4031bfdc90036..d24f62df7d2393b87201b0e12c8b41bf472e7ec4 100644 (file)
@@ -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;
     }
 
index 4a11e07d39c34a7781881fbaded86d5b57259f6c..86b61ee225d32fac3a5c0887f8211733ba93e361 100644 (file)
@@ -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 ();
index dbc6c47502aeeee7dfbc7cebbe622f6a02bc5a5a..958357ff26253d47449666160ce25ed717d21c10 100644 (file)
@@ -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;
index cea95adca7642a1078289a804752194207358563..7cac982c0b45c52dcb2ebc0bafcd71578ded2778 100644 (file)
@@ -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;
index 9a07498f78df44096fba2bf171ea56ccf85d2262..58d7f29536711ae22eb7a8a41faeeafcbe042279 100644 (file)
@@ -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;
index 7cb6335de538fc32e6db83e4014204f94121bae3..fc90013800995a0ce50da4f7ff11d89d6c4549ab 100644 (file)
@@ -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 ();
index 7e342327ea3b7876c8dfb75435a9855fe889df48..241311693100a117e69d6d22b087289dcfac8597 100644 (file)
@@ -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 ();
index 77a2878c97a4dcb708bc5da383bd12eed99301e8..420c8020212ef1802425a7a6e21fd342edc3f48a 100644 (file)
@@ -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 ();
index eca70837167cd3e27818e9329fd09c7c26887dd3..ccbd94751879ad4a29d2d3bf43bd98d5b9d77f7e 100644 (file)
@@ -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 ();
             }
       }
index 4de08b4f54736f248cc6fc2d58233137303bd4d3..d2eab70e0d3313612d2c3ea1c67056351222019a 100644 (file)
@@ -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;
             }
 
index 4e89fe123a2703bd0bd9ff61a1d8440855fe22a9..d60c5dc1ddc787392bc7553863e9137e60700eff 100644 (file)
@@ -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;
             }
 
index 93ebb774090c3432eba71611f98e72910a0c2880..e4d50ce573d25377c2854d98272dac7b82dbf243 100644 (file)
@@ -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;
             }
 
index c1b1acf6cfdbe404a10777a408cae9c3cd24cdd9..672b8d3d6db049d3d253fef169b1d589dcd524cb 100644 (file)
@@ -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;
             }
 
index 23261eb4f69db24f6bbca2b3271bf628bc2a107d..3cf00b8c8e7e0f63427fb040a54531bab8e8bd6d 100644 (file)
@@ -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 ();
index 7745179056e8b615c0743aa9927ae9f8ddaea98d..6299420a9a78924f9d63d3ef665312b3e55d67dc 100644 (file)
@@ -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 ();
index 1b47188f1b81fc73424599547a0f9556b5c3eb3b..c0b1b41b5505e9ee0049bc7c624b7f82679aef26 100644 (file)
@@ -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 ();
             }
       }
index 84eb720cd6f45a903f72dd2c87c11e1be3c78512..05f499e7114ca0bd8539da60f8df022cb81e7886 100644 (file)
@@ -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;
             }
 
index 94ea3e4a447582f893a3ec9d23904e7cc4d842b2..4ded88c3384e349d2845105077642d8cfe024d4d 100644 (file)
@@ -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.;
index 6c601200e6684e106e68357db7c6909c54343607..ef3e17b1552c03f2cfce5579b649b1f728b8065b 100644 (file)
@@ -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.;
index 7649ecb009a3d31ba11ca72d056b9f65430e5263..bd92dc2c9c8cc23adebb8c64835d0b9904d66f16 100644 (file)
@@ -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.;
index 68269260279c726f91cd34860e10bed6eb4f7d99..0c950d7cd1d546dc5a8a1d4ce8c6c24e6e859c05 100644 (file)
@@ -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;
index 34febf212669bc6b01379cb4cb7405eb8689e711..3f5ce10a7176654d6f8f9eb6f19481e84c2a51d3 100644 (file)
@@ -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;
index 802565b7c74fd942c502a9e45d3e2b453c4de4c1..d33074117c2ad39411da82d7a83325dc0d492ffd 100644 (file)
@@ -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);
index 1d4e91eca8eab70e99d6f2fceec8a1d755de93be..5ec674841ce85f7f25a0616d88088c67b27c0efa 100644 (file)
@@ -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);
index 0b75bfa0488e3fd6db8634e19e44115d4182fe17..0a616f8d4d0fa9cab3f5f968d60e27a2b18663b3 100644 (file)
@@ -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);
index 471f68f55942366f5833d2538e34c872e91307f6..6e10fe988b5f3e58c6c25aa7543dbe2e8deb2f0c 100644 (file)
 // 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
index 764c49bd5084434271c340664b1d6d1ea205aa9c..562adb36d05d350fe10ae59680930f5b43e1d7fe 100644 (file)
@@ -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;
index 23adcbb7bd6b1cd00216ffe6f58429b12e23f75e..831139c5bffb0af16b58eb6a0da1cf956a9370bd 100644 (file)
@@ -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;
index c0cb3860b41741710c8db9602f0d74e817064f4c..54662507fd2568d545f4f53dc7d35e7290e8a290 100644 (file)
@@ -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);
index 6acb1f7c1b86d3a3137635ccc2173921b53a1751..d536357df789b09bbe7622f74c96527c8e903dc6 100644 (file)
@@ -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);
index e462b5154d4e1a586c2f5a5c0865eadd2b0744da..8e1760ca9691afaf8d3c956a46582beaa7baa9b9 100644 (file)
@@ -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);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.