ThreadedRadialAverageLoop Class Reference

RadialAverage threaded loop. More...

#include <ThreadedRadialAverageLoop.h>

Public Member Functions
	ThreadedRadialAverageLoop (RadialAverage &)
	ThreadedRadialAverageLoop (const ThreadedRadialAverageLoop &x, Threads::split split)
	Splitting constructor.
virtual	~ThreadedRadialAverageLoop ()
	dummy virtual destructor
void	operator() (const QPDataRange &range)
	parens operator with the code that is executed in threads
virtual void	join (const ThreadedRadialAverageLoop &)
	thread join method

Protected Attributes
RadialAverage &	_radavg
	rasterizer to manage the sample data
THREAD_ID	_tid
	ID number of the current thread.

Detailed Description

RadialAverage threaded loop.

Definition at line 22 of file ThreadedRadialAverageLoop.h.

Constructor & Destructor Documentation

ThreadedRadialAverageLoop() [1/2]

ThreadedRadialAverageLoop::ThreadedRadialAverageLoop

(

RadialAverage &

green

)

Definition at line 25 of file ThreadedRadialAverageLoop.C.

: _radavg(green) {}

ThreadedRadialAverageLoop() [2/2]

ThreadedRadialAverageLoop::ThreadedRadialAverageLoop	(	const ThreadedRadialAverageLoop &	x,
		Threads::split	split
	)

Splitting constructor.

Definition at line 28 of file ThreadedRadialAverageLoop.C.

  : _radavg(x._radavg)
{
}

~ThreadedRadialAverageLoop()

virtual ThreadedRadialAverageLoop::~ThreadedRadialAverageLoop ( )

inlinevirtual

dummy virtual destructor

Definition at line 31 of file ThreadedRadialAverageLoop.h.

{}

Member Function Documentation

join()

virtual void ThreadedRadialAverageLoop::join ( const ThreadedRadialAverageLoop &  )

inlinevirtual

thread join method

Definition at line 37 of file ThreadedRadialAverageLoop.h.

{}

operator()()

void ThreadedRadialAverageLoop::operator()

(

const QPDataRange &

range

)

parens operator with the code that is executed in threads

Definition at line 35 of file ThreadedRadialAverageLoop.C.

{
  // fetch data from parent
  const auto radius = _radavg._radius;
  const auto & qp_data = _radavg._qp_data;
  const auto & kd_tree = _radavg._kd_tree;
  const auto & weights_type = _radavg._weights_type;
 
  // tree search data structures
  std::vector<nanoflann::ResultItem<std::size_t, Real>> ret_matches;
  nanoflann::SearchParameters search_params;
 
  // result map entry
  const auto end_it = _radavg._average.end();
  auto it = end_it;
 
  // iterate over qp range
  for (auto && local_qp : qpdata_range)
  {
    // Look up result map iterator only if we enter a new element. this saves a bunch
    // of map lookups because same element entries are consecutive in the qp_data vector.
    if (it == end_it || it->first != local_qp._elem_id)
      it = _radavg._average.find(local_qp._elem_id);
 
    // initialize result entry
    mooseAssert(it != end_it, "Current element id not found in result set.");
    auto & sum = it->second[local_qp._qp];
    sum = 0.0;
 
    ret_matches.clear();
    std::size_t n_result =
        kd_tree->radiusSearch(&(local_qp._q_point(0)), radius * radius, ret_matches, search_params);
    Real total_vol = 0.0;
    Real weight = 1.0;
    for (std::size_t j = 0; j < n_result; ++j)
    {
      const auto & other_qp = qp_data[ret_matches[j].first];
      switch (weights_type)
      {
        case RadialAverage::WeightsType::CONSTANT:
          break;
 
        case RadialAverage::WeightsType::LINEAR:
          weight = radius - std::sqrt(ret_matches[j].second);
          break;
 
        case RadialAverage::WeightsType::COSINE:
          weight = std::cos(std::sqrt(ret_matches[j].second) / radius * libMesh::pi) + 1.0;
          break;
      }
 
      sum += other_qp._value * other_qp._volume * weight;
      total_vol += other_qp._volume * weight;
    }
    sum /= total_vol;
  }
}

Member Data Documentation

_radavg

RadialAverage& ThreadedRadialAverageLoop::_radavg

protected

rasterizer to manage the sample data

Definition at line 41 of file ThreadedRadialAverageLoop.h.

Referenced by operator()().

_tid

THREAD_ID ThreadedRadialAverageLoop::_tid

protected

ID number of the current thread.

Definition at line 44 of file ThreadedRadialAverageLoop.h.

---

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

• include/userobjects/ThreadedRadialAverageLoop.h
• src/userobjects/ThreadedRadialAverageLoop.C