Fixed some parallel compilation problems and added FETI-DP interface.

AMDiS/CMakeLists.txt

@@ -199,6 +199,7 @@ if(ENABLE_PARALLEL_DOMAIN)
                 ${SOURCE_DIR}/parallel/ParMetisPartitioner.cc
 		${SOURCE_DIR}/parallel/PetscProblemStat.cc
 		${SOURCE_DIR}/parallel/PetscSolver.cc
+		${SOURCE_DIR}/parallel/PetscSolverFeti.cc
 		${SOURCE_DIR}/parallel/PetscSolverGlobalMatrix.cc
 		${SOURCE_DIR}/parallel/PetscSolverSchur.cc
 		${SOURCE_DIR}/parallel/StdMpi.cc

AMDiS/src/parallel/ParallelProblemStatBase.h

@@ -52,10 +52,6 @@ namespace AMDiS {
     MeshDistributor *meshDistributor;     
   };
 
-#ifdef HAVE_PARALLEL_DOMAIN_AMDIS
-  typedef ProblemStatSeq ProblemStat;
-#endif
-
 }
 
 #endif

AMDiS/src/parallel/PetscProblemStat.h

@@ -27,6 +27,7 @@
 #include "Global.h"
 #include "parallel/ParallelProblemStatBase.h"
 #include "parallel/PetscSolver.h"
+#include "parallel/PetscSolverFeti.h"
 #include "parallel/PetscSolverGlobalMatrix.h"
 #include "parallel/PetscSolverSchur.h"
 
@@ -46,7 +47,13 @@ namespace AMDiS {
 #ifdef HAVE_PETSC_DEV
 	petscSolver = new PetscSolverSchur();
 #else
-	ERROR_EXIT("Petsc schur complement solver is only supported when petsc-dev is used!\n");
+	ERROR_EXIT("PETSc schur complement solver is only supported when petsc-dev is used!\n");
+#endif
+      } else if (name == "petsc-feti") {
+#ifdef HAVE_PETSC_DEV
+	petscSolver = new PetscSolverFeti();
+#else
+	ERROR_EXIT("PETSc FETI-DP solver is only supported when petsc-dev is used!\n");
 #endif
       } else if (name == "petsc" || name == "") {
 	petscSolver = new PetscSolverGlobalMatrix();
@@ -75,6 +82,10 @@ namespace AMDiS {
 
   typedef PetscProblemStat ParallelProblemStat;
 
+#ifdef HAVE_PARALLEL_DOMAIN_AMDIS
+  typedef PetscProblemStat ProblemStat;
+#endif
+
 } // namespace AMDiS
 
 #endif

AMDiS/src/parallel/PetscSolverFeti.cc

+//
+// Software License for AMDiS
+//
+// Copyright (c) 2010 Dresden University of Technology 
+// All rights reserved.
+// Authors: Simon Vey, Thomas Witkowski et al.
+//
+// This file is part of AMDiS
+//
+// See also license.opensource.txt in the distribution.
+
+
+#include "parallel/PetscSolverFeti.h"
+#include "parallel/StdMpi.h"
+#include "parallel/MpiHelper.h"
+
+namespace AMDiS {
+
+  using namespace std;
+
+
+#ifdef HAVE_PETSC_DEV 
+  void PetscSolverFeti::fillPetscMatrix(Matrix<DOFMatrix*> *mat, SystemVector *vec)
+  {
+    FUNCNAME("PetscSolverFeti::fillPetscMatrix()");   
+  }
+
+
+  void PetscSolverFeti::solvePetscMatrix(SystemVector &vec, AdaptInfo *adaptInfo)
+  {
+    FUNCNAME("PetscSolverFeti::solvePetscMatrix()");
+  }
+#endif
+
+}

AMDiS/src/parallel/PetscSolverFeti.h

+// ============================================================================
+// ==                                                                        ==
+// == AMDiS - Adaptive multidimensional simulations                          ==
+// ==                                                                        ==
+// ==  http://www.amdis-fem.org                                              ==
+// ==                                                                        ==
+// ============================================================================
+//
+// Software License for AMDiS
+//
+// Copyright (c) 2010 Dresden University of Technology 
+// All rights reserved.
+// Authors: Simon Vey, Thomas Witkowski et al.
+//
+// This file is part of AMDiS
+//
+// See also license.opensource.txt in the distribution.
+
+
+
+/** \file PetscSolverFeti.h */
+
+#include "parallel/PetscSolver.h"
+
+#ifndef AMDIS_PETSC_SOLVER_FETI_H
+#define AMDIS_PETSC_SOLVER_FETI_H
+
+namespace AMDiS {
+
+  using namespace std;
+
+
+#ifdef HAVE_PETSC_DEV
+  class PetscSolverFeti : public PetscSolver
+  {
+  public:
+    PetscSolverFeti() 
+      : PetscSolver()
+    {}
+
+    void fillPetscMatrix(Matrix<DOFMatrix*> *mat, SystemVector *vec);
+
+    void solvePetscMatrix(SystemVector &vec, AdaptInfo *adaptInfo);
+  };
+#endif
+
+}
+
+#endif