Add nullspace support on the level of PETSc solver object.

AMDiS/src/io/MacroInfo.cc

@@ -167,7 +167,8 @@ namespace AMDiS {
       TEST_EXIT(i_key >= 0)
 	("macro file %s must not contain key %s on line %d\n",
 	 filename.c_str(), key, line_no);
-      TEST_EXIT(!key_def[i_key])("key %s defined second time on line %d in file %s\n");
+      TEST_EXIT(!key_def[i_key])
+	("key %s defined second time on line %d in file %s\n");
 
       sort_key[n_keys++] = i_key;
       key_def[i_key] = true;

AMDiS/src/parallel/PetscSolver.h

@@ -127,6 +127,19 @@ namespace AMDiS {
       nullspace.push_back(vec);
     }
 
+    /// Sets the nullspace to be constant for some specific components.
+    void setConstantNullspace(vector<int> &components)
+    {
+      constNullspaceComponent = components;
+    }
+
+    /// Sets the nullspace to be constant for a specific component.
+    void setConstantNullspace(int component)
+    {
+      constNullspaceComponent.clear();
+      constNullspaceComponent.push_back(component);
+    }
+
     inline bool isCoarseSpace(const FiniteElemSpace *feSpace, 
 			      DegreeOfFreedom dof)
     {