Fixed bug with dirichlet values.

AMDiS/src/DOFMatrix.cc

@@ -537,7 +537,7 @@ namespace AMDiS {
     inserter_type &ins = *inserter;  
     for (std::set<int>::iterator it = dirichletDofs.begin(); 
 	 it != dirichletDofs.end(); ++it)
-	ins[*it][*it] = 1.0;
+      ins[*it][*it] = 1.0;
   }
 
 

AMDiS/src/DirichletBC.cc

@@ -81,4 +81,13 @@ namespace AMDiS {
     }
   }
 
+
+  void DirichletBC::initVector(DOFVectorBase<double>* vec)
+  {
+    FUNCNAME("DirichletBC::initVector()");
+
+    if (dynamic_cast<DOFVector<double>*>(vec))
+      dynamic_cast<DOFVector<double>*>(vec)->getDirichletValues().clear();
+  }
+
 }

AMDiS/src/DirichletBC.h

@@ -68,6 +68,9 @@ namespace AMDiS {
 			       const BoundaryType* localBound,
 			       int nBasFcts);
 
+    ///
+    void initVector(DOFVectorBase<double>*);
+
     /// Implementation of BoundaryCondition::boundResidual().
     double boundResidual(ElInfo*, 
 			 DOFMatrix *,

AMDiS/src/parallel/PetscSolverGlobalMatrix.cc

@@ -11,6 +11,7 @@
 
 #include <mpi.h>
 #include "DirichletBC.h"
+#include "DOFVector.h"
 #include "parallel/PetscSolverGlobalMatrix.h"
 #include "parallel/StdMpi.h"
 #include "parallel/MpiHelper.h"
@@ -750,7 +751,7 @@ namespace AMDiS {
     for (dofIt.reset(); !dofIt.end(); ++dofIt) {
       
       DegreeOfFreedom dof = dofIt.getDOFIndex();
-      
+
       if (rankOnly && !(*interiorMap)[rowComp].isRankDof(dof))
 	continue;