So

AMDiS/src/parallel/PetscSolverFeti.cc

@@ -240,9 +240,13 @@ namespace AMDiS {
 
     if (enableStokesMode) {
       MSG("WARNING: MAKE THIS MORE GENERAL!!!!!\n");
-      subdomain->setCoarseSpaceDofMapping(&primalDofMap, 0);
+      subdomain->setCoarseSpaceDofMapping(&interfaceDofMap, 0);
+      //      subdomain->setCoarseSpaceDofMapping(&primalDofMap, 0);
       subdomain->setCoarseSpaceDofMapping(&primalDofMap, 1);
-      subdomain->setCoarseSpaceDofMapping(&interfaceDofMap, 2);
+      subdomain->setCoarseSpaceDofMapping(&primalDofMap, 2);
+      subdomain->setCoarseSpaceDofMapping(&primalDofMap, 3);
+      subdomain->setCoarseSpaceDofMapping(&primalDofMap, 4);
+
     } else {	  
       subdomain->setCoarseSpaceDofMapping(&primalDofMap); 
     }
@@ -1556,9 +1560,11 @@ namespace AMDiS {
 	      valuesLocal.push_back(value(*icursor));
 	    }
 
-	    int rowIndex = dualDofMap.getLocalMatIndex(rowComponent, *cursor);		
-	    MatSetValues(mat_duals_duals, 1, &rowIndex, colsLocal.size(),
-			 &(colsLocal[0]), &(valuesLocal[0]), INSERT_VALUES);
+	    {
+	      int rowIndex = dualDofMap.getLocalMatIndex(rowComponent, *cursor);		
+	      MatSetValues(mat_duals_duals, 1, &rowIndex, colsLocal.size(),
+			   &(colsLocal[0]), &(valuesLocal[0]), INSERT_VALUES);
+	    }
 	    break;
 
 

AMDiS/src/parallel/PetscSolverFetiOperators.cc

@@ -250,6 +250,7 @@ namespace AMDiS {
       yvec = y_lagrange;
 
       VecCopy(x_interface, y_interface);
+      //      VecScale(y_interface, 0.2519);
     } else {
       xvec = x;
       yvec = y;