initial commit.

Author: mdh266

CMakeLists.txt

@@ -0,0 +1,50 @@
+##
+#  CMake script for the step-42 tutorial program:
+##
+
+# Set the name of the project and target:
+SET(TARGET "LDGPoisson")
+
+# Declare all source files the target consists of:
+SET(TARGET_SRC
+  ${TARGET}.cc
+  Functions.cc
+  )
+
+# Define the output that should be cleaned:
+SET(CLEAN_UP_FILES *.vtu *.pvtu *.visit)
+
+# Usually, you will not need to modify anything beyond this point...
+
+CMAKE_MINIMUM_REQUIRED(VERSION 2.8.8)
+
+FIND_PACKAGE(deal.II 8.0 QUIET
+  HINTS ${deal.II_DIR} ${DEAL_II_DIR} ../ ../../ $ENV{DEAL_II_DIR}
+  )
+IF(NOT ${deal.II_FOUND})
+  MESSAGE(FATAL_ERROR "\n"
+    "*** Could not locate deal.II. ***\n\n"
+    "You may want to either pass a flag -DDEAL_II_DIR=/path/to/deal.II to cmake\n"
+    "or set an environment variable \"DEAL_II_DIR\" that contains this path."
+    )
+ENDIF()
+
+#
+# Are all dependencies fullfilled?
+#
+IF( NOT DEAL_II_WITH_MPI OR
+    NOT DEAL_II_WITH_P4EST OR
+    NOT DEAL_II_WITH_TRILINOS )
+  MESSAGE(FATAL_ERROR "
+Error! The deal.II library found at ${DEAL_II_PATH} was not configured with
+    DEAL_II_WITH_MPI = ON
+    DEAL_II_WITH_P4EST = ON
+    DEAL_II_WITH_TRILINOS = ON
+One or all of these are OFF in your installation but are required for this tutorial step."
+    )
+ENDIF()
+
+DEAL_II_INITIALIZE_CACHED_VARIABLES()
+PROJECT(${TARGET})
+
+DEAL_II_INVOKE_AUTOPILOT()

Functions.cc

@@ -0,0 +1,86 @@
+// @sect3{Functions.cc}
+// In this file we keep right hand side function, Dirichlet boundary
+// conditions and solution to our Poisson equation problem.  Since
+// these classes and functions have been discussed extensively in
+// the deal.ii tutorials we won't discuss them any further.
+#include <deal.II/base/function.h>
+#include <deal.II/base/tensor_function.h>
+#include <deal.II/lac/vector.h>
+#include <cmath>
+
+using namespace dealii;
+
+template <int dim>
+class RightHandSide : public Function<dim>
+{
+public:
+    RightHandSide() : Function<dim>(1)
+    {}
+
+    virtual double value(const Point<dim> &p,
+                         const unsigned int component = 0 ) const;
+};
+
+template <int dim>
+class DirichletBoundaryValues : public Function<dim>
+{
+public:
+    DirichletBoundaryValues() : Function<dim>(1)
+    {}
+
+    virtual double value(const Point<dim> &p,
+                         const unsigned int component = 0 ) const;
+};
+
+template<int dim>
+class TrueSolution : public Function<dim>
+{
+public:
+    TrueSolution() : Function<dim>(dim+1)
+    {}
+
+    virtual void vector_value(const Point<dim> & p,
+                              Vector<double> &valuess) const;
+};
+
+template <int dim>
+double
+RightHandSide<dim>::
+value(const Point<dim> &p,
+      const unsigned int ) const
+{
+    const double x = p[0];
+    const double y = p[1];
+    return 4*M_PI*M_PI*(cos(2*M_PI*y) - sin(2*M_PI*x));
+
+}
+
+template <int dim>
+double
+DirichletBoundaryValues<dim>::
+value(const Point<dim> &p,
+      const unsigned int ) const
+{
+    const double x = p[0];
+    const double y = p[1];
+    return cos(2*M_PI*y) -sin(2*M_PI*x) - x;
+}
+
+
+template <int dim>
+void
+TrueSolution<dim>::
+vector_value(const Point<dim> &p,
+             Vector<double> &values) const
+{
+    Assert(values.size() == dim+1,
+           ExcDimensionMismatch(values.size(), dim+1) );
+
+    double x = p[0];
+    double y = p[1];
+
+    values(0) = 1 + 2*M_PI*cos(2*M_PI*x);
+    values(1) = 2*M_PI*sin(2*M_PI*y);
+
+    values(2) = cos(2*M_PI*y) - sin(2*M_PI*x) - x;
+}