move to cpp: TimeDiscretizationBuilder

NumLib/ODESolver/TimeDiscretizationBuilder.cpp

+/**
+ * \copyright
+ * Copyright (c) 2012-2016, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#include "TimeDiscretizationBuilder.h"
+
+#include "BaseLib/Error.h"
+
+namespace NumLib
+{
+std::unique_ptr<TimeDiscretization> createTimeDiscretization(
+    BaseLib::ConfigTree const& config)
+{
+    using T = std::unique_ptr<TimeDiscretization>;
+
+    //! \ogs_file_param{process__time_discretization__type}
+    auto const type = config.getConfigParameter<std::string>("type");
+
+    if (type == "BackwardEuler")
+    {
+        using ConcreteTD = BackwardEuler;
+        return T(new ConcreteTD);
+    }
+    else if (type == "ForwardEuler")
+    {
+        using ConcreteTD = ForwardEuler;
+        return T(new ConcreteTD);
+    }
+    else if (type == "CrankNicolson")
+    {
+        //! \ogs_file_param{process__time_discretization__CrankNicolson__theta}
+        auto const theta = config.getConfigParameter<double>("theta");
+        using ConcreteTD = CrankNicolson;
+        return T(new ConcreteTD(theta));
+    }
+    else if (type == "BackwardDifferentiationFormula")
+    {
+        //! \ogs_file_param{process__time_discretization__BackwardDifferentiationFormula__order}
+        auto const order = config.getConfigParameter<unsigned>("order");
+        using ConcreteTD = BackwardDifferentiationFormula;
+        return T(new ConcreteTD(order));
+    }
+    else
+    {
+        OGS_FATAL("Unrecognized time discretization type `%s'", type.c_str());
+    }
+}
+}

NumLib/ODESolver/TimeDiscretizationBuilder.h

@@ -19,42 +19,7 @@ namespace NumLib
 {
 
 std::unique_ptr<TimeDiscretization> createTimeDiscretization(
-    BaseLib::ConfigTree const& config)
-{
-    using T = std::unique_ptr<TimeDiscretization>;
-
-    //! \ogs_file_param{process__time_discretization__type}
-    auto const type = config.getConfigParameter<std::string>("type");
-
-    if (type == "BackwardEuler")
-    {
-        using ConcreteTD = BackwardEuler;
-        return T(new ConcreteTD);
-    }
-    else if (type == "ForwardEuler")
-    {
-        using ConcreteTD = ForwardEuler;
-        return T(new ConcreteTD);
-    }
-    else if (type == "CrankNicolson")
-    {
-        //! \ogs_file_param{process__time_discretization__CrankNicolson__theta}
-        auto const theta = config.getConfigParameter<double>("theta");
-        using ConcreteTD = CrankNicolson;
-        return T(new ConcreteTD(theta));
-    }
-    else if (type == "BackwardDifferentiationFormula")
-    {
-        //! \ogs_file_param{process__time_discretization__BackwardDifferentiationFormula__order}
-        auto const order = config.getConfigParameter<unsigned>("order");
-        using ConcreteTD = BackwardDifferentiationFormula;
-        return T(new ConcreteTD(order));
-    }
-    else
-    {
-        OGS_FATAL("Unrecognized time discretization type `%s'", type.c_str());
-    }
-}
+    BaseLib::ConfigTree const& config);
 }
 
 #endif  // NUMLIB_TIMEDISCRETIZATION_BUILDER_H