ASPECT
rk_4.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_4_h
22 #define _aspect_particle_integrator_rk_4_h
23 
25 
26 namespace aspect
27 {
28  namespace Particle
29  {
30  namespace Integrator
31  {
39  template <int dim>
40  class RK4 : public Interface<dim>, public SimulatorAccess<dim>
41  {
42  public:
43  RK4();
44 
49  virtual
50  void
51  initialize () override;
52 
71  void
72  local_integrate_step(const typename ParticleHandler<dim>::particle_iterator &begin_particle,
73  const typename ParticleHandler<dim>::particle_iterator &end_particle,
74  const std::vector<Tensor<1,dim>> &old_velocities,
75  const std::vector<Tensor<1,dim>> &velocities,
76  const double dt) override;
77 
88  bool new_integration_step() override;
89 
101  static const unsigned int n_integrator_properties = 4*dim;
102 
103  private:
108  unsigned int integrator_substep;
109 
113  std::array<unsigned int,4> property_index_k;
114  };
115  }
116  }
117 }
118 
119 #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
std::array< unsigned int, 4 > property_index_k
Definition: rk_4.h:113
static const unsigned int n_integrator_properties
Definition: rk_4.h:101
unsigned int integrator_substep
Definition: rk_4.h:108
virtual void initialize() override