Prescribed Solution Post-processor

Author: lhy11009

benchmarks/annulus/plugin/annulus.cc

@@ -663,6 +663,7 @@ namespace aspect
         }
     };
 
+    // todo_post
     /**
       * A postprocessor that visualizes the analytical solution.
       */

include/aspect/postprocess/visualization/prescribed_solution.h

@@ -0,0 +1,84 @@
+/*
+  Copyright (C) 2015 - 2021 by the authors of the ASPECT code.
+
+  This file is part of ASPECT.
+
+  ASPECT is free software; you can redistribute it and/or modify
+  it under the terms of the GNU General Public License as published by
+  the Free Software Foundation; either version 2, or (at your option)
+  any later version.
+
+  ASPECT is distributed in the hope that it will be useful,
+  but WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+  GNU General Public License for more details.
+
+  You should have received a copy of the GNU General Public License
+  along with ASPECT; see the file LICENSE.  If not see
+  <http://www.gnu.org/licenses/>.
+*/
+
+
+#ifndef _aspect_postprocess_visualization_prescribed_solution_h
+#define _aspect_postprocess_visualization_prescribed_solution_h
+
+#include <aspect/postprocess/visualization.h>
+#include <aspect/simulator_access.h>
+
+#include <deal.II/numerics/data_out.h>
+
+
+namespace aspect
+{
+  namespace Postprocess
+  {
+    namespace VisualizationPostprocessors
+    {
+      // todo:move the class to here later
+      /**
+      * A class derived from DataPostprocessor that outputs prescribed solution
+      * values and indicators
+      */
+      template <int dim>
+      class PrescribedSolutionPostprocessor
+        : public DataPostprocessor<dim>,
+          public SimulatorAccess<dim>,
+          public Interface<dim>
+      {
+        public:
+          PrescribedSolutionPostprocessor ();
+
+          std::vector<std::string>
+          get_names () const override;
+
+          std::vector<DataComponentInterpretation::DataComponentInterpretation>
+          get_data_component_interpretation () const override;
+
+          UpdateFlags
+          get_needed_update_flags () const override;
+
+          void
+          evaluate_vector_field(const DataPostprocessorInputs::Vector<dim> &input_data,
+                                std::vector<Vector<double>> &computed_quantities) const  override;
+
+          /**
+           * Declare the parameters this class takes through input files.
+           */
+          static
+          void
+          declare_parameters (ParameterHandler &prm);
+
+          /**
+           * Read the parameters this class declares from the parameter file.
+           */
+          void
+          parse_parameters (ParameterHandler &prm) override;
+
+        private:
+          std::vector<std::string> model_names;
+      };
+    }
+  }
+}
+
+#endif

include/aspect/prescribed_solution/interface.h

@@ -71,6 +71,7 @@ namespace aspect
      *
      * @ingroup PrescribedSolution
      */
+    // todo_post
     template <int dim>
     class Manager : public Plugins::ManagerBase<Interface<dim>>, public SimulatorAccess<dim>
     {

include/aspect/simulator_access.h

@@ -152,6 +152,11 @@ namespace aspect
     template <int dim> class Manager;
   }
 
+  namespace PrescribedSolution
+  {
+    template <int dim> class  Manager;
+  }
+
   /**
    * SimulatorAccess is a base class for different plugins like postprocessors.
    * It provides access to the various variables of the main class that
@@ -1020,6 +1025,13 @@ namespace aspect
       const StokesMatrixFreeHandler<dim> &
       get_stokes_matrix_free () const;
 
+      /**
+       * Return a reference to the PrescribedSolution::Manager that manages the
+       * Prescribed solution plugins.
+       */
+      const PrescribedSolution::Manager<dim> &
+      get_prescribed_solution () const;
+
       /**
        * Compute the angular momentum and other rotation properties
        * of the velocities in the given solution vector.