www.mooseframework.org
Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
SamplerBase Class Reference

Base class for VectorPostprocessors that need to do "sampling" of values in the domain. More...

#include <SamplerBase.h>

Inheritance diagram for SamplerBase:
[legend]

Public Member Functions

 SamplerBase (const InputParameters &parameters, VectorPostprocessor *vpp, const libMesh::Parallel::Communicator &comm)
 
virtual ~SamplerBase ()=default
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

void setupVariables (const std::vector< std::string > &variable_names)
 You MUST call this in the constructor of the child class and pass down the name of the variables. More...
 
virtual void addSample (const Point &p, const Real &id, const std::vector< Real > &values)
 Call this with the value of every variable at each point you want to sample at. More...
 
virtual void initialize ()
 Initialize the datastructures. More...
 
virtual void finalize ()
 Finalize the values. More...
 
virtual void threadJoin (const SamplerBase &y)
 Join the values. More...
 

Protected Attributes

const InputParameters_sampler_params
 The child params. More...
 
VectorPostprocessor_vpp
 The child VectorPostprocessor. More...
 
const libMesh::Parallel::Communicator_comm
 The communicator of the child. More...
 
std::vector< std::string > _variable_names
 The variable names. More...
 
const unsigned int _sort_by
 What to sort by. More...
 
VectorPostprocessorValue_x
 x coordinate of the points More...
 
VectorPostprocessorValue_y
 y coordinate of the points More...
 
VectorPostprocessorValue_z
 x coordinate of the points More...
 
VectorPostprocessorValue_id
 The node ID of each point. More...
 
std::vector< VectorPostprocessorValue * > _values
 

Detailed Description

Base class for VectorPostprocessors that need to do "sampling" of values in the domain.

Definition at line 36 of file SamplerBase.h.

Constructor & Destructor Documentation

◆ SamplerBase()

SamplerBase::SamplerBase ( const InputParameters parameters,
VectorPostprocessor vpp,
const libMesh::Parallel::Communicator comm 
)
Parameters
parametersThe parameters for the object
vppA pointer to the child object
commThe communicator of the child

Definition at line 36 of file SamplerBase.C.

39  : _sampler_params(parameters),
40  _vpp(vpp),
41  _comm(comm),
42  _sort_by(parameters.get<MooseEnum>("sort_by")),
43  _x(vpp->declareVector("x")),
44  _y(vpp->declareVector("y")),
45  _z(vpp->declareVector("z")),
46  _id(vpp->declareVector("id"))
47 {
48 }
std::vector< std::pair< R1, R2 > > get(const std::string &param1, const std::string &param2) const
Combine two vector parameters into a single vector of pairs.
VectorPostprocessorValue & _y
y coordinate of the points
Definition: SamplerBase.h:111
const unsigned int _sort_by
What to sort by.
Definition: SamplerBase.h:106
VectorPostprocessor * _vpp
The child VectorPostprocessor.
Definition: SamplerBase.h:97
This is a "smart" enum class intended to replace many of the shortcomings in the C++ enum type It sho...
Definition: MooseEnum.h:31
VectorPostprocessorValue & declareVector(const std::string &vector_name)
Register a new vector to fill up.
VectorPostprocessorValue & _x
x coordinate of the points
Definition: SamplerBase.h:109
const InputParameters & _sampler_params
The child params.
Definition: SamplerBase.h:94
VectorPostprocessorValue & _id
The node ID of each point.
Definition: SamplerBase.h:116
VectorPostprocessorValue & _z
x coordinate of the points
Definition: SamplerBase.h:113
const libMesh::Parallel::Communicator & _comm
The communicator of the child.
Definition: SamplerBase.h:100

◆ ~SamplerBase()

virtual SamplerBase::~SamplerBase ( )
virtualdefault

Member Function Documentation

◆ addSample()

void SamplerBase::addSample ( const Point &  p,
const Real id,
const std::vector< Real > &  values 
)
protectedvirtual

Call this with the value of every variable at each point you want to sample at.

Parameters
pThe point where you took the sample
idThis can either be an actual ID or a distance or anything else you want
valuesThe value of each variable

Definition at line 61 of file SamplerBase.C.

Referenced by SideValueSampler::execute(), LineFunctionSampler::execute(), ElementValueSampler::execute(), NodalValueSampler::execute(), and PointSamplerBase::finalize().

62 {
63  _x.push_back(p(0));
64  _y.push_back(p(1));
65  _z.push_back(p(2));
66 
67  _id.push_back(id);
68 
69  mooseAssert(values.size() == _variable_names.size(), "Mismatch of variable names to vector size");
70  for (MooseIndex(values) i = 0; i < values.size(); ++i)
71  _values[i]->emplace_back(values[i]);
72 }
std::vector< std::string > _variable_names
The variable names.
Definition: SamplerBase.h:103
VectorPostprocessorValue & _y
y coordinate of the points
Definition: SamplerBase.h:111
std::vector< VectorPostprocessorValue * > _values
Definition: SamplerBase.h:118
VectorPostprocessorValue & _x
x coordinate of the points
Definition: SamplerBase.h:109
VectorPostprocessorValue & _id
The node ID of each point.
Definition: SamplerBase.h:116
VectorPostprocessorValue & _z
x coordinate of the points
Definition: SamplerBase.h:113

◆ finalize()

void SamplerBase::finalize ( )
protectedvirtual

Finalize the values.

YOU MUST CALL THIS DURING finalize() in the child class!

We have several vectors that all need to be processed in the same way. Rather than enumerate them all, let's just create a vector of pointers and work on them that way.

We now have one sorted vector. The remaining vectors need to be sorted according to that vector. We'll need a temporary vector to hold values as we map them according to the sorted indices. After that, we'll swap the vector contents with the sorted vector to get the values back into the original vector.

Reimplemented in LineMaterialSamplerBase< T >, LineMaterialSamplerBase< Real >, PointSamplerBase, ElementValueSampler, NodalValueSampler, LineFunctionSampler, and SideValueSampler.

Definition at line 91 of file SamplerBase.C.

Referenced by SideValueSampler::finalize(), LineFunctionSampler::finalize(), NodalValueSampler::finalize(), ElementValueSampler::finalize(), PointSamplerBase::finalize(), and LineMaterialSamplerBase< Real >::finalize().

92 {
98  constexpr auto NUM_ID_VECTORS = 4;
99 
100  std::vector<VectorPostprocessorValue *> vec_ptrs;
101  vec_ptrs.reserve(_values.size() + NUM_ID_VECTORS);
102  // Initialize the pointer vector with the position and ID vectors
103  vec_ptrs = {{&_x, &_y, &_z, &_id}};
104  // Now extend the vector by all the remaining values vector before processing
105  vec_ptrs.insert(vec_ptrs.end(), _values.begin(), _values.end());
106 
107  // Gather up each of the partial vectors
108  for (auto vec_ptr : vec_ptrs)
109  _comm.allgather(*vec_ptr, /* identical buffer lengths = */ false);
110 
111  // Now create an index vector by using an indirect sort
112  std::vector<std::size_t> sorted_indices;
113  Moose::indirectSort(vec_ptrs[_sort_by]->begin(), vec_ptrs[_sort_by]->end(), sorted_indices);
114 
122  // This vector is used as temp storage to sort each of the remaining vectors according to the
123  // first
124  auto vector_length = sorted_indices.size();
125  VectorPostprocessorValue tmp_vector(vector_length);
126 
127 #ifndef NDEBUG
128  for (const auto vec_ptr : vec_ptrs)
129  if (vec_ptr->size() != vector_length)
130  mooseError("Vector length mismatch");
131 #endif
132 
133  // Sort each of the vectors using the same sorted indices
134  for (auto & vec_ptr : vec_ptrs)
135  {
136  for (MooseIndex(sorted_indices) j = 0; j < sorted_indices.size(); ++j)
137  tmp_vector[j] = (*vec_ptr)[sorted_indices[j]];
138 
139  // Swap vector storage with sorted vector
140  vec_ptr->swap(tmp_vector);
141  }
142 }
void allgather(const T &send_data, std::vector< T, A > &recv_data) const
void mooseError(Args &&... args)
Emit an error message with the given stringified, concatenated args and terminate the application...
Definition: MooseError.h:284
void indirectSort(RandomAccessIterator beg, RandomAccessIterator end, std::vector< size_t > &b)
Definition: IndirectSort.h:68
VectorPostprocessorValue & _y
y coordinate of the points
Definition: SamplerBase.h:111
const unsigned int _sort_by
What to sort by.
Definition: SamplerBase.h:106
std::vector< VectorPostprocessorValue * > _values
Definition: SamplerBase.h:118
VectorPostprocessorValue & _x
x coordinate of the points
Definition: SamplerBase.h:109
std::vector< Real > VectorPostprocessorValue
Definition: MooseTypes.h:192
VectorPostprocessorValue & _id
The node ID of each point.
Definition: SamplerBase.h:116
VectorPostprocessorValue & _z
x coordinate of the points
Definition: SamplerBase.h:113
const libMesh::Parallel::Communicator & _comm
The communicator of the child.
Definition: SamplerBase.h:100

◆ initialize()

void SamplerBase::initialize ( )
protectedvirtual

Initialize the datastructures.

YOU MUST CALL THIS DURING initialize() in the child class!

Reimplemented in LineMaterialSamplerBase< T >, LineMaterialSamplerBase< Real >, PointSamplerBase, PointVariableSamplerBase, PositionsFunctorValueSampler, ElementValueSampler, NodalValueSampler, LineFunctionSampler, PointValueSampler, and SideValueSampler.

Definition at line 75 of file SamplerBase.C.

Referenced by SideValueSampler::initialize(), LineFunctionSampler::initialize(), NodalValueSampler::initialize(), ElementValueSampler::initialize(), PointVariableSamplerBase::initialize(), PointSamplerBase::initialize(), and LineMaterialSamplerBase< Real >::initialize().

76 {
77  // Don't reset the vectors if we want to retain history
78  if (_vpp->containsCompleteHistory() && _comm.rank() == 0)
79  return;
80 
81  _x.clear();
82  _y.clear();
83  _z.clear();
84  _id.clear();
85 
86  std::for_each(
87  _values.begin(), _values.end(), [](VectorPostprocessorValue * vec) { vec->clear(); });
88 }
processor_id_type rank() const
VectorPostprocessorValue & _y
y coordinate of the points
Definition: SamplerBase.h:111
std::vector< VectorPostprocessorValue * > _values
Definition: SamplerBase.h:118
bool containsCompleteHistory() const
Return whether or not this VectorPostprocessor contains complete history.
VectorPostprocessor * _vpp
The child VectorPostprocessor.
Definition: SamplerBase.h:97
VectorPostprocessorValue & _x
x coordinate of the points
Definition: SamplerBase.h:109
std::vector< Real > VectorPostprocessorValue
Definition: MooseTypes.h:192
VectorPostprocessorValue & _id
The node ID of each point.
Definition: SamplerBase.h:116
VectorPostprocessorValue & _z
x coordinate of the points
Definition: SamplerBase.h:113
const libMesh::Parallel::Communicator & _comm
The communicator of the child.
Definition: SamplerBase.h:100

◆ setupVariables()

void SamplerBase::setupVariables ( const std::vector< std::string > &  variable_names)
protected

You MUST call this in the constructor of the child class and pass down the name of the variables.

Parameters
variable_namesThe names of the variables. Note: The order of the variables sets the order of the values for addSample()

Definition at line 51 of file SamplerBase.C.

Referenced by ElementValueSampler::ElementValueSampler(), LineFunctionSampler::LineFunctionSampler(), LineMaterialSamplerBase< Real >::LineMaterialSamplerBase(), NodalValueSampler::NodalValueSampler(), PointVariableSamplerBase::PointVariableSamplerBase(), PositionsFunctorValueSampler::PositionsFunctorValueSampler(), and SideValueSampler::SideValueSampler().

52 {
53  _variable_names = variable_names;
54  _values.reserve(variable_names.size());
55 
56  for (const auto & variable_name : variable_names)
57  _values.push_back(&_vpp->declareVector(variable_name));
58 }
std::vector< std::string > _variable_names
The variable names.
Definition: SamplerBase.h:103
std::vector< VectorPostprocessorValue * > _values
Definition: SamplerBase.h:118
VectorPostprocessor * _vpp
The child VectorPostprocessor.
Definition: SamplerBase.h:97
VectorPostprocessorValue & declareVector(const std::string &vector_name)
Register a new vector to fill up.

◆ threadJoin()

void SamplerBase::threadJoin ( const SamplerBase y)
protectedvirtual

Join the values.

YOU MUST CALL THIS DURING threadJoin() in the child class!

Parameters
yYou must cast the UserObject to your child class type first then you can pass it in here.

Definition at line 145 of file SamplerBase.C.

Referenced by SideValueSampler::threadJoin(), ElementValueSampler::threadJoin(), and NodalValueSampler::threadJoin().

146 {
147  _x.insert(_x.end(), y._x.begin(), y._x.end());
148  _y.insert(_y.end(), y._y.begin(), y._y.end());
149  _z.insert(_z.end(), y._z.begin(), y._z.end());
150 
151  _id.insert(_id.end(), y._id.begin(), y._id.end());
152 
153  for (MooseIndex(_variable_names) i = 0; i < _variable_names.size(); i++)
154  _values[i]->insert(_values[i]->end(), y._values[i]->begin(), y._values[i]->end());
155 }
std::vector< std::string > _variable_names
The variable names.
Definition: SamplerBase.h:103
VectorPostprocessorValue & _y
y coordinate of the points
Definition: SamplerBase.h:111
std::vector< VectorPostprocessorValue * > _values
Definition: SamplerBase.h:118
VectorPostprocessorValue & _x
x coordinate of the points
Definition: SamplerBase.h:109
VectorPostprocessorValue & _id
The node ID of each point.
Definition: SamplerBase.h:116
VectorPostprocessorValue & _z
x coordinate of the points
Definition: SamplerBase.h:113

◆ validParams()

InputParameters SamplerBase::validParams ( )
static

Definition at line 22 of file SamplerBase.C.

Referenced by SideValueSampler::validParams(), LineFunctionSampler::validParams(), NodalValueSampler::validParams(), ElementValueSampler::validParams(), PointVariableSamplerBase::validParams(), PointSamplerBase::validParams(), and LineMaterialSamplerBase< Real >::validParams().

23 {
25 
26  MooseEnum sort_options("x y z id");
27  params.addRequiredParam<MooseEnum>("sort_by", sort_options, "What to sort the samples by");
28 
29  // The value from this VPP is naturally already on every processor
30  // TODO: Make this not the case! See #11415
31  params.set<bool>("_auto_broadcast") = false;
32 
33  return params;
34 }
T & set(const std::string &name, bool quiet_mode=false)
Returns a writable reference to the named parameters.
The main MOOSE class responsible for handling user-defined parameters in almost every MOOSE system...
void addRequiredParam(const std::string &name, const std::string &doc_string)
This method adds a parameter and documentation string to the InputParameters object that will be extr...
InputParameters emptyInputParameters()
This is a "smart" enum class intended to replace many of the shortcomings in the C++ enum type It sho...
Definition: MooseEnum.h:31

Member Data Documentation

◆ _comm

const libMesh::Parallel::Communicator& SamplerBase::_comm
protected

The communicator of the child.

Definition at line 100 of file SamplerBase.h.

Referenced by PointSamplerBase::finalize(), finalize(), and initialize().

◆ _id

VectorPostprocessorValue& SamplerBase::_id
protected

The node ID of each point.

Definition at line 116 of file SamplerBase.h.

Referenced by addSample(), finalize(), LineValueSampler::getValue(), initialize(), and threadJoin().

◆ _sampler_params

const InputParameters& SamplerBase::_sampler_params
protected

The child params.

Definition at line 94 of file SamplerBase.h.

◆ _sort_by

const unsigned int SamplerBase::_sort_by
protected

What to sort by.

Definition at line 106 of file SamplerBase.h.

Referenced by finalize(), and LineValueSampler::getValue().

◆ _values

std::vector<VectorPostprocessorValue *> SamplerBase::_values
protected

◆ _variable_names

std::vector<std::string> SamplerBase::_variable_names
protected

The variable names.

Definition at line 103 of file SamplerBase.h.

Referenced by addSample(), setupVariables(), and threadJoin().

◆ _vpp

VectorPostprocessor* SamplerBase::_vpp
protected

The child VectorPostprocessor.

Definition at line 97 of file SamplerBase.h.

Referenced by initialize(), and setupVariables().

◆ _x

VectorPostprocessorValue& SamplerBase::_x
protected

x coordinate of the points

Definition at line 109 of file SamplerBase.h.

Referenced by addSample(), finalize(), initialize(), and threadJoin().

◆ _y

VectorPostprocessorValue& SamplerBase::_y
protected

y coordinate of the points

Definition at line 111 of file SamplerBase.h.

Referenced by addSample(), finalize(), initialize(), and threadJoin().

◆ _z

VectorPostprocessorValue& SamplerBase::_z
protected

x coordinate of the points

Definition at line 113 of file SamplerBase.h.

Referenced by addSample(), finalize(), initialize(), and threadJoin().


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