SSPRK104.hxx

Go to the documentation of this file.

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *\
 * Copyright (c) 2025, Davide Stocco and Enrico Bertolazzi.                                      *
 *                                                                                               *
 * The Sandals project is distributed under the BSD 2-Clause License.                            *
 *                                                                                               *
 * Davide Stocco                                                               Enrico Bertolazzi *
 * University of Trento                                                     University of Trento *
 * e-mail: davide.stocco@unitn.it                             e-mail: enrico.bertolazzi@unitn.it *
\* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
 
#pragma once
 
#ifndef SANDALS_SSPRK104_HXX
#define SANDALS_SSPRK104_HXX
 
namespace Sandals
{
  class SSPRK104Tableau : public Tableau<10>
  {
  public:
    using Tableau<10>::Type;
    using Tableau<10>::Vector;
    using Tableau<10>::Matrix;

    SSPRK104Tableau() {
      this->name  = "SSPRK104";
      this->type  = Type::ERK;
      this->order = 4;
      this->A << 0.0,      0.0,      0.0,      0.0,      0.0,      0.0,     0.0,     0.0,     0.0,     0.0,
                 1.0/6.0,  0.0,      0.0,      0.0,      0.0,      0.0,     0.0,     0.0,     0.0,     0.0,
                 1.0/6.0,  1.0/6.0,  0.0,      0.0,      0.0,      0.0,     0.0,     0.0,     0.0,     0.0,
                 1.0/6.0,  1.0/6.0,  1.0/6.0,  0.0,      0.0,      0.0,     0.0,     0.0,     0.0,     0.0,
                 1.0/6.0,  1.0/6.0,  1.0/6.0,  1.0/6.0,  0.0,      0.0,     0.0,     0.0,     0.0,     0.0,
                 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 0.0,     0.0,     0.0,     0.0,     0.0,
                 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/6.0, 0.0,     0.0,     0.0,     0.0,
                 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/6.0, 1.0/6.0, 0.0,     0.0,     0.0,
                 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/6.0, 1.0/6.0, 1.0/6.0, 0.0,     0.0,
                 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/15.0, 1.0/6.0, 1.0/6.0, 1.0/6.0, 1.0/6.0, 0.0;
      this->b << 1.0/10.0, 1.0/10.0, 1.0/10.0, 1.0/10.0, 1.0/10.0, 1.0/10.0, 1.0/10.0, 1.0/10.0, 1.0/10.0, 1.0/10.0;
      this->c << 0.0,      1.0/6.0,  1.0/3.0,  1.0/2.0,  2.0/3.0,  1.0/3.0, 1.0/2.0, 2.0/3.0, 5.0/6.0, 1.0;
    }
  }; // class SSPRK104Tableau

  template <Integer N, Integer M = 0>
  class SSPRK104 : public RungeKutta<10, N, M>
  {
  public:
    using System = typename Implicit<N, M>::Pointer; 

    SSPRK104() : RungeKutta<10, N, M>(SSPRK104Tableau()) {}

    SSPRK104(System t_system) : RungeKutta<10, N, M>(SSPRK104Tableau(), t_system) {}
 
  }; // class SSPRK104
 
} // namespace Sandals
 
#endif // SANDALS_SSPRK104_HXX