ASPECT
stokes_residual.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2015 - 2019 by the authors of the ASPECT code.
3 
4  This file is part of ASPECT.
5 
6  ASPECT is free software; you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation; either version 2, or (at your option)
9  any later version.
10 
11  ASPECT is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU General Public License for more details.
15 
16  You should have received a copy of the GNU General Public License
17  along with ASPECT; see the file LICENSE. If not see
18  <http://www.gnu.org/licenses/>.
19 */
20 
21 
22 #ifndef _aspect_postprocess_stokes_residual_h
23 #define _aspect_postprocess_stokes_residual_h
24 
27 
28 #include <deal.II/base/data_out_base.h>
29 
30 
31 namespace aspect
32 {
33  namespace Postprocess
34  {
35 
41  template <int dim>
42  class StokesResidual : public Interface<dim>, public ::aspect::SimulatorAccess<dim>
43  {
44  public:
48  StokesResidual ();
49 
53  void initialize() override;
54 
58  std::pair<std::string,std::string>
59  execute (TableHandler &statistics) override;
60 
64  void save (std::map<std::string, std::string> &status_strings) const override;
65 
69  void load (const std::map<std::string, std::string> &status_strings) override;
70 
75  template <class Archive>
76  void serialize (Archive &ar, const unsigned int version);
77 
78  private:
82  struct DataPoint
83  {
84  double time;
85  unsigned int solve_index;
86  std::vector<double> values;
87 
88  template <class Archive>
89  void serialize (Archive &ar, const unsigned int version);
90  };
91 
95  void stokes_solver_callback (const SolverControl &solver_control_cheap,
96  const SolverControl &solver_control_expensive);
97 
101  std::vector<DataPoint> entries;
102 
103  };
104  }
105 }
106 
107 
108 #endif
void serialize(Archive &ar, const unsigned int version)
void stokes_solver_callback(const SolverControl &solver_control_cheap, const SolverControl &solver_control_expensive)
std::vector< DataPoint > entries
void serialize(Archive &ar, const unsigned int version)
std::pair< std::string, std::string > execute(TableHandler &statistics) override
void save(std::map< std::string, std::string > &status_strings) const override
void load(const std::map< std::string, std::string > &status_strings) override