]> https://gitweb.dealii.org/ - dealii.git/commitdiff
fix a few broken tests 596/head
authorTimo Heister <timo.heister@gmail.com>
Tue, 24 Feb 2015 17:51:50 +0000 (12:51 -0500)
committerTimo Heister <timo.heister@gmail.com>
Tue, 24 Feb 2015 22:55:47 +0000 (17:55 -0500)
Looks like somebody messed up the VectorOperation in compress() in many
places.

95 files changed:
tests/mpi/petsc_distribute_01_block.cc
tests/petsc/01.cc
tests/petsc/05.cc
tests/petsc/06.cc
tests/petsc/07.cc
tests/petsc/08.cc
tests/petsc/09.cc
tests/petsc/10.cc
tests/petsc/11.cc
tests/petsc/12.cc
tests/petsc/17.cc
tests/petsc/19.cc
tests/petsc/20.cc
tests/petsc/21.cc
tests/petsc/22.cc
tests/petsc/23.cc
tests/petsc/24.cc
tests/petsc/25.cc
tests/petsc/26.cc
tests/petsc/27.cc
tests/petsc/28.cc
tests/petsc/31.cc
tests/petsc/32.cc
tests/petsc/33.cc
tests/petsc/34.cc
tests/petsc/35.cc
tests/petsc/36.cc
tests/petsc/37.cc
tests/petsc/38.cc
tests/petsc/39.cc
tests/petsc/40.cc
tests/petsc/41.cc
tests/petsc/42.cc
tests/petsc/43.cc
tests/petsc/45.cc
tests/petsc/46.cc
tests/petsc/47.cc
tests/petsc/48.cc
tests/petsc/51.cc
tests/petsc/52.cc
tests/petsc/55.cc
tests/petsc/56.cc
tests/petsc/64.cc
tests/petsc/deal_solver_01.cc
tests/petsc/deal_solver_02.cc
tests/petsc/deal_solver_03.cc
tests/petsc/deal_solver_04.cc
tests/petsc/deal_solver_05.cc
tests/petsc/full_matrix_00.cc
tests/petsc/full_matrix_01.cc
tests/petsc/full_matrix_05.cc
tests/petsc/full_matrix_06.cc
tests/petsc/full_matrix_08.cc
tests/petsc/full_matrix_09.cc
tests/petsc/full_matrix_10.cc
tests/petsc/full_matrix_iterator_01.cc
tests/petsc/full_matrix_vector_01.cc
tests/petsc/full_matrix_vector_02.cc
tests/petsc/full_matrix_vector_03.cc
tests/petsc/full_matrix_vector_05.cc
tests/petsc/full_matrix_vector_06.cc
tests/petsc/full_matrix_vector_07.cc
tests/petsc/solver_01.cc
tests/petsc/solver_02.cc
tests/petsc/solver_03.cc
tests/petsc/solver_03_mf.cc
tests/petsc/solver_03_precondition_boomeramg.cc
tests/petsc/solver_03_precondition_boomeramg_symmetric.cc
tests/petsc/solver_03_precondition_eisenstat.cc
tests/petsc/solver_03_precondition_icc.cc
tests/petsc/solver_03_precondition_lu.cc
tests/petsc/solver_03_precondition_parasails.cc
tests/petsc/solver_03_precondition_sor.cc
tests/petsc/solver_03_precondition_ssor.cc
tests/petsc/solver_04.cc
tests/petsc/solver_05.cc
tests/petsc/solver_06.cc
tests/petsc/solver_07.cc
tests/petsc/solver_08.cc
tests/petsc/solver_09.cc
tests/petsc/solver_10.cc
tests/petsc/solver_11.cc
tests/petsc/solver_12.cc
tests/petsc/solver_13.cc
tests/petsc/sparse_direct_mumps.cc
tests/petsc/sparse_matrix_01.cc
tests/petsc/sparse_matrix_02.cc
tests/petsc/sparse_matrix_iterator_01.cc
tests/petsc/sparse_matrix_vector_01.cc
tests/petsc/sparse_matrix_vector_02.cc
tests/petsc/sparse_matrix_vector_03.cc
tests/petsc/sparse_matrix_vector_04.cc
tests/petsc/sparse_matrix_vector_05.cc
tests/petsc/sparse_matrix_vector_06.cc
tests/petsc/sparse_matrix_vector_07.cc

index 631a1d9dab8dae27c9da032f8d7b271077fc0686..43fafab6b55bda19e96584fb6febb97ca4f3bd06 100644 (file)
@@ -52,7 +52,7 @@ void test()
     vec.block(0)(i) = i;
   for (unsigned int i=vec.block(1).local_range().first; i<vec.block(1).local_range().second; ++i)
     vec.block(1)(i) = i;
-  vec.compress(VectorOperation::add);
+  vec.compress(VectorOperation::insert);
 
   // verify correctness so far
   {
index a4de9a246abaf73a269d4f32d126e810d5afe074..3e29ed4b4526b5589c13fe18b98837e8baebae74 100644 (file)
@@ -32,7 +32,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then make sure we retrieve the same ones
   for (unsigned int i=0; i<m.m(); ++i)
index 27bb210372ecf1d7707ea4aaa214d0da6a35e9d9..07330cc82e5b766fa0567b981e332260299e4c3f 100644 (file)
@@ -37,7 +37,7 @@ void test (PETScWrappers::SparseMatrix &m)
           ++counter;
         }
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   deallog << m.n_nonzero_elements() << std::endl;
   Assert (m.n_nonzero_elements() == counter,
index aefef366b54f2524a758b627deddf89e3800a10f..2c9d927e7c52e94ef6875784018fd669ef274fb8 100644 (file)
@@ -32,7 +32,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // l1-norm (max col-sum)
index c61e90ecb574caabb081d63e14a2054c62ab5450..90a9b745dba4b7ed0ef6fd5c32e8627b28949344 100644 (file)
@@ -32,7 +32,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // linfty-norm (max row-sum)
index 2fb1dfa25345259b64e8365d508e2a1ab849fea8..6fe6ec05efe4bc314b3efe6d388013222a3a6449 100644 (file)
@@ -37,7 +37,7 @@ void test (PETScWrappers::SparseMatrix &m)
         }
   norm = std::sqrt(norm);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // l2-norm (max row-sum)
index fb6b7a6e14ea0ac9c1d8205a24f271634b7d0b50..4914c887142dd3ea9a6fd559d2f311c5775c8391 100644 (file)
@@ -31,7 +31,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then multiply everything by 1.25 and
   // make sure we retrieve the values we
index c0bc8b65efd4a5cd4fb2e247407d969bc86ede92..eeba560d9e0212f79c996b1b7d63f8d137f41982 100644 (file)
@@ -31,7 +31,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then divide everything by 4/3 and
   // make sure we retrieve the values we
index 4cc621ed2d2694c1f02a18cc3c828f5a5dd2740f..0e8c572e932de5f4a8bca62593fbdcadc17e01a2 100644 (file)
@@ -29,7 +29,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); i+=1+i)
     v(i) = i;
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   Assert (v.size() == 100, ExcInternalError());
 
index 74c51d4e7dd73a96ddc789eac302d4c51426838c..021e3e5e8a3de664145aba6a40d71387a9152bd5 100644 (file)
@@ -36,7 +36,7 @@ void test (PETScWrappers::Vector &v)
       pattern[i] = true;
     }
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // check that they are ok, and this time
   // all of them
index 3aa9d9951a786e730ee7b8971dc533d3a1fcfab4..59f9ed383424bb414c0671876fe1191e9b010f01 100644 (file)
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm += std::fabs(1.*i);
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   Assert (v.l1_norm() == norm, ExcInternalError());
index f04a700336e90e9129ec154a4e9df01c68d25d3e..2f801f6488f43db4ba7d800ecf142858d5c75779 100644 (file)
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm = std::max(norm,fabs(i));
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   Assert (v.linfty_norm() == norm, ExcInternalError());
index 1610b28891b26c691aaea0309d9af914dcad775e..6c845ce1cec7c6f01d21466433dbe371690ed166 100644 (file)
@@ -35,7 +35,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       pattern[i] = true;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // multiply v with 5/4
   v *= 5./4.;
index 036ed63cf4317df75f15812f0b4226531ea937db..eee63e666b8120ccceafa519eebe33fd29e4e43b 100644 (file)
@@ -35,7 +35,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       pattern[i] = true;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // multiply v with 3/4
   v /= 4./3.;
index e6f2fc698210431b9b26525666f52a0f95eff1e5..821a1ff3bc11258f4ab4a1ea22b66fab15f24719 100644 (file)
@@ -35,8 +35,8 @@ void test (PETScWrappers::Vector &v,
       v(i) = i;
     else
       w(i) = i;
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // make sure the scalar product is zero
   Assert (v*w == 0, ExcInternalError());
index ef2935852227632014d46b1d1665e5445e9688ed..25f843160b321e8d00eead8ea76c2f7f111e0a88 100644 (file)
@@ -42,8 +42,8 @@ void test (PETScWrappers::Vector &v,
         }
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // make sure the scalar product is zero
   Assert (v*w == product, ExcInternalError());
index 3bbd795d0ef4e6c0013c38ef98126e9f0533ed00..47ade91b2a0f968a6ddf125a6624a10ec877a41b 100644 (file)
@@ -32,7 +32,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then clear it again and make sure the
   // vector is really empty
index 993f4f4e0ad912cfe9401c7704d80989f8e933f1..aa6715587b12c54092009c7ff90c2c068e291491 100644 (file)
@@ -30,7 +30,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then clear it again and make sure the
   // vector is really empty
index 25309df9afbfadf53e01afa4bee39c2494bf59ab..3ded09884877e917ce4e6d02535d113668cea8b3 100644 (file)
@@ -31,7 +31,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   const unsigned int sz = v.size();
   v = 2;
index 3c3632c06f3d74e4681f5f992e1a3eddcb7656b3..b773cb4933f500703be968e1a9e55e5279bd25fe 100644 (file)
@@ -30,7 +30,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then copy it
   PETScWrappers::Vector w (v.size());
index 5e96a0a06d88300f5232b76a989618ee9be44d77..130b4c73e59f2524828cc1ec9424f2dec6b5259d 100644 (file)
@@ -31,7 +31,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then copy it to a vector of different
   // size
index 875f12311c1479adf960c0e7b3ab54284fd3ce81..0a1bcb31aa4d6939fa61eba23417ab03325ba4c0 100644 (file)
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       norm += std::fabs(1.*i)*std::fabs(1.*i);
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   const double eps=typeid(PetscScalar)==typeid(double) ? 1e-14 : 1e-5;
index d490ac0f3e8d1379005653f2157ce754df13f4d8..a3727120061ac347cafdee1da0379cf34a132ad4 100644 (file)
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       sum += i;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   const double eps=typeid(PetscScalar)==typeid(double) ? 1e-14 : 1e-5;
index f367eb5594d3192624c607d592379b12a9e7a350..62e686c7482d5395c3dab13ad3aa19b3f1397535 100644 (file)
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       sum += i*i*i;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then check the norm
   const double eps=typeid(PetscScalar)==typeid(double) ? 1e-14 : 1e-5;
index 707bbe2b0aaf8d9ee5a310b5aae95699350650e3..3d685b81dd7dac5e48d495c0065b460743c1be43 100644 (file)
@@ -33,7 +33,7 @@ void test (PETScWrappers::Vector &v)
       v(i) = i;
       sum += i*i*i;
     }
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // set them to zero again
   v = 0;
index 036130efdd9699b7722af0b65c21fef3f402ee02..d3a1b0758a2a822ba9f2446a619cd75a557a33a6 100644 (file)
@@ -36,8 +36,8 @@ void test (PETScWrappers::Vector &v,
         w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v += w;
 
index ff457604bb03edd4526d906e746e249fea878830..718aaf7721dbdb8b6e6a579ffd09407c54a21199 100644 (file)
@@ -36,8 +36,8 @@ void test (PETScWrappers::Vector &v,
         w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v -= w;
 
index 91d11537e4e0322c0af67482a73a3aa8eb43fd29..bf8522ae1ea3da6b1c459b39297fb74c992715b3 100644 (file)
@@ -29,7 +29,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   v.add (1.);
 
index dc24186477365fe19d4620d62891bd4bd6bdf764..2817de384ad01390ecdde496108d965fbce9b4f3 100644 (file)
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.add (w);
 
index 13ba1db9ed7296236e1b7a25a3306127af626e79..3f85e85ab0ff5a813e7362bbb5801291c25364de 100644 (file)
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.add (2, w);
 
index 483edb86886fc98b134e71eb647ef558d303ffa5..3918118b0ea13c576da87016459e829cf9ca1ca1 100644 (file)
@@ -35,9 +35,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
-  x.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
+  x.compress (VectorOperation::insert);
 
   v.add (2, w, 3, x);
 
index 32d369d7ae98cdea956ff8bdf0bb64e2082594b4..5c9eac41b121781cba02fd3500ad15b02c499244 100644 (file)
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.sadd (2, w);
 
index 319960ff5dc0ccddf69ef374f937803cec714c4a..fe35044bb92ae2073fefba166901f6da7124a5f7 100644 (file)
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.sadd (3, 2, w);
 
index 7fa4fb03e641b6a542f3364d4c4d9dcc6a298636..1e137386bc661327cc017e00e3a1c46a8a2110e5 100644 (file)
@@ -35,9 +35,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
-  x.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
+  x.compress (VectorOperation::insert);
 
   v.sadd (1.5, 2, w, 3, x);
 
index 2feef69fc02a9425729f467665adda0dda03b096..b5d70c399aeeafe85d6a3f8d5320a45c73e234e4 100644 (file)
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.scale (w);
 
index 5bf4056acbbe4ca112923fc9a6d2ee073983e370..a24193216b7b6e0000e265ca2a07d750528a5942 100644 (file)
@@ -33,8 +33,8 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   v.equ (2, w);
 
index 0763eaa0e8e6dfaaf82957d91615dcf7ed624345..c9ccac82b19055eed4b32ea4a37b7dd0084d7bb8 100644 (file)
@@ -35,9 +35,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
-  x.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
+  x.compress (VectorOperation::insert);
 
   v.equ (2, w, 3, x);
 
index eea0b002a9f82dfac3397bd98d78927294a0b1fd..97e658e7b66c2ca6f51932704bbe39acbc851795 100644 (file)
@@ -35,9 +35,9 @@ void test (PETScWrappers::Vector &v,
       x(i) = i+2.;
     }
 
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
-  x.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
+  x.compress (VectorOperation::insert);
 
   v.ratio (w, x);
 
index eb98b24895e4e496259b782e57ddf2d944039b34..2eb2ccea39a75f3fd6e33cf86be137d2f41bad9c 100644 (file)
@@ -30,11 +30,10 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     if (i%3 == 0)
       v(i) = i+1.;
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // then copy it
   PETScWrappers::Vector w (v);
-  w.compress (VectorOperation::add);
 
   // make sure they're equal
   deallog << v *w << ' ' << v.l2_norm() * w.l2_norm()
index 591944405a93bbe4048e79e4b357c375f0445546..1491e969b7ca049abb9e95c6aee2038934473067 100644 (file)
@@ -33,7 +33,7 @@ void test (PETScWrappers::SparseMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then make sure we retrieve the same ones
   for (unsigned int i=0; i<m.m(); ++i)
index 0013f4da53135467ffbdfd1c273f5fa6f6455e0c..7949d49cf7da75612b0cbf9875f52df377c6d791 100644 (file)
@@ -38,7 +38,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); i+=1+i)
     v(i) *= 2;
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // check that they are ok, and this time
   // all of them
index 34b6eda15915ecdee684f1ba07491404bea946bf..7bd342ad935728857806f6812e8848c21d6b10b3 100644 (file)
@@ -38,7 +38,7 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); i+=1+i)
     v(i) /= 2;
 
-  v.compress (VectorOperation::add);
+  v.compress (VectorOperation::insert);
 
   // check that they are ok, and this time
   // all of them
index 0f17ad5f6d49cc591f08b2c19598bee74c90ca7f..59e804ea6835ac76a40838d2ba1cfdac3244a2c0 100644 (file)
@@ -34,8 +34,8 @@ template<typename MatrixType>
 void test (MatrixType &m)
 {
   m.add(0,0,1);
+  m.compress(VectorOperation::add);
   m = 0;
-  m.compress();
 
   Assert(fabs(m.frobenius_norm())<1e-15, ExcInternalError());
 
index 51fe0842809e1158dedae9d9764d4f0ff8dadc28..c9ad7e4c877709a38b98355b951a834879201ede 100644 (file)
@@ -83,13 +83,11 @@ int main(int argc, char **argv)
     FDMatrix testproblem(size, size);
     PETScWrappers::SparseMatrix  A(dim, dim, 5);
     testproblem.five_point(A);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverCG<PETScWrappers::Vector> solver(control,mem);
index 0ef008f3e702ce3524a90c14ee1d1fc442bed4f8..07ca880f3493d77acf8e31185444cef0bff29c07 100644 (file)
@@ -84,13 +84,11 @@ int main(int argc, char **argv)
     FDMatrix testproblem(size, size);
     PETScWrappers::SparseMatrix  A(dim, dim, 5);
     testproblem.five_point(A);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverBicgstab<PETScWrappers::Vector> solver(control,mem);
index 3edb73d612c1b81a4dd1b12b40289c5c86711765..27daaebc5f8a86563065bae7f242b6875fbcc173 100644 (file)
@@ -89,9 +89,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverGMRES<PETScWrappers::Vector> solver(control,mem);
index d3f932664b862913476646b782e1097fb71ec9cb..99c19b0b6e3729d3b28753b671ea07451bc420db 100644 (file)
@@ -88,9 +88,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverMinRes<PETScWrappers::Vector> solver(control,mem);
index ce225df08d4eb8f6c64913b72ad1c63d6ea3b06c..72bd9bdade5e77d7e979c5ee772b5870a7243ab8 100644 (file)
@@ -87,9 +87,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     GrowingVectorMemory<PETScWrappers::Vector> mem;
     SolverQMRS<PETScWrappers::Vector> solver(control,mem);
index 824d1fbe6afbb680d10bf63dc749918cd9c4aadc..c9914e21272986d585e9b7a474b17e0c11465f64 100644 (file)
@@ -34,7 +34,7 @@ void test (PETScWrappers::FullMatrix &m)
     for (unsigned int j=0; j<m.n(); ++j)
       m.set (i, j, i+2*j);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // things we know
   Assert (m.m() == 3, ExcInternalError());
@@ -57,7 +57,7 @@ void test (PETScWrappers::FullMatrix &m)
   for (unsigned int i=0; i<m.m(); ++i)
     for (unsigned int j=0; j<m.n(); ++j)
       m.set (i, j, j+2*i);
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
   
   // things we know
   Assert (m.m() == 5, ExcInternalError());
index a19fb62e18a138607436f272af6dabc7f15f04fc..f3975ca74057259f7454f4bfe095b9498543fbc5 100644 (file)
@@ -32,7 +32,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then make sure we retrieve the same ones
   for (unsigned int i=0; i<m.m(); ++i)
index 776152ba56b303e628001395c20917f138a242d8..d3c45b6947502372d4f92ed69ee660c01441c4ae 100644 (file)
@@ -37,7 +37,7 @@ void test (PETScWrappers::FullMatrix &m)
           ++counter;
         }
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // check n_nonzero_elements then
   // output a dummy number.
index a7b4a2b3fc66bdecb478ce75d546a27da10bee64..76bc6c60d308d8e2e100b87e1d5f9e2bf0d8a823 100644 (file)
@@ -32,7 +32,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // l1-norm (max col-sum)
index d6f0fdcbe5e2ab59e5d0ccf6463ed03e1da1e516..3ef0233d99eb497116723fade8683316e18fbb49 100644 (file)
@@ -37,7 +37,7 @@ void test (PETScWrappers::FullMatrix &m)
         }
   norm = std::sqrt(norm);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // compare against the exact value of the
   // l2-norm (max row-sum)
index 744fd1cc6b3b3f425700503402c63dd4bf2079ef..513ad58f0931212d5f0863ed352e424c7e2ed47d 100644 (file)
@@ -31,7 +31,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then multiply everything by 1.25 and
   // make sure we retrieve the values we
index 2f1088d211cef01e233ff1af90fdca5a79ce258c..d0b87e0d3d3a6f0d9881c867c085e8dce7988567 100644 (file)
@@ -31,7 +31,7 @@ void test (PETScWrappers::FullMatrix &m)
       if ((i+2*j+1) % 3 == 0)
         m.set (i,j, i*j*.5+.5);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   // then divide everything by 4/3 and
   // make sure we retrieve the values we
index 21fdf5414e5e254f69752cc67195f32884a07c50..108d4c13f421ead9f9f88750add9df859a169a3b 100644 (file)
@@ -29,7 +29,7 @@ void test ()
   m.set (0,0,1);
   m.set (1,1,2);
   m.set (1,2,3);
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
   PETScWrappers::FullMatrix::const_iterator i = m.begin();
   deallog << i->row() << ' ' << i->column() << ' ' << i->value() << std::endl;
   ++i;
index 9685ee438245cb48e999cc462ed70d2335470ed8..ef9a9a83048d4410061defecfad45824984e4785 100644 (file)
@@ -36,9 +36,8 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // w:=Mv
   m.vmult (w,v);
index 12d2454b6d1d62a694b6f4ea099c32edbdd5d3db..b63166e148ea9d2c3542ea39216dcbc98b2ae0c0 100644 (file)
@@ -36,9 +36,8 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // w:=Mv
   m.Tvmult (w,v);
index a9c15194310b9862f0b3b6427ad189500216e737..5abc75fa6d67c21d834a891be885b95d857d01f3 100644 (file)
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // w:=Mv
   m.vmult_add (w,v);
index d5dd8bc86f018437bcbfb1a57f35cec6f71ac22a..27982b5dadc48fd673000f1120e48ee8a80ed44a 100644 (file)
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // <w,Mv>
   const PetscScalar s = m.matrix_scalar_product (w,v);
index a93f05b9d09c68b07d6fe4406d9dc4dd2a4538d1..5a70894fcc1ead97ae2e9f30e5880948d5da949b 100644 (file)
@@ -35,8 +35,8 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // <w,Mv>
   const PetscScalar s = m.matrix_norm_square (v);
index a8226aa8b02e1a40c03697e7c2a856d3e88ba09b..001829635e841968b6b1369257a5920240ba281b 100644 (file)
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // x=w-Mv
   const double s = m.residual (x, v, w);
index 0294903e5504430c0d99d7502e6e14a27d509e40..0a48809690730ff92b217fe2e699588c337e098b 100644 (file)
@@ -81,9 +81,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverRichardson solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index 372fe4acc783f417db76f9c317f3c6a4db61ec98..58de09df2e75a6f9f342cda9b3191a7f6edc0f42 100644 (file)
@@ -84,9 +84,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverChebychev solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index c8ba190a513893bfb4e21a0d92d789cdbbdbc902..4c839e3c80fe785b814886a6fcaf13b19aca0241 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index ad3b8eae4c649c2c879f290505921cdabf62da7c..f651d2bc45b6f0830f0235146cd6c6b137278888 100644 (file)
@@ -80,9 +80,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionNone preconditioner(A);
index 3952dd18c1a11d88818c693ed499a88dc837ae29..404b096793f10ffa8a14d67b427bfe306d1bfa4a 100644 (file)
@@ -83,9 +83,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionBoomerAMG preconditioner(A);
index c301083b4ac221c5777591aa7bc5317685d94af7..ce04b203bd85a53b08d91ba8120dc6c04877365c 100644 (file)
@@ -84,9 +84,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionBoomerAMG preconditioner(A, true);
index 8528015206656887c446a62396a01dfa57bec690..5747cf0e356669857dff8b259ba8bd80f563eda0 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionEisenstat preconditioner(A);
index 6f3fa7708a7d60b2fe1c37df2dfb787f5fc1e55a..7a04bc77708e8f8bd03fa29b6b2c0cb8b57a10b8 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionICC preconditioner(A);
index 3bb0a590b430b9cbfb44cab5e946c893c226af05..b7b88a0f378c0638aa08131a9485de4e77440719 100644 (file)
@@ -89,9 +89,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionLU preconditioner(A);
index def46b1520c94de76e733c315b72acc31e34359f..269fdaae8fae321d42d7ab37834835d16069edd1 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionParaSails preconditioner(A);
index ceab998ba09fa01ace5252ddaeaafe4a59c74e96..91edc6de6d7214632cb65333a544f62a2e48298f 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionSOR preconditioner(A);
index 5f47ab3d907a79072965f92063abf966fcbe92e5..9856f14bb78d4029899141c66787b35e29760212 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCG solver(control);
     PETScWrappers::PreconditionSSOR preconditioner(A);
index 895845c74e9f0a34fc7622ece07b62ad9178de1e..0f57a698b6b9816b23c05ce0ed13b64c268b4965 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverBiCG solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index c88c9e4786f80bd2197bc1ad69ac25855a763252..d3bcd3d20a1f3c33c2a0f873fea71fe0ba80b179 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverGMRES solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index 1a81d8e151858b89d1531c10dfe44995b0f33aaa..95c65ac3aff3c2cbe9628fc8a701dbbfe0cc23ee 100644 (file)
@@ -85,9 +85,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverBicgstab solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index 6d677ade0535424b090081fb1cfe3e5be82a7d66..0a568fc9407a56dee04915cd25c9d225599c4f62 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCGS solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index 4b749fd338f9d5836f4a88a70a13b94f5c276381..87286242f8997e2878f00d0caaf5c32aa548cae2 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverTFQMR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index 08520c11ba25da46db5774e91ff8cc0458dc4e40..01e44da87853f1d7c0f2c2a6136ed454abc48682 100644 (file)
@@ -86,9 +86,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverTCQMR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index 5fc7cd7efc8c7f530f5b294374f025d956da464a..3197a839fe8a3e88950470e10fc5e969a1c998e5 100644 (file)
@@ -82,9 +82,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverCR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index c22613b52897fc486c3525d2ee2b5c13352c3359..5f8a6c81d9947cb26aa4cfe749ab054ef5a60899 100644 (file)
@@ -86,9 +86,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverLSQR solver(control);
     PETScWrappers::PreconditionJacobi preconditioner(A);
index b0c4032d39fdeda14e3613276c96db1552277f10..79a6e9d43e0625684e249950e079cf8dc801cb1a 100644 (file)
@@ -75,9 +75,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  f(dim);
     PETScWrappers::Vector  u(dim);
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     PETScWrappers::SolverPreOnly solver(control);
 
index 50909db7cd45ab328cca144645b203a0ebf405a5..5deb2efdd823fd4aa202e53ca5dfad1ed245dd45 100644 (file)
@@ -73,9 +73,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  u(dim);
     u = 0.;
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
     PETScWrappers::Vector t=u;
 
     PETScWrappers::SolverPreOnly solver(control);
index 4670fb127e8abe043d5d93743af796da99cbbe36..3432a94fc6b35bb826c574616cca85fec0e801ba 100644 (file)
@@ -57,9 +57,7 @@ int main(int argc, char **argv)
     PETScWrappers::Vector  u(dim);
     u = 0.;
     f = 1.;
-    A.compress (VectorOperation::add);
-    f.compress (VectorOperation::add);
-    u.compress (VectorOperation::add);
+    A.compress (VectorOperation::insert);
 
     SolverControl cn;
     PETScWrappers::SparseDirectMUMPS solver(cn);
index 3d39fb34b2856a8cb437b06c09248a8a3318e630..41b4d0d97a8f44460929f9db97e0bd517a3354df 100644 (file)
@@ -33,13 +33,13 @@ void test ()
     for (unsigned int j=0; j<=i; ++j)
       m.set (i,j, i+2*j);
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   PETScWrappers::SparseMatrix m2(s,s,s);
   m2.set(0,1,5.0);
   for (unsigned int i=0; i<m2.m(); ++i)
     m2.set(i,i,1.0+ i);
-  m2.compress (VectorOperation::add);
+  m2.compress (VectorOperation::insert);
 
   // we now print the matrix m,
   // print after adding m2, and then subtract m2 again
index 92329cfc3ea523c607b74ee58fc5f2c921b21e96..acab595b6ffccb51ec2a90eff73284637aa1dce5 100644 (file)
@@ -34,7 +34,7 @@ void test ()
       m.set (i,j, i+2*j);
 
 
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
 
   {
 
index e9091f7e8cf21451bf00e4901b5f21ad6bd708ca..740bc0d846527d43126cfa06a691d3fb68632cba 100644 (file)
@@ -29,7 +29,7 @@ void test ()
   m.set (0,0,1);
   m.set (1,1,2);
   m.set (1,2,3);
-  m.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
   PETScWrappers::SparseMatrix::const_iterator i = m.begin();
   deallog << i->row() << ' ' << i->column() << ' ' << i->value() << std::endl;
   ++i;
index e510513ac2c97d23baddf355f7761cdc6f688d08..d54aa90c13f2cdc784b0b38fc1c6f26368709689 100644 (file)
@@ -36,9 +36,8 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // w:=Mv
   m.vmult (w,v);
index bc98ad15b8c79dc26e72b7a9920fba212a29471d..39d9cc574407c4f0db3e36e6db3a99bea5f1edf6 100644 (file)
@@ -36,9 +36,8 @@ void test (PETScWrappers::Vector &v,
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // w:=Mv
   m.Tvmult (w,v);
index 9f2df38acd8b834dab75c7007abbf79ecf03227d..ec264536bace1a3d415914d1c30a2ff4ba109d1a 100644 (file)
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // w:=Mv
   m.vmult_add (w,v);
index 287cdd6a1e4fc310ed8d37cb7979b00ddf30c398..ce746f37eb3ba04ed1af221529f15e506747d773 100644 (file)
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // w:=Mv
   m.Tvmult_add (w,v);
index 49fa71e3ac6a43e27c7e0ccb0b0590aa20e0eb53..d08b2642f1179923013d53d0c714ae774acedfc1 100644 (file)
@@ -39,9 +39,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // <w,Mv>
   const PetscScalar s = m.matrix_scalar_product (w,v);
index d589fe11fd360b27a9a263baf0b3a019097c45ff..a84f749f095f4f0f83bf374a62e5271bb9ad6266 100644 (file)
@@ -35,8 +35,8 @@ void test (PETScWrappers::Vector &v)
   for (unsigned int i=0; i<v.size(); ++i)
     v(i) = i;
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
 
   // <w,Mv>
   const PetscScalar s = m.matrix_norm_square (v);
index 26041ff6bc0f3abf70868a19e5201d000ea98a19..35fc56c8705016c2f0c0e55933a3547ab5dcd30a 100644 (file)
@@ -40,9 +40,9 @@ void test (PETScWrappers::Vector &v,
       w(i) = i+1;
     }
 
-  m.compress (VectorOperation::add);
-  v.compress (VectorOperation::add);
-  w.compress (VectorOperation::add);
+  m.compress (VectorOperation::insert);
+  v.compress (VectorOperation::insert);
+  w.compress (VectorOperation::insert);
 
   // x=w-Mv
   const double s = m.residual (x, v, w);

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.