20 #ifndef LIBMESH_ADJOINT_RESIDUAL_ERROR_ESTIMATOR_H 21 #define LIBMESH_ADJOINT_RESIDUAL_ERROR_ESTIMATOR_H 24 #include "libmesh/error_estimator.h" 25 #include "libmesh/qoi_set.h" 119 bool estimate_parent_error =
false)
override;
144 #endif // LIBMESH_ADJOINT_RESIDUAL_ERROR_ESTIMATOR_H std::unique_ptr< ErrorEstimator > & dual_error_estimator()
Access to the "subestimator" (default: PatchRecovery) to use on the dual/adjoint solution.
This class implements a goal oriented error indicator, by weighting residual-based estimates from the...
Data structure for specifying which Quantities of Interest should be calculated in an adjoint or a pa...
The ErrorVector is a specialization of the StatisticsVector for error data computed on a finite eleme...
ErrorEstimatorType
Defines an enum for the different types of error estimators which are available.
QoISet & qoi_set()
Access to the QoISet (default: weight all QoIs equally) to use when computing errors.
QoISet _qoi_set
A QoISet to handle cases with multiple QoIs available.
The libMesh namespace provides an interface to certain functionality in the library.
AdjointResidualErrorEstimator()
Constructor.
std::string error_plot_suffix
To aid in investigating error estimator behavior, set this string to a suffix with which to plot (pre...
Manages consistently variables, degrees of freedom, and coefficient vectors.
This class holds functions that will estimate the error in a finite element solution on a given mesh...
virtual ~AdjointResidualErrorEstimator()=default
const QoISet & qoi_set() const
Access to the QoISet (default: weight all QoIs equally) to use when computing errors.
std::unique_ptr< ErrorEstimator > _dual_error_estimator
An error estimator for the adjoint problem.
std::unique_ptr< ErrorEstimator > & primal_error_estimator()
Access to the "subestimator" (default: PatchRecovery) to use on the primal/forward solution...
std::unique_ptr< ErrorEstimator > _primal_error_estimator
An error estimator for the forward problem.
AdjointResidualErrorEstimator & operator=(const AdjointResidualErrorEstimator &)=delete
virtual void estimate_error(const System &system, ErrorVector &error_per_cell, const NumericVector< Number > *solution_vector=nullptr, bool estimate_parent_error=false) override
Compute the adjoint-weighted error on each element and place it in the error_per_cell vector...
virtual ErrorEstimatorType type() const override