libMesh::QoISet Class Reference

Data structure for specifying which Quantities of Interest should be calculated in an adjoint or a parameter sensitivity calculation. More...

#include <qoi_set.h>

Classes
class	iterator

Public Member Functions
	QoISet ()
	Empty constructor: "calculate all QoIs in the System". More...
	QoISet (const System &sys)
	Default constructor: "calculate all QoIs in the System", "give every QoI weight 1.0". More...
	QoISet (std::vector< bool > indices)
	Constructor-from-vector-of-bool: "calculate the QoIs for which \p indices[q] is true". More...
	QoISet (const std::vector< unsigned int > &indices)
	Constructor-from-vector: "calculate the listed QoIs", "give every QoI weight 1.0". More...
void	clear ()
	Resets to "calculate all QoIs, give every QoI weight 1.0". More...
std::size_t	size (const System &sys) const
void	add_indices (const std::vector< unsigned int > &indices)
	Add this indices to the set to be calculated. More...
void	add_index (std::size_t)
	Add this index to the set to be calculated. More...
void	remove_indices (const std::vector< unsigned int > &indices)
	Remove these indices from the set to be calculated. More...
void	remove_index (std::size_t)
	Remove this index from the set to be calculated. More...
void	set_weight (std::size_t, Real)
	Set the weight for this index. More...
Real	weight (std::size_t) const
	Get the weight for this index (default 1.0) More...
bool	has_index (std::size_t) const
	Return whether or not this index is in the set to be calculated. More...
iterator	begin () const
	Return an iterator pointing to the first index in the set. More...

Private Attributes
std::vector< bool >	_indices
	Interpret _indices.empty() to mean "calculate all indices". More...
std::vector< Real >	_weights
	Interpret _weights.size() <= i to mean "weight i = 1.0". More...

Detailed Description

Data structure for specifying which Quantities of Interest should be calculated in an adjoint or a parameter sensitivity calculation.

Author

Roy Stogner

Date

2009 Used to specify quantities of interest in a simulation.

Definition at line 45 of file qoi_set.h.

Constructor & Destructor Documentation

QoISet() [1/4]

libMesh::QoISet::QoISet ( )

inline

Empty constructor: "calculate all QoIs in the System".

No further changes to this special QoISet should be made; it doesn't even know how many QoIs your system has, it just knows to instruct a function to use all of them.

Definition at line 100 of file qoi_set.h.

: _indices(), _weights() {}

QoISet() [2/4]

libMesh::QoISet::QoISet ( const System &  sys )

explicit

Default constructor: "calculate all QoIs in the System", "give every QoI weight 1.0".

Definition at line 31 of file qoi_set.C.

: _indices(sys.n_qois(), true) {}

QoISet() [3/4]

libMesh::QoISet::QoISet ( std::vector< bool >  indices )

inlineexplicit

Constructor-from-vector-of-bool: "calculate the QoIs for which \p indices[q] is true".

Definition at line 114 of file qoi_set.h.

                                  :
    _indices(std::move(indices)), _weights() {}

QoISet() [4/4]

libMesh::QoISet::QoISet ( const std::vector< unsigned int > &  indices )

inlineexplicit

Constructor-from-vector: "calculate the listed QoIs", "give every QoI weight 1.0".

Definition at line 195 of file qoi_set.h.

References add_indices().

                                                      :
  _indices(), _weights()
{
  this->add_indices(indices);
}

Member Function Documentation

add_index()

void libMesh::QoISet::add_index ( std::size_t  i )

inline

Add this index to the set to be calculated.

Definition at line 204 of file qoi_set.h.

References _indices.

{
  if (i >= _indices.size())
    _indices.resize(i+1, true);
  _indices[i] = true;
}

add_indices()

void libMesh::QoISet::add_indices

(

const std::vector< unsigned int > &

indices

)

Add this indices to the set to be calculated.

Definition at line 46 of file qoi_set.C.

References _indices.

Referenced by main(), and QoISet().

{
  unsigned int max_size = 0;
  for (const auto & i : indices)
    max_size = std::max(max_size, i + 1);

  _indices.resize(max_size);

  for (const auto & i : indices)
    _indices[i] = true;
}

begin()

iterator libMesh::QoISet::begin ( ) const

inline

Return an iterator pointing to the first index in the set.

Definition at line 173 of file qoi_set.h.

References _indices.

{ return iterator(0, _indices); }

clear()

void libMesh::QoISet::clear ( )

inline

Resets to "calculate all QoIs, give every QoI weight 1.0".

Definition at line 127 of file qoi_set.h.

References _indices, and _weights.

{ _indices.clear(); _weights.clear(); }

has_index()

bool libMesh::QoISet::has_index ( std::size_t  i ) const

inline

Return whether or not this index is in the set to be calculated.

Definition at line 224 of file qoi_set.h.

References _indices.

Referenced by libMesh::ImplicitSystem::adjoint_qoi_parameter_sensitivity(), libMesh::ImplicitSystem::adjoint_solve(), libMesh::SensitivityData::allocate_data(), libMesh::SensitivityData::allocate_hessian_data(), libMesh::ExplicitSystem::assemble_qoi(), libMesh::FEMSystem::assemble_qoi(), libMesh::ExplicitSystem::assemble_qoi_derivative(), libMesh::FEMSystem::assemble_qoi_derivative(), libMesh::AdjointRefinementEstimator::estimate_error(), libMesh::AdjointResidualErrorEstimator::estimate_error(), libMesh::ImplicitSystem::forward_qoi_parameter_sensitivity(), libMesh::Euler2Solver::integrate_adjoint_refinement_error_estimate(), libMesh::TwostepTimeSolver::integrate_adjoint_refinement_error_estimate(), libMesh::EulerSolver::integrate_adjoint_refinement_error_estimate(), libMesh::SteadySolver::integrate_adjoint_refinement_error_estimate(), libMesh::ImplicitSystem::qoi_parameter_hessian(), libMesh::ImplicitSystem::qoi_parameter_hessian_vector_product(), size(), and libMesh::ImplicitSystem::weighted_sensitivity_adjoint_solve().

{
  return (_indices.size() <= i || _indices[i]);
}

remove_index()

void libMesh::QoISet::remove_index ( std::size_t  i )

inline

Remove this index from the set to be calculated.

Definition at line 214 of file qoi_set.h.

References _indices.

{
  if (i >= _indices.size())
    _indices.resize(i+1, true);
  _indices[i] = false;
}

remove_indices()

void libMesh::QoISet::remove_indices ( const std::vector< unsigned int > &  indices )

inline

Remove these indices from the set to be calculated.

Definition at line 61 of file qoi_set.C.

References _indices.

{
  for (const auto & i : indices)
    _indices[i] = false;
}

set_weight()

void libMesh::QoISet::set_weight ( std::size_t  i, Real  w  )

inline

Set the weight for this index.

Definition at line 232 of file qoi_set.h.

References _weights.

Referenced by main().

{
  if (_weights.size() <= i)
    _weights.resize(i+1, 1.0);

  _weights[i] = w;
}

size()

std::size_t libMesh::QoISet::size

(

const System &

sys

)

const

Returns

The number of QoIs that would be computed for the System sys

Definition at line 35 of file qoi_set.C.

References has_index(), libMesh::make_range(), and libMesh::System::n_qois().

Referenced by libMesh::TwostepTimeSolver::integrate_adjoint_sensitivity(), libMesh::UnsteadySolver::integrate_adjoint_sensitivity(), and libMesh::System::qoi_parameter_sensitivity().

{
  std::size_t qoi_count = 0;
  for (auto i : make_range(sys.n_qois()))
    if (this->has_index(i))
      qoi_count++;
  return qoi_count;
}

weight()

Real libMesh::QoISet::weight ( std::size_t  i ) const

inline

Get the weight for this index (default 1.0)

Definition at line 243 of file qoi_set.h.

References _weights.

Referenced by libMesh::AdjointRefinementEstimator::estimate_error(), and libMesh::AdjointResidualErrorEstimator::estimate_error().

{
  if (_weights.size() <= i)
    return 1.0;
  return _weights[i];
}

Member Data Documentation

_indices

std::vector<bool> libMesh::QoISet::_indices

private

Interpret _indices.empty() to mean "calculate all indices".

Definition at line 179 of file qoi_set.h.

Referenced by add_index(), add_indices(), begin(), clear(), has_index(), remove_index(), and remove_indices().

_weights

std::vector<Real> libMesh::QoISet::_weights

private

Interpret _weights.size() <= i to mean "weight i = 1.0".

Definition at line 184 of file qoi_set.h.

Referenced by clear(), set_weight(), and weight().

---

The documentation for this class was generated from the following files:

• include/systems/qoi_set.h
• src/systems/qoi_set.C