ASPECT
rk_2.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2015 - 2021 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 #ifndef _aspect_particle_integrator_rk_2_h
22 #define _aspect_particle_integrator_rk_2_h
23 
25 
27 
28 
29 namespace aspect
30 {
31  namespace Particle
32  {
33  namespace Integrator
34  {
41  template <int dim>
42  class RK2 : public Interface<dim>, public SimulatorAccess<dim>
43  {
44  public:
45  RK2();
46 
51  virtual
52  void
53  initialize () override;
54 
72  void
73  local_integrate_step(const typename ParticleHandler<dim>::particle_iterator &begin_particle,
74  const typename ParticleHandler<dim>::particle_iterator &end_particle,
75  const std::vector<Tensor<1,dim>> &old_velocities,
76  const std::vector<Tensor<1,dim>> &velocities,
77  const double dt) override;
78 
89  bool new_integration_step() override;
90 
94  static
95  void
96  declare_parameters (ParameterHandler &prm);
97 
101  virtual
102  void
103  parse_parameters (ParameterHandler &prm) override;
104 
116  static const unsigned int n_integrator_properties = dim;
117 
118  private:
123  unsigned int integrator_substep;
124 
129 
135  };
136 
137  }
138  }
139 }
140 
141 #endif
void local_integrate_step(const typename ParticleHandler< dim >::particle_iterator &begin_particle, const typename ParticleHandler< dim >::particle_iterator &end_particle, const std::vector< Tensor< 1, dim >> &old_velocities, const std::vector< Tensor< 1, dim >> &velocities, const double dt) override
unsigned int property_index_old_location
Definition: rk_2.h:128
unsigned int integrator_substep
Definition: rk_2.h:123
static void declare_parameters(ParameterHandler &prm)
virtual void initialize() override
static const unsigned int n_integrator_properties
Definition: rk_2.h:116
virtual void parse_parameters(ParameterHandler &prm) override