Merge pull request #28 from nicola-giuliani/bug_fix_gcc

Bug fix gcc

Author: luca-heltai

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_fma.cc

@@ -723,7 +723,7 @@ void BEMFMA<dim>::direct_integrals()
 
                             //MAGARI USARE FEVALUES CON IL DOFHANDLER CRETINO DISCONTINUO E IL MAPPING bem_fma
                             Point<dim> D;
-                            double s;
+                            double s=0.;
 
                             const Tensor<1, dim> R =  this->quadPoints.at(cell)[*pos] - support_points[nodeIndex];
                             LaplaceKernel::kernels(R, D, s);
@@ -800,7 +800,7 @@ void BEMFMA<dim>::direct_integrals()
                         for (unsigned int q=0; q<singular_quadrature->size(); ++q)
                           {
                             Point<dim> D;
-                            double s;
+                            double s=0.;
 
                             const Tensor<1, dim> R = singular_q_points[q] - support_points[nodeIndex];
                             LaplaceKernel::kernels(R, D, s);
@@ -978,7 +978,7 @@ void BEMFMA<dim>::direct_integrals()
                     // here we compute the distance R between the node and the quad point
                     const Tensor<1,dim> R =  quadPoints[cell][*pos] - support_points[nodeIndex];
                     Point<dim> D;
-                    double s;
+                    double s=0.;
 
                     LaplaceKernel::kernels(R, D, s);
 

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);
 }