]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Added an MPI_InitFinalize to the files, so that at least some tests work now.
authorkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 6 Nov 2008 10:21:30 +0000 (10:21 +0000)
committerkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 6 Nov 2008 10:21:30 +0000 (10:21 +0000)
git-svn-id: https://svn.dealii.org/trunk@17484 0785d39b-7218-0410-832d-ea1e28bc413d

111 files changed:
tests/trilinos/01.cc
tests/trilinos/02.cc
tests/trilinos/03.cc
tests/trilinos/03a.cc
tests/trilinos/03b.cc
tests/trilinos/04.cc
tests/trilinos/05.cc
tests/trilinos/06.cc
tests/trilinos/07.cc
tests/trilinos/08.cc
tests/trilinos/09.cc
tests/trilinos/10.cc
tests/trilinos/11.cc
tests/trilinos/12.cc
tests/trilinos/13.cc
tests/trilinos/14.cc
tests/trilinos/15.cc
tests/trilinos/16.cc
tests/trilinos/17.cc
tests/trilinos/18.cc
tests/trilinos/19.cc
tests/trilinos/20.cc
tests/trilinos/21.cc
tests/trilinos/22.cc
tests/trilinos/23.cc
tests/trilinos/24.cc
tests/trilinos/25.cc
tests/trilinos/26.cc
tests/trilinos/27.cc
tests/trilinos/28.cc
tests/trilinos/29.cc
tests/trilinos/30.cc
tests/trilinos/31.cc
tests/trilinos/32.cc
tests/trilinos/33.cc
tests/trilinos/34.cc
tests/trilinos/35.cc
tests/trilinos/36.cc
tests/trilinos/37.cc
tests/trilinos/38.cc
tests/trilinos/39.cc
tests/trilinos/40.cc
tests/trilinos/41.cc
tests/trilinos/42.cc
tests/trilinos/43.cc
tests/trilinos/44.cc
tests/trilinos/45.cc
tests/trilinos/46.cc
tests/trilinos/47.cc
tests/trilinos/48.cc
tests/trilinos/49.cc
tests/trilinos/50.cc
tests/trilinos/51.cc
tests/trilinos/52.cc
tests/trilinos/53.cc
tests/trilinos/54.cc
tests/trilinos/55.cc
tests/trilinos/56.cc
tests/trilinos/57.cc
tests/trilinos/58.cc
tests/trilinos/59.cc
tests/trilinos/60.cc
tests/trilinos/61.cc
tests/trilinos/62.cc
tests/trilinos/63.cc
tests/trilinos/64.cc
tests/trilinos/65.cc
tests/trilinos/66.cc
tests/trilinos/67.cc
tests/trilinos/68.cc
tests/trilinos/69.cc
tests/trilinos/70.cc
tests/trilinos/block_vector_iterator_01.cc
tests/trilinos/block_vector_iterator_02.cc
tests/trilinos/block_vector_iterator_03.cc
tests/trilinos/deal_solver_01.cc
tests/trilinos/deal_solver_02.cc
tests/trilinos/deal_solver_03.cc
tests/trilinos/deal_solver_04.cc
tests/trilinos/deal_solver_05.cc
tests/trilinos/deal_solver_06.cc
tests/trilinos/parallel_sparse_matrix_01.cc
tests/trilinos/slowness_01.cc
tests/trilinos/slowness_02.cc
tests/trilinos/slowness_03.cc
tests/trilinos/slowness_04.cc
tests/trilinos/solver_01.cc
tests/trilinos/solver_02.cc
tests/trilinos/solver_03.cc
tests/trilinos/solver_04.cc
tests/trilinos/solver_05.cc
tests/trilinos/solver_06.cc
tests/trilinos/solver_07.cc
tests/trilinos/solver_08.cc
tests/trilinos/solver_09.cc
tests/trilinos/solver_10.cc
tests/trilinos/solver_11.cc
tests/trilinos/sparse_matrix_iterator_01.cc
tests/trilinos/sparse_matrix_vector_01.cc
tests/trilinos/sparse_matrix_vector_02.cc
tests/trilinos/sparse_matrix_vector_03.cc
tests/trilinos/sparse_matrix_vector_04.cc
tests/trilinos/sparse_matrix_vector_05.cc
tests/trilinos/sparse_matrix_vector_06.cc
tests/trilinos/sparse_matrix_vector_07.cc
tests/trilinos/vector_assign_01.cc
tests/trilinos/vector_assign_02.cc
tests/trilinos/vector_equality_1.cc
tests/trilinos/vector_equality_2.cc
tests/trilinos/vector_equality_3.cc
tests/trilinos/vector_equality_4.cc

index 02a0f5df50b3882d6c8ff861e849a4989c0211d7..8c5ac4294ed77b47de24b1f68b149a6722d2ca9d 100644 (file)
 //----------------------------  trilinos_01.cc  ---------------------------
 
 
-// check setting elements in a petsc matrix using
+// check setting elements in a trilinos matrix using
 // TrilinosWrappers::SparseMatrix::set()
 
-#include "../tests.h"
-#include <lac/trilinos_sparse_matrix.h>    
+#include "../tests.h" 
+#include <base/utilities.h>
+#include <lac/trilinos_sparse_matrix.h>
 #include <fstream>
 #include <iostream>
 
@@ -30,7 +31,7 @@ void test (TrilinosWrappers::SparseMatrix &m)
         m.set (i,j, i*j*.5+.5);
 
   m.compress ();
-  
+
                                    // then make sure we retrieve the same ones
   for (unsigned int i=0; i<m.m(); ++i)
     for (unsigned int j=0; j<m.m(); ++j)
@@ -41,7 +42,6 @@ void test (TrilinosWrappers::SparseMatrix &m)
         }
       else
         {
-          Assert (m(i,j) == 0, ExcInternalError());
           Assert (m.el(i,j) == 0, ExcInternalError());
         }
 
@@ -55,7 +55,9 @@ int main (int argc,char **argv)
   std::ofstream logfile("01/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
 
   try
     {
index cc6d0949b939a51e135d824aeb1fa944468215fa..504c39b2a5231bd9fd61e69108a7bdd0e9b93705 100644 (file)
@@ -15,7 +15,8 @@
 // check setting elements in a petsc matrix using
 // TrilinosWrappers::SparseMatrix::add()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -41,7 +42,6 @@ void test (TrilinosWrappers::SparseMatrix &m)
         }
       else
         {
-          Assert (m(i,j) == 0, ExcInternalError());
           Assert (m.el(i,j) == 0, ExcInternalError());
         }
 
@@ -55,7 +55,9 @@ int main (int argc,char **argv)
   std::ofstream logfile("02/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
 
   try
     {
index c4b11986a770c56c0dad65da0dcbef634c810647..21e7a224c189ba44d2992bed2d9540cc20e17cd3 100644 (file)
@@ -16,7 +16,8 @@
 // intermixed. this poses PETSc some problems, since one has to flush some
 // buffer in between these two types of operations
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -47,7 +48,6 @@ void test (TrilinosWrappers::SparseMatrix &m)
         }
       else
         {
-          Assert (m(i,j) == 0, ExcInternalError());
           Assert (m.el(i,j) == 0, ExcInternalError());
         }
 
@@ -61,7 +61,9 @@ int main (int argc,char **argv)
   std::ofstream logfile("03/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
 
   try
     {
index b9f3369e807bfc3cacacd29c2036f12176d752ea..c87845486755aa889805f8f6c7e765a0dd9c2ad6 100644 (file)
@@ -19,7 +19,8 @@
 // in contrast to trilinos_03, we set and add the same elements here twice, to
 // get double the original value
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -32,6 +33,8 @@ void test (TrilinosWrappers::SparseMatrix &m)
     for (unsigned int j=0; j<m.m(); ++j)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
+
+  m.compress();
   
                                    // then add the same elements again
   for (unsigned int i=0; i<m.m(); ++i)
@@ -51,7 +54,6 @@ void test (TrilinosWrappers::SparseMatrix &m)
         }
       else
         {
-          Assert (m(i,j) == 0, ExcInternalError());
           Assert (m.el(i,j) == 0, ExcInternalError());
         }
 
@@ -65,7 +67,9 @@ int main (int argc,char **argv)
   std::ofstream logfile("03a/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
 
   try
     {
index d6a444701042d03f7599ad009686ddceba0bcc39..baa62d34df735e21f60a2ca59380b8c83d9e4b19 100644 (file)
@@ -19,7 +19,8 @@
 // in contrast to trilinos_03a, we set and add the same elements here twice, then
 // overwrite them again to get the original value back
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -69,7 +70,9 @@ int main (int argc,char **argv)
   std::ofstream logfile("03b/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
 
   try
     {
index bf9bf622cf18d0c201a137a6bafd3684ac8a14ea..671ad9f60dbfce298f7690269616a275ddaa2791 100644 (file)
@@ -14,7 +14,8 @@
 
 // check querying matrix sizes
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -35,7 +36,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("04/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 903d56e47985f5e750be3ca6ec132f94df608f6a..f0f64e80acf420925677c39a1b8c8206d22cd71f 100644 (file)
@@ -15,7 +15,8 @@
 // check querying the number of nonzero elements in
 // TrilinosWrappers::SparseMatrix
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -50,7 +51,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("05/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index cfc5530b87e960cf7ea823413d6f4d349a3eb107..8ffba06838e879df717a9a75472fa77ef04f8155 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::SparseMatrix::l1_norm
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -46,7 +47,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("06/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index cf1d666b39b09d6c05b83e2b269ad757bb9640ed..3bfe0e6debe6beca1e7d5ae76471fe3d3f1dacd0 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::SparseMatrix::linfty_norm
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -46,7 +47,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("07/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 2d2f9c535e82181b3bc6b733728a9dfaefc2400f..c10af21adae2c361a4a0d923a7e1a1759b4bc6d2 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::SparseMatrix::frobenius_norm
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -51,7 +52,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("08/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 79e628d01ddcd1f4c1a48113b888479983bc1248..985ed7df15627b57761bd8b5cb5e6af281ab0e0d 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::SparseMatrix::operator *=
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -58,7 +59,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("09/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 43d2af75fd20e99ad7f2aa7198519b2c70917338..7e6cba76f1b15a7cb1459c27d72d18622a2f4dde 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::SparseMatrix::operator /=
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -58,7 +59,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("10/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index d906d06a7108898fcf3452fc47a0127e6f4b0d90..b6a059dfe7f137b25ad8f35c5c4fe2813f0d3bc1 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::size()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -40,7 +41,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("11/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 50414f14dbfc5db55a44fa26e8fae337b075fee3..a931c891e68ff03cfe7283d17869f5318274a53e 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator() in set-mode
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -53,7 +54,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("12/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index ba362f52749bc8c358547b444033d46e8b62cd76..692730b23e5db8254aec34355b2759c50dd85169 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator() in add-mode
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -53,7 +54,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("13/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 39f68945527ff7b0f7a31561a8412876f30f4421..a1f3ffbb1f8c961873c805f911c5b7c049fedc4b 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator() in set/add-mode alternatingly
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -59,7 +60,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("14/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 19d8e8af9d07a68bc4f206b6d94cff40e2e0ea7f..5565564b0b91e07a3ff8e236b147c360a673209e 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator() in set/add-mode alternatingly, but
 // writing to the same elements
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -56,7 +57,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("15/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 7d869b2bc8eee0ae64d82e3f84ff2a7f01d9a92e..211487fa510619d7e8cf61235caeac07271273bb 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator() in set/add-mode alternatingly, but
 // writing and overwriting the same elements
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -58,7 +59,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("16/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index f73f95f0a17d21be2b2b6b1b40b682b3360e8a92..768efbe78cb658ddee14e5e7f001917e3f579d4b 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::l1_norm()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -45,7 +46,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("17/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index b22380f1ba6006d801b12b5fcadec3f6c28c3845..c67b5bc6e759b8e3ecb44df6c5160c370bd49f06 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::l2_norm()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -46,7 +47,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("18/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index ea2c7e5399244e9ddde339169cf80bda4ad5ec45..2b135fe6684c7ed240686806f54a22bb50e247c2 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::linfty_norm()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -45,7 +46,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("19/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index e416571023b7a51a5fe30914a0e5e75841d4b33c..b7d2e87103a1cce7ca68af1e82fdcc6c763f2bbc 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator *=
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -54,7 +55,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("20/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 1e819ec72646947017f2cdc5f8009918824c76ca..0fe7b8ab1bb55e2a2fa77c11e8831310632b3f2a 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator /=
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -54,7 +55,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("21/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 0d0f1778fe758fd0e3987e96afc69b356d7ca53c..18baad7134f8879a49108a509b1b2f13b233bf60 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator*(Vector) on two vectors that are
 // orthogonal
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -48,7 +49,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("22/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 547972b4d11c679186a895055cefb2857949dc04..cc72557858e3a09872152afb8896b9000ee018a7 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator*(Vector) on two vectors that are
 // not orthogonal
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -55,7 +56,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("23/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index d987aa00eea33d4afb20f9ca1ba21c0c0fef232e..b9c316cbe87cdd7b3eac58abfa2005d473907305 100644 (file)
@@ -16,7 +16,8 @@
 // function has since been removed, so we test for v=0 instead, although that
 // may be covered by one of the other tests
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -48,7 +49,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("24/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 1dc46c11aa3733a9a92273cd7ba6fddbe1c5825c..c757c0780abc1416141385954d5e3c2ad6f2fbbd 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator = (TrilinosScalar) with setting to zero
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -46,7 +47,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("25/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index d41cdce143ced1a3d5b9ebc8f66b84c860c4c4bd..6586959a84782bd01c0ea46d61a1d987cedc8294 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator = (TrilinosScalar) with setting to a
 // nonzero value
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -45,7 +46,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("26/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index bc9746a300bfc9df60f1410eadcef94150044c32..0ab53efa1e7b43ab7e23acec7478cfe89808a196 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator = (Vector) 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -50,7 +51,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("27/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 2abc9660ccd02f86fa0e436af6161853ef310ffd..06e5cae8ad8eae65948273fed1857908b02467e4 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator = (Vector), except that we don't
 // resize the vector to be copied to beforehand
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -52,7 +53,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("28/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 9b39d4fdd74c5aac6893eeb5a9ad3c5bc5b5109c..e74640cae44c3d8c877c3978057be5f1f8aad1fe 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::reinit(fast)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -37,7 +38,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("29/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index d40004a6aa8cbe565862a19499c510c38031f93a..6a350b5019251b1dffc9bacb0554db636fb66015 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::reinit(!fast)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -45,7 +46,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("30/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index d0a1dcb2935d7d122e43b8a856e0a0a3a7e008a7..e04fdb939c73c24e3708910434f90a9be32853c8 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::l2_norm()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -47,7 +48,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("31/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 16fc8917c4ab2c309040ca7b626d26635ee9b782..600030cd732f352e81b72ce0f05efbe8808b7115 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::mean_value()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -47,7 +48,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("32/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 76b87d81058e9f387d9e4b9de6628f00f8494ed8..c84482ffa7dde9b76ee662fb9490d8b242193157 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::lp_norm(3)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -48,7 +49,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("33/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index d73c5aab83eddc1312530c1889f9f40c966fd7d6..bdd678f11751ff48680cbb00494a9a658fe16340 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::all_zero
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -48,7 +49,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("34/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 2f26fc9f513f2ebaa5ceb515c0fd5ffad89a37f9..32e78b216c3926fef0709dc3c92884a6fe07bb3d 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator+=(Vector) 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -64,7 +65,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("35/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 4a2fa81a13905ce24876a1a257b6a2ab737cba11..47150c47d917fca7ef166b556f3044ac15525094 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator-=(Vector) 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -64,7 +65,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("36/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 007980b4b949256a942e870796dcd9185e72aee7..741816565706ef6db00bc0692d317a9295c6f3cc 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::add (scalar)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -44,7 +45,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("37/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 1f0b6eb93aa4d95780f926ced630fd69bd197ac2..a260dd0318288eb4433016abe6f1e0b121eed3e9 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::add(Vector) 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -52,7 +53,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("38/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index dd92d462f770b5d818cd69ae9f0332e06b97b74d..5102866dd46b37a2f65fb36da3eefad8dbe77aa9 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::add(scalar, Vector) 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -52,7 +53,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("39/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 391ddd23ad66c613824a82a033f8f1e2c54c5043..e2ce534b526ae608429162ea24f2e9964ba2a6eb 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::add(s,V,s,V)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -56,7 +57,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("40/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index cd8a05ae32ad1bebe685ac00efc54f4f5ce6323a..370afb6db0970eb19bc50a48938052a9ccbd1291 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::sadd(s, Vector) 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -52,7 +53,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("41/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 261179f05fa221d880da6a49598d99a9b5db950a..5ee5942bbb8049ed7763498718f6b5fcac668bce 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::sadd(scalar, scalar, Vector) 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -52,7 +53,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("42/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index afe7d29d19c8d8d6da3f34077681e86796242c59..436b3df948691bee7be2f37f458c069c0751e10b 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::sadd(s,s,V,s,V)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -56,7 +57,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("43/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 28f88bb55ec925defb9e5131399ee87993d41e5a..db0a53b90732fa97db5808f3f30410ee23f027fb 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::sadd(s,s,V,s,V,s,V)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -60,7 +61,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("44/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 5eba7bc0229a37e22eff2a2ebb5869d849499408..649a7c2b56813cc89caf2d0bf59cb3c73e944b0a 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::scale
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -52,7 +53,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("45/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index fdf1ba7f991b7bd93f032f46021e066c63da7789..12fc510f231e0a2552dcd4b8f965f0411346617f 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::equ (s,V)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -52,7 +53,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("46/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 32eb4177744f4407fbf07aa4d9a3dbcb5b0bfb49..7bbfce11bcec4228505386d751097a9629d02030 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::equ (s,V,s,V)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -56,7 +57,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("47/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 30c07869ceaf3e3a59a08ecb4ed4b37e9da13352..2c0e95568a1489b85186b164f9f0e4a771eff74a 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::ratio
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -58,7 +59,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("48/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 4730d1fa7dce913092f2eeaffd6c9988a6ace6b4..7fda29f65f84ab02f2b2e20f6bae1e434e63cde2 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::operator = (Vector<PetscVector>)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/vector.h>
 #include <fstream>
@@ -49,7 +50,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("49/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index e0c0a2c39d4453e8da73d88b72328b5657e1606c..8179b0f25f3742b6c603ad5e7a1c6e85c97cbb3d 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::operator = (Vector<T>) with T!=TrilinosScalar
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/vector.h>
 #include <fstream>
@@ -65,7 +66,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("50/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 84153d9b3f77a20ec89c8169d76ad6c1328c6154..e33f0773223f1beb61ea43c833b791d37396828d 100644 (file)
@@ -14,7 +14,8 @@
 
 // check copy constructor TrilinosWrappers::Vector::Vector(Vector) 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -50,7 +51,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("51/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 5fe1e5d21b045dcadc47b905a935cffb864d1afb..1abbf3f57c09cf6abfa4c7c93ad143e912d52a81 100644 (file)
@@ -16,7 +16,8 @@
 // TrilinosWrappers::SparseMatrix::set(). like trilinos_01, but use a different
 // constructor for the sparse matrix
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -56,7 +57,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("52/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index ef71214cbf15b078b696328d03740ad3850228df..08677580442415bb493dc683c55a8906383446ff 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator() in set, and later in add-mode
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -55,7 +56,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("53/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 8699955b8167091fd6fdfeec07476677cdaa0485..12d88e785e1fc79598d73d3b2579838a3228dd13 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator() in set, and later in -= mode
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -55,7 +56,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("54/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 23831adf340d5b92af62ab0fc125d1807252098e..bea5d6e9c07fecbe7747c52829051d337d9f3447 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator() in set, and later in *= mode
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -55,7 +56,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("55/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 4eae42d5aae2ff215810ae777308b942ccafb396..e93cab77af41a10bfeec824136a4c0d254f5d753 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator() in set, and later in /= mode
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -55,7 +56,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("56/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 0180745d34e1d0984765dea27f321941d83b2d3c..3166c77e427c4f7a80e106acde6a35b931fa4fdb 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::is_non_zero
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -53,7 +54,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("57/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 6290a9c541b279327a7381803f75d9686ce70510..001bd3b81c39ea0626370d82305db669157375fe 100644 (file)
@@ -14,7 +14,8 @@
 
 // check ::Vector (const TrilinosWrappers::Vector &) copy constructor
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/vector.h>
 
@@ -55,7 +56,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("58/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 7b130dd4b3deae43efaaa516772763658693c15c..4bd71dabe357c6eb1e43d42e3767fd49b1180c91 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector (const ::Vector &) copy constructor
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/vector.h>
 
@@ -58,7 +59,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("59/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index f94de694a12346a841094f51abf73129f87c669d..01e81aafb5862cd3824170e1159a4aafa447e004 100644 (file)
@@ -14,7 +14,8 @@
 
 // check ::Vector::operator = (const TrilinosWrappers::Vector &)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/vector.h>
 
@@ -55,7 +56,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("60/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 01f833cb723abd1cdf1756bce54f74f74ddce21d..248f15beb06068821d04d89b6dba28e79ce9448b 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::Vector::operator = (const ::Vector &)
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/vector.h>
 
@@ -58,7 +59,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("61/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index cbb026efb2e08dc4cb5b2774147708ba7f2837bf..5946d73333a97f2771cee9382137318b370c2c43 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::MatrixBase::clear ()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/vector.h>
 
@@ -43,7 +44,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("62/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 3ba5b3031a4594d995d2d3c697b885d3897d591f..04db7ad9cb375d574513bd9e6d30400ff8d80cbe 100644 (file)
@@ -12,7 +12,8 @@
 //----------------------------  trilinos_63.cc  ---------------------------
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/vector.h>
 
@@ -41,7 +42,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("63/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index df2236546aadb694f775a8f7463235744a58e2df..41c1f1f417a25a6f8cc36011b6bdc0673630b1ff 100644 (file)
@@ -17,7 +17,8 @@
 // TrilinosWrappers::MatrixBase::operator=
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/vector.h>
@@ -46,7 +47,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("64/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 5b264213bb3ed21d83d8cf5569e008c8f5dc0a97..862c98050efe808da33ffda95a421222b073903c 100644 (file)
@@ -15,7 +15,8 @@
 // This test used to fail after upgrading to petsc 2.2.1
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 
 #include <fstream>
@@ -39,7 +40,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("65/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 58b3651c2a46881ebf8d3b0ef4382005acd7f6c0..e90259a05640be4e47107e201d18eda7c919133d 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::MatrixBase::clear_row ()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/vector.h>
 
@@ -97,7 +98,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("66/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 2afaaa5eb17d7e5a6f32c4eaed82a1e15553b35e..13fce014cd96ea9b51cf60cd3c364e60cd10923b 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::MatrixBase::clear_rows ()
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/vector.h>
 
@@ -103,7 +104,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("67/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 3abb4927a07e2f17a63117c5930ac999ccd87c99..0f37d9809690602f367dab3e0f1b743bbe943126 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::MatrixBase::clear_row () with used second argument
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/vector.h>
 
@@ -101,7 +102,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("68/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index a34547709d50d2867b4dc56b4a575efc87521db9..8fb6a26755bb930351344bdd0183aeb79c6fb98a 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosWrappers::MatrixBase::clear_rows () with used second argument
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/vector.h>
 
@@ -107,7 +108,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("69/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index c425a11839a99e8f6545dedd4af550c3cb70ea80..8e503d06a517a8d32a0d0d3f2267108b048a6d83 100644 (file)
@@ -14,7 +14,8 @@
 
 // check TrilinosScalar
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <base/logstream.h>
 #include <lac/trilinos_vector.h>
 
@@ -25,7 +26,10 @@ int main ()
   std::ofstream logfile("70/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   if (typeid(TrilinosScalar)==typeid(double))
     deallog << "double" << std::endl;
index 9f9395bf80a652bd190b4e2c926227a73fc14b40..2e4018f033c11fc9aca1007ab623ad64628b518f 100644 (file)
@@ -14,7 +14,8 @@
 
 // make sure that block vector iterator allows reading and writing correctly
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_block_vector.h>
 #include <fstream>
 #include <iostream>
@@ -80,7 +81,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("block_vector_iterator_01/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 5bb002554885b86358eef3de2c5472301b2b4b9d..6fce3e14713c7385e3c5fd4cf932fba9023305b9 100644 (file)
@@ -14,7 +14,8 @@
 
 // like _01, except that we use operator[] instead of operator*
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_block_vector.h>
 #include <fstream>
 #include <iostream>
@@ -74,7 +75,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("block_vector_iterator_02/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 07fdbfc8b3603c114fee164895f813600456ceb9..94789b532d5b15c9495657f9d4e3c7502f5d1f3f 100644 (file)
@@ -14,7 +14,8 @@
 // this test is an adaptation of lac/block_vector_iterator for PETSc block
 // vectors
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <base/logstream.h>
 #include <lac/trilinos_block_vector.h>
 #include <fstream>
@@ -325,7 +326,10 @@ int main (int argc,char **argv)
   logfile.precision(3);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   try
     {
index 4469d40e5e3a97f4782ea6690c72845c3645ea71..5b97fd080999278aa797eab282d2a951f5969e76 100644 (file)
@@ -14,7 +14,8 @@
 // test the CG solver using the Trilinos matrix and vector classes
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -63,7 +64,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   {
     SolverControl control(100, 1.e-3);
index 0fd4d5d9fbcb51d1902ed1c01fe514c2880af50c..327d5e1cf07de44d59f2ffdbc0a0953db50085b6 100644 (file)
@@ -15,7 +15,8 @@
 
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -65,7 +66,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index 73c626761bda03164691f2f58b74f17f33298c16..e9a3fe4f6f058d5bd8893c0d74a71ea68e362c13 100644 (file)
@@ -15,7 +15,8 @@
 
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -65,7 +66,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index 651101cbeff4865338116c932a019256b7d6ecef..d045c1255c0a12d83b2d0e82fa5d90a9a52046a6 100644 (file)
@@ -14,7 +14,8 @@
 // test the MINRES solver using the Trilinos matrix and vector classes
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -64,7 +65,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index 204530fb417d437b6e7064fb1fdd339cbc5fd877..8c8c0b5359461767e6986c7aa7462251f61b2c43 100644 (file)
@@ -14,7 +14,8 @@
 // test the QMRS solver using the Trilinos matrix and vector classes
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -63,7 +64,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index dd742d9dc6526750669c1a87fc4f35f5db02bbbf..3ff908ecd6cf8ee0e7cdb0672746271b0f20513c 100644 (file)
@@ -14,7 +14,8 @@
 // test the Richardson solver using the Trilinos matrix and vector classes
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -62,7 +63,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   {
     const unsigned int size = 32;
index 9f3f3cda0ec89fec4a751b9c00a09624a1b0d35a..0b80cc8cf231a8bfc3fe59038a2a769aeeff45fe 100644 (file)
@@ -23,7 +23,8 @@
 // and with the option -log_info. All the output should say that no additional
 // malloc calls have been performed
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/compressed_sparsity_pattern.h>
 #include <fstream>
@@ -118,7 +119,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("sparse_matrix_01/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 292efb50c37a034d474d8b2dbcd7edef294e934b..496e98e2db5f1fad7766d1ffb60a9ec637467cce 100644 (file)
@@ -19,7 +19,8 @@
 //
 // the tests build the 5-point stencil matrix for a uniform grid of size N*N
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/sparse_matrix.h>
 #include <lac/vector.h>
 #include <lac/trilinos_sparse_matrix.h>
index d21205e42839c295ac57eb66ce82f3cc004e1c53..9156f13e71268107619f4ad94aebb6939225076a 100644 (file)
@@ -19,7 +19,8 @@
 //
 // the tests build the 5-point stencil matrix for a uniform grid of size N*N
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/sparse_matrix.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/trilinos_vector.h>
index 8ed90df490711f5c15ee4842a2a4826b62e05e61..9a1cdf49341c571ba6aefb81593ded6996aeaa92 100644 (file)
@@ -19,7 +19,8 @@
 //
 // the tests build the 5-point stencil matrix for a uniform grid of size N*N
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/sparse_matrix.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/trilinos_vector.h>
index f4def4cc4b6c174114b9bb42d1c403569b8d21f4..135cbc4b5dbb8278279c1bb4357b1b3939565d5e 100644 (file)
@@ -25,7 +25,8 @@
 // matrix in a consecutive fashion, but rather according to the order of
 // degrees of freedom in the sequence of cells that we traverse
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/sparse_matrix.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <lac/trilinos_vector.h>
index 22e1ed7c5b9125fa19ddb75f7d305810d0b7b583..ac6774aac6b73304fa5290e258cefc2557db2b29 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos Richardson solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -60,7 +61,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index 613d89aea6e9d22ff6349bba0a14d748461ea19a..ec6cedce9803f64bbce3a70e11bb2dfcf788ca85 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos Chebychev solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -59,7 +60,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   {  
     SolverControl control(100, 1.e-3);
index f4e135edde03a1a65bc3cd81636783dee0422b79..5fa013f82ba1a741e2a371f09d47b08967484a7d 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos CG solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -58,7 +59,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   {  
     SolverControl control(100, 1.e-3);
index 04b965a56b9605ec88672f11a8f583ca7d7e2724..aae751de12a1653f90ff5dcc76659e98da0a4f4b 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos BiCG solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -58,7 +59,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   {  
     SolverControl control(100, 1.e-3);
index c97a80a672178fa347c03077021345e90b4ca5ac..d4d28b01d5f49509fa7f6fd210842137258bef75 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos GMRES solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -58,7 +59,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   {  
     SolverControl control(100, 1.e-3);
index 570093ebda1933987838d9da4287c644758df312..5fe5bcd4891da00c57c171e8281018cb87105419 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos Bicgstab solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -58,7 +59,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index 0a89e3e4a68c99338663edca0920222713ef346d..74f4362acf441c6e6ca3154884fc53ea0b24f4cc 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos CGS solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -58,7 +59,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   {
     SolverControl control(100, 1.e-3);
index fbb2adf667ad174170c5271d4a412aa2e918c734..71f522757489f29073a6cb33ad1cca2981d9635e 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos TFQMR solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -58,7 +59,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index b4e47661d21789b2c83bf22321ca4ada1309274b..4e2995179d74a806ff4703eaf778bb18fad7c213 100644 (file)
@@ -18,7 +18,8 @@
 
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -62,7 +63,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index 66d59cf2dcff267faa056b7223bedf1851daafbe..cb7a6ac4aecf67a1a814aa7a076f2eb8b8d304bc 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos CR solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -58,7 +59,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index 46cc98200b74c185337507dc433190eba87c1cff..f68de16c1530674d3b17aa4961e27e9e8083e264 100644 (file)
@@ -14,7 +14,8 @@
 // test the Trilinos LSQR solver
 
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include "../lac/testmatrix.h"
 #include <cmath>
 #include <fstream>
@@ -59,7 +60,10 @@ int main(int argc, char **argv)
   logfile.precision(4);
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
   
   {
     SolverControl control(100, 1.e-3);
index d8e3c29e6b915b0ceee0c00feebb03146797ba47..341cd27bbd799f841ac388d3cb0f4a6202a16f2f 100644 (file)
@@ -14,7 +14,8 @@
 
 // test TrilinosWrappers::MatrixBase::const_iterator
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_sparse_matrix.h>    
 #include <fstream>
 #include <iostream>
@@ -44,7 +45,10 @@ int main (int argc,char **argv)
   std::ofstream logfile("sparse_matrix_iterator_01/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index bc961f0feaa526fd01605241d3911c8f1f321c00..cd8308bcc6a486d1645f3603ea6c1d11cc9ea622 100644 (file)
@@ -14,7 +14,8 @@
 
 // check SparseMatrix::vmult
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <fstream>
@@ -61,7 +62,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("sparse_matrix_vector_01/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index a557d076abc1262b9a333d68737f7bf4eaaecc09..9f5b17c20a89f1f5759451ca6035575bd0fc1267 100644 (file)
@@ -14,7 +14,8 @@
 
 // check SparseMatrix::Tvmult
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <fstream>
@@ -61,7 +62,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("sparse_matrix_vector_02/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 4a4e4670e4724be421f2d5f9e60eac2b7d46e370..8f6b858194a7ad39035881e9010ed8fe7407081d 100644 (file)
@@ -14,7 +14,8 @@
 
 // check SparseMatrix::vmult_add
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <fstream>
@@ -64,7 +65,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("sparse_matrix_vector_03/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 48b8454bfa74ebe3828b70e0ab5cc0f794e836d4..edf3b7e0341a5ece711b17d59078923ff61a9755 100644 (file)
@@ -14,7 +14,8 @@
 
 // check SparseMatrix::Tvmult_add
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <fstream>
@@ -64,7 +65,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("sparse_matrix_vector_04/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 93ea33ec2925b97a780c145ad6865347c63dd547..bbea123f50f5f9687b213b20f7b0c06fc8f17a6a 100644 (file)
@@ -14,7 +14,8 @@
 
 // check SparseMatrix::matrix_scalar_product
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <fstream>
@@ -67,7 +68,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("sparse_matrix_vector_05/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 72b98ce2650344d75ee8c91f301edf9952424409..842f1bc760dab09fe2750556ff2dd76961a0ae51 100644 (file)
@@ -14,7 +14,8 @@
 
 // check SparseMatrix::matrix_norm_square
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <fstream>
@@ -59,7 +60,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("sparse_matrix_vector_06/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index f0cf939a91e96df16963e3e7cbfde45fe1b6b253..b860e4f89f4baf2702c86c8be642be24c3899f1c 100644 (file)
@@ -14,7 +14,8 @@
 
 // check SparseMatrix::matrix_norm_square
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>
 #include <lac/trilinos_sparse_matrix.h>
 #include <fstream>
@@ -69,7 +70,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("sparse_matrix_vector_07/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index a5b6d1dad9570bee90b3f773b7319c02d8ff883e..b77a142d843437360b7a3df1adfb4eefd0d066d8 100644 (file)
@@ -20,7 +20,8 @@
 //
 // this was fixed 2004-04-05, and this test checks that it works
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -51,7 +52,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("vector_assign_01/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 9ffa3b52e51831365ca0c3e38c6581b571826140..ac49b02033e25d4c54d3d1f64ffc51886b4b874c 100644 (file)
@@ -19,7 +19,8 @@
 // argument to the user-defined operator+=. This is not exciting, but since I
 // wrote the test to make sure it works this way, let's keep it then...
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -50,7 +51,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("vector_assign_02/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 3089b094d0bc4bec3b55f655a4d6aff85dfb9343..40ab8f0755478e5a384beb33b55af841ec804d47 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator==(TrilinosWrappers::Vector) for vectors that are not
 // equal
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -46,7 +47,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("vector_equality_1/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 3fb91d6d948bd6bf0a33cc205ee2a94aee2fc0d0..000ca72fcfc5d2813fcf0fc5ca02be5fda4f161f 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator==(TrilinosWrappers::Vector) for vectors that are
 // equal
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -48,7 +49,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("vector_equality_2/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index 9b6fccc11023aaa8c278054183d31e110cb928e7..89172f50433eb3c7eedd5c49aca3f8b9fd7d3d27 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator!=(TrilinosWrappers::Vector) for vectors that are not
 // equal
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -46,7 +47,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("vector_equality_3/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {
index c2bbf3995d6aae29243ed4c46f81d5703de38a8f..46c1a826a849b160178356aadbf052b716320611 100644 (file)
@@ -15,7 +15,8 @@
 // check TrilinosWrappers::Vector::operator!=(TrilinosWrappers::Vector) for vectors that are
 // equal
 
-#include "../tests.h"
+#include "../tests.h" 
+#include <base/utilities.h>
 #include <lac/trilinos_vector.h>    
 #include <fstream>
 #include <iostream>
@@ -48,7 +49,10 @@ int main (int argc, char **argv)
   std::ofstream logfile("vector_equality_4/output");
   deallog.attach(logfile);
   deallog.depth_console(0);
-  deallog.threshold_double(1.e-10);
+  deallog.threshold_double(1.e-10); 
+
+  Utilities::System::MPI_InitFinalize mpi_initialization (argc, argv);
+
 
   try
     {

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.