libMesh
Loading...
Searching...
No Matches
first_order_unsteady_solver.h
Go to the documentation of this file.
1// The libMesh Finite Element Library.
2// Copyright (C) 2002-2026 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3
4// This library is free software; you can redistribute it and/or
5// modify it under the terms of the GNU Lesser General Public
6// License as published by the Free Software Foundation; either
7// version 2.1 of the License, or (at your option) any later version.
8
9// This library is distributed in the hope that it will be useful,
10// but WITHOUT ANY WARRANTY; without even the implied warranty of
11// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12// Lesser General Public License for more details.
13
14// You should have received a copy of the GNU Lesser General Public
15// License along with this library; if not, write to the Free Software
16// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17
18#ifndef LIBMESH_FIRST_ORDER_UNSTEADY_SOLVER_H
19#define LIBMESH_FIRST_ORDER_UNSTEADY_SOLVER_H
20
21#include "libmesh/unsteady_solver.h"
22
23namespace libMesh
24{
75{
76public:
81 explicit
84
88 virtual ~FirstOrderUnsteadySolver () = default;
89
90 virtual unsigned int time_order() const override
91 { return 1; }
92
93 virtual void integrate_qoi_timestep() override = 0;
94
95#ifdef LIBMESH_ENABLE_AMR
96 virtual void integrate_adjoint_refinement_error_estimate(AdjointRefinementEstimator & adjoint_refinement_error_estimator, ErrorVector & QoI_elementwise_error) override = 0;
97#endif // LIBMESH_ENABLE_AMR
98
99protected:
100
106 void prepare_accel(DiffContext & context);
107
116
117};
118
119} // end namespace libMesh
120
121#endif // LIBMESH_FIRST_ORDER_UNSTEADY_SOLVER_H
This class implements a "brute force" goal-oriented error estimator which computes an estimate of err...
This class provides all data required for a physics package (e.g.
This class provides a specific system class.
Definition diff_system.h:57
The ErrorVector is a specialization of the StatisticsVector for error data computed on a finite eleme...
Generic class from which first order UnsteadySolvers should subclass.
bool compute_second_order_eqns(bool compute_jacobian, DiffContext &c)
If there are second order variables, then we need to compute their residual equations and correspondi...
virtual void integrate_adjoint_refinement_error_estimate(AdjointRefinementEstimator &adjoint_refinement_error_estimator, ErrorVector &QoI_elementwise_error) override=0
A method to compute the adjoint refinement error estimate at the current timestep.
void prepare_accel(DiffContext &context)
If there are second order variables in the system, then we also prepare the accel for those variables...
virtual void integrate_qoi_timestep() override=0
A method to integrate the system::QoI functionals.
virtual ~FirstOrderUnsteadySolver()=default
Destructor.
virtual unsigned int time_order() const override
This is a generic class that defines a solver to handle time integration of DifferentiableSystems.
The libMesh namespace provides an interface to certain functionality in the library.
void compute_jacobian(const NumericVector< Number > &, SparseMatrix< Number > &J, NonlinearImplicitSystem &system)
Definition assembly.C:315