Fixed problem with new KSP solver interface.

AMDiS/src/parallel/ParallelDofMapping.cc

@@ -213,7 +213,7 @@ namespace AMDiS {
       levelData(NULL),
       dofComm(NULL),
       globalMapping(true),
-      needMatIndexFromGlobal(false),
+      needMatIndexFromGlobal(true),
       nRankDofs(1),
       nLocalDofs(1),
       nOverallDofs(1),

AMDiS/src/parallel/PetscSolverGlobalMatrix.cc

@@ -605,7 +605,7 @@ namespace AMDiS {
   }
 
 
-  void PetscSolverGlobalMatrix::initSolver(KSP ksp)
+  void PetscSolverGlobalMatrix::initSolver(KSP &ksp)
   {
     FUNCNAME("PetscSolverGlobalMatrix::initSolver()");
 

AMDiS/src/parallel/PetscSolverGlobalMatrix.h

@@ -88,7 +88,7 @@ namespace AMDiS {
      */
     void createFieldSplit(PC pc, string splitName, vector<int> &components);
 
-    virtual void initSolver(KSP ksp);
+    virtual void initSolver(KSP &ksp);
 
     virtual void exitSolver(KSP ksp);
 

AMDiS/src/parallel/PetscSolverNavierStokes.cc

@@ -18,7 +18,7 @@ namespace AMDiS {
   using namespace std;
 
   
-  void PetscSolverNavierStokes::initSolver(KSP ksp)
+  void PetscSolverNavierStokes::initSolver(KSP &ksp)
   {
     FUNCNAME("PetscSolverNavierStokes::initSolver()");
 

AMDiS/src/parallel/PetscSolverNavierStokes.h

@@ -37,7 +37,7 @@ namespace AMDiS {
     {}
 
   protected:
-    void initSolver(KSP ksp);
+    void initSolver(KSP &ksp);
 
     void initPreconditioner(PC pc);
   };