Fixed a bug when using openmp parallelization for assemlbing matrices.

AMDiS/src/DOFMatrix.cc

@@ -18,7 +18,6 @@
 namespace AMDiS {
 
   using namespace mtl;
-  using boost::lexical_cast;
 
   DOFMatrix *DOFMatrix::traversePtr = NULL;
 

AMDiS/src/Operator.cc

@@ -6,13 +6,8 @@
 #include "Quadrature.h"
 #include "OpenMP.h"
 
-#include "boost/lexical_cast.hpp"
-#include <fstream>
-
 namespace AMDiS {
 
-  using boost::lexical_cast;
-
   const Flag OperatorTerm::PW_CONST = 1;
   const Flag OperatorTerm::SYMMETRIC = 2;
 

AMDiS/src/SecondOrderAssembler.cc

@@ -188,13 +188,13 @@ namespace AMDiS {
 	tmpLALt[myRank][j] = new DimMat<double>(dim, NO_INIT);      
 
 #ifdef _OPENMP
-#pragma omp critical (dofIndexAccess)
+#pragma omp critical
 #endif
       {   
-      psiFast = updateFastQuadrature(psiFast, owner->getRowFESpace()->getBasisFcts(), 
-				     INIT_GRD_PHI);
-      phiFast = updateFastQuadrature(phiFast, owner->getRowFESpace()->getBasisFcts(), 
-				     INIT_GRD_PHI);
+        psiFast = updateFastQuadrature(psiFast, owner->getRowFESpace()->getBasisFcts(), 
+  				       INIT_GRD_PHI);
+        phiFast = updateFastQuadrature(phiFast, owner->getRowFESpace()->getBasisFcts(), 
+				       INIT_GRD_PHI);
       }
       firstCall = false;
     }