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

#include <SamplerBase.h>

Inheritance diagram for SamplerBase:

Public Member Functions
	SamplerBase (const InputParameters &parameters, VectorPostprocessor *vpp, const libMesh::Parallel::Communicator &comm)

virtual	~SamplerBase ()=default

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 41 of file SamplerBase.h.

Constructor & Destructor Documentation

◆ SamplerBase()

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

Parameters

parameters	The parameters for the object
vpp	A pointer to the child object
comm	The communicator of the child

Definition at line 35 of file SamplerBase.C.

   : _sampler_params(parameters),
     _vpp(vpp),
     _comm(comm),
     _sort_by(parameters.get<MooseEnum>("sort_by")),
     _x(vpp->declareVector("x")),
     _y(vpp->declareVector("y")),
     _z(vpp->declareVector("z")),
     _id(vpp->declareVector("id"))
 {
 }

◆ ~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

p	The point where you took the sample
id	This can either be an actual ID or a distance or anything else you want
values	The value of each variable

Definition at line 60 of file SamplerBase.C.

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

 {
   _x.push_back(p(0));
   _y.push_back(p(1));
   _z.push_back(p(2));
 
   _id.push_back(id);
 
   mooseAssert(values.size() == _variable_names.size(), "Mismatch of variable names to vector size");
   for (auto i = beginIndex(values); i < values.size(); ++i)
     _values[i]->emplace_back(values[i]);
 }

◆ 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, NodalValueSampler, ElementValueSampler, LineFunctionSampler, and SideValueSampler.

Definition at line 90 of file SamplerBase.C.

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

 {
   constexpr auto NUM_ID_VECTORS = 4;
 
   std::vector<VectorPostprocessorValue *> vec_ptrs;
   vec_ptrs.reserve(_values.size() + NUM_ID_VECTORS);
   // Initialize the pointer vector with the position and ID vectors
   vec_ptrs = {{&_x, &_y, &_z, &_id}};
   // Now extend the vector by all the remaining values vector before processing
   vec_ptrs.insert(vec_ptrs.end(), _values.begin(), _values.end());
 
   // Gather up each of the partial vectors
   for (auto vec_ptr : vec_ptrs)
     _comm.allgather(*vec_ptr, /* identical buffer lengths = */ false);
 
   // Now create an index vector by using an indirect sort
   std::vector<std::size_t> sorted_indices;
   Moose::indirectSort(vec_ptrs[_sort_by]->begin(), vec_ptrs[_sort_by]->end(), sorted_indices);
 
   // This vector is used as temp storage to sort each of the remaining vectors according to the
   // first
   auto vector_length = sorted_indices.size();
   VectorPostprocessorValue tmp_vector(vector_length);
 
 #ifndef NDEBUG
   for (const auto vec_ptr : vec_ptrs)
     if (vec_ptr->size() != vector_length)
       mooseError("Vector length mismatch");
 #endif
 
   // Sort each of the vectors using the same sorted indices
   for (auto i = beginIndex(vec_ptrs); i < vec_ptrs.size(); ++i)
   {
     for (auto j = beginIndex(sorted_indices); j < sorted_indices.size(); ++j)
       tmp_vector[j] = (*vec_ptrs[i])[sorted_indices[j]];
 
     // Swap vector storage with sorted vector
     vec_ptrs[i]->swap(tmp_vector);
   }
 }

◆ 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, NodalValueSampler, ElementValueSampler, LineFunctionSampler, and SideValueSampler.

Definition at line 74 of file SamplerBase.C.

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

 {
   // Don't reset the vectors if we want to retain history
   if (_vpp->containsCompleteHistory() && _comm.rank() == 0)
     return;
 
   _x.clear();
   _y.clear();
   _z.clear();
   _id.clear();
 
   std::for_each(
       _values.begin(), _values.end(), [](VectorPostprocessorValue * vec) { vec->clear(); });
 }

◆ 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_names The names of the variables. Note: The order of the variables sets the order of the values for addSample()

Definition at line 50 of file SamplerBase.C.

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

 {
   _variable_names = variable_names;
   _values.reserve(variable_names.size());
 
   for (const auto & variable_name : variable_names)
     _values.push_back(&_vpp->declareVector(variable_name));
 }

◆ threadJoin()

void SamplerBase::threadJoin ( const SamplerBase & y )

protectedvirtual

Join the values.

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

Parameters

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

Definition at line 144 of file SamplerBase.C.

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

 {
   _x.insert(_x.end(), y._x.begin(), y._x.end());
   _y.insert(_y.end(), y._y.begin(), y._y.end());
   _z.insert(_z.end(), y._z.begin(), y._z.end());
 
   _id.insert(_id.end(), y._id.begin(), y._id.end());
 
   for (unsigned int i = 0; i < _variable_names.size(); i++)
     _values[i]->insert(_values[i]->end(), y._values[i]->begin(), y._values[i]->end());
 }