ExplicitTimeIntegrator.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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"
#include "MeshChangedInterface.h"

#include "libmesh/linear_solver.h"
#include "libmesh/preconditioner.h"
#include "libmesh/numeric_vector.h"

#include "LumpedPreconditioner.h"

class ExplicitTimeIntegrator : public TimeIntegrator, public MeshChangedInterface
{
public:
  static InputParameters validParams();

  ExplicitTimeIntegrator(const InputParameters & parameters);

  virtual void initialSetup() override;
  virtual void init() override;
  virtual void preSolve() override;
  virtual void meshChanged() override;
  virtual bool isExplicit() const override { return true; }

protected:
  enum SolveType
  {
    CONSISTENT,
    LUMPED,
    LUMP_PRECONDITIONED
  };

  virtual bool performExplicitSolve(SparseMatrix<Number> & mass_matrix);

  bool solveLinearSystem(SparseMatrix<Number> & mass_matrix);

  bool checkLinearConvergence();

  virtual TagID massMatrixTagID() const { return _Ke_time_tag; }

  MooseEnum _solve_type;

  NumericVector<Real> * _explicit_residual;

  NumericVector<Real> * _solution_update;

  NumericVector<Real> * _mass_matrix_diag_inverted;

  NumericVector<Real> * _ones;

  TagID _Ke_time_tag;

  std::unique_ptr<libMesh::LinearSolver<Number>> _linear_solver;

  std::unique_ptr<LumpedPreconditioner> _preconditioner;

  Real _current_time;
};