Add some more debug functions.

AMDiS/src/ArhReader.cc

@@ -46,7 +46,7 @@ namespace AMDiS {
       file.read(reinterpret_cast<char*>(&nStructureCodes), 4);
       file.read(reinterpret_cast<char*>(&codeSize), 4);
 
-      vector<unsigned long> structureCode(nStructureCodes);
+      vector<unsigned long int> structureCode(nStructureCodes);
       file.read(reinterpret_cast<char*>(&(structureCode[0])), 8 * nStructureCodes);
 
       MeshStructure elementStructure;

AMDiS/src/Debug.cc

@@ -557,6 +557,28 @@ namespace AMDiS {
     }
 
 
+    void exportDofVectorByCoords(const DOFVector<double>* vec, 
+				 std::string filename)
+    {
+      DOFVector<WorldVector<double> > coords(vec->getFeSpace(), "dofCoords");
+      vec->getFeSpace()->getMesh()->getDofIndexCoords(vec->getFeSpace(), coords);
+      int dim = vec->getFeSpace()->getMesh()->getDim();
+
+      std::ofstream file;
+      file.open(filename.c_str());
+      file << vec->getUsedSize() << "\n";
+
+      DOFIterator<WorldVector<double> > it(&coords, USED_DOFS);
+      for (it.reset(); !it.end(); ++it) {
+        for (int i = 0; i < dim; i++)
+	  file << (*it)[i] << " ";
+        file << (*vec)[it.getDOFIndex()] << "\n";
+      }
+    
+      file.close();
+    }
+
+
     void createSortedDofs(Mesh *mesh, ElementIdxToDofs &elMap)
     {
       FUNCNAME("debug::dbgCreateElementMap()");

AMDiS/src/Debug.h

@@ -141,6 +141,9 @@ namespace AMDiS {
 
     int getLocalNeighbourIndex(Mesh *mesh, int elIndex, int neighIndex);
 
+    void exportDofVectorByCoords(const DOFVector<double>* vec, 
+				 std::string filename);
+
     /** \brief
      * Traverse a mesh and store for each element all its vertex DOFs in local sorted 
      * order (by values).

AMDiS/src/SecondOrderTerm.cc

@@ -194,7 +194,7 @@ namespace AMDiS {
   {
     const DimVec<WorldVector<double> > &grdLambda = elInfo->getGrdLambda();
     for (int iq = 0; iq < nPoints; iq++)
-      l1lt(grdLambda, *(LALt[iq]), (*f)(vecAtQPs[iq]));
+      l1lt(grdLambda, *(LALt[iq]), (*f)(vecAtQPs[iq]));    
   }
 
   void VecAtQP_SOT::eval(int nPoints,

AMDiS/src/ZeroOrderAssembler.cc

@@ -68,9 +68,9 @@ namespace AMDiS {
       newAssembler = new StandardZOA(op, assembler, quad);
     } else {
       if (pwConst)
-  	newAssembler = new PrecalcZOA(op, assembler, quad);
+	newAssembler = new PrecalcZOA(op, assembler, quad);
       else
-  	newAssembler = new FastQuadZOA(op, assembler, quad);      
+	newAssembler = new FastQuadZOA(op, assembler, quad);      
     }
 
     subAssemblers->push_back(newAssembler);
@@ -130,7 +130,7 @@ namespace AMDiS {
 	for (int i = 0; i < nRow; i++) {
 	  double psival = (*(psi->getPhi(i)))(quadrature->getLambda(iq));
 	  for (int j = 0; j < nCol; j++)
-	    mat[i][j] += quadrature->getWeight(iq) * c[iq] * psival * phival[j];
+	    mat[i][j] += quadrature->getWeight(iq) * c[iq] * psival * phival[j];	  
 	}
       }
     }