LStableDirk3.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "TimeIntegrator.h"

class LStableDirk3;

template <>
InputParameters validParams<LStableDirk3>();

class LStableDirk3 : public TimeIntegrator
{
public:
  LStableDirk3(const InputParameters & parameters);

  virtual int order() override { return 3; }
  virtual void computeTimeDerivatives() override;
  virtual void computeADTimeDerivatives(DualReal & ad_u_dot,
                                        const dof_id_type & dof) const override;
  virtual void solve() override;
  virtual void postResidual(NumericVector<Number> & residual) override;

protected:
  template <typename T, typename T2>
  void computeTimeDerivativeHelper(T & u_dot, const T2 & u_old) const;

  // Indicates the current stage.
  unsigned int _stage;

  // Store pointers to the various stage residuals
  NumericVector<Number> * _stage_residuals[3];

  // The parameter of the method, set at construction time and cannot be changed.
  const Real _gamma; // 0.4358665215084589

  // Butcher tableau "C" parameters derived from _gamma
  // 0.4358665215084589, 0.7179332607542295, 1.0000000000000000
  Real _c[3];

  // Butcher tableau "A" values derived from _gamma.  We only use the
  // lower triangle of this.
  // 0.4358665215084589
  // 0.2820667392457705,  0.4358665215084589
  // 1.2084966491760099, -0.6443631706844688, 0.4358665215084589
  Real _a[3][3];
};

template <typename T, typename T2>
void
LStableDirk3::computeTimeDerivativeHelper(T & u_dot, const T2 & u_old) const
{
  u_dot -= u_old;
  u_dot *= 1. / _dt;
}