cleaned compress

Author: nicola-giuliani

include/constrained_matrix.h

@@ -96,7 +96,7 @@ void ConstrainedOperator<VEC,MATRIX>::vmult(VEC &dst, const VEC &src) const
   //        std::cout<<i<<" ("<<this_mpi_process<<")  "<<loc_src(i)<<std::endl;
 
   matrix.vmult(dst, src);
-  dst.compress(VectorOperation::insert);
+  // dst.compress(VectorOperation::insert);
   // IndexSet dummy(dst.locally_owned_elements());
 
   for (unsigned int i=0; i<src.size(); ++i)

source/bem_problem.cc

@@ -1172,6 +1172,7 @@ void BEMProblem<dim>::vmult(TrilinosWrappers::MPI::Vector &dst, const TrilinosWr
   // one
   if (!have_dirichlet_bc)
     vector_shift(dst, -dst.l2_norm());
+  dst.compress(VectorOperation::add);
 }