That's it for today.

AMDiS/src/ProblemStat.cc

@@ -870,6 +870,17 @@ namespace AMDiS {
     INFO(info, 8)("buildAfterCoarsen needed %.5f seconds\n", 
 		  TIME_USED(first, clock()));
 #endif
+
+#if 0
+    VtkWriter::writeFile(rhs, "rhs0.vtu");
+    for (int i = 0; i < nComponents; i++) {
+      double s = rhs->getDOFVector(i)->sum();
+      s /= -static_cast<double>(rhs->getDOFVector(i)->getUsedSize());
+      MSG("MOVE RHS: %f\n", s);
+      add(*(rhs->getDOFVector(i)), s, *(rhs->getDOFVector(i)));
+    }
+    VtkWriter::writeFile(rhs, "rhs1.vtu");
+#endif
   }
 
 

AMDiS/src/parallel/FeSpaceMapping.cc

@@ -41,6 +41,8 @@ namespace AMDiS {
 
     if (overlap)
       computeNonLocalIndices();
+
+    computeMatIndex(0);
   }
 
 
@@ -89,6 +91,52 @@ namespace AMDiS {
   }
 
 
+  void GlobalDofMap::computeMatIndex(int offset)
+  {
+    FUNCNAME("GlobalDofMap::computeMatIndex()");
+
+    map<DegreeOfFreedom, int> dofToMatIndex;
+
+    for (DofMapping::iterator it = dofMap.begin(); it != dofMap.end(); ++it) {
+      if (!nonRankDofs.count(it->first)) {
+	int globalMatIndex = it->second - rStartDofs + offset;
+	dofToMatIndex[it->first] = globalMatIndex;
+      }
+    }
+
+    if (!overlap)
+      return;
+
+    TEST_EXIT_DBG(sendDofs != NULL && recvDofs != NULL)
+      ("No communicator given!\n");
+
+    StdMpi<vector<DegreeOfFreedom> > stdMpi(*mpiComm);
+    for (DofComm::Iterator it(*sendDofs, feSpace); !it.end(); it.nextRank()) {
+      vector<DegreeOfFreedom> sendGlobalDofs;
+
+      for (; !it.endDofIter(); it.nextDof())
+	if (dofMap.count(it.getDofIndex()))
+	  sendGlobalDofs.push_back(dofToMatIndex[it.getDofIndex()]);
+
+      stdMpi.send(it.getRank(), sendGlobalDofs);
+    }
+
+    for (DofComm::Iterator it(*recvDofs, feSpace); !it.end(); it.nextRank())
+      stdMpi.recv(it.getRank());
+
+    stdMpi.startCommunication();
+
+    {
+      int i = 0;
+      for (DofComm::Iterator it(*recvDofs, feSpace); !it.end(); it.nextRank()) 
+	for (; !it.endDofIter(); it.nextDof())
+	  if (dofMap.count(it.getDofIndex()))
+	    dofToMatIndex[it.getDofIndex()] = 
+	      stdMpi.getRecvData(it.getRank())[i++];
+    }
+  }
+
+
   void GlobalDofMap::print()
   {
     FUNCNAME("GlobalDofMap::print()");

AMDiS/src/parallel/FeSpaceMapping.h

@@ -46,6 +46,8 @@ namespace AMDiS {
     GlobalDofMap(MPI::Intracomm* m)
       : mpiComm(m),
 	feSpace(NULL),
+	sendDofs(NULL),
+	recvDofs(NULL),
 	nRankDofs(0),
 	nOverallDofs(0),
 	rStartDofs(0),
@@ -102,6 +104,8 @@ namespace AMDiS {
 
     void computeNonLocalIndices();
 
+    void computeMatIndex(int offset);
+
     void print();
 
     void setFeSpace(const FiniteElemSpace *fe)

AMDiS/src/parallel/MeshDistributor.cc

@@ -31,6 +31,7 @@
 #include "parallel/DofComm.h"
 #include "io/ElementFileWriter.h"
 #include "io/MacroInfo.h"
+#include "io/MacroWriter.h"
 #include "io/VtkWriter.h"
 #include "Mesh.h"
 #include "Traverse.h"

AMDiS/src/parallel/PetscSolverFeti.cc

@@ -260,8 +260,7 @@ namespace AMDiS {
     primalDofMap[feSpace].update();
 
     MSG("nRankPrimals = %d   nOverallPrimals = %d\n", 
-	primalDofMap[feSpace].nRankDofs, 
-	primalDofMap[feSpace].nOverallDofs);
+	primalDofMap[feSpace].nRankDofs, primalDofMap[feSpace].nOverallDofs);
 
     TEST_EXIT_DBG(primals.size() == primalDofMap[feSpace].size())
       ("Number of primals %d, but number of global primals on this rank is %d!\n",
@@ -433,8 +432,7 @@ namespace AMDiS {
     
     localDofMap[feSpace].update(false);
 
-    TEST_EXIT_DBG(nLocalInterior + 
-		  primalDofMap[feSpace].size() + 
+    TEST_EXIT_DBG(nLocalInterior + primalDofMap[feSpace].size() + 
 		  dualDofMap[feSpace].size() == 
 		  static_cast<unsigned int>(admin->getUsedDofs()))
       ("Should not happen!\n");

AMDiS/src/parallel/PetscSolverGlobalMatrix.cc

@@ -211,7 +211,6 @@ namespace AMDiS {
     // Print iteration counter and residual norm of the solution.
     printSolutionInfo(adaptInfo);
 
-
     // === Destroy PETSc's variables. ===
 
     VecDestroy(&petscRhsVec);