https://mooseframework.inl.gov
Classes | Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Static Protected Member Functions | Protected Attributes | Static Protected Attributes | Private Member Functions | Private Attributes | List of all members
RayTracingStudy Class Referenceabstract

Base class for Ray tracing studies that will generate Rays and then propagate all of them to termination. More...

#include <RayTracingStudy.h>

Inheritance diagram for RayTracingStudy:
[legend]

Classes

class  AcquireMoveDuringTraceKey
 Key that is used for restricting access to moveRayToBufferDuringTrace() and acquireRayDuringTrace(). More...
 
class  AcquireRayInternalKey
 Key that is used for restricting access to acquireRayInternal(). More...
 

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 RayTracingStudy (const InputParameters &parameters)
 
virtual void initialSetup () override
 
virtual void residualSetup () override
 
virtual void jacobianSetup () override
 
virtual void meshChanged () override
 
virtual void timestepSetup () override
 
virtual void initialize () override
 
virtual void finalize () override
 
virtual void execute () override
 Executes the study (generates and propagates Rays) More...
 
virtual void segmentSubdomainSetup (const SubdomainID subdomain, const THREAD_ID tid, const RayID ray_id)
 Setup for on subdomain change or subdomain AND ray change during ray tracing. More...
 
virtual void reinitSegment (const Elem *elem, const Point &start, const Point &end, const Real length, THREAD_ID tid)
 Reinitialize objects for a Ray segment for ray tracing. More...
 
virtual void postOnSegment (const THREAD_ID tid, const std::shared_ptr< Ray > &ray)
 Called at the end of a Ray segment. More...
 
virtual void preTrace (const THREAD_ID, const std::shared_ptr< Ray > &)
 Called at the beginning of a trace for a ray. More...
 
void executeStudy ()
 Method for executing the study so that it can be called out of the standard UO execute() More...
 
unsigned long long int endingProcessorCrossings () const
 Total number of processor crossings for Rays that finished on this processor. More...
 
unsigned int endingMaxProcessorCrossings () const
 Max number of total processor crossings for Rays that finished on this processor. More...
 
unsigned long long int totalProcessorCrossings () const
 Total number of processor crossings. More...
 
unsigned int maxProcessorCrossings () const
 Max number of processor crossings for all Rays. More...
 
unsigned long long int endingIntersections () const
 Total number of Ray/element intersections for Rays that finished on this processor. More...
 
unsigned int endingMaxIntersections () const
 Max number of intersections for Rays that finished on this processor. More...
 
unsigned long long int totalIntersections () const
 Total number of Ray/element intersections. More...
 
unsigned int maxIntersections () const
 Max number of intersections for a Ray. More...
 
unsigned int maxTrajectoryChanges () const
 Max number of trajectory changes for a Ray. More...
 
Real endingDistance () const
 Total amount of distance traveled by the rays that end on this processor. More...
 
Real totalDistance () const
 Total distance traveled by all Rays. More...
 
unsigned long long int localTraceRayResult (const int result) const
 
const ParallelRayStudyparallelRayStudy () const
 
Real rayMaxDistance () const
 Max distance any Ray can travel. More...
 
Real executionTime ()
 Duration for execute() in seconds. More...
 
Real executionTimeNano ()
 Duration for execute() in nanoseconds. More...
 
Real generationTime () const
 Duration for creation of all Rays in seconds. More...
 
Real propagationTime () const
 Duration for creation of all Rays in seconds. More...
 
bool tolerateFailure () const
 Whether or not to tolerate failure. More...
 
bool bankRaysOnCompletion () const
 Whether or not to bank Rays on completion. More...
 
bool rayDependentSubdomainSetup () const
 Whether or not to use Ray dependent subdomain setup. More...
 
RayDataIndex registerRayData (const std::string &name)
 Register a value to be filled in the data on a Ray with a given name. More...
 
std::vector< RayDataIndexregisterRayData (const std::vector< std::string > &names)
 Register values to be filled in the data on a Ray with a given name. More...
 
RayDataIndex getRayDataIndex (const std::string &name, const bool graceful=false) const
 Gets the index associated with a registered value in the Ray data. More...
 
std::vector< RayDataIndexgetRayDataIndices (const std::vector< std::string > &names, const bool graceful=false) const
 Gets the indices associated with registered values in the Ray data. More...
 
const std::string & getRayDataName (const RayDataIndex index) const
 Gets the name associated with a registered value in the Ray data. More...
 
std::vector< std::string > getRayDataNames (const std::vector< RayDataIndex > &indices) const
 Gets the names associated with registered values in the Ray data. More...
 
std::size_t rayDataSize () const
 The registered size of values in the Ray data. More...
 
bool hasRayData () const
 Whether or not any Ray data are registered. More...
 
const std::vector< std::string > & rayDataNames () const
 The Ray data names. More...
 
RayDataIndex registerRayAuxData (const std::string &name)
 Register a value to be filled in the aux data on a Ray with a given name. More...
 
std::vector< RayDataIndexregisterRayAuxData (const std::vector< std::string > &names)
 Register values to be filled in the aux data on a Ray with a given name. More...
 
RayDataIndex getRayAuxDataIndex (const std::string &name, const bool graceful=false) const
 Gets the index associated with a registered value in the Ray aux data. More...
 
std::vector< RayDataIndexgetRayAuxDataIndices (const std::vector< std::string > &names, const bool graceful=false) const
 Gets the indices associated with registered values in the Ray aux data. More...
 
const std::string & getRayAuxDataName (const RayDataIndex index) const
 Gets the name associated with a registered value in the Ray aux data. More...
 
std::vector< std::string > getRayAuxDataNames (const std::vector< RayDataIndex > &indices) const
 Gets the names associated with registered values in the Ray aux data. More...
 
std::size_t rayAuxDataSize () const
 The registered size of values in the Ray aux data. More...
 
bool hasRayAuxData () const
 Whether or not any Ray aux data are registered. More...
 
const std::vector< std::string > & rayAuxDataNames () const
 The Ray aux data names. More...
 
bool hasRayKernels (const THREAD_ID tid)
 Whether or not there are currently any active RayKernel objects. More...
 
void getRayKernels (std::vector< RayKernelBase *> &result, SubdomainID id, THREAD_ID tid)
 Fills the active RayKernels associated with this study and a block into result. More...
 
template<typename T >
void getRayKernels (std::vector< T *> &result, THREAD_ID tid)
 Fills the active RayKernels associated with this study into result. More...
 
void getRayKernels (std::vector< RayKernelBase *> &result, SubdomainID id, THREAD_ID tid, RayID ray_id)
 Fills the active RayKernels associeted with this study, block, and potentially Ray into result. More...
 
void getRayBCs (std::vector< RayBoundaryConditionBase *> &result, BoundaryID id, THREAD_ID tid)
 Fills the active RayBCs associated with this study and a boundary into result. More...
 
template<typename T >
void getRayBCs (std::vector< T *> &result, const std::vector< BoundaryID > &ids, THREAD_ID tid)
 Fills the active RayBCs associated with this study and boundaries result. More...
 
template<typename T >
void getRayBCs (std::vector< T *> &result, THREAD_ID tid)
 Fills the active RayBCs associated with this study into result. More...
 
virtual void getRayBCs (std::vector< RayBoundaryConditionBase *> &result, const std::vector< TraceRayBndElement > &bnd_elems, THREAD_ID tid, RayID ray_id)
 Fills the active RayBCs associated with thie study, boundary elements, and potentially Ray into result. More...
 
const std::vector< RayKernelBase * > & currentRayKernels (THREAD_ID tid) const
 Gets the current RayKernels for a thread, which are set in segmentSubdomainSetup() More...
 
const BoundingBox & boundingBox () const
 Get the nodal bounding box for the domain. More...
 
const BoundingBox & looseBoundingBox () const
 Get the loose nodal bounding box for the domain. More...
 
Real domainMaxLength () const
 Get the inflated maximum length across the domain. More...
 
Real totalVolume () const
 Get the current total volume of the domain. More...
 
bool isRectangularDomain () const
 Whether or not the domain is rectangular (if it is prefectly encompassed by its bounding box) More...
 
bool hasInternalSidesets () const
 Whether or not the local mesh has internal sidesets that have RayBCs on them. More...
 
const std::vector< std::vector< BoundaryID > > & getInternalSidesets (const Elem *elem) const
 Get the internal sidesets (that have RayBC(s)) for each side for a given element. More...
 
const std::set< BoundaryID > & getInternalSidesets () const
 Gets the internal sidesets (that have RayBCs) within the local domain. More...
 
bool sideIsNonPlanar (const Elem *elem, const unsigned short s) const
 Whether or not the side on elem elem is non-planar. More...
 
bool hasSameLevelActiveElems () const
 Whether or not the mesh has active elements of the same level. More...
 
void moveRayToBufferDuringTrace (std::shared_ptr< Ray > &ray, const THREAD_ID tid, const AcquireMoveDuringTraceKey &)
 INTERNAL method for moving a Ray into the buffer during tracing. More...
 
MeshBase & meshBase () const
 Access to the libMesh MeshBase. More...
 
MooseMeshmesh ()
 
virtual const PointgetSideNormal (const Elem *elem, const unsigned short side, const THREAD_ID tid)
 Get the outward normal for a given element side. More...
 
virtual const Point * getElemNormals (const Elem *, const THREAD_ID)
 Gets the outward normals for a given element. More...
 
RayData getBankedRayData (const RayID ray_id, const RayDataIndex index) const
 Gets the data value for a banked ray with a given ID. More...
 
RayData getBankedRayAuxData (const RayID ray_id, const RayDataIndex index) const
 Gets the data value for a banked ray with a given ID. More...
 
RayID registeredRayID (const std::string &name, const bool graceful=false) const
 Gets the ID of a registered ray. More...
 
const std::string & registeredRayName (const RayID ray_id) const
 Gets the name of a registered ray. More...
 
bool useRayRegistration () const
 Whether or not ray registration is being used. More...
 
bool dataOnCacheTraces () const
 Whether or not to store the Ray data on the cached Ray traces. More...
 
bool auxDataOnCacheTraces () const
 Whether or not to store the Ray aux data on the cached Ray traces. More...
 
bool segmentsOnCacheTraces () const
 Whether or not to cache individual element segments when _cache_traces = true. More...
 
virtual bool shouldCacheTrace (const std::shared_ptr< Ray > &) const
 Virtual that allows for selection in if a Ray should be cached or not (only used when _cache_traces). More...
 
TraceDatainitThreadedCachedTrace (const std::shared_ptr< Ray > &ray, THREAD_ID tid)
 Initialize a Ray in the threaded cached trace map to be filled with segments. More...
 
const std::vector< TraceData > & getCachedTraces () const
 Get the cached trace data structure. More...
 
Real subdomainHmax (const SubdomainID subdomain_id) const
 Get the cached hmax for all elements in a subdomain. More...
 
virtual void onCompleteRay (const std::shared_ptr< Ray > &ray)
 Entry point for acting on a ray when it is completed (shouldContinue() == false) More...
 
void verifyUniqueRayIDs (const std::vector< std::shared_ptr< Ray >>::const_iterator begin, const std::vector< std::shared_ptr< Ray >>::const_iterator end, const bool global, const std::string &error_suffix) const
 Verifies that the Rays in the given range have unique Ray IDs. More...
 
void verifyUniqueRays (const std::vector< std::shared_ptr< Ray >>::const_iterator begin, const std::vector< std::shared_ptr< Ray >>::const_iterator end, const std::string &error_suffix)
 Verifies that the Rays in the given range are unique. More...
 
bool currentlyPropagating () const
 Whether or not the study is propagating (tracing Rays) More...
 
bool currentlyGenerating () const
 Whether or not the study is generating. More...
 
bool verifyRays () const
 Whether or not to verify if Rays have valid information before being traced. More...
 
bool verifyTraceIntersections () const
 Whether or not trace verification is enabled in devel/dbg modes. More...
 
bool sideIsIncoming (const Elem *const elem, const unsigned short side, const Point &direction, const THREAD_ID tid)
 Whether or not side is incoming on element elem in direction direction. More...
 
bool warnNonPlanar () const
 Whether or not to produce a warning when interacting with a non-planar mesh. More...
 
ParallelStudy< std::shared_ptr< Ray >, Ray > * parallelStudy ()
 The underlying parallel study: used for the context for calling the packed range routines. More...
 
const libMesh::ElemelemSide (const libMesh::Elem &elem, const unsigned int s, const THREAD_ID tid=0)
 Get an element's side pointer without excessive memory allocation. More...
 
virtual Real spatialValue (const Point &) const
 
virtual const std::vector< Point > spatialPoints () const
 
void setPrimaryThreadCopy (UserObject *primary)
 
UserObjectprimaryThreadCopy ()
 
SubProblemgetSubProblem () const
 
bool shouldDuplicateInitialExecution () const
 
void gatherSum (T &value)
 
void gatherMax (T &value)
 
void gatherMin (T &value)
 
void gatherProxyValueMax (T1 &proxy, T2 &value)
 
void gatherProxyValueMin (T1 &proxy, T2 &value)
 
std::set< UserObjectName > getDependObjects () const
 
const std::set< std::string > & getRequestedItems () override
 
const std::set< std::string > & getSuppliedItems () override
 
unsigned int systemNumber () const
 
virtual bool needThreadedCopy () const
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
bool isKokkosObject () const
 
MooseAppgetMooseApp () const
 
const std::string & type () const
 
const std::string & name () const
 
std::string typeAndName () const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
MooseObjectName uniqueName () const
 
const InputParametersparameters () const
 
const hit::Node * getHitNode () const
 
bool hasBase () const
 
const std::string & getBase () const
 
const TgetParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const TqueryParam (const std::string &name) const
 
const TgetRenamedParam (const std::string &old_name, const std::string &new_name) const
 
T getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool haveParameter (const std::string &name) const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &name) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
std::string messagePrefix (const bool hit_prefix=true) const
 
std::string errorPrefix (const std::string &) const
 
void mooseError (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecatedNoTrace (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
void callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
virtual void customSetup (const ExecFlagType &)
 
const ExecFlagEnumgetExecuteOnEnum () const
 
UserObjectName getUserObjectName (const std::string &param_name) const
 
const TgetUserObject (const std::string &param_name, bool is_dependency=true) const
 
const TgetUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const UserObjectBasegetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectBasegetUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const std::vector< MooseVariableScalar *> & getCoupledMooseScalarVars ()
 
const std::set< TagID > & getScalarVariableCoupleableVectorTags () const
 
const std::set< TagID > & getScalarVariableCoupleableMatrixTags () const
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyByName (const std::string &prop_name_in)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOldByName (const std::string &prop_name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOlderByName (const std::string &prop_name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialProperty (const std::string &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOld (const std::string &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOlder (const std::string &name)
 
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > > getBlockMaterialProperty (const MaterialPropertyName &name)
 
std::pair< Moose::Kokkos::MaterialProperty< T, dimension >, std::set< SubdomainID > > getKokkosBlockMaterialProperty (const MaterialPropertyName &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty ()
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialPropertyByName (const std::string &prop_name)
 
const MaterialProperty< T > & getZeroMaterialProperty (Ts... args)
 
std::set< SubdomainIDgetMaterialPropertyBlocks (const std::string &name)
 
std::vector< SubdomainName > getMaterialPropertyBlockNames (const std::string &name)
 
std::set< BoundaryIDgetMaterialPropertyBoundaryIDs (const std::string &name)
 
std::vector< BoundaryName > getMaterialPropertyBoundaryNames (const std::string &name)
 
void checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)
 
std::unordered_map< SubdomainID, std::vector< MaterialBase *> > buildRequiredMaterials (bool allow_stateful=true)
 
void statefulPropertiesAllowed (bool)
 
virtual bool getMaterialPropertyCalled () const
 
virtual const std::unordered_set< unsigned int > & getMatPropDependencies () const
 
virtual void resolveOptionalProperties ()
 
const GenericMaterialProperty< T, is_ad > & getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)
 
bool isImplicit ()
 
Moose::StateArg determineState () const
 
virtual void threadJoin (const UserObject &) override
 
virtual void threadJoin (const UserObject &) override
 
virtual void subdomainSetup () override
 
virtual void subdomainSetup () override
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
const GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false)
 
bool hasMaterialProperty (const std::string &name)
 
bool hasMaterialProperty (const std::string &name)
 
bool hasMaterialPropertyByName (const std::string &name)
 
bool hasMaterialPropertyByName (const std::string &name)
 
bool hasADMaterialProperty (const std::string &name)
 
bool hasADMaterialProperty (const std::string &name)
 
bool hasADMaterialPropertyByName (const std::string &name)
 
bool hasADMaterialPropertyByName (const std::string &name)
 
bool hasKokkosMaterialProperty (const std::string &name)
 
bool hasKokkosMaterialProperty (const std::string &name)
 
bool hasKokkosMaterialPropertyByName (const std::string &name)
 
bool hasKokkosMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
const FunctiongetFunction (const std::string &name) const
 
const FunctiongetFunctionByName (const FunctionName &name) const
 
bool hasFunction (const std::string &param_name) const
 
bool hasFunctionByName (const FunctionName &name) const
 
Moose::Kokkos::Function getKokkosFunction (const std::string &name) const
 
const TgetKokkosFunction (const std::string &name) const
 
Moose::Kokkos::Function getKokkosFunctionByName (const FunctionName &name) const
 
const TgetKokkosFunctionByName (const FunctionName &name) const
 
bool hasKokkosFunction (const std::string &param_name) const
 
bool hasKokkosFunctionByName (const FunctionName &name) const
 
bool isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessor (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessorByName (const PostprocessorName &name) const
 
std::size_t coupledPostprocessors (const std::string &param_name) const
 
const PostprocessorName & getPostprocessorName (const std::string &param_name, const unsigned int index=0) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name) const
 
const VectorPostprocessorName & getVectorPostprocessorName (const std::string &param_name) const
 
TgetSampler (const std::string &name)
 
SamplergetSampler (const std::string &name)
 
TgetSamplerByName (const SamplerName &name)
 
SamplergetSamplerByName (const SamplerName &name)
 
virtual void meshDisplaced ()
 
PerfGraphperfGraph ()
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const TgetDistribution (const std::string &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const TgetDistribution (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const TgetDistributionByName (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const TgetDistributionByName (const std::string &name) const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 
std::shared_ptr< RayacquireRayDuringTrace (const THREAD_ID tid, const AcquireMoveDuringTraceKey &)
 INTERNAL methods for acquiring a Ray during a trace in RayKernels and RayBCs. More...
 
std::shared_ptr< RayacquireRayInternal (const RayID id, const std::size_t data_size, const std::size_t aux_data_size, const bool reset, const AcquireRayInternalKey &)
 
TraceRaytraceRay (const THREAD_ID tid)
 Gets the threaded TraceRay object for tid. More...
 
const TraceRaytraceRay (const THREAD_ID tid) const
 

Static Public Member Functions

static InputParameters validParams ()
 
static void callMooseError (MooseApp *const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node *node, const bool show_trace=true)
 
static void sort (typename std::vector< T > &vector)
 
static void sortDFS (typename std::vector< T > &vector)
 
static void cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header, NameFunc &&name_func)
 
static void cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)
 

Public Attributes

 usingCombinedWarningSolutionWarnings
 
const ConsoleStream _console
 

Static Public Attributes

static const std::string type_param
 
static const std::string name_param
 
static const std::string unique_name_param
 
static const std::string app_param
 
static const std::string moose_base_param
 
static const std::string kokkos_object_param
 
static constexpr PropertyValue::id_type default_property_id
 
static constexpr PropertyValue::id_type zero_property_id
 
static constexpr auto SYSTEM
 
static constexpr auto NAME
 

Protected Member Functions

virtual void generateRays ()=0
 Subclasses should override this to determine how to generate Rays. More...
 
virtual void preExecuteStudy ()
 Entry point before study execution. More...
 
virtual void postExecuteStudy ()
 Entry point after study execution. More...
 
Real computeTotalVolume ()
 Helper function for computing the total domain volume. More...
 
std::vector< RayKernelBase * > & currentRayKernelsWrite (THREAD_ID tid)
 Gets the writeable current RayKernels for a thread. More...
 
void reserveRayBuffer (const std::size_t size)
 Reserve size entires in the Ray buffer. More...
 
bool sameLevelActiveElems () const
 Determine whether or not the mesh currently has active elements that are all the same level. More...
 
virtual void buildSegmentQuadrature (const Point &start, const Point &end, const Real length, std::vector< Point > &points, std::vector< Real > &weights) const
 Builds quadrature points for a given segment using the _segment_qrule. More...
 
const std::vector< std::shared_ptr< Ray > > & rayBank () const
 Get the Ray bank. More...
 
std::shared_ptr< RaygetBankedRay (const RayID ray_id) const
 Gets the Ray with the ID ray_id from the Ray bank. More...
 
void resetUniqueRayIDs ()
 Resets the generation of unique RayIDs via generateUniqueRayID() to the beginning of the range. More...
 
void resetReplicatedRayIDs ()
 Resets the generation of unique replicated RayIDs accessed via generateReplicatedRayID(). More...
 
std::vector< RayTracingObject * > getRayTracingObjects ()
 Gets all of the currently active RayTracingObjects. More...
 
virtual RayID generateUniqueRayID (const THREAD_ID tid)
 Generates a unique RayID to be used for a Ray. More...
 
RayID generateReplicatedRayID ()
 Generates a Ray ID that is replicated across all processors. More...
 
void moveRayToBuffer (std::shared_ptr< Ray > &ray)
 Moves a ray to the buffer to be traced during generateRays(). More...
 
void moveRaysToBuffer (std::vector< std::shared_ptr< Ray >> &rays)
 Moves rays to the buffer to be traced during generateRays(). More...
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &name) const override
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override
 
virtual void addUserObjectDependencyHelper (const UserObjectBase &uo) const override
 
void addReporterDependencyHelper (const ReporterName &reporter_name) override
 
void flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
 
InvalidSolutionID registerInvalidSolutionInternal (const std::string &message, const bool warning) const
 
const ReporterContextBasegetReporterContextBaseByName (const ReporterName &reporter_name) const
 
const ReporterNamegetReporterName (const std::string &param_name) const
 
TdeclareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< TdeclareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const TgetRestartableData (const std::string &data_name) const
 
TdeclareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
TdeclareRecoverableData (const std::string &data_name, Args &&... args)
 
TdeclareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
TdeclareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
 
std::string restartableName (const std::string &data_name) const
 
const TgetMeshProperty (const std::string &data_name, const std::string &prefix)
 
const TgetMeshProperty (const std::string &data_name)
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name) const
 
bool hasMeshProperty (const std::string &data_name) const
 
std::string meshPropertyName (const std::string &data_name) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const
 
std::string timedSectionName (const std::string &section_name) const
 
bool isCoupledScalar (const std::string &var_name, unsigned int i=0) const
 
unsigned int coupledScalarComponents (const std::string &var_name) const
 
unsigned int coupledScalar (const std::string &var_name, unsigned int comp=0) const
 
libMesh::Order coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< is_ad > & coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const VariableValuecoupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalargetScalarVar (const std::string &var_name, unsigned int comp) const
 
virtual void checkMaterialProperty (const std::string &name, const unsigned int state)
 
virtual void getKokkosMaterialPropertyHook (const std::string &, const unsigned int)
 
void markMatPropRequested (const std::string &)
 
MaterialPropertyName getMaterialPropertyName (const std::string &name) const
 
void checkExecutionStage ()
 
const TgetReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const TgetReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const TgetReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const TgetReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 
std::shared_ptr< RayacquireRay ()
 User APIs for constructing Rays within the RayTracingStudy. More...
 
std::shared_ptr< RayacquireUnsizedRay ()
 Acquire a Ray from the pool of Rays within generateRays(), without resizing the data (sizes the data to zero). More...
 
std::shared_ptr< RayacquireReplicatedRay ()
 Acquire a Ray from the pool of Rays within generateRays() in a replicated fashion. More...
 
std::shared_ptr< RayacquireCopiedRay (const Ray &ray)
 Acquires a Ray that that is copied from another Ray within generateRays(). More...
 
std::shared_ptr< RayacquireRegisteredRay (const std::string &name)
 Acquires a Ray with a given name within generateRays(). More...
 

Static Protected Member Functions

static std::string meshPropertyName (const std::string &data_name, const std::string &prefix)
 

Protected Attributes

MooseMesh_mesh
 The Mesh. More...
 
const Parallel::Communicator & _comm
 The Communicator. More...
 
const processor_id_type _pid
 The rank of this processor (this actually takes time to lookup - so just do it once) More...
 
const bool _ray_kernel_coverage_check
 Whether or not to perform coverage checks on RayKernels. More...
 
const bool _warn_non_planar
 Whether not to warn if non-planar faces are found. More...
 
const bool _use_ray_registration
 Whether or not to use Ray registration. More...
 
const bool _use_internal_sidesets
 Whether or not to use the internal sidesets in ray tracing. More...
 
const bool _tolerate_failure
 Whether or not to tolerate a Ray Tracing failure. More...
 
const bool _bank_rays_on_completion
 Whether or not to bank rays on completion. More...
 
const bool _ray_dependent_subdomain_setup
 Whether or not subdomain setup is dependent on the Ray. More...
 
bool _always_cache_traces
 Whether or not to cache traces on every trace execution. More...
 
const bool _data_on_cache_traces
 Whether or not to store the Ray data on the cache traces. More...
 
const bool _aux_data_on_cache_traces
 Whether or not to store the Ray aux data on the cache traces. More...
 
const bool _segments_on_cache_traces
 Whether or not to cache individual element segments when caching. More...
 
const Real _ray_max_distance
 Max distance a Ray can travel before being killed (can change) More...
 
const bool _verify_rays
 Whether or not to verify if Rays have valid information before being traced. More...
 
const bool _verify_trace_intersections
 Whether or not to verify the trace intersections in devel and dbg modes. More...
 
const Moose::CoordinateSystemType_coord_sys
 
const THREAD_ID _tid
 
SubProblem_subproblem
 
FEProblemBase_fe_problem
 
SystemBase_sys
 
Assembly_assembly
 
const bool _duplicate_initial_execution
 
std::set< std::string > _depend_uo
 
const bool & _enabled
 
MooseApp_app
 
Factory_factory
 
ActionFactory_action_factory
 
const std::string & _type
 
const std::string & _name
 
const InputParameters_pars
 
const ExecFlagEnum_execute_enum
 
const ExecFlagType_current_execute_flag
 
MooseApp_restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase_mci_feproblem
 
FEProblemBase_mdi_feproblem
 
MooseApp_pg_moose_app
 
const std::string _prefix
 
FEProblemBase_sc_fe_problem
 
const THREAD_ID _sc_tid
 
const Real_real_zero
 
const VariableValue_scalar_zero
 
const Point & _point_zero
 
const InputParameters_mi_params
 
const std::string _mi_name
 
const MooseObjectName _mi_moose_object_name
 
FEProblemBase_mi_feproblem
 
SubProblem_mi_subproblem
 
const THREAD_ID _mi_tid
 
const bool _is_kokkos_object
 
const Moose::MaterialDataType _material_data_type
 
MaterialData_material_data
 
bool _stateful_allowed
 
bool _get_material_property_called
 
std::vector< std::unique_ptr< PropertyValue > > _default_properties
 
std::unordered_set< unsigned int_material_property_dependencies
 
const MaterialPropertyName _get_suffix
 
const bool _use_interpolated_state
 
const InputParameters_ti_params
 
FEProblemBase_ti_feproblem
 
bool _is_implicit
 
Real_t
 
const Real_t_old
 
int_t_step
 
Real_dt
 
Real_dt_old
 
bool _is_transient
 
const Parallel::Communicator & _communicator
 

Static Protected Attributes

static const std::string _interpolated_old
 
static const std::string _interpolated_older
 

Private Member Functions

void coverageChecks ()
 Perform coverage checks (coverage of RayMaterials and RayKernels, if enabled) More...
 
void dependencyChecks ()
 Perform checks to see if the listed dependencies in the RayTracingObjects exist. More...
 
void verifyDependenciesExist (const std::vector< RayTracingObject *> &rtos)
 Verifies that the dependencies exist for a set of RayTracingObjects. More...
 
void traceableMeshChecks ()
 Check for if all of the element types in the mesh are supported by ray tracing. More...
 
void periodicBoundaryChecks ()
 Check for overlapping PeriodicRayBC boundaries and check for cases in which ghosting may not be sufficient with distributed mesh for overlapping periodic boundaries. More...
 
void internalSidesetSetup ()
 Does the setup for internal sidesets. More...
 
void nonPlanarSideSetup ()
 Sets up the caching of whether or not each element side is non-planar, which is stored in _non_planar_sides. More...
 
void localElemIndexSetup ()
 Sets up the _elem_index_helper, which is used for obtaining a contiguous index for all elements that this processor knows about. More...
 
void registeredRaySetup ()
 Sets up the maps from Ray to associated RayTracingObjects if _use_ray_registration. More...
 
void zeroAuxVariables ()
 Zero the AuxVariables that the registered AuxRayKernels contribute to. More...
 
void subdomainHMaxSetup ()
 Caches the hmax for all elements in each subdomain. More...
 
RayDataIndex registerRayDataInternal (const std::string &name, const bool aux)
 Internal method for registering Ray data or Ray aux data with a name. More...
 
std::vector< RayDataIndexregisterRayDataInternal (const std::vector< std::string > &names, const bool aux)
 Internal method for registering Ray data or Ray aux data with names. More...
 
RayDataIndex getRayDataIndexInternal (const std::string &name, const bool aux, const bool graceful) const
 Internal method for getting the index of Ray data or Ray aux data. More...
 
std::vector< RayDataIndexgetRayDataIndicesInternal (const std::vector< std::string > &names, const bool aux, const bool graceful) const
 Internal method for getting the indicies of Ray data or Ray aux data. More...
 
const std::string & getRayDataNameInternal (const RayDataIndex index, const bool aux) const
 Internal method for getting the name of Ray data or Ray aux data. More...
 
RayData getBankedRayDataInternal (const RayID ray_id, const RayDataIndex index, const bool aux) const
 Internal method for getting the value (replicated across all processors) in a Ray's data or aux data from the Ray banks. More...
 
RayID registerRay (const std::string &name)
 Registers a Ray with a given name. More...
 

Private Attributes

std::vector< libMesh::ElemSideBuilder_threaded_elem_side_builders
 Threaded helpers for building element sides without extraneous allocation. More...
 
std::unordered_map< std::string, RayDataIndex_ray_data_map
 The map from Ray data names to index. More...
 
std::unordered_map< std::string, RayDataIndex_ray_aux_data_map
 The map from Ray aux data names to index. More...
 
std::vector< std::string > _ray_data_names
 The names for each Ray data entry. More...
 
std::vector< std::string > _ray_aux_data_names
 The names for each Ray aux data entry. More...
 
std::unordered_map< std::string, RayID > & _registered_ray_map
 Map from registered Ray name to ID. More...
 
std::vector< std::string > & _reverse_registered_ray_map
 Map from registered Ray ID to name. More...
 
std::vector< TraceData_cached_traces
 Storage for the cached traces. More...
 
std::vector< std::vector< TraceData > > _threaded_cached_traces
 The threaded storage for cached traces. More...
 
std::vector< std::size_t > _num_cached
 Number of currently cached objects for Jacobian/residual for each thread. More...
 
std::set< BoundaryID_internal_sidesets
 The BoundaryIDs on the local mesh that have internal RayBCs. More...
 
std::vector< std::vector< std::vector< BoundaryID > > > _internal_sidesets_map
 Internal sideset data, if internal sidesets exist (indexed with getLocalElemIndex()) More...
 
bool _has_non_planar_sides
 Whether or not the local mesh has elements with non-planar sides. More...
 
std::vector< std::vector< unsigned short > > _non_planar_sides
 Non planar side data, which is for quick checking if an elem side is non-planar We use unsigned short here to avoid a std::vector<bool>; 0 = false, otherwise true. More...
 
bool _has_same_level_active_elems
 Whether or not the mesh has active elements of the same level. More...
 
libMesh::BoundingBox _b_box
 Nodal bounding box for the domain. More...
 
libMesh::BoundingBox _loose_b_box
 Loose nodal bounding box for the domain. More...
 
Real _domain_max_length
 An inflated max distance for the domain. More...
 
Real _total_volume
 The total volume of the domain. More...
 
std::vector< TheWarehouse::QueryCache< AttribSubdomains > > _threaded_cache_ray_kernel
 Threaded cached subdomain query for RayKernelBase objects pertaining to this study. More...
 
std::vector< TheWarehouse::QueryCache< AttribBoundaries > > _threaded_cache_ray_bc
 Threaded cached boundary query for RayBC objects pertaining to this study. More...
 
std::vector< std::vector< std::set< const RayTracingObject * > > > _threaded_ray_object_registration
 Threaded storage for all of the RayTracingObjects associated with a single Ray. More...
 
std::vector< std::vector< RayKernelBase * > > _threaded_current_ray_kernels
 The current RayKernel objects for each thread. More...
 
std::vector< std::shared_ptr< TraceRay > > _threaded_trace_ray
 The TraceRay objects for each thread (they do the physical tracing) More...
 
std::vector< std::unique_ptr< libMesh::FEBase > > _threaded_fe_face
 Face FE used for computing face normals for each thread. More...
 
std::vector< std::unique_ptr< libMesh::QBase > > _threaded_q_face
 Face quadrature used for computing face normals for each thread. More...
 
std::vector< std::unordered_map< std::pair< const Elem *, unsigned short >, Point > > _threaded_cached_normals
 Threaded cache for side normals that have been computed already during tracing. More...
 
std::vector< std::shared_ptr< Ray > > _ray_bank
 Cumulative Ray bank - stored only when _bank_rays_on_completion. More...
 
std::vector< RayID_threaded_next_ray_id
 Storage for the next available unique RayID, obtained via generateUniqueRayID() More...
 
RayID _replicated_next_ray_id
 Storage for the next available replicated RayID, obtained via generateReplicatedRayID() More...
 
const std::unique_ptr< ParallelRayStudy_parallel_ray_study
 The study that used is to actually execute (trace) the Rays. More...
 
std::unique_ptr< libMesh::QBase_segment_qrule
 Quadrature rule for laying points across a 1D ray segment. More...
 
unsigned long long int _ending_processor_crossings
 Total number of processor crossings for Rays that finished on this processor. More...
 
unsigned int _ending_max_processor_crossings
 Max number of total processor crossings for Rays that finished on this processor. More...
 
unsigned long long int _total_processor_crossings
 Total number of processor crossings. More...
 
unsigned int _max_processor_crossings
 Max number of processor crossings for all Rays. More...
 
unsigned long long int _ending_intersections
 Total number of Ray/element intersections for Rays that finished on this processor. More...
 
unsigned int _ending_max_intersections
 Max number of intersections for Rays that finished on this processor. More...
 
unsigned int _ending_max_trajectory_changes
 Max number of trajectory changes for Rays that finished on this processor. More...
 
unsigned long long int _total_intersections
 Total number of Ray/element intersections. More...
 
unsigned int _max_intersections
 Max number of intersections for a single Ray. More...
 
unsigned int _max_trajectory_changes
 Max number of trajectory changes for a single Ray. More...
 
Real _ending_distance
 Total distance traveled by Rays that end on this processor. More...
 
Real _total_distance
 Total distance traveled by all Rays. More...
 
std::vector< unsigned long long int_local_trace_ray_results
 Cumulative results on this processor from the threaded TraceRay objects. More...
 
std::unordered_map< SubdomainID, Real_subdomain_hmax
 The cached hmax for all elements in a subdomain. More...
 
bool _called_initial_setup
 Whether or not we've called initial setup - used to stop from late registration. More...
 
ElemIndexHelper _elem_index_helper
 Helper for defining a local contiguous index for each element. More...
 
Threads::spin_mutex _spin_mutex
 Spin mutex object for locks. More...
 
std::chrono::steady_clock::time_point _execution_start_time
 Timing. More...
 
std::chrono::steady_clock::duration _execution_time
 
std::chrono::steady_clock::duration _generation_time
 
std::chrono::steady_clock::duration _propagation_time
 

Detailed Description

Base class for Ray tracing studies that will generate Rays and then propagate all of them to termination.

Subclasses must override generateRays()

Definition at line 40 of file RayTracingStudy.h.

Constructor & Destructor Documentation

◆ RayTracingStudy()

RayTracingStudy::RayTracingStudy ( const InputParameters parameters)

Definition at line 130 of file RayTracingStudy.C.

133  _comm(_mesh.comm()),
134  _pid(_comm.rank()),
135 
136  _ray_kernel_coverage_check(getParam<bool>("ray_kernel_coverage_check")),
137  _warn_non_planar(getParam<bool>("warn_non_planar")),
138  _use_ray_registration(getParam<bool>("_use_ray_registration")),
139  _use_internal_sidesets(getParam<bool>("use_internal_sidesets")),
140  _tolerate_failure(getParam<bool>("tolerate_failure")),
141  _bank_rays_on_completion(getParam<bool>("_bank_rays_on_completion")),
142  _ray_dependent_subdomain_setup(getParam<bool>("_ray_dependent_subdomain_setup")),
143 
144  _always_cache_traces(getParam<bool>("always_cache_traces")),
145  _data_on_cache_traces(getParam<bool>("data_on_cache_traces")),
146  _aux_data_on_cache_traces(getParam<bool>("aux_data_on_cache_traces")),
147  _segments_on_cache_traces(getParam<bool>("segments_on_cache_traces")),
148  _ray_max_distance(getParam<Real>("ray_distance")),
149  _verify_rays(getParam<bool>("verify_rays")),
150 #ifndef NDEBUG
151  _verify_trace_intersections(getParam<bool>("verify_trace_intersections")),
152 #endif
153 
155 
157  declareRestartableData<std::unordered_map<std::string, RayID>>("registered_ray_map")),
159  declareRestartableData<std::vector<std::string>>("reverse_registered_ray_map")),
160 
162 
164 
165  _has_non_planar_sides(true),
167 
168  _b_box(MeshTools::create_nodal_bounding_box(_mesh.getMesh())),
169  _domain_max_length(1.01 * (_b_box.max() - _b_box.min()).norm()),
171 
181 
182  _parallel_ray_study(std::make_unique<ParallelRayStudy>(*this, _threaded_trace_ray)),
183 
185 
186  _called_initial_setup(false),
187 
188  _elem_index_helper(_mesh.getMesh(), name() + "_elem_index")
189 {
190  // Initialize a tracing object for each thread
191  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
192  {
193  // Initialize a tracing object for each thread
194  _threaded_trace_ray[tid] = std::make_shared<TraceRay>(*this, tid);
195 
196  // Setup the face FEs for normal computation on the fly
197  _threaded_fe_face[tid] = FEBase::build(_mesh.dimension(), FEType(CONSTANT, MONOMIAL));
198  _threaded_q_face[tid] = QBase::build(QGAUSS, _mesh.dimension() - 1, CONSTANT);
199  _threaded_fe_face[tid]->attach_quadrature_rule(_threaded_q_face[tid].get());
200  _threaded_fe_face[tid]->get_normals();
201  }
202 
203  // Evaluating on residual and Jacobian evaluation
204  if (_execute_enum.isValueSet(EXEC_PRE_KERNELS))
205  {
206  if (!getParam<bool>("allow_other_flags_with_prekernels") && _execute_enum.size() > 1)
207  paramError("execute_on",
208  "PRE_KERNELS cannot be mixed with any other execution flag.\nThat is, you cannot "
209  "currently "
210  "mix RayKernels that contribute to the Jacobian/residual with those that do not.");
211 
212  if (_app.useEigenvalue())
213  mooseError("Execution on residual and Jacobian evaluation (execute_on = PRE_KERNELS)\n",
214  "is not supported for an eigenvalue solve.");
215  }
216 
219 
220  // Scale the bounding box for loose checking
222  _loose_b_box.scale(TOLERANCE * TOLERANCE);
223 }
const bool _tolerate_failure
Whether or not to tolerate a Ray Tracing failure.
std::unique_ptr< FEGenericBase< Real > > build(const unsigned int dim, const FEType &fet)
std::vector< std::vector< std::set< const RayTracingObject * > > > _threaded_ray_object_registration
Threaded storage for all of the RayTracingObjects associated with a single Ray.
std::vector< std::string > & _reverse_registered_ray_map
Map from registered Ray ID to name.
unsigned int n_threads()
bool _has_same_level_active_elems
Whether or not the mesh has active elements of the same level.
bool sameLevelActiveElems() const
Determine whether or not the mesh currently has active elements that are all the same level...
const bool _segments_on_cache_traces
Whether or not to cache individual element segments when caching.
void paramError(const std::string &param, Args... args) const
MooseMesh & _mesh
The Mesh.
std::vector< std::vector< TraceData > > _threaded_cached_traces
The threaded storage for cached traces.
const bool _use_internal_sidesets
Whether or not to use the internal sidesets in ray tracing.
const Real _ray_max_distance
Max distance a Ray can travel before being killed (can change)
const InputParameters & parameters() const
libMesh::BoundingBox _loose_b_box
Loose nodal bounding box for the domain.
unsigned int size() const
processor_id_type rank() const
const Parallel::Communicator & comm() const
const bool _verify_trace_intersections
Whether or not to verify the trace intersections in devel and dbg modes.
T & declareRestartableData(const std::string &data_name, Args &&... args)
std::vector< std::vector< RayKernelBase * > > _threaded_current_ray_kernels
The current RayKernel objects for each thread.
std::vector< std::shared_ptr< TraceRay > > _threaded_trace_ray
The TraceRay objects for each thread (they do the physical tracing)
std::vector< RayID > _threaded_next_ray_id
Storage for the next available unique RayID, obtained via generateUniqueRayID()
FEProblemBase & _fe_problem
const bool _ray_dependent_subdomain_setup
Whether or not subdomain setup is dependent on the Ray.
CONSTANT
const std::string & name() const
bool & useEigenvalue()
void resetReplicatedRayIDs()
Resets the generation of unique replicated RayIDs accessed via generateReplicatedRayID().
const bool _data_on_cache_traces
Whether or not to store the Ray data on the cache traces.
MeshBase & getMesh()
const bool _bank_rays_on_completion
Whether or not to bank rays on completion.
bool _always_cache_traces
Whether or not to cache traces on every trace execution.
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.
virtual unsigned int dimension() const
const Point & min() const
std::vector< TheWarehouse::QueryCache< AttribSubdomains > > _threaded_cache_ray_kernel
Threaded cached subdomain query for RayKernelBase objects pertaining to this study.
std::vector< std::unique_ptr< libMesh::FEBase > > _threaded_fe_face
Face FE used for computing face normals for each thread.
std::vector< TheWarehouse::QueryCache< AttribBoundaries > > _threaded_cache_ray_bc
Threaded cached boundary query for RayBC objects pertaining to this study.
bool isValueSet(const std::string &value) const
const bool _warn_non_planar
Whether not to warn if non-planar faces are found.
libMesh::BoundingBox _b_box
Nodal bounding box for the domain.
const ExecFlagEnum & _execute_enum
ElemIndexHelper _elem_index_helper
Helper for defining a local contiguous index for each element.
std::vector< std::unordered_map< std::pair< const Elem *, unsigned short >, Point > > _threaded_cached_normals
Threaded cache for side normals that have been computed already during tracing.
Real _domain_max_length
An inflated max distance for the domain.
Real _total_volume
The total volume of the domain.
const bool _aux_data_on_cache_traces
Whether or not to store the Ray aux data on the cache traces.
bool _has_non_planar_sides
Whether or not the local mesh has elements with non-planar sides.
const bool _verify_rays
Whether or not to verify if Rays have valid information before being traced.
auto norm(const T &a)
std::vector< std::size_t > _num_cached
Number of currently cached objects for Jacobian/residual for each thread.
bool _called_initial_setup
Whether or not we&#39;ve called initial setup - used to stop from late registration.
virtual MooseMesh & mesh() override
GeneralUserObject(const InputParameters &parameters)
void mooseError(Args &&... args) const
const processor_id_type _pid
The rank of this processor (this actually takes time to lookup - so just do it once) ...
std::unordered_map< std::string, RayID > & _registered_ray_map
Map from registered Ray name to ID.
const Parallel::Communicator & _comm
The Communicator.
const bool _use_ray_registration
Whether or not to use Ray registration.
const bool _ray_kernel_coverage_check
Whether or not to perform coverage checks on RayKernels.
const Point & max() const
std::vector< unsigned long long int > _local_trace_ray_results
Cumulative results on this processor from the threaded TraceRay objects.
void scale(const Real factor)
std::vector< libMesh::ElemSideBuilder > _threaded_elem_side_builders
Threaded helpers for building element sides without extraneous allocation.
void resetUniqueRayIDs()
Resets the generation of unique RayIDs via generateUniqueRayID() to the beginning of the range...
std::vector< std::unique_ptr< libMesh::QBase > > _threaded_q_face
Face quadrature used for computing face normals for each thread.
unsigned int THREAD_ID
Real computeTotalVolume()
Helper function for computing the total domain volume.

Member Function Documentation

◆ acquireCopiedRay()

std::shared_ptr< Ray > RayTracingStudy::acquireCopiedRay ( const Ray ray)
protected

Acquires a Ray that that is copied from another Ray within generateRays().

All of the information is copied except for the counters (intersections, processor crossings, etc), which are reset.

Definition at line 1835 of file RayTracingStudy.C.

Referenced by RepeatableRayStudyBaseTest::defineRays(), RayTracingStudyTest::generateRays(), and RepeatableRayStudyBase::generateRays().

1836 {
1837  mooseAssert(currentlyGenerating(), "Can only use during generateRays()");
1838  return _parallel_ray_study->acquireParallelData(
1839  /* tid = */ 0, &ray, Ray::ConstructRayKey());
1840 }
bool currentlyGenerating() const
Whether or not the study is generating.
Class that is used as a parameter to the public constructors/reset methods.
Definition: Ray.h:99
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.

◆ acquireRay()

std::shared_ptr< Ray > RayTracingStudy::acquireRay ( )
protected

User APIs for constructing Rays within the RayTracingStudy.

Rays can ONLY be constructed by users within the RayTracingStudy via the following methods. Acquire a Ray from the pool of Rays within generateRays().

A unique ID is generated and assigned to the acquired Ray. The data and aux data sizes are set according to the sizes required by the RayTracingStudy.

Definition at line 1771 of file RayTracingStudy.C.

Referenced by LotsOfRaysRayStudy::defineRay(), StationaryRayStudyTest::generateRays(), TestRay::generateRays(), RayTracingStudyTest::generateRays(), TestReuseRaysStudy::generateRays(), TestTransientRaysStudy::generateRays(), and ViewFactorRayStudy::generateRays().

1772 {
1773  mooseAssert(currentlyGenerating(), "Can only use during generateRays()");
1774 
1775  return _parallel_ray_study->acquireParallelData(
1776  /* tid = */ 0,
1777  this,
1778  generateUniqueRayID(/* tid = */ 0),
1779  rayDataSize(),
1780  rayAuxDataSize(),
1781  /* reset = */ true,
1783 }
bool currentlyGenerating() const
Whether or not the study is generating.
Class that is used as a parameter to the public constructors/reset methods.
Definition: Ray.h:99
std::size_t rayAuxDataSize() const
The registered size of values in the Ray aux data.
virtual RayID generateUniqueRayID(const THREAD_ID tid)
Generates a unique RayID to be used for a Ray.
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.
std::size_t rayDataSize() const
The registered size of values in the Ray data.

◆ acquireRayDuringTrace()

std::shared_ptr< Ray > RayTracingStudy::acquireRayDuringTrace ( const THREAD_ID  tid,
const AcquireMoveDuringTraceKey  
)

INTERNAL methods for acquiring a Ray during a trace in RayKernels and RayBCs.

You should not use these APIs directly. If you wish to acquire a Ray during generation during generateRays()), use the protected RayTracingStudy::acquire{}Ray() methods. If you wish to acquire a Ray during propagation in RayKernels and RayBC, use the protected RayKernelBase::acquireRay() and RayBoundaryConditionBase::acquireRay(), respectively.

Definition at line 1843 of file RayTracingStudy.C.

Referenced by RayBoundaryConditionBase::acquireRay(), and RayKernelBase::acquireRay().

1844 {
1845  mooseAssert(currentlyPropagating(), "Can only use during propagation");
1846  return _parallel_ray_study->acquireParallelData(tid,
1847  this,
1848  generateUniqueRayID(tid),
1849  rayDataSize(),
1850  rayAuxDataSize(),
1851  /* reset = */ true,
1853 }
Class that is used as a parameter to the public constructors/reset methods.
Definition: Ray.h:99
std::size_t rayAuxDataSize() const
The registered size of values in the Ray aux data.
virtual RayID generateUniqueRayID(const THREAD_ID tid)
Generates a unique RayID to be used for a Ray.
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)
std::size_t rayDataSize() const
The registered size of values in the Ray data.

◆ acquireRayInternal()

std::shared_ptr<Ray> RayTracingStudy::acquireRayInternal ( const RayID  id,
const std::size_t  data_size,
const std::size_t  aux_data_size,
const bool  reset,
const AcquireRayInternalKey  
)
inline

Definition at line 489 of file RayTracingStudy.h.

Referenced by dataLoad(), and libMesh::Parallel::Packing< std::shared_ptr< Ray > >::unpack().

494  {
495  return _parallel_ray_study->acquireParallelData(
496  0, this, id, data_size, aux_data_size, reset, Ray::ConstructRayKey());
497  }
Class that is used as a parameter to the public constructors/reset methods.
Definition: Ray.h:99
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.

◆ acquireRegisteredRay()

std::shared_ptr< Ray > RayTracingStudy::acquireRegisteredRay ( const std::string &  name)
protected

Acquires a Ray with a given name within generateRays().

Used when ray registration is enabled, that is, the private paramater '_use_ray_registration' == true.

This method must be called on all processors at the same time with the same name. This method can only be called on thread 0, which is why there is no thread argument.

Definition at line 1816 of file RayTracingStudy.C.

Referenced by RepeatableRayStudyBaseTest::defineRays(), RepeatableRayStudy::defineRays(), and RayTracingStudyTest::generateRays().

1817 {
1818  mooseAssert(currentlyGenerating(), "Can only use during generateRays()");
1819 
1820  // Either register a Ray or get an already registered Ray id
1821  const RayID id = registerRay(name);
1822 
1823  // Acquire a Ray with the properly sized data initialized to zero
1824  return _parallel_ray_study->acquireParallelData(
1825  /* tid = */ 0,
1826  this,
1827  id,
1828  rayDataSize(),
1829  rayAuxDataSize(),
1830  /* reset = */ true,
1832 }
bool currentlyGenerating() const
Whether or not the study is generating.
RayID registerRay(const std::string &name)
Registers a Ray with a given name.
unsigned long int RayID
Type for a Ray&#39;s ID.
Definition: Ray.h:44
Class that is used as a parameter to the public constructors/reset methods.
Definition: Ray.h:99
std::size_t rayAuxDataSize() const
The registered size of values in the Ray aux data.
const std::string & name() const
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.
std::size_t rayDataSize() const
The registered size of values in the Ray data.

◆ acquireReplicatedRay()

std::shared_ptr< Ray > RayTracingStudy::acquireReplicatedRay ( )
protected

Acquire a Ray from the pool of Rays within generateRays() in a replicated fashion.

That is, this method must be called on all processors at the same time and the ID of the resulting Ray is the same across all processors.

The data and aux data sizes are set according to the sizes required by the RayTracingStudy.

Definition at line 1800 of file RayTracingStudy.C.

Referenced by ConeRayStudy::defineRays(), and TestPICRayStudy::generateRays().

1801 {
1802  mooseAssert(currentlyGenerating(), "Can only use during generateRays()");
1803  libmesh_parallel_only(comm());
1804 
1805  return _parallel_ray_study->acquireParallelData(
1806  /* tid = */ 0,
1807  this,
1809  rayDataSize(),
1810  rayAuxDataSize(),
1811  /* reset = */ true,
1813 }
bool currentlyGenerating() const
Whether or not the study is generating.
Class that is used as a parameter to the public constructors/reset methods.
Definition: Ray.h:99
std::size_t rayAuxDataSize() const
The registered size of values in the Ray aux data.
const Parallel::Communicator & comm() const
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.
RayID generateReplicatedRayID()
Generates a Ray ID that is replicated across all processors.
std::size_t rayDataSize() const
The registered size of values in the Ray data.

◆ acquireUnsizedRay()

std::shared_ptr< Ray > RayTracingStudy::acquireUnsizedRay ( )
protected

Acquire a Ray from the pool of Rays within generateRays(), without resizing the data (sizes the data to zero).

A unique ID is generated and assigned to the acquired Ray.

Definition at line 1786 of file RayTracingStudy.C.

Referenced by LotsOfRaysRayStudy::defineRay().

1787 {
1788  mooseAssert(currentlyGenerating(), "Can only use during generateRays()");
1789 
1790  return _parallel_ray_study->acquireParallelData(/* tid = */ 0,
1791  this,
1792  generateUniqueRayID(/* tid = */ 0),
1793  /* data_size = */ 0,
1794  /* aux_data_size = */ 0,
1795  /* reset = */ true,
1797 }
bool currentlyGenerating() const
Whether or not the study is generating.
Class that is used as a parameter to the public constructors/reset methods.
Definition: Ray.h:99
virtual RayID generateUniqueRayID(const THREAD_ID tid)
Generates a unique RayID to be used for a Ray.
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.

◆ auxDataOnCacheTraces()

bool RayTracingStudy::auxDataOnCacheTraces ( ) const
inline

Whether or not to store the Ray aux data on the cached Ray traces.

Definition at line 573 of file RayTracingStudy.h.

Referenced by TraceRay::continueTraceOffProcessor(), TraceRay::onCompleteTrace(), TraceRay::onContinueTrace(), TraceRay::onTrajectoryChanged(), RayTracingMeshOutput::RayTracingMeshOutput(), and TraceRay::trace().

573 { return _aux_data_on_cache_traces; }
const bool _aux_data_on_cache_traces
Whether or not to store the Ray aux data on the cache traces.

◆ bankRaysOnCompletion()

bool RayTracingStudy::bankRaysOnCompletion ( ) const
inline

Whether or not to bank Rays on completion.

Definition at line 214 of file RayTracingStudy.h.

Referenced by RayIntegralValue::initialize(), and RayDataValue::RayDataValue().

214 { return _bank_rays_on_completion; }
const bool _bank_rays_on_completion
Whether or not to bank rays on completion.

◆ boundingBox()

const BoundingBox& RayTracingStudy::boundingBox ( ) const
inline

Get the nodal bounding box for the domain.

Definition at line 422 of file RayTracingStudy.h.

Referenced by TestRayDataStudy::auxDataValue(), TestRayDataStudy::dataValue(), LotsOfRaysRayStudy::defineRays(), and TraceRay::trace().

422 { return _b_box; }
libMesh::BoundingBox _b_box
Nodal bounding box for the domain.

◆ buildSegmentQuadrature()

void RayTracingStudy::buildSegmentQuadrature ( const Point start,
const Point end,
const Real  length,
std::vector< Point > &  points,
std::vector< Real > &  weights 
) const
protectedvirtual

Builds quadrature points for a given segment using the _segment_qrule.

Parameters
startStart point of the segment
endEnd point of the segment
lengthThe lengh of the start -> end segment
pointsPoints to fill into (should be sized ahead of time)
weightsWeights to fill into (should be sized ahead of time)

Definition at line 814 of file RayTracingStudy.C.

Referenced by reinitSegment().

819 {
820  points.resize(_segment_qrule->n_points());
821  weights.resize(_segment_qrule->n_points());
822 
823  const Point diff = end - start;
824  const Point sum = end + start;
825  mooseAssert(MooseUtils::absoluteFuzzyEqual(length, diff.norm()), "Invalid length");
826 
827  // The standard quadrature rule should be on x = [-1, 1]
828  // To scale the points, you...
829  // - Scale to size of the segment in 3D
830  // initial_scaled_qp = x_qp * 0.5 * (end - start) = 0.5 * x_qp * diff
831  // - Shift quadrature midpoint to segment midpoint
832  // final_qp = initial_scaled_qp + 0.5 * (end - start) = initial_scaled_qp + 0.5 * sum
833  // = 0.5 * (x_qp * diff + sum)
834  for (unsigned int qp = 0; qp < _segment_qrule->n_points(); ++qp)
835  {
836  points[qp] = 0.5 * (_segment_qrule->qp(qp)(0) * diff + sum);
837  weights[qp] = 0.5 * _segment_qrule->w(qp) * length;
838  }
839 }
auto norm() const
std::unique_ptr< libMesh::QBase > _segment_qrule
Quadrature rule for laying points across a 1D ray segment.

◆ computeTotalVolume()

Real RayTracingStudy::computeTotalVolume ( )
protected

Helper function for computing the total domain volume.

Definition at line 1491 of file RayTracingStudy.C.

1492 {
1493  Real volume = 0;
1494  for (const auto & elem : *_mesh.getActiveLocalElementRange())
1495  volume += elem->volume();
1496  _communicator.sum(volume);
1497  return volume;
1498 }
libMesh::ConstElemRange * getActiveLocalElementRange()
MooseMesh & _mesh
The Mesh.
const Parallel::Communicator & _communicator
Real volume(const MeshBase &mesh, unsigned int dim=libMesh::invalid_uint)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ coverageChecks()

void RayTracingStudy::coverageChecks ( )
private

Perform coverage checks (coverage of RayMaterials and RayKernels, if enabled)

Definition at line 319 of file RayTracingStudy.C.

Referenced by initialSetup().

320 {
321  // Check for coverage of RayKernels on domain
323  {
324  std::vector<RayKernelBase *> ray_kernels;
325  getRayKernels(ray_kernels, 0);
326 
327  std::set<SubdomainID> ray_kernel_blocks;
328  for (const auto & rk : ray_kernels)
329  ray_kernel_blocks.insert(rk->blockIDs().begin(), rk->blockIDs().end());
330 
331  std::set<SubdomainID> missing;
332  std::set_difference(_mesh.meshSubdomains().begin(),
333  _mesh.meshSubdomains().end(),
334  ray_kernel_blocks.begin(),
335  ray_kernel_blocks.end(),
336  std::inserter(missing, missing.begin()));
337 
338  if (!missing.empty() && !ray_kernel_blocks.count(Moose::ANY_BLOCK_ID))
339  {
340  std::ostringstream error;
341  error << "Subdomains { ";
342  std::copy(missing.begin(), missing.end(), std::ostream_iterator<SubdomainID>(error, " "));
343  error << "} do not have RayKernels defined!";
344 
345  mooseError(error.str());
346  }
347  }
348 }
void getRayKernels(std::vector< RayKernelBase *> &result, SubdomainID id, THREAD_ID tid)
Fills the active RayKernels associated with this study and a block into result.
MooseMesh & _mesh
The Mesh.
const SubdomainID ANY_BLOCK_ID
void mooseError(Args &&... args) const
const bool _ray_kernel_coverage_check
Whether or not to perform coverage checks on RayKernels.
const std::set< SubdomainID > & meshSubdomains() const

◆ currentlyGenerating()

bool RayTracingStudy::currentlyGenerating ( ) const
inline

Whether or not the study is generating.

Definition at line 641 of file RayTracingStudy.h.

Referenced by acquireCopiedRay(), acquireRay(), acquireRegisteredRay(), acquireReplicatedRay(), acquireUnsizedRay(), dataLoad(), moveRaysToBuffer(), moveRayToBuffer(), rayBank(), reserveRayBuffer(), Ray::resetCounters(), resetReplicatedRayIDs(), and resetUniqueRayIDs().

641 { return _parallel_ray_study->currentlyPreExecuting(); }
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.

◆ currentlyPropagating()

bool RayTracingStudy::currentlyPropagating ( ) const
inline

◆ currentRayKernels()

const std::vector<RayKernelBase *>& RayTracingStudy::currentRayKernels ( THREAD_ID  tid) const
inline

Gets the current RayKernels for a thread, which are set in segmentSubdomainSetup()

Definition at line 414 of file RayTracingStudy.h.

Referenced by TraceRay::onCompleteTrace(), TraceRay::onSegment(), TraceRay::onSubdomainChanged(), reinitSegment(), and TraceRay::trace().

415  {
416  return _threaded_current_ray_kernels[tid];
417  }
std::vector< std::vector< RayKernelBase * > > _threaded_current_ray_kernels
The current RayKernel objects for each thread.

◆ currentRayKernelsWrite()

std::vector<RayKernelBase *>& RayTracingStudy::currentRayKernelsWrite ( THREAD_ID  tid)
inlineprotected

Gets the writeable current RayKernels for a thread.

Allows for other ray studies to fill the current ray kernels in a custom manner

Definition at line 721 of file RayTracingStudy.h.

722  {
723  return _threaded_current_ray_kernels[tid];
724  }
std::vector< std::vector< RayKernelBase * > > _threaded_current_ray_kernels
The current RayKernel objects for each thread.

◆ dataOnCacheTraces()

bool RayTracingStudy::dataOnCacheTraces ( ) const
inline

Whether or not to store the Ray data on the cached Ray traces.

Definition at line 569 of file RayTracingStudy.h.

Referenced by TraceRay::continueTraceOffProcessor(), TraceRay::onCompleteTrace(), TraceRay::onContinueTrace(), TraceRay::onTrajectoryChanged(), RayTracingMeshOutput::RayTracingMeshOutput(), and TraceRay::trace().

569 { return _data_on_cache_traces; }
const bool _data_on_cache_traces
Whether or not to store the Ray data on the cache traces.

◆ dependencyChecks()

void RayTracingStudy::dependencyChecks ( )
private

Perform checks to see if the listed dependencies in the RayTracingObjects exist.

Definition at line 351 of file RayTracingStudy.C.

Referenced by initialSetup().

352 {
353  std::vector<RayTracingObject *> ray_tracing_objects;
354 
355  getRayKernels(ray_tracing_objects, 0);
356  verifyDependenciesExist(ray_tracing_objects);
357 
358  getRayBCs(ray_tracing_objects, 0);
359  verifyDependenciesExist(ray_tracing_objects);
360 }
void getRayKernels(std::vector< RayKernelBase *> &result, SubdomainID id, THREAD_ID tid)
Fills the active RayKernels associated with this study and a block into result.
void verifyDependenciesExist(const std::vector< RayTracingObject *> &rtos)
Verifies that the dependencies exist for a set of RayTracingObjects.
void getRayBCs(std::vector< RayBoundaryConditionBase *> &result, BoundaryID id, THREAD_ID tid)
Fills the active RayBCs associated with this study and a boundary into result.

◆ domainMaxLength()

Real RayTracingStudy::domainMaxLength ( ) const
inline

Get the inflated maximum length across the domain.

Definition at line 433 of file RayTracingStudy.h.

Referenced by TraceRay::exitsElem(), and TraceRay::trace().

433 { return _domain_max_length; }
Real _domain_max_length
An inflated max distance for the domain.

◆ elemSide()

const libMesh::Elem& RayTracingStudy::elemSide ( const libMesh::Elem elem,
const unsigned int  s,
const THREAD_ID  tid = 0 
)
inline

Get an element's side pointer without excessive memory allocation.

Parameters
elemThe element to build a side for
sThe side to build
tidThe thread id
Returns
A pointer to the side element

Definition at line 689 of file RayTracingStudy.h.

Referenced by LotsOfRaysRayStudy::defineRays(), StationaryRayStudyTest::generateRays(), nonPlanarSideSetup(), and ClaimRays::postClaimRay().

690  {
691  return _threaded_elem_side_builders[tid](elem, s);
692  }
std::vector< libMesh::ElemSideBuilder > _threaded_elem_side_builders
Threaded helpers for building element sides without extraneous allocation.

◆ endingDistance()

Real RayTracingStudy::endingDistance ( ) const
inline

Total amount of distance traveled by the rays that end on this processor.

Definition at line 168 of file RayTracingStudy.h.

168 { return _ending_distance; }
Real _ending_distance
Total distance traveled by Rays that end on this processor.

◆ endingIntersections()

unsigned long long int RayTracingStudy::endingIntersections ( ) const
inline

Total number of Ray/element intersections for Rays that finished on this processor.

Definition at line 147 of file RayTracingStudy.h.

147 { return _ending_intersections; }
unsigned long long int _ending_intersections
Total number of Ray/element intersections for Rays that finished on this processor.

◆ endingMaxIntersections()

unsigned int RayTracingStudy::endingMaxIntersections ( ) const
inline

Max number of intersections for Rays that finished on this processor.

Definition at line 151 of file RayTracingStudy.h.

151 { return _ending_max_intersections; }
unsigned int _ending_max_intersections
Max number of intersections for Rays that finished on this processor.

◆ endingMaxProcessorCrossings()

unsigned int RayTracingStudy::endingMaxProcessorCrossings ( ) const
inline

Max number of total processor crossings for Rays that finished on this processor.

Definition at line 134 of file RayTracingStudy.h.

unsigned int _ending_max_processor_crossings
Max number of total processor crossings for Rays that finished on this processor. ...

◆ endingProcessorCrossings()

unsigned long long int RayTracingStudy::endingProcessorCrossings ( ) const
inline

Total number of processor crossings for Rays that finished on this processor.

Definition at line 130 of file RayTracingStudy.h.

130 { return _ending_processor_crossings; }
unsigned long long int _ending_processor_crossings
Total number of processor crossings for Rays that finished on this processor.

◆ execute()

void RayTracingStudy::execute ( )
overridevirtual

Executes the study (generates and propagates Rays)

Implements GeneralUserObject.

Reimplemented in SingleTraceLineSourceTest.

Definition at line 313 of file RayTracingStudy.C.

Referenced by SingleTraceLineSourceTest::execute().

314 {
315  executeStudy();
316 }
void executeStudy()
Method for executing the study so that it can be called out of the standard UO execute() ...

◆ executeStudy()

void RayTracingStudy::executeStudy ( )

Method for executing the study so that it can be called out of the standard UO execute()

Definition at line 875 of file RayTracingStudy.C.

Referenced by execute().

876 {
877  TIME_SECTION("executeStudy", 2, "Executing Study");
878 
879  mooseAssert(_called_initial_setup, "Initial setup not called");
880 
881  // Reset ray start/complete timers
883  _max_intersections = 0;
885 
886  // Reset physical tracing stats
887  for (auto & val : _local_trace_ray_results)
888  val = 0;
889 
890  // Reset crossing and intersection
897  _ending_distance = 0;
898  _total_distance = 0;
899 
900  // Zero the AuxVariables that our AuxRayKernels contribute to before they accumulate
902 
903  preExecuteStudy();
904  for (auto & rto : getRayTracingObjects())
905  rto->preExecuteStudy();
906 
907  _ray_bank.clear();
908 
909  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
910  {
911  _threaded_trace_ray[tid]->preExecute();
912  _threaded_cached_normals[tid].clear();
913  }
914 
916  _execution_start_time = std::chrono::steady_clock::now();
917 
918  _parallel_ray_study->preExecute();
919 
920  {
921  {
922  auto generation_start_time = std::chrono::steady_clock::now();
923 
924  TIME_SECTION("generateRays", 2, "Generating Rays");
925 
926  generateRays();
927 
928  _generation_time = std::chrono::steady_clock::now() - generation_start_time;
929  }
930 
931  // At this point, nobody is working so this is good time to make sure
932  // Rays are unique across all processors in the working buffer
933  if (verifyRays())
934  {
935  verifyUniqueRays(_parallel_ray_study->workBuffer().begin(),
936  _parallel_ray_study->workBuffer().end(),
937  /* error_suffix = */ "after generateRays()");
938 
939  verifyUniqueRayIDs(_parallel_ray_study->workBuffer().begin(),
940  _parallel_ray_study->workBuffer().end(),
941  /* global = */ true,
942  /* error_suffix = */ "after generateRays()");
943  }
944 
946 
947  {
948  TIME_SECTION("propagateRays", 2, "Propagating Rays");
949 
950  const auto propagation_start_time = std::chrono::steady_clock::now();
951 
952  _parallel_ray_study->execute();
953 
954  _propagation_time = std::chrono::steady_clock::now() - propagation_start_time;
955  }
956  }
957 
958  _execution_time = std::chrono::steady_clock::now() - _execution_start_time;
959 
960  if (verifyRays())
961  {
962  verifyUniqueRays(_parallel_ray_study->workBuffer().begin(),
963  _parallel_ray_study->workBuffer().end(),
964  /* error_suffix = */ "after tracing completed");
965 
966 #ifndef NDEBUG
967  // Outside of debug, _ray_bank always holds all of the Rays that have ended on this processor
968  // We can use this as a global point to check for unique IDs for every Ray that has traced
970  _ray_bank.end(),
971  /* global = */ true,
972  /* error_suffix = */ "after tracing completed");
973 #endif
974  }
975 
976  // Update counters from the threaded trace objects
977  for (const auto & tr : _threaded_trace_ray)
978  for (std::size_t i = 0; i < _local_trace_ray_results.size(); ++i)
979  _local_trace_ray_results[i] += tr->results()[i];
980 
981  // Update local ending counters
988  // ...and communicate the global values
995 
996  // Throw a warning with the number of failed (tolerated) traces
997  if (_tolerate_failure)
998  {
1000  _communicator.sum(failures);
1001  if (failures)
1002  mooseWarning(
1003  type(), " '", name(), "': ", failures, " ray tracing failures were tolerated.\n");
1004  }
1005 
1006  // Clear the current RayKernels
1007  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
1009 
1010  // Move the threaded cache trace information into the full cached trace vector
1011  // Here, we only clear the cached vectors so that we might not have to
1012  // reallocate on future traces
1013  std::size_t num_entries = 0;
1014  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
1015  num_entries += _threaded_cached_traces[tid].size();
1016  _cached_traces.clear();
1017  _cached_traces.reserve(num_entries);
1018  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
1019  {
1020  for (const auto & entry : _threaded_cached_traces[tid])
1021  _cached_traces.emplace_back(std::move(entry));
1022  _threaded_cached_traces[tid].clear();
1023  }
1024 
1025  // Add any stragglers that contribute to the Jacobian or residual
1026  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
1027  if (_num_cached[tid] != 0)
1028  {
1029  mooseAssert(_fe_problem.currentlyComputingJacobian() ||
1031  "Should not have cached values without Jacobian/residual computation");
1032 
1035  else
1037 
1038  _num_cached[tid] = 0;
1039  }
1040 
1041  // AuxRayKernels may have modified AuxVariables
1044 
1045  // Clear FE
1046  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
1047  {
1050  }
1051 
1052  postExecuteStudy();
1053  for (auto & rto : getRayTracingObjects())
1054  rto->postExecuteStudy();
1055 }
const bool _tolerate_failure
Whether or not to tolerate a Ray Tracing failure.
unsigned long long int _total_intersections
Total number of Ray/element intersections.
unsigned long long int _ending_intersections
Total number of Ray/element intersections for Rays that finished on this processor.
unsigned int _max_trajectory_changes
Max number of trajectory changes for a single Ray.
void verifyUniqueRayIDs(const std::vector< std::shared_ptr< Ray >>::const_iterator begin, const std::vector< std::shared_ptr< Ray >>::const_iterator end, const bool global, const std::string &error_suffix) const
Verifies that the Rays in the given range have unique Ray IDs.
unsigned int n_threads()
NumericVector< Number > & solution()
std::vector< std::vector< TraceData > > _threaded_cached_traces
The threaded storage for cached traces.
Real _ending_distance
Total distance traveled by Rays that end on this processor.
const Parallel::Communicator & _communicator
void clearActiveMaterialProperties(const THREAD_ID tid)
std::vector< std::vector< RayKernelBase * > > _threaded_current_ray_kernels
The current RayKernel objects for each thread.
std::vector< TraceData > _cached_traces
Storage for the cached traces.
std::vector< std::shared_ptr< TraceRay > > _threaded_trace_ray
The TraceRay objects for each thread (they do the physical tracing)
FEProblemBase & _fe_problem
unsigned long long int _total_processor_crossings
Total number of processor crossings.
bool verifyRays() const
Whether or not to verify if Rays have valid information before being traced.
const std::string & name() const
unsigned int _ending_max_intersections
Max number of intersections for Rays that finished on this processor.
sideset clear()
virtual void postExecuteStudy()
Entry point after study execution.
unsigned int _max_processor_crossings
Max number of processor crossings for all Rays.
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.
std::vector< RayTracingObject * > getRayTracingObjects()
Gets all of the currently active RayTracingObjects.
const bool & currentlyComputingResidual() const
unsigned int _max_intersections
Max number of intersections for a single Ray.
unsigned int _ending_max_trajectory_changes
Max number of trajectory changes for Rays that finished on this processor.
unsigned int _ending_max_processor_crossings
Max number of total processor crossings for Rays that finished on this processor. ...
void zeroAuxVariables()
Zero the AuxVariables that the registered AuxRayKernels contribute to.
const std::string & type() const
std::chrono::steady_clock::time_point _execution_start_time
Timing.
virtual void preExecuteStudy()
Entry point before study execution.
AuxiliarySystem & getAuxiliarySystem()
virtual void generateRays()=0
Subclasses should override this to determine how to generate Rays.
virtual void close()=0
unsigned long long int _ending_processor_crossings
Total number of processor crossings for Rays that finished on this processor.
std::vector< std::unordered_map< std::pair< const Elem *, unsigned short >, Point > > _threaded_cached_normals
Threaded cache for side normals that have been computed already during tracing.
void max(const T &r, T &o, Request &req) const
std::vector< std::size_t > _num_cached
Number of currently cached objects for Jacobian/residual for each thread.
std::chrono::steady_clock::duration _propagation_time
void mooseWarning(Args &&... args) const
bool _called_initial_setup
Whether or not we&#39;ve called initial setup - used to stop from late registration.
void verifyUniqueRays(const std::vector< std::shared_ptr< Ray >>::const_iterator begin, const std::vector< std::shared_ptr< Ray >>::const_iterator end, const std::string &error_suffix)
Verifies that the Rays in the given range are unique.
Real _total_distance
Total distance traveled by all Rays.
const bool & currentlyComputingJacobian() const
std::vector< std::shared_ptr< Ray > > _ray_bank
Cumulative Ray bank - stored only when _bank_rays_on_completion.
std::vector< unsigned long long int > _local_trace_ray_results
Cumulative results on this processor from the threaded TraceRay objects.
std::chrono::steady_clock::duration _generation_time
void registeredRaySetup()
Sets up the maps from Ray to associated RayTracingObjects if _use_ray_registration.
virtual void addCachedResidual(const THREAD_ID tid) override
std::chrono::steady_clock::duration _execution_time
virtual void clearActiveElementalMooseVariables(const THREAD_ID tid) override
unsigned int THREAD_ID
virtual void addCachedJacobian(const THREAD_ID tid) override

◆ executionTime()

Real RayTracingStudy::executionTime ( )
inline

Duration for execute() in seconds.

Definition at line 189 of file RayTracingStudy.h.

189 { return std::chrono::duration<Real>(_execution_time).count(); }
std::chrono::steady_clock::duration _execution_time

◆ executionTimeNano()

Real RayTracingStudy::executionTimeNano ( )
inline

Duration for execute() in nanoseconds.

Definition at line 193 of file RayTracingStudy.h.

194  {
195  return std::chrono::duration<Real, std::nano>(_execution_time).count();
196  }
std::chrono::steady_clock::duration _execution_time

◆ finalize()

virtual void RayTracingStudy::finalize ( )
inlineoverridevirtual

Implements GeneralUserObject.

Definition at line 81 of file RayTracingStudy.h.

81 {}

◆ generateRays()

virtual void RayTracingStudy::generateRays ( )
protectedpure virtual

Subclasses should override this to determine how to generate Rays.

This will be called within execute() and makes up the "generation phase" of the algorithm.

Implemented in ViewFactorRayStudy, RayTracingStudyWithRegistrationTest, RayTracingStudyNoBankingTest, RepeatableRayStudyBase, TestTransientRaysStudy, TestReuseRaysStudy, TestPICRayStudy, RayTracingStudyTest, StationaryRayStudyTest, and TestRay.

Referenced by executeStudy().

◆ generateReplicatedRayID()

RayID RayTracingStudy::generateReplicatedRayID ( )
protected

Generates a Ray ID that is replicated across all processors.

Definition at line 1754 of file RayTracingStudy.C.

Referenced by acquireReplicatedRay().

1755 {
1756  return _replicated_next_ray_id++;
1757 }
RayID _replicated_next_ray_id
Storage for the next available replicated RayID, obtained via generateReplicatedRayID() ...

◆ generateUniqueRayID()

RayID RayTracingStudy::generateUniqueRayID ( const THREAD_ID  tid)
protectedvirtual

Generates a unique RayID to be used for a Ray.

This is used internally when acquiring new Rays.

Definition at line 1729 of file RayTracingStudy.C.

Referenced by acquireRay(), acquireRayDuringTrace(), and acquireUnsizedRay().

1730 {
1731  // Get the current ID to return
1732  const auto id = _threaded_next_ray_id[tid];
1733 
1734  // Advance so that the next call has the correct ID
1736 
1737  return id;
1738 }
unsigned long int RayID
Type for a Ray&#39;s ID.
Definition: Ray.h:44
unsigned int n_threads()
std::vector< RayID > _threaded_next_ray_id
Storage for the next available unique RayID, obtained via generateUniqueRayID()
processor_id_type n_processors() const

◆ generationTime()

Real RayTracingStudy::generationTime ( ) const
inline

Duration for creation of all Rays in seconds.

Definition at line 200 of file RayTracingStudy.h.

Referenced by PerProcessorRayTracingResultsVectorPostprocessor::execute().

200 { return std::chrono::duration<Real>(_generation_time).count(); }
std::chrono::steady_clock::duration _generation_time

◆ getBankedRay()

std::shared_ptr< Ray > RayTracingStudy::getBankedRay ( const RayID  ray_id) const
protected

Gets the Ray with the ID ray_id from the Ray bank.

If the Ray with ray_id is not found across all processors, this will error.

This will ONLY return a valid Ray (not a null shared_ptr) on the processor that has the Ray.

Definition at line 1381 of file RayTracingStudy.C.

Referenced by getBankedRayDataInternal().

1382 {
1383  // This is only a linear search - can be improved on with a map in the future
1384  // if this is used on a larger scale
1385  std::shared_ptr<Ray> ray;
1386  for (const std::shared_ptr<Ray> & possible_ray : rayBank())
1387  if (possible_ray->id() == ray_id)
1388  {
1389  ray = possible_ray;
1390  break;
1391  }
1392 
1393  // Make sure one and only one processor has the Ray
1394  unsigned int have_ray = ray ? 1 : 0;
1395  _communicator.sum(have_ray);
1396  if (have_ray == 0)
1397  mooseError("Could not find a Ray with the ID ", ray_id, " in the Ray banks.");
1398 
1399  // This should never happen... but let's make sure
1400  mooseAssert(have_ray == 1, "Multiple rays with the same ID were found in the Ray banks");
1401 
1402  return ray;
1403 }
const std::vector< std::shared_ptr< Ray > > & rayBank() const
Get the Ray bank.
const Parallel::Communicator & _communicator
void mooseError(Args &&... args) const

◆ getBankedRayAuxData()

RayData RayTracingStudy::getBankedRayAuxData ( const RayID  ray_id,
const RayDataIndex  index 
) const

Gets the data value for a banked ray with a given ID.

This will return the value replicated across all processors

Definition at line 1425 of file RayTracingStudy.C.

Referenced by RayDataValue::getValue().

1426 {
1427  return getBankedRayDataInternal(ray_id, index, /* aux = */ true);
1428 }
RayData getBankedRayDataInternal(const RayID ray_id, const RayDataIndex index, const bool aux) const
Internal method for getting the value (replicated across all processors) in a Ray&#39;s data or aux data ...

◆ getBankedRayData()

RayData RayTracingStudy::getBankedRayData ( const RayID  ray_id,
const RayDataIndex  index 
) const

Gets the data value for a banked ray with a given ID.

This will return the value replicated across all processors

Definition at line 1419 of file RayTracingStudy.C.

Referenced by RayDataValue::getValue(), and RayIntegralValue::getValue().

1420 {
1421  return getBankedRayDataInternal(ray_id, index, /* aux = */ false);
1422 }
RayData getBankedRayDataInternal(const RayID ray_id, const RayDataIndex index, const bool aux) const
Internal method for getting the value (replicated across all processors) in a Ray&#39;s data or aux data ...

◆ getBankedRayDataInternal()

RayData RayTracingStudy::getBankedRayDataInternal ( const RayID  ray_id,
const RayDataIndex  index,
const bool  aux 
) const
private

Internal method for getting the value (replicated across all processors) in a Ray's data or aux data from the Ray banks.

Definition at line 1406 of file RayTracingStudy.C.

Referenced by getBankedRayAuxData(), and getBankedRayData().

1409 {
1410  // Will be a nullptr shared_ptr if this processor doesn't own the Ray
1411  const std::shared_ptr<Ray> ray = getBankedRay(ray_id);
1412 
1413  Real value = ray ? (aux ? ray->auxData(index) : ray->data(index)) : 0;
1414  _communicator.sum(value);
1415  return value;
1416 }
const Parallel::Communicator & _communicator
std::shared_ptr< Ray > getBankedRay(const RayID ray_id) const
Gets the Ray with the ID ray_id from the Ray bank.
Real value(unsigned n, unsigned alpha, unsigned beta, Real x)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ getCachedTraces()

const std::vector<TraceData>& RayTracingStudy::getCachedTraces ( ) const
inline

Get the cached trace data structure.

Definition at line 595 of file RayTracingStudy.h.

Referenced by RayTracingMeshOutput::buildIDMap(), RayTracingMeshOutput::buildSegmentMesh(), RayTracingMeshOutput::fillFields(), and RayTracingMeshOutput::output().

595 { return _cached_traces; }
std::vector< TraceData > _cached_traces
Storage for the cached traces.

◆ getElemNormals()

virtual const Point* RayTracingStudy::getElemNormals ( const Elem *  ,
const THREAD_ID   
)
inlinevirtual

Gets the outward normals for a given element.

Returns a pointer to the normal for the zeroth side.

Reimplemented in BackfaceCullingStudyTest.

Definition at line 531 of file RayTracingStudy.h.

Referenced by TraceRay::moveThroughNeighbor(), RayTracingStudyTest::RayTracingStudyTest(), and TraceRay::trace().

532  {
533  mooseError("Unimplemented element normal caching in ", type(), "::getElemNormals()");
534  }
const std::string & type() const
void mooseError(Args &&... args) const

◆ getInternalSidesets() [1/2]

const std::vector<std::vector<BoundaryID> >& RayTracingStudy::getInternalSidesets ( const Elem *  elem) const

Get the internal sidesets (that have RayBC(s)) for each side for a given element.

This will be empty if the elem does not have any internal sidesets that have RayBC(s)

Referenced by TraceRay::applyOnInternalBoundary().

◆ getInternalSidesets() [2/2]

const std::set<BoundaryID>& RayTracingStudy::getInternalSidesets ( ) const
inline

Gets the internal sidesets (that have RayBCs) within the local domain.

Definition at line 459 of file RayTracingStudy.h.

Referenced by ViewFactorRayStudy::initialSetup().

459 { return _internal_sidesets; }
std::set< BoundaryID > _internal_sidesets
The BoundaryIDs on the local mesh that have internal RayBCs.

◆ getRayAuxDataIndex()

RayDataIndex RayTracingStudy::getRayAuxDataIndex ( const std::string &  name,
const bool  graceful = false 
) const

Gets the index associated with a registered value in the Ray aux data.

Parameters
nameThe value name to get the index of
gracefulWhether or not to exit gracefully if none is found (return INVALID_RAY_DATA_INDEX)
Returns
The index for the value

Definition at line 1214 of file RayTracingStudy.C.

Referenced by RayDataValue::initialize(), RayTracingStudyTest::RayTracingStudyTest(), and RayTracingMeshOutput::setupEquationSystem().

1216 {
1217  return getRayDataIndexInternal(name, /* aux = */ true, graceful);
1218 }
const std::string & name() const
RayDataIndex getRayDataIndexInternal(const std::string &name, const bool aux, const bool graceful) const
Internal method for getting the index of Ray data or Ray aux data.

◆ getRayAuxDataIndices()

std::vector< RayDataIndex > RayTracingStudy::getRayAuxDataIndices ( const std::vector< std::string > &  names,
const bool  graceful = false 
) const

Gets the indices associated with registered values in the Ray aux data.

Parameters
namesThe value names to get the indices of
gracefulWhether or not to exit gracefully if none is found (index is filled with INVALID_RAY_DATA_INDEX)
Returns
The indices for the values

Definition at line 1221 of file RayTracingStudy.C.

1223 {
1224  return getRayDataIndicesInternal(names, /* aux = */ true, graceful);
1225 }
std::vector< RayDataIndex > getRayDataIndicesInternal(const std::vector< std::string > &names, const bool aux, const bool graceful) const
Internal method for getting the indicies of Ray data or Ray aux data.

◆ getRayAuxDataName()

const std::string & RayTracingStudy::getRayAuxDataName ( const RayDataIndex  index) const

Gets the name associated with a registered value in the Ray aux data.

Parameters
indexThe index to get the name of
Returns
The name associated with index

Definition at line 1228 of file RayTracingStudy.C.

Referenced by Ray::getInfo().

1229 {
1230  return getRayDataNameInternal(index, /* aux = */ true);
1231 }
const std::string & getRayDataNameInternal(const RayDataIndex index, const bool aux) const
Internal method for getting the name of Ray data or Ray aux data.

◆ getRayAuxDataNames()

std::vector<std::string> RayTracingStudy::getRayAuxDataNames ( const std::vector< RayDataIndex > &  indices) const

Gets the names associated with registered values in the Ray aux data.

Parameters
indicesThe indices to get the names of
Returns
The associated names

◆ getRayBCs() [1/4]

void RayTracingStudy::getRayBCs ( std::vector< RayBoundaryConditionBase *> &  result,
BoundaryID  id,
THREAD_ID  tid 
)

Fills the active RayBCs associated with this study and a boundary into result.

Definition at line 1292 of file RayTracingStudy.C.

Referenced by dependencyChecks(), getRayBCs(), ViewFactorRayStudy::initialSetup(), internalSidesetSetup(), TraceRay::onBoundary(), periodicBoundaryChecks(), and RayTracingStudyTest::RayTracingStudyTest().

1295 {
1296  // If the cache doesn't have any attributes yet, it means that we haven't set
1297  // the conditions yet. We do this so that it can be generated on the fly on first use.
1298  if (!_threaded_cache_ray_bc[tid].numAttribs())
1299  {
1300  if (!_called_initial_setup)
1301  mooseError("Should not call getRayBCs() before initialSetup()");
1302 
1303  auto query = _fe_problem.theWarehouse()
1304  .query()
1305  .condition<AttribRayTracingStudy>(this)
1306  .condition<AttribSystem>("RayBoundaryCondition")
1307  .condition<AttribThread>(tid);
1308  _threaded_cache_ray_bc[tid] = query.clone();
1309  }
1310 
1311  _threaded_cache_ray_bc[tid].queryInto(result, std::make_tuple(id, false));
1312 }
Attribute for the RayTracingStudy a RayTracingObject is associated with.
FEProblemBase & _fe_problem
TheWarehouse & theWarehouse() const
std::vector< TheWarehouse::QueryCache< AttribBoundaries > > _threaded_cache_ray_bc
Threaded cached boundary query for RayBC objects pertaining to this study.
query_obj query
Query query()
bool _called_initial_setup
Whether or not we&#39;ve called initial setup - used to stop from late registration.
void mooseError(Args &&... args) const

◆ getRayBCs() [2/4]

template<typename T >
void RayTracingStudy::getRayBCs ( std::vector< T *> &  result,
const std::vector< BoundaryID > &  ids,
THREAD_ID  tid 
)
inline

Fills the active RayBCs associated with this study and boundaries result.

Definition at line 376 of file RayTracingStudy.h.

377  {
379  .query()
380  .condition<AttribRayTracingStudy>(this)
381  .condition<AttribSystem>("RayBoundaryCondition")
382  .condition<AttribBoundaries>(ids)
383  .condition<AttribThread>(tid)
384  .queryInto(result);
385  }
Attribute for the RayTracingStudy a RayTracingObject is associated with.
FEProblemBase & _fe_problem
TheWarehouse & theWarehouse() const
Query query()

◆ getRayBCs() [3/4]

template<typename T >
void RayTracingStudy::getRayBCs ( std::vector< T *> &  result,
THREAD_ID  tid 
)
inline

Fills the active RayBCs associated with this study into result.

Definition at line 390 of file RayTracingStudy.h.

391  {
393  .query()
394  .condition<AttribRayTracingStudy>(this)
395  .condition<AttribSystem>("RayBoundaryCondition")
396  .condition<AttribThread>(tid)
397  .queryInto(result);
398  }
Attribute for the RayTracingStudy a RayTracingObject is associated with.
FEProblemBase & _fe_problem
TheWarehouse & theWarehouse() const
Query query()

◆ getRayBCs() [4/4]

void RayTracingStudy::getRayBCs ( std::vector< RayBoundaryConditionBase *> &  result,
const std::vector< TraceRayBndElement > &  bnd_elems,
THREAD_ID  tid,
RayID  ray_id 
)
virtual

Fills the active RayBCs associated with thie study, boundary elements, and potentially Ray into result.

This is purposely virtual because it allows derived studies to optimize the retrieval of RayBCs during the trace in TraceRay.

Definition at line 1315 of file RayTracingStudy.C.

1319 {
1320  // No Ray registration: no need to sift through objects
1321  if (!_use_ray_registration)
1322  {
1323  if (bnd_elems.size() == 1)
1324  getRayBCs(result, bnd_elems[0].bnd_id, tid);
1325  else
1326  {
1327  std::vector<BoundaryID> bnd_ids(bnd_elems.size());
1328  for (MooseIndex(bnd_elems.size()) i = 0; i < bnd_elems.size(); ++i)
1329  bnd_ids[i] = bnd_elems[i].bnd_id;
1330  getRayBCs(result, bnd_ids, tid);
1331  }
1332  }
1333  // Has Ray registration: only pick the objects associated with ray_id
1334  else
1335  {
1336  // Get all of the RayBCs on these boundaries
1337  std::vector<RayBoundaryConditionBase *> rbcs;
1338  if (bnd_elems.size() == 1)
1339  getRayBCs(rbcs, bnd_elems[0].bnd_id, tid);
1340  else
1341  {
1342  std::vector<BoundaryID> bnd_ids(bnd_elems.size());
1343  for (MooseIndex(bnd_elems.size()) i = 0; i < bnd_elems.size(); ++i)
1344  bnd_ids[i] = bnd_elems[i].bnd_id;
1345  getRayBCs(rbcs, bnd_ids, tid);
1346  }
1347 
1348  // The RayTracingObjects associated with this ray
1349  mooseAssert(ray_id < _threaded_ray_object_registration[tid].size(), "Not in registration");
1350  const auto & ray_id_rtos = _threaded_ray_object_registration[tid][ray_id];
1351 
1352  // The result is the union of all of the kernels and the objects associated with this Ray
1353  result.clear();
1354  for (auto rbc : rbcs)
1355  if (ray_id_rtos.count(rbc))
1356  result.push_back(rbc);
1357  }
1358 }
std::vector< std::vector< std::set< const RayTracingObject * > > > _threaded_ray_object_registration
Threaded storage for all of the RayTracingObjects associated with a single Ray.
void getRayBCs(std::vector< RayBoundaryConditionBase *> &result, BoundaryID id, THREAD_ID tid)
Fills the active RayBCs associated with this study and a boundary into result.
const bool _use_ray_registration
Whether or not to use Ray registration.

◆ getRayDataIndex()

RayDataIndex RayTracingStudy::getRayDataIndex ( const std::string &  name,
const bool  graceful = false 
) const

Gets the index associated with a registered value in the Ray data.

Parameters
nameThe value name to get the index of
gracefulWhether or not to exit gracefully if none is found (return INVALID_RAY_DATA_INDEX)
Returns
The index for the value

Definition at line 1183 of file RayTracingStudy.C.

Referenced by RayDataValue::initialize(), RayIntegralValue::initialize(), RayTracingStudyTest::RayTracingStudyTest(), and RayTracingMeshOutput::setupEquationSystem().

1184 {
1185  return getRayDataIndexInternal(name, /* aux = */ false, graceful);
1186 }
const std::string & name() const
RayDataIndex getRayDataIndexInternal(const std::string &name, const bool aux, const bool graceful) const
Internal method for getting the index of Ray data or Ray aux data.

◆ getRayDataIndexInternal()

RayDataIndex RayTracingStudy::getRayDataIndexInternal ( const std::string &  name,
const bool  aux,
const bool  graceful 
) const
private

Internal method for getting the index of Ray data or Ray aux data.

Definition at line 1119 of file RayTracingStudy.C.

Referenced by getRayAuxDataIndex(), getRayDataIndex(), and getRayDataIndicesInternal().

1122 {
1123  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
1124 
1125  const auto & map = aux ? _ray_aux_data_map : _ray_data_map;
1126  const auto find = map.find(name);
1127  if (find != map.end())
1128  return find->second;
1129 
1130  if (graceful)
1132 
1133  const auto & other_map = aux ? _ray_data_map : _ray_aux_data_map;
1134  if (other_map.find(name) != other_map.end())
1135  mooseError("Ray data with name '",
1136  name,
1137  "' was not found.\n\n",
1138  "However, Ray ",
1139  (aux ? "non-aux" : "aux"),
1140  " data with said name was found.\n",
1141  "Did you mean to use ",
1142  (aux ? "getRayDataIndex()/getRayDataIndices()?"
1143  : "getRayAuxDataIndex()/getRayAuxDataIndices()"),
1144  "?");
1145 
1146  mooseError("Unknown Ray ", (aux ? "aux " : ""), "data with name ", name);
1147 }
KOKKOS_INLINE_FUNCTION const T * find(const T &target, const T *const begin, const T *const end)
static const RayDataIndex INVALID_RAY_DATA_INDEX
Invalid index into a Ray&#39;s data.
Definition: Ray.h:212
if(subdm)
Threads::spin_mutex _spin_mutex
Spin mutex object for locks.
const std::string & name() const
std::unordered_map< std::string, RayDataIndex > _ray_aux_data_map
The map from Ray aux data names to index.
std::unordered_map< std::string, RayDataIndex > _ray_data_map
The map from Ray data names to index.
void mooseError(Args &&... args) const

◆ getRayDataIndices()

std::vector< RayDataIndex > RayTracingStudy::getRayDataIndices ( const std::vector< std::string > &  names,
const bool  graceful = false 
) const

Gets the indices associated with registered values in the Ray data.

Parameters
namesThe value names to get the indices of
gracefulWhether or not to exit gracefully if none is found (index is filled with INVALID_RAY_DATA_INDEX)
Returns
The indices for the values

Definition at line 1189 of file RayTracingStudy.C.

1191 {
1192  return getRayDataIndicesInternal(names, /* aux = */ false, graceful);
1193 }
std::vector< RayDataIndex > getRayDataIndicesInternal(const std::vector< std::string > &names, const bool aux, const bool graceful) const
Internal method for getting the indicies of Ray data or Ray aux data.

◆ getRayDataIndicesInternal()

std::vector< RayDataIndex > RayTracingStudy::getRayDataIndicesInternal ( const std::vector< std::string > &  names,
const bool  aux,
const bool  graceful 
) const
private

Internal method for getting the indicies of Ray data or Ray aux data.

Definition at line 1150 of file RayTracingStudy.C.

Referenced by getRayAuxDataIndices(), and getRayDataIndices().

1153 {
1154  std::vector<RayDataIndex> indices(names.size());
1155  for (std::size_t i = 0; i < names.size(); ++i)
1156  indices[i] = getRayDataIndexInternal(names[i], aux, graceful);
1157  return indices;
1158 }
RayDataIndex getRayDataIndexInternal(const std::string &name, const bool aux, const bool graceful) const
Internal method for getting the index of Ray data or Ray aux data.

◆ getRayDataName()

const std::string & RayTracingStudy::getRayDataName ( const RayDataIndex  index) const

Gets the name associated with a registered value in the Ray data.

Parameters
indexThe index to get the name of
Returns
The name associated with index

Definition at line 1196 of file RayTracingStudy.C.

Referenced by Ray::getInfo(), and RayTracingStudyTest::RayTracingStudyTest().

1197 {
1198  return getRayDataNameInternal(index, /* aux = */ false);
1199 }
const std::string & getRayDataNameInternal(const RayDataIndex index, const bool aux) const
Internal method for getting the name of Ray data or Ray aux data.

◆ getRayDataNameInternal()

const std::string & RayTracingStudy::getRayDataNameInternal ( const RayDataIndex  index,
const bool  aux 
) const
private

Internal method for getting the name of Ray data or Ray aux data.

Definition at line 1161 of file RayTracingStudy.C.

Referenced by getRayAuxDataName(), and getRayDataName().

1162 {
1163  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
1164 
1165  if ((aux ? rayAuxDataSize() : rayDataSize()) < index)
1166  mooseError("Unknown Ray ", aux ? "aux " : "", "data with index ", index);
1167  return aux ? _ray_aux_data_names[index] : _ray_data_names[index];
1168 }
std::size_t rayAuxDataSize() const
The registered size of values in the Ray aux data.
std::vector< std::string > _ray_aux_data_names
The names for each Ray aux data entry.
Threads::spin_mutex _spin_mutex
Spin mutex object for locks.
void mooseError(Args &&... args) const
std::vector< std::string > _ray_data_names
The names for each Ray data entry.
std::size_t rayDataSize() const
The registered size of values in the Ray data.

◆ getRayDataNames()

std::vector<std::string> RayTracingStudy::getRayDataNames ( const std::vector< RayDataIndex > &  indices) const

Gets the names associated with registered values in the Ray data.

Parameters
indicesThe indices to get the names of
Returns
The associated names

◆ getRayKernels() [1/3]

void RayTracingStudy::getRayKernels ( std::vector< RayKernelBase *> &  result,
SubdomainID  id,
THREAD_ID  tid 
)

Fills the active RayKernels associated with this study and a block into result.

Definition at line 1242 of file RayTracingStudy.C.

Referenced by coverageChecks(), dependencyChecks(), getRayKernels(), hasRayKernels(), initialSetup(), RayTracingStudyTest::RayTracingStudyTest(), segmentSubdomainSetup(), and zeroAuxVariables().

1243 {
1244  // If the cache doesn't have any attributes yet, it means that we haven't set
1245  // the conditions yet. We do this so that it can be generated on the fly on first use.
1246  if (!_threaded_cache_ray_kernel[tid].numAttribs())
1247  {
1248  if (!_called_initial_setup)
1249  mooseError("Should not call getRayKernels() before initialSetup()");
1250 
1251  auto query = _fe_problem.theWarehouse()
1252  .query()
1253  .condition<AttribRayTracingStudy>(this)
1254  .condition<AttribSystem>("RayKernel")
1255  .condition<AttribThread>(tid);
1256  _threaded_cache_ray_kernel[tid] = query.clone();
1257  }
1258 
1259  _threaded_cache_ray_kernel[tid].queryInto(result, id);
1260 }
Attribute for the RayTracingStudy a RayTracingObject is associated with.
FEProblemBase & _fe_problem
TheWarehouse & theWarehouse() const
std::vector< TheWarehouse::QueryCache< AttribSubdomains > > _threaded_cache_ray_kernel
Threaded cached subdomain query for RayKernelBase objects pertaining to this study.
query_obj query
Query query()
bool _called_initial_setup
Whether or not we&#39;ve called initial setup - used to stop from late registration.
void mooseError(Args &&... args) const

◆ getRayKernels() [2/3]

template<typename T >
void RayTracingStudy::getRayKernels ( std::vector< T *> &  result,
THREAD_ID  tid 
)
inline

Fills the active RayKernels associated with this study into result.

Definition at line 354 of file RayTracingStudy.h.

355  {
357  .query()
358  .condition<AttribRayTracingStudy>(this)
359  .condition<AttribSystem>("RayKernel")
360  .condition<AttribThread>(tid)
361  .queryInto(result);
362  }
Attribute for the RayTracingStudy a RayTracingObject is associated with.
FEProblemBase & _fe_problem
TheWarehouse & theWarehouse() const
Query query()

◆ getRayKernels() [3/3]

void RayTracingStudy::getRayKernels ( std::vector< RayKernelBase *> &  result,
SubdomainID  id,
THREAD_ID  tid,
RayID  ray_id 
)

Fills the active RayKernels associeted with this study, block, and potentially Ray into result.

Definition at line 1263 of file RayTracingStudy.C.

1267 {
1268  // No Ray registration: no need to sift through objects
1269  if (!_use_ray_registration)
1270  {
1271  getRayKernels(result, id, tid);
1272  }
1273  // Has Ray registration: only pick the objects associated with ray_id
1274  else
1275  {
1276  // Get all of the kernels on this block
1277  std::vector<RayKernelBase *> rkbs;
1278  getRayKernels(rkbs, id, tid);
1279 
1280  // The RayTracingObjects associated with this ray
1281  const auto & ray_id_rtos = _threaded_ray_object_registration[tid][ray_id];
1282 
1283  // The result is the union of all of the kernels and the objects associated with this Ray
1284  result.clear();
1285  for (auto rkb : rkbs)
1286  if (ray_id_rtos.count(rkb))
1287  result.push_back(rkb);
1288  }
1289 }
void getRayKernels(std::vector< RayKernelBase *> &result, SubdomainID id, THREAD_ID tid)
Fills the active RayKernels associated with this study and a block into result.
std::vector< std::vector< std::set< const RayTracingObject * > > > _threaded_ray_object_registration
Threaded storage for all of the RayTracingObjects associated with a single Ray.
const bool _use_ray_registration
Whether or not to use Ray registration.

◆ getRayTracingObjects()

std::vector< RayTracingObject * > RayTracingStudy::getRayTracingObjects ( )
protected

Gets all of the currently active RayTracingObjects.

Definition at line 1361 of file RayTracingStudy.C.

Referenced by executeStudy(), initialSetup(), jacobianSetup(), registeredRaySetup(), residualSetup(), and timestepSetup().

1362 {
1363  std::vector<RayTracingObject *> result;
1364  _fe_problem.theWarehouse().query().condition<AttribRayTracingStudy>(this).queryInto(result);
1365  return result;
1366 }
Attribute for the RayTracingStudy a RayTracingObject is associated with.
FEProblemBase & _fe_problem
TheWarehouse & theWarehouse() const
Query query()

◆ getSideNormal()

const Point & RayTracingStudy::getSideNormal ( const Elem elem,
const unsigned short  side,
const THREAD_ID  tid 
)
virtual

Get the outward normal for a given element side.

Definition at line 1654 of file RayTracingStudy.C.

Referenced by LotsOfRaysRayStudy::defineRays(), BackfaceCullingStudyTest::generateNormals(), ViewFactorRayStudy::generateRays(), TraceRay::getPointNeighbors(), ReflectRayBC::onBoundary(), sideIsIncoming(), and TraceRay::trace().

1655 {
1656  std::unordered_map<std::pair<const Elem *, unsigned short>, Point> & cache =
1658 
1659  // See if we've already cached this side normal
1660  const auto elem_side_pair = std::make_pair(elem, side);
1661  const auto search = cache.find(elem_side_pair);
1662 
1663  // Haven't cached this side normal: compute it and then cache it
1664  if (search == cache.end())
1665  {
1666  _threaded_fe_face[tid]->reinit(elem, side);
1667  const auto & normal = _threaded_fe_face[tid]->get_normals()[0];
1668  cache.emplace(elem_side_pair, normal);
1669  return normal;
1670  }
1671 
1672  // Have cached this side normal: simply return it
1673  return search->second;
1674 }
std::vector< std::unique_ptr< libMesh::FEBase > > _threaded_fe_face
Face FE used for computing face normals for each thread.
std::vector< std::unordered_map< std::pair< const Elem *, unsigned short >, Point > > _threaded_cached_normals
Threaded cache for side normals that have been computed already during tracing.

◆ hasInternalSidesets()

bool RayTracingStudy::hasInternalSidesets ( ) const
inline

Whether or not the local mesh has internal sidesets that have RayBCs on them.

NOTE: if useInternalSidesets() == false, this will be false even if the mesh does have internal sidesets

Definition at line 449 of file RayTracingStudy.h.

Referenced by TraceRay::applyOnInternalBoundary(), and TraceRay::trace().

449 { return _internal_sidesets.size(); }
std::set< BoundaryID > _internal_sidesets
The BoundaryIDs on the local mesh that have internal RayBCs.

◆ hasRayAuxData()

bool RayTracingStudy::hasRayAuxData ( ) const
inline

Whether or not any Ray aux data are registered.

Definition at line 336 of file RayTracingStudy.h.

336 { return _ray_aux_data_names.size(); }
std::vector< std::string > _ray_aux_data_names
The names for each Ray aux data entry.

◆ hasRayData()

bool RayTracingStudy::hasRayData ( ) const
inline

Whether or not any Ray data are registered.

Definition at line 276 of file RayTracingStudy.h.

Referenced by RayTracingMeshOutput::fillFields().

276 { return _ray_data_names.size(); }
std::vector< std::string > _ray_data_names
The names for each Ray data entry.

◆ hasRayKernels()

bool RayTracingStudy::hasRayKernels ( const THREAD_ID  tid)

Whether or not there are currently any active RayKernel objects.

Definition at line 1234 of file RayTracingStudy.C.

Referenced by ViewFactorRayStudy::initialSetup(), and TraceRay::preExecute().

1235 {
1236  std::vector<RayKernelBase *> result;
1237  getRayKernels(result, tid);
1238  return result.size();
1239 }
void getRayKernels(std::vector< RayKernelBase *> &result, SubdomainID id, THREAD_ID tid)
Fills the active RayKernels associated with this study and a block into result.

◆ hasSameLevelActiveElems()

bool RayTracingStudy::hasSameLevelActiveElems ( ) const
inline

Whether or not the mesh has active elements of the same level.

Use this over sameLevelActiveElems(), which is for internally setting _has_same_level_active_elems

Definition at line 475 of file RayTracingStudy.h.

Referenced by TraceRay::applyOnInternalBoundary(), and TraceRay::trace().

bool _has_same_level_active_elems
Whether or not the mesh has active elements of the same level.

◆ initialize()

virtual void RayTracingStudy::initialize ( )
inlineoverridevirtual

Implements GeneralUserObject.

Definition at line 80 of file RayTracingStudy.h.

80 {}

◆ initialSetup()

void RayTracingStudy::initialSetup ( )
overridevirtual

Reimplemented from GeneralUserObject.

Reimplemented in BackfaceCullingStudyTest, and ViewFactorRayStudy.

Definition at line 226 of file RayTracingStudy.C.

Referenced by ViewFactorRayStudy::initialSetup(), and BackfaceCullingStudyTest::initialSetup().

227 {
228  // Keep track of initialSetup call to avoid registration of various things
229  _called_initial_setup = true;
230 
231  // Sets up a local index for each elem this proc knows about
233 
234  // Check for RayKernel coverage
235  coverageChecks();
236 
237  // Make sure the dependencies exist, if any
239 
240  // Check for traceable element types
242 
243  // Check for sane periodic boundaries
245 
246  // Setup for internal sidesets
248 
250 
251  // Setup approximate hmax for each subdomain
253 
254  // Call initial setup on all of the objects
255  for (auto & rto : getRayTracingObjects())
256  rto->initialSetup();
257 
258  // Check for proper exec flags with RayKernels
259  std::vector<RayKernelBase *> ray_kernels;
260  getRayKernels(ray_kernels, 0);
261  for (const auto & rkb : ray_kernels)
262  if (dynamic_cast<RayKernel *>(rkb) && !_execute_enum.isValueSet(EXEC_PRE_KERNELS))
263  mooseError("This study has RayKernel objects that contribute to residuals and Jacobians.",
264  "\nIn this case, the study must use the execute_on = PRE_KERNELS");
265 
266  // Build 1D quadrature rule for along a segment
268  QBase::build(QGAUSS, 1, _fe_problem.getSystemBase(_sys.number()).getMinQuadratureOrder());
269 }
void getRayKernels(std::vector< RayKernelBase *> &result, SubdomainID id, THREAD_ID tid)
Fills the active RayKernels associated with this study and a block into result.
void internalSidesetSetup()
Does the setup for internal sidesets.
void periodicBoundaryChecks()
Check for overlapping PeriodicRayBC boundaries and check for cases in which ghosting may not be suffi...
void traceableMeshChecks()
Check for if all of the element types in the mesh are supported by ray tracing.
FEProblemBase & _fe_problem
void nonPlanarSideSetup()
Sets up the caching of whether or not each element side is non-planar, which is stored in _non_planar...
void coverageChecks()
Perform coverage checks (coverage of RayMaterials and RayKernels, if enabled)
std::vector< RayTracingObject * > getRayTracingObjects()
Gets all of the currently active RayTracingObjects.
void dependencyChecks()
Perform checks to see if the listed dependencies in the RayTracingObjects exist.
bool isValueSet(const std::string &value) const
unsigned int number() const
const ExecFlagEnum & _execute_enum
SystemBase & _sys
virtual const SystemBase & getSystemBase(const unsigned int sys_num) const
bool _called_initial_setup
Whether or not we&#39;ve called initial setup - used to stop from late registration.
void mooseError(Args &&... args) const
void localElemIndexSetup()
Sets up the _elem_index_helper, which is used for obtaining a contiguous index for all elements that ...
void subdomainHMaxSetup()
Caches the hmax for all elements in each subdomain.
std::unique_ptr< libMesh::QBase > _segment_qrule
Quadrature rule for laying points across a 1D ray segment.

◆ initThreadedCachedTrace()

TraceData & RayTracingStudy::initThreadedCachedTrace ( const std::shared_ptr< Ray > &  ray,
THREAD_ID  tid 
)

Initialize a Ray in the threaded cached trace map to be filled with segments.

Definition at line 1513 of file RayTracingStudy.C.

Referenced by TraceRay::trace().

1514 {
1515  mooseAssert(shouldCacheTrace(ray), "Not caching trace");
1516  mooseAssert(currentlyPropagating(), "Should only use while tracing");
1517 
1518  _threaded_cached_traces[tid].emplace_back(ray);
1519  return _threaded_cached_traces[tid].back();
1520 }
std::vector< std::vector< TraceData > > _threaded_cached_traces
The threaded storage for cached traces.
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)
virtual bool shouldCacheTrace(const std::shared_ptr< Ray > &) const
Virtual that allows for selection in if a Ray should be cached or not (only used when _cache_traces)...

◆ internalSidesetSetup()

void RayTracingStudy::internalSidesetSetup ( )
private

Does the setup for internal sidesets.

This includes:

  • Setting if the local mesh has internal sidesets that have RayBCs and on which boundary (stored in _internal_sidesets)
  • Sets up the internal sideset mapping that maps element sides to the BoundaryIDs on said side that have internal sidesets that have RayBCs (stored in _internal_sidesets_map)

If _use_internal_sidesets == false, this will still check to make sure that the user does not have internal sidesets with RayBCs on them, and will report an error if so.

Definition at line 516 of file RayTracingStudy.C.

Referenced by initialSetup(), and meshChanged().

517 {
518  // Even if we have _use_internal_sidesets == false, we will make sure the user didn't add RayBCs
519  // on internal boundaries
520 
521  // Clear the data structures and size the map based on the elements that we know about
522  _internal_sidesets.clear();
523  _internal_sidesets_map.clear();
525 
526  // First, we are going to store all elements with internal sidesets (if any) that have active
527  // RayBCs on them as elem -> vector of (side, vector of boundary ids)
528  for (const auto & bnd_elem : *_mesh.getBoundaryElementRange())
529  {
530  Elem * elem = bnd_elem->_elem;
531  const unsigned int side = bnd_elem->_side;
532  const auto bnd_id = bnd_elem->_bnd_id;
533 
534  // Not internal
535  const Elem * const neighbor = elem->neighbor_ptr(side);
536  if (!neighbor || neighbor == remote_elem)
537  continue;
538 
539  // No RayBCs on this sideset
540  std::vector<RayBoundaryConditionBase *> result;
541  getRayBCs(result, bnd_id, 0);
542  if (result.empty())
543  continue;
544 
545  if (neighbor->subdomain_id() == elem->subdomain_id())
546  mooseError("RayBCs exist on internal sidesets that are not bounded by a different",
547  "\nsubdomain on each side.",
548  "\n\nIn order to use RayBCs on internal sidesets, said sidesets must have",
549  "\na different subdomain on each side.");
550 
551  // Mark that this boundary is an internal sideset with RayBC(s)
552  _internal_sidesets.insert(bnd_id);
553 
554  // Get elem's entry in the internal sidset data structure
555  const auto index = _elem_index_helper.getIndex(elem);
556  auto & entry = _internal_sidesets_map[index];
557 
558  // Initialize this elem's sides if they have not been already
559  if (entry.empty())
560  entry.resize(elem->n_sides(), std::vector<BoundaryID>());
561 
562  // Add the internal boundary to the side entry
563  entry[side].push_back(bnd_id);
564  }
565 
567  mooseError("RayBCs are defined on internal sidesets, but the study is not set to use ",
568  "internal sidesets during tracing.",
569  "\n\nSet the parameter use_internal_sidesets = true to enable this capability.");
570 }
MooseMesh & _mesh
The Mesh.
const bool _use_internal_sidesets
Whether or not to use the internal sidesets in ray tracing.
libMesh::dof_id_type maxIndex() const
Gets the maximum index generated using this object.
std::set< BoundaryID > _internal_sidesets
The BoundaryIDs on the local mesh that have internal RayBCs.
libMesh::dof_id_type getIndex(const libMesh::Elem *elem) const
Get the index associated with the element elem.
ElemIndexHelper _elem_index_helper
Helper for defining a local contiguous index for each element.
virtual unsigned int n_sides() const=0
const Elem * neighbor_ptr(unsigned int i) const
subdomain_id_type subdomain_id() const
void getRayBCs(std::vector< RayBoundaryConditionBase *> &result, BoundaryID id, THREAD_ID tid)
Fills the active RayBCs associated with this study and a boundary into result.
void mooseError(Args &&... args) const
libMesh::StoredRange< MooseMesh::const_bnd_elem_iterator, const BndElement *> * getBoundaryElementRange()
std::vector< std::vector< std::vector< BoundaryID > > > _internal_sidesets_map
Internal sideset data, if internal sidesets exist (indexed with getLocalElemIndex()) ...

◆ isRectangularDomain()

bool RayTracingStudy::isRectangularDomain ( ) const

Whether or not the domain is rectangular (if it is prefectly encompassed by its bounding box)

Definition at line 1705 of file RayTracingStudy.C.

Referenced by LotsOfRaysRayStudy::LotsOfRaysRayStudy(), and TraceRay::preExecute().

1706 {
1707  Real bbox_volume = 1;
1708  for (unsigned int d = 0; d < _mesh.dimension(); ++d)
1709  bbox_volume *= std::abs(_b_box.max()(d) - _b_box.min()(d));
1710 
1711  return MooseUtils::absoluteFuzzyEqual(bbox_volume, totalVolume(), TOLERANCE);
1712 }
Real totalVolume() const
Get the current total volume of the domain.
MooseMesh & _mesh
The Mesh.
virtual unsigned int dimension() const
const Point & min() const
libMesh::BoundingBox _b_box
Nodal bounding box for the domain.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const Point & max() const

◆ jacobianSetup()

void RayTracingStudy::jacobianSetup ( )
overridevirtual

Reimplemented from GeneralUserObject.

Definition at line 282 of file RayTracingStudy.C.

283 {
284  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
285  mooseAssert(_num_cached[tid] == 0, "Cached residuals/Jacobians not empty");
286 
287  for (auto & rto : getRayTracingObjects())
288  rto->jacobianSetup();
289 }
unsigned int n_threads()
std::vector< RayTracingObject * > getRayTracingObjects()
Gets all of the currently active RayTracingObjects.
std::vector< std::size_t > _num_cached
Number of currently cached objects for Jacobian/residual for each thread.
unsigned int THREAD_ID

◆ localElemIndexSetup()

void RayTracingStudy::localElemIndexSetup ( )
private

Sets up the _elem_index_helper, which is used for obtaining a contiguous index for all elements that this processor knows about.

Said index is used for quickly accessing things like internal sidesets and checking if an element side is non-planar during tracing.

Definition at line 508 of file RayTracingStudy.C.

Referenced by initialSetup(), and meshChanged().

509 {
510  // TODO: We could probably minimize this to local active elements followed by
511  // boundary point neighbors, but if using distibuted mesh it really shouldn't matter
512  _elem_index_helper.initialize(_mesh.getMesh().active_element_ptr_range());
513 }
MooseMesh & _mesh
The Mesh.
void initialize(const libMesh::SimpleRange< libMesh::MeshBase::element_iterator > elems)
Initializes the indices in a contiguous manner for the given element range.
MeshBase & getMesh()
ElemIndexHelper _elem_index_helper
Helper for defining a local contiguous index for each element.

◆ localTraceRayResult()

unsigned long long int RayTracingStudy::localTraceRayResult ( const int  result) const
inline

Definition at line 174 of file RayTracingStudy.h.

Referenced by PerProcessorRayTracingResultsVectorPostprocessor::execute().

175  {
176  return _local_trace_ray_results[result];
177  }
std::vector< unsigned long long int > _local_trace_ray_results
Cumulative results on this processor from the threaded TraceRay objects.

◆ looseBoundingBox()

const BoundingBox& RayTracingStudy::looseBoundingBox ( ) const
inline

Get the loose nodal bounding box for the domain.

BoundingBox::contains_point does not have a tolerance, so we need to stretch the box a little bit for some contains_point checks.

Definition at line 429 of file RayTracingStudy.h.

Referenced by Ray::setStart(), and Ray::setStartingEndPoint().

429 { return _loose_b_box; }
libMesh::BoundingBox _loose_b_box
Loose nodal bounding box for the domain.

◆ maxIntersections()

unsigned int RayTracingStudy::maxIntersections ( ) const
inline

Max number of intersections for a Ray.

Definition at line 159 of file RayTracingStudy.h.

159 { return _max_intersections; }
unsigned int _max_intersections
Max number of intersections for a single Ray.

◆ maxProcessorCrossings()

unsigned int RayTracingStudy::maxProcessorCrossings ( ) const
inline

Max number of processor crossings for all Rays.

Definition at line 142 of file RayTracingStudy.h.

Referenced by RayTracingStudyResult::getValue().

142 { return _max_processor_crossings; }
unsigned int _max_processor_crossings
Max number of processor crossings for all Rays.

◆ maxTrajectoryChanges()

unsigned int RayTracingStudy::maxTrajectoryChanges ( ) const
inline

Max number of trajectory changes for a Ray.

Definition at line 163 of file RayTracingStudy.h.

163 { return _max_trajectory_changes; }
unsigned int _max_trajectory_changes
Max number of trajectory changes for a single Ray.

◆ mesh()

MooseMesh& RayTracingStudy::mesh ( )
inline
Returns
A reference to the MooseMesh associated with this study.

Definition at line 520 of file RayTracingStudy.h.

520 { return _mesh; }
MooseMesh & _mesh
The Mesh.

◆ meshBase()

MeshBase& RayTracingStudy::meshBase ( ) const
inline

◆ meshChanged()

void RayTracingStudy::meshChanged ( )
overridevirtual

Reimplemented from GeneralUserObject.

Reimplemented in RepeatableRayStudyBase, and BackfaceCullingStudyTest.

Definition at line 299 of file RayTracingStudy.C.

Referenced by RepeatableRayStudyBase::meshChanged().

300 {
305 
307 
308  for (const auto & trace_ray : _threaded_trace_ray)
309  trace_ray->meshChanged();
310 }
void internalSidesetSetup()
Does the setup for internal sidesets.
bool _has_same_level_active_elems
Whether or not the mesh has active elements of the same level.
bool sameLevelActiveElems() const
Determine whether or not the mesh currently has active elements that are all the same level...
std::vector< std::shared_ptr< TraceRay > > _threaded_trace_ray
The TraceRay objects for each thread (they do the physical tracing)
void nonPlanarSideSetup()
Sets up the caching of whether or not each element side is non-planar, which is stored in _non_planar...
void localElemIndexSetup()
Sets up the _elem_index_helper, which is used for obtaining a contiguous index for all elements that ...
void subdomainHMaxSetup()
Caches the hmax for all elements in each subdomain.

◆ moveRaysToBuffer()

void RayTracingStudy::moveRaysToBuffer ( std::vector< std::shared_ptr< Ray >> &  rays)
protected

Moves rays to the buffer to be traced during generateRays().

This moves the rays into the buffer (with std::move), therefore all valid rays in rays will be nullptr after this call.

Definition at line 1619 of file RayTracingStudy.C.

Referenced by TestPICRayStudy::generateRays(), TestReuseRaysStudy::generateRays(), and TestTransientRaysStudy::generateRays().

1620 {
1621  mooseAssert(currentlyGenerating(), "Can only use while generating");
1622 #ifndef NDEBUG
1623  for (const std::shared_ptr<Ray> & ray : rays)
1624  {
1625  mooseAssert(ray, "Null ray");
1626  mooseAssert(ray->shouldContinue(), "Ray is not continuing");
1627  }
1628 #endif
1629 
1630  _parallel_ray_study->moveWorkToBuffer(rays, /* tid = */ 0);
1631 }
bool currentlyGenerating() const
Whether or not the study is generating.
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.

◆ moveRayToBuffer()

void RayTracingStudy::moveRayToBuffer ( std::shared_ptr< Ray > &  ray)
protected

Moves a ray to the buffer to be traced during generateRays().

This moves the Ray into the buffer (with std::move), therefore ray will be nullptr after this call.

Definition at line 1609 of file RayTracingStudy.C.

Referenced by RayTracingStudyTest::generateRays(), StationaryRayStudyTest::generateRays(), TestReuseRaysStudy::generateRays(), TestTransientRaysStudy::generateRays(), RepeatableRayStudyBase::generateRays(), and ViewFactorRayStudy::generateRays().

1610 {
1611  mooseAssert(currentlyGenerating(), "Can only use while generating");
1612  mooseAssert(ray, "Null ray");
1613  mooseAssert(ray->shouldContinue(), "Ray is not continuing");
1614 
1615  _parallel_ray_study->moveWorkToBuffer(ray, /* tid = */ 0);
1616 }
bool currentlyGenerating() const
Whether or not the study is generating.
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.

◆ moveRayToBufferDuringTrace()

void RayTracingStudy::moveRayToBufferDuringTrace ( std::shared_ptr< Ray > &  ray,
const THREAD_ID  tid,
const AcquireMoveDuringTraceKey  
)

INTERNAL method for moving a Ray into the buffer during tracing.

You should not and cannot use this API. It is protected by the AcquireMoveDuringTraceKey, which is only constructable by RayKernelBase and RayBoundaryConditionBase.

Definition at line 1634 of file RayTracingStudy.C.

Referenced by RayBoundaryConditionBase::moveRayToBuffer(), and RayKernelBase::moveRayToBuffer().

1637 {
1638  mooseAssert(ray, "Null ray");
1639  mooseAssert(currentlyPropagating(), "Can only use while tracing");
1640 
1641  _parallel_ray_study->moveWorkToBuffer(ray, tid);
1642 }
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)

◆ nonPlanarSideSetup()

void RayTracingStudy::nonPlanarSideSetup ( )
private

Sets up the caching of whether or not each element side is non-planar, which is stored in _non_planar_sides.

This is done because this check is made within TraceRay, and determining whether or not a side is non planar is quite expensive and we don't want to call it mid-trace.

Definition at line 573 of file RayTracingStudy.C.

Referenced by initialSetup(), and meshChanged().

574 {
575  _has_non_planar_sides = false;
576  bool warned = !_warn_non_planar;
577 
578  // Nothing to do here for 2D or 1D
579  if (_mesh.dimension() != 3)
580  return;
581 
582  // Clear the data structure and size it based on the elements that we know about
583  _non_planar_sides.clear();
585 
586  for (const Elem * elem : _mesh.getMesh().active_element_ptr_range())
587  {
588  const auto index = _elem_index_helper.getIndex(elem);
589  auto & entry = _non_planar_sides[index];
590  entry.resize(elem->n_sides(), 0);
591 
592  for (const auto s : elem->side_index_range())
593  {
594  const auto & side = elemSide(*elem, s);
595  if (side.n_vertices() < 4)
596  continue;
597 
598  if (!side.has_affine_map())
599  {
600  entry[s] = 1;
601  _has_non_planar_sides = true;
602 
603  if (!warned)
604  {
605  mooseWarning("The mesh contains non-planar faces.\n\n",
606  "Ray tracing on non-planar faces is an approximation and may fail.\n\n",
607  "Use at your own risk! You can disable this warning by setting the\n",
608  "parameter 'warn_non_planar' to false.");
609  warned = true;
610  }
611  }
612  }
613  }
614 }
MooseMesh & _mesh
The Mesh.
libMesh::dof_id_type maxIndex() const
Gets the maximum index generated using this object.
std::vector< std::vector< unsigned short > > _non_planar_sides
Non planar side data, which is for quick checking if an elem side is non-planar We use unsigned short...
MeshBase & getMesh()
virtual unsigned int dimension() const
const bool _warn_non_planar
Whether not to warn if non-planar faces are found.
libMesh::dof_id_type getIndex(const libMesh::Elem *elem) const
Get the index associated with the element elem.
ElemIndexHelper _elem_index_helper
Helper for defining a local contiguous index for each element.
bool _has_non_planar_sides
Whether or not the local mesh has elements with non-planar sides.
const libMesh::Elem & elemSide(const libMesh::Elem &elem, const unsigned int s, const THREAD_ID tid=0)
Get an element&#39;s side pointer without excessive memory allocation.
void mooseWarning(Args &&... args) const

◆ onCompleteRay()

void RayTracingStudy::onCompleteRay ( const std::shared_ptr< Ray > &  ray)
virtual

Entry point for acting on a ray when it is completed (shouldContinue() == false)

Reimplemented in TestRayDataStudy.

Definition at line 1058 of file RayTracingStudy.C.

Referenced by TestRayDataStudy::onCompleteRay(), and ParallelRayStudy::postExecuteChunk().

1059 {
1060  mooseAssert(currentlyPropagating(), "Should only be called during Ray propagation");
1061 
1062  _ending_processor_crossings += ray->processorCrossings();
1064  std::max(_ending_max_processor_crossings, ray->processorCrossings());
1065  _ending_intersections += ray->intersections();
1066  _ending_max_intersections = std::max(_ending_max_intersections, ray->intersections());
1068  std::max(_ending_max_trajectory_changes, ray->trajectoryChanges());
1069  _ending_distance += ray->distance();
1070 
1071 #ifdef NDEBUG
1072  // In non-opt modes, we will always bank the Rays for debugging
1074 #endif
1075  _ray_bank.emplace_back(ray);
1076 }
unsigned long long int _ending_intersections
Total number of Ray/element intersections for Rays that finished on this processor.
Real _ending_distance
Total distance traveled by Rays that end on this processor.
const bool _bank_rays_on_completion
Whether or not to bank rays on completion.
unsigned int _ending_max_intersections
Max number of intersections for Rays that finished on this processor.
unsigned int _ending_max_trajectory_changes
Max number of trajectory changes for Rays that finished on this processor.
unsigned int _ending_max_processor_crossings
Max number of total processor crossings for Rays that finished on this processor. ...
unsigned long long int _ending_processor_crossings
Total number of processor crossings for Rays that finished on this processor.
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)
std::vector< std::shared_ptr< Ray > > _ray_bank
Cumulative Ray bank - stored only when _bank_rays_on_completion.

◆ parallelRayStudy()

const ParallelRayStudy& RayTracingStudy::parallelRayStudy ( ) const
inline

Definition at line 179 of file RayTracingStudy.h.

Referenced by PerProcessorRayTracingResultsVectorPostprocessor::execute(), and RayTracingStudyResult::getValue().

179 { return *_parallel_ray_study; }
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.

◆ parallelStudy()

ParallelStudy<std::shared_ptr<Ray>, Ray>* RayTracingStudy::parallelStudy ( )
inline

The underlying parallel study: used for the context for calling the packed range routines.

Definition at line 678 of file RayTracingStudy.h.

Referenced by RayTracingStudyTest::generateRays(), and RepeatableRayStudyBase::verifyReplicatedRays().

678 { return _parallel_ray_study.get(); }
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.

◆ periodicBoundaryChecks()

void RayTracingStudy::periodicBoundaryChecks ( )
private

Check for overlapping PeriodicRayBC boundaries and check for cases in which ghosting may not be sufficient with distributed mesh for overlapping periodic boundaries.

Definition at line 401 of file RayTracingStudy.C.

Referenced by initialSetup().

402 {
403  // Collect the PeriodicRayBCs
404  std::vector<const RayBoundaryConditionBase *> rbc_ptrs;
405  getRayBCs(rbc_ptrs, 0);
406  std::vector<const PeriodicRayBC *> prbc_ptrs;
407  for (const auto rbc_ptr : rbc_ptrs)
408  if (const auto prbc_ptr = dynamic_cast<const PeriodicRayBC *>(rbc_ptr))
409  prbc_ptrs.push_back(prbc_ptr);
410  if (prbc_ptrs.empty())
411  return;
412 
413  // Collect each of the periodic boundaries
414  std::map<boundary_id_type,
415  std::tuple<const PeriodicRayBC *,
417  std::unordered_set<dof_id_type>>>
418  boundary_map;
419  for (const auto prbc_ptr : prbc_ptrs)
420  {
421  for (const auto & [bid, pb] : prbc_ptr->getPeriodicBoundaries())
422  {
423  const auto [it, inserted] = boundary_map.emplace(
424  std::piecewise_construct,
425  std::tuple{bid},
426  std::forward_as_tuple(prbc_ptr, pb.get(), std::unordered_set<dof_id_type>()));
427  if (!inserted)
428  prbc_ptr->mooseError("The periodic boundary '",
430  "' has been defined in both ",
431  prbc_ptr->typeAndName(),
432  " and ",
433  std::get<0>(it->second)->typeAndName());
434  }
435  }
436 
437  // Because we don't have ghosting setup correctly yet, we need to check if any
438  // of the periodic boundaries are neighbors with distributed mesh. If the
439  // mesh is replicated, we don't need to check this. See #31280.
440  if (comm().size() == 1 || !_mesh.isDistributedMesh())
441  return;
442 
443  // Collect all of the nodes that are on each periodic boundary
444  const auto & sideset_map = _mesh.getMesh().get_boundary_info().get_sideset_map();
445  for (const auto & [elem, side_bid_pair] : sideset_map)
446  {
447  const auto [side, bid] = side_bid_pair;
448  if (auto it = boundary_map.find(bid); it != boundary_map.end())
449  for (const auto n : elem->nodes_on_side(side))
450  std::get<2>(it->second).insert(elem->node_ref(n).id());
451  }
452 
453  // Distributed meshes have distributed boundary information, so sync
454  for (auto & bid_tuple_pair : boundary_map)
455  comm().set_union(std::get<2>(bid_tuple_pair.second));
456 
457  // Check for periodic boundaries that share nodes
458  std::map<std::pair<boundary_id_type, boundary_id_type>,
459  std::pair<const PeriodicRayBC *, const PeriodicRayBC *>>
460  warn_boundaries;
461  for (auto it = boundary_map.begin(); it != boundary_map.end(); ++it)
462  {
463  const auto & [bid, tup] = *it;
464  const auto [prbc_ptr, pb, node_ids] = tup;
465 
466  for (auto other_it = std::next(it); other_it != boundary_map.end(); ++other_it)
467  {
468  const auto & [other_bid, other_tup] = *other_it;
469 
470  // Don't check against boundaries that are paried together
471  if (pb->pairedboundary == other_bid)
472  continue;
473 
474  const auto other_prbc_ptr = std::get<0>(other_tup);
475  const auto & other_node_ids = std::get<2>(other_tup);
476  for (const auto node_id : node_ids)
477  if (other_node_ids.count(node_id))
478  {
479  if (!warn_boundaries.count(std::make_pair(other_bid, bid)))
480  warn_boundaries.emplace(std::make_pair(bid, other_bid),
481  std::make_pair(prbc_ptr, other_prbc_ptr));
482  break;
483  }
484  }
485  }
486 
487  if (warn_boundaries.size())
488  {
489  std::ostringstream oss;
490  oss << warn_boundaries.size()
491  << " ray tracing periodic boundaries were found to be neighbors:\n\n";
492  for (const auto & [bids_pair, prbc_ptrs_pair] : warn_boundaries)
493  {
494  const auto [bid, paired_bid] = bids_pair;
495  const auto [prbc_ptr, paired_prbc_ptr] = prbc_ptrs_pair;
496  oss << " '" << _mesh.getBoundaryString(bid) << "' (in " << prbc_ptr->typeAndName()
497  << ") <-> '" << _mesh.getBoundaryString(paired_bid) << "' (in "
498  << paired_prbc_ptr->typeAndName() << ")\n";
499  }
500  oss << "\nThe periodic propagation of rays at points where two or more periodic"
501  << "\nboundaries meet is not fully supported with a distributed mesh."
502  << "\n\nIf you encounter trace failures, you should use a replicated mesh.";
503  mooseWarning(oss.str());
504  }
505 }
RayBC that enforces periodic boundaries.
Definition: PeriodicRayBC.h:17
MooseMesh & _mesh
The Mesh.
const Parallel::Communicator & comm() const
const BoundaryInfo & get_boundary_info() const
processor_id_type size() const
int8_t boundary_id_type
MeshBase & getMesh()
std::string typeAndName() const
std::string getBoundaryString(const BoundaryID boundary_id) const
const std::multimap< const Elem *, std::pair< unsigned short int, boundary_id_type > > & get_sideset_map() const
void getRayBCs(std::vector< RayBoundaryConditionBase *> &result, BoundaryID id, THREAD_ID tid)
Fills the active RayBCs associated with this study and a boundary into result.
void mooseWarning(Args &&... args) const
virtual bool isDistributedMesh() const
void set_union(T &data, const unsigned int root_id) const

◆ postExecuteStudy()

virtual void RayTracingStudy::postExecuteStudy ( )
inlineprotectedvirtual

Entry point after study execution.

Reimplemented in ViewFactorRayStudy, TestTransientRaysStudy, TestReuseRaysStudy, TestPICRayStudy, RayTracingStudyTest, and TestRayLots.

Definition at line 709 of file RayTracingStudy.h.

Referenced by executeStudy().

709 {}

◆ postOnSegment()

void RayTracingStudy::postOnSegment ( const THREAD_ID  tid,
const std::shared_ptr< Ray > &  ray 
)
virtual

Called at the end of a Ray segment.

Parameters
tidThread id
rayThe ray

Definition at line 842 of file RayTracingStudy.C.

Referenced by TraceRay::trace().

843 {
844  mooseAssert(currentlyPropagating(), "Should not call while not propagating");
846  mooseAssert(_num_cached[tid] == 0,
847  "Values should only be cached when computing Jacobian/residual");
848 
849  // Fill into cached Jacobian/residuals if necessary
851  {
853 
854  if (++_num_cached[tid] == 20)
855  {
856  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
858  _num_cached[tid] = 0;
859  }
860  }
862  {
864 
865  if (++_num_cached[tid] == 20)
866  {
867  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
869  _num_cached[tid] = 0;
870  }
871  }
872 }
FEProblemBase & _fe_problem
Threads::spin_mutex _spin_mutex
Spin mutex object for locks.
const bool & currentlyComputingResidual() const
virtual void cacheResidual(const THREAD_ID tid) override
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)
std::vector< std::size_t > _num_cached
Number of currently cached objects for Jacobian/residual for each thread.
const bool & currentlyComputingJacobian() const
virtual void cacheJacobian(const THREAD_ID tid) override
virtual void addCachedResidual(const THREAD_ID tid) override
virtual void addCachedJacobian(const THREAD_ID tid) override

◆ preExecuteStudy()

virtual void RayTracingStudy::preExecuteStudy ( )
inlineprotectedvirtual

Entry point before study execution.

Reimplemented in ViewFactorRayStudy.

Definition at line 705 of file RayTracingStudy.h.

Referenced by executeStudy().

705 {}

◆ preTrace()

virtual void RayTracingStudy::preTrace ( const THREAD_ID  ,
const std::shared_ptr< Ray > &   
)
inlinevirtual

Called at the beginning of a trace for a ray.

Parameters
tidThread id
rayThe ray

Definition at line 120 of file RayTracingStudy.h.

Referenced by TraceRay::trace().

120 {}

◆ propagationTime()

Real RayTracingStudy::propagationTime ( ) const
inline

Duration for creation of all Rays in seconds.

Definition at line 204 of file RayTracingStudy.h.

Referenced by PerProcessorRayTracingResultsVectorPostprocessor::execute().

204 { return std::chrono::duration<Real>(_propagation_time).count(); }
std::chrono::steady_clock::duration _propagation_time

◆ rayAuxDataNames()

const std::vector<std::string>& RayTracingStudy::rayAuxDataNames ( ) const
inline

The Ray aux data names.

Definition at line 340 of file RayTracingStudy.h.

Referenced by RayTracingMeshOutput::setupEquationSystem().

340 { return _ray_aux_data_names; }
std::vector< std::string > _ray_aux_data_names
The names for each Ray aux data entry.

◆ rayAuxDataSize()

std::size_t RayTracingStudy::rayAuxDataSize ( ) const
inline

The registered size of values in the Ray aux data.

Definition at line 332 of file RayTracingStudy.h.

Referenced by acquireRay(), acquireRayDuringTrace(), acquireRegisteredRay(), acquireReplicatedRay(), Ray::auxData(), and getRayDataNameInternal().

332 { return _ray_aux_data_names.size(); }
std::vector< std::string > _ray_aux_data_names
The names for each Ray aux data entry.

◆ rayBank()

const std::vector< std::shared_ptr< Ray > > & RayTracingStudy::rayBank ( ) const
protected

Get the Ray bank.

This is the bank of Rays that have completed on this processor after an execution of the study.

This is only available when the private parameter _bank_rays_on_completion is set to true.

Definition at line 1369 of file RayTracingStudy.C.

Referenced by RayTracingStudyTest::generateRays(), getBankedRay(), TestRayLots::postExecuteStudy(), RayTracingStudyTest::postExecuteStudy(), TestPICRayStudy::postExecuteStudy(), TestReuseRaysStudy::postExecuteStudy(), TestTransientRaysStudy::postExecuteStudy(), and RayTracingStudyNoBankingTest::RayTracingStudyNoBankingTest().

1370 {
1372  mooseError("The Ray bank is not available because the private parameter "
1373  "'_bank_rays_on_completion' is set to false.");
1375  mooseError("Cannot get the Ray bank during generation or propagation.");
1376 
1377  return _ray_bank;
1378 }
bool currentlyGenerating() const
Whether or not the study is generating.
const bool _bank_rays_on_completion
Whether or not to bank rays on completion.
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)
void mooseError(Args &&... args) const
std::vector< std::shared_ptr< Ray > > _ray_bank
Cumulative Ray bank - stored only when _bank_rays_on_completion.

◆ rayDataNames()

const std::vector<std::string>& RayTracingStudy::rayDataNames ( ) const
inline

The Ray data names.

Definition at line 280 of file RayTracingStudy.h.

Referenced by RayTracingExodus::outputMesh(), RayTracingNemesis::outputMesh(), and RayTracingMeshOutput::setupEquationSystem().

280 { return _ray_data_names; }
std::vector< std::string > _ray_data_names
The names for each Ray data entry.

◆ rayDataSize()

std::size_t RayTracingStudy::rayDataSize ( ) const
inline

The registered size of values in the Ray data.

Definition at line 272 of file RayTracingStudy.h.

Referenced by acquireRay(), acquireRayDuringTrace(), acquireRegisteredRay(), acquireReplicatedRay(), Ray::data(), and getRayDataNameInternal().

272 { return _ray_data_names.size(); }
std::vector< std::string > _ray_data_names
The names for each Ray data entry.

◆ rayDependentSubdomainSetup()

bool RayTracingStudy::rayDependentSubdomainSetup ( ) const
inline

Whether or not to use Ray dependent subdomain setup.

Definition at line 219 of file RayTracingStudy.h.

Referenced by TraceRay::trace().

const bool _ray_dependent_subdomain_setup
Whether or not subdomain setup is dependent on the Ray.

◆ rayMaxDistance()

Real RayTracingStudy::rayMaxDistance ( ) const
inline

Max distance any Ray can travel.

Definition at line 184 of file RayTracingStudy.h.

Referenced by TraceRay::trace().

184 { return _ray_max_distance; }
const Real _ray_max_distance
Max distance a Ray can travel before being killed (can change)

◆ registeredRayID()

RayID RayTracingStudy::registeredRayID ( const std::string &  name,
const bool  graceful = false 
) const

Gets the ID of a registered ray.

Parameters
nameThe name of said ray
gracefulWhether or not to exit gracefully if none is found (with invalid_id)

Definition at line 1455 of file RayTracingStudy.C.

Referenced by RayDataValue::getValue(), RayIntegralValue::initialize(), RayTracingStudyTest::RayTracingStudyTest(), RayTracingStudyWithRegistrationTest::RayTracingStudyWithRegistrationTest(), and registeredRaySetup().

1456 {
1457  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
1458 
1459  if (!_use_ray_registration)
1460  mooseError("Should not use registeredRayID() with Ray registration disabled");
1461 
1462  const auto search = _registered_ray_map.find(name);
1463  if (search != _registered_ray_map.end())
1464  return search->second;
1465 
1466  if (graceful)
1467  return Ray::INVALID_RAY_ID;
1468 
1469  mooseError("Attempted to obtain ID of registered Ray ",
1470  name,
1471  ", but a Ray with said name is not registered.");
1472 }
if(subdm)
Threads::spin_mutex _spin_mutex
Spin mutex object for locks.
const std::string & name() const
void mooseError(Args &&... args) const
std::unordered_map< std::string, RayID > & _registered_ray_map
Map from registered Ray name to ID.
const bool _use_ray_registration
Whether or not to use Ray registration.
static const RayID INVALID_RAY_ID
Invalid Ray ID.
Definition: Ray.h:214

◆ registeredRayName()

const std::string & RayTracingStudy::registeredRayName ( const RayID  ray_id) const

Gets the name of a registered ray.

Parameters
ray_idThe ID of said ray

Definition at line 1475 of file RayTracingStudy.C.

Referenced by Ray::getInfo(), RayTracingStudyTest::RayTracingStudyTest(), and RayTracingStudyWithRegistrationTest::RayTracingStudyWithRegistrationTest().

1476 {
1477  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
1478 
1479  if (!_use_ray_registration)
1480  mooseError("Should not use registeredRayName() with Ray registration disabled");
1481 
1482  if (_reverse_registered_ray_map.size() > ray_id)
1483  return _reverse_registered_ray_map[ray_id];
1484 
1485  mooseError("Attempted to obtain name of registered Ray with ID ",
1486  ray_id,
1487  ", but a Ray with said ID is not registered.");
1488 }
std::vector< std::string > & _reverse_registered_ray_map
Map from registered Ray ID to name.
Threads::spin_mutex _spin_mutex
Spin mutex object for locks.
void mooseError(Args &&... args) const
const bool _use_ray_registration
Whether or not to use Ray registration.

◆ registeredRaySetup()

void RayTracingStudy::registeredRaySetup ( )
private

Sets up the maps from Ray to associated RayTracingObjects if _use_ray_registration.

If ray registration is disabled, this makes sure no RayTracingObjects provided rays.

Definition at line 671 of file RayTracingStudy.C.

Referenced by executeStudy().

672 {
673  // First, clear the objects associated with each Ray on each thread
674  const auto num_rays = _registered_ray_map.size();
675  for (auto & entry : _threaded_ray_object_registration)
676  {
677  entry.clear();
678  entry.resize(num_rays);
679  }
680 
681  const auto rtos = getRayTracingObjects();
682 
684  {
685  // All of the registered ray names - used when a RayTracingObject did not specify
686  // any Rays so it should be associated with all Rays.
687  std::vector<std::string> all_ray_names;
688  all_ray_names.reserve(_registered_ray_map.size());
689  for (const auto & pair : _registered_ray_map)
690  all_ray_names.push_back(pair.first);
691 
692  for (auto & rto : rtos)
693  {
694  // The Ray names associated with this RayTracingObject
695  const auto & ray_names = rto->parameters().get<std::vector<std::string>>("rays");
696  // The registration for RayTracingObjects for the thread rto is on
697  const auto tid = rto->parameters().get<THREAD_ID>("_tid");
698  auto & registration = _threaded_ray_object_registration[tid];
699 
700  // Register each Ray for this object in the registration
701  for (const auto & ray_name : (ray_names.empty() ? all_ray_names : ray_names))
702  {
703  const auto id = registeredRayID(ray_name, /* graceful = */ true);
704  if (ray_names.size() && id == Ray::INVALID_RAY_ID)
705  rto->paramError(
706  "rays", "Supplied ray '", ray_name, "' is not a registered Ray in ", typeAndName());
707  registration[id].insert(rto);
708  }
709  }
710  }
711  // Not using Ray registration
712  else
713  {
714  for (const auto & rto : rtos)
715  if (rto->parameters().get<std::vector<std::string>>("rays").size())
716  rto->paramError(
717  "rays",
718  "Rays cannot be supplied when the study does not require Ray registration.\n\n",
719  type(),
720  " does not require Ray registration.");
721  }
722 }
std::vector< std::vector< std::set< const RayTracingObject * > > > _threaded_ray_object_registration
Threaded storage for all of the RayTracingObjects associated with a single Ray.
std::vector< RayTracingObject * > getRayTracingObjects()
Gets all of the currently active RayTracingObjects.
const std::string & type() const
std::string typeAndName() const
std::unordered_map< std::string, RayID > & _registered_ray_map
Map from registered Ray name to ID.
RayID registeredRayID(const std::string &name, const bool graceful=false) const
Gets the ID of a registered ray.
const bool _use_ray_registration
Whether or not to use Ray registration.
static const RayID INVALID_RAY_ID
Invalid Ray ID.
Definition: Ray.h:214
unsigned int THREAD_ID

◆ registerRay()

RayID RayTracingStudy::registerRay ( const std::string &  name)
private

Registers a Ray with a given name.

This is for internal use only and is an internal map for replicated Ray IDs -> names for easy access for studies that have a small number of rays to generate.

Parameters
nameThe name to register
Returns
The allocated ID for the registered ray

Definition at line 1431 of file RayTracingStudy.C.

Referenced by acquireRegisteredRay().

1432 {
1433  libmesh_parallel_only(comm());
1434 
1435  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
1436 
1437  if (!_use_ray_registration)
1438  mooseError("Cannot use registerRay() with Ray registration disabled");
1439 
1440  // This is parallel only for now. We could likely stagger the ID building like we do with
1441  // the unique IDs, but it would require a sync point which isn't there right now
1442  libmesh_parallel_only(comm());
1443 
1444  const auto & it = _registered_ray_map.find(name);
1445  if (it != _registered_ray_map.end())
1446  return it->second;
1447 
1448  const auto id = _reverse_registered_ray_map.size();
1449  _registered_ray_map.emplace(name, id);
1450  _reverse_registered_ray_map.push_back(name);
1451  return id;
1452 }
std::vector< std::string > & _reverse_registered_ray_map
Map from registered Ray ID to name.
const Parallel::Communicator & comm() const
Threads::spin_mutex _spin_mutex
Spin mutex object for locks.
const std::string & name() const
void mooseError(Args &&... args) const
std::unordered_map< std::string, RayID > & _registered_ray_map
Map from registered Ray name to ID.
const bool _use_ray_registration
Whether or not to use Ray registration.

◆ registerRayAuxData() [1/2]

RayDataIndex RayTracingStudy::registerRayAuxData ( const std::string &  name)

Register a value to be filled in the aux data on a Ray with a given name.

Parameters
nameThe name to register
Returns
The Ray aux data index for the registered value

Note that this does not actually allocate the storage. It is simply a registration system to keep track of the values stored in the Ray aux data.

Definition at line 1202 of file RayTracingStudy.C.

Referenced by RayTracingStudyTest::RayTracingStudyTest(), StationaryRayStudyTest::StationaryRayStudyTest(), TestRayDataStudy::TestRayDataStudy(), and TestRayLots::TestRayLots().

1203 {
1204  return registerRayDataInternal(name, /* aux = */ true);
1205 }
RayDataIndex registerRayDataInternal(const std::string &name, const bool aux)
Internal method for registering Ray data or Ray aux data with a name.
const std::string & name() const

◆ registerRayAuxData() [2/2]

std::vector< RayDataIndex > RayTracingStudy::registerRayAuxData ( const std::vector< std::string > &  names)

Register values to be filled in the aux data on a Ray with a given name.

Parameters
namesThe names to register
Returns
The Ray aux data indices for the registered values

Note that this does not actually allocate the storage. It is simply a registration system to keep track of the values stored in the Ray aux data.

Definition at line 1208 of file RayTracingStudy.C.

1209 {
1210  return registerRayDataInternal(names, /* aux = */ true);
1211 }
RayDataIndex registerRayDataInternal(const std::string &name, const bool aux)
Internal method for registering Ray data or Ray aux data with a name.

◆ registerRayData() [1/2]

RayDataIndex RayTracingStudy::registerRayData ( const std::string &  name)

Register a value to be filled in the data on a Ray with a given name.

Parameters
nameThe name to register
Returns
The Ray data index for the registered value

Note that this does not actually allocate the storage. It is simply a registration system to keep track of the values stored in the Ray data.

Definition at line 1171 of file RayTracingStudy.C.

Referenced by RayTracingStudyTest::generateRays(), RayTracingStudyTest::RayTracingStudyTest(), StationaryRayStudyTest::StationaryRayStudyTest(), TestRayDataStudy::TestRayDataStudy(), and TestRayLots::TestRayLots().

1172 {
1173  return registerRayDataInternal(name, /* aux = */ false);
1174 }
RayDataIndex registerRayDataInternal(const std::string &name, const bool aux)
Internal method for registering Ray data or Ray aux data with a name.
const std::string & name() const

◆ registerRayData() [2/2]

std::vector< RayDataIndex > RayTracingStudy::registerRayData ( const std::vector< std::string > &  names)

Register values to be filled in the data on a Ray with a given name.

Parameters
namesThe names to register
Returns
The Ray data indices for the registered value

Note that this does not actually allocate the storage. It is simply a registration system to keep track of the values stored in the Ray data.

Definition at line 1177 of file RayTracingStudy.C.

1178 {
1179  return registerRayDataInternal(names, /* aux = */ false);
1180 }
RayDataIndex registerRayDataInternal(const std::string &name, const bool aux)
Internal method for registering Ray data or Ray aux data with a name.

◆ registerRayDataInternal() [1/2]

RayDataIndex RayTracingStudy::registerRayDataInternal ( const std::string &  name,
const bool  aux 
)
private

Internal method for registering Ray data or Ray aux data with a name.

Definition at line 1079 of file RayTracingStudy.C.

Referenced by registerRayAuxData(), registerRayData(), and registerRayDataInternal().

1080 {
1081  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
1082 
1084  mooseError("Cannot register Ray ", (aux ? "aux " : ""), "data after initialSetup()");
1085 
1086  auto & map = aux ? _ray_aux_data_map : _ray_data_map;
1087  const auto find = map.find(name);
1088  if (find != map.end())
1089  return find->second;
1090 
1091  auto & other_map = aux ? _ray_data_map : _ray_aux_data_map;
1092  if (other_map.find(name) != other_map.end())
1093  mooseError("Cannot register Ray aux data with name ",
1094  name,
1095  " because Ray ",
1096  (aux ? "(non-aux)" : "aux"),
1097  " data already exists with said name.");
1098 
1099  // Add into the name -> index map
1100  map.emplace(name, map.size());
1101 
1102  // Add into the index -> names vector
1103  auto & vector = aux ? _ray_aux_data_names : _ray_data_names;
1104  vector.push_back(name);
1105 
1106  return map.size() - 1;
1107 }
KOKKOS_INLINE_FUNCTION const T * find(const T &target, const T *const begin, const T *const end)
if(subdm)
std::vector< std::string > _ray_aux_data_names
The names for each Ray aux data entry.
Threads::spin_mutex _spin_mutex
Spin mutex object for locks.
const std::string & name() const
std::unordered_map< std::string, RayDataIndex > _ray_aux_data_map
The map from Ray aux data names to index.
std::unordered_map< std::string, RayDataIndex > _ray_data_map
The map from Ray data names to index.
bool _called_initial_setup
Whether or not we&#39;ve called initial setup - used to stop from late registration.
void mooseError(Args &&... args) const
std::vector< std::string > _ray_data_names
The names for each Ray data entry.

◆ registerRayDataInternal() [2/2]

std::vector< RayDataIndex > RayTracingStudy::registerRayDataInternal ( const std::vector< std::string > &  names,
const bool  aux 
)
private

Internal method for registering Ray data or Ray aux data with names.

Definition at line 1110 of file RayTracingStudy.C.

1111 {
1112  std::vector<RayDataIndex> indices(names.size());
1113  for (std::size_t i = 0; i < names.size(); ++i)
1114  indices[i] = registerRayDataInternal(names[i], aux);
1115  return indices;
1116 }
RayDataIndex registerRayDataInternal(const std::string &name, const bool aux)
Internal method for registering Ray data or Ray aux data with a name.

◆ reinitSegment()

void RayTracingStudy::reinitSegment ( const Elem elem,
const Point start,
const Point end,
const Real  length,
THREAD_ID  tid 
)
virtual

Reinitialize objects for a Ray segment for ray tracing.

Parameters
elemThe elem the segment is in
startStart point of the segment
endEnd point of the segment
lengthThe length of the start -> end segment
tidThread id

Definition at line 778 of file RayTracingStudy.C.

Referenced by TraceRay::onSegment().

780 {
781  mooseAssert(MooseUtils::absoluteFuzzyEqual((start - end).norm(), length), "Invalid length");
782  mooseAssert(currentlyPropagating(), "Should not call while not propagating");
783 
785 
786  // If we have any variables or material properties that are active, we definitely need to reinit
787  bool reinit = _fe_problem.hasActiveElementalMooseVariables(tid) ||
789  // If not, make sure that the RayKernels have not requested a reinit (this could happen when a
790  // RayKernel doesn't have variables or materials but still does an integration and needs qps)
791  if (!reinit)
792  for (const RayKernelBase * rk : currentRayKernels(tid))
793  if (rk->needSegmentReinit())
794  {
795  reinit = true;
796  break;
797  }
798 
799  if (reinit)
800  {
801  _fe_problem.prepare(elem, tid);
802 
803  std::vector<Point> points;
804  std::vector<Real> weights;
805  buildSegmentQuadrature(start, end, length, points, weights);
806  _fe_problem.reinitElemPhys(elem, points, tid);
807  _fe_problem.assembly(tid, _sys.number()).modifyArbitraryWeights(weights);
808 
810  }
811 }
virtual void prepare(const Elem *elem, const THREAD_ID tid) override
Base object for the RayKernel syntax.
Definition: RayKernelBase.h:27
const std::vector< RayKernelBase * > & currentRayKernels(THREAD_ID tid) const
Gets the current RayKernels for a thread, which are set in segmentSubdomainSetup() ...
virtual Assembly & assembly(const THREAD_ID tid, const unsigned int sys_num) override
FEProblemBase & _fe_problem
unsigned int number() const
virtual void setCurrentSubdomainID(const Elem *elem, const THREAD_ID tid) override
SystemBase & _sys
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)
subdomain_id_type subdomain_id() const
void reinitMaterials(SubdomainID blk_id, const THREAD_ID tid, bool swap_stateful=true)
auto norm(const T &a)
virtual void buildSegmentQuadrature(const Point &start, const Point &end, const Real length, std::vector< Point > &points, std::vector< Real > &weights) const
Builds quadrature points for a given segment using the _segment_qrule.
virtual bool hasActiveElementalMooseVariables(const THREAD_ID tid) const
bool hasActiveMaterialProperties(const THREAD_ID tid) const
virtual void reinitElemPhys(const Elem *elem, const std::vector< Point > &phys_points_in_elem, const THREAD_ID tid) override

◆ reserveRayBuffer()

void RayTracingStudy::reserveRayBuffer ( const std::size_t  size)
protected

Reserve size entires in the Ray buffer.

This can only be used within generateRays() and should be used when possible to reserve entires in the buffer before adding Rays via moveRay(s)ToBuffer.

Definition at line 1645 of file RayTracingStudy.C.

Referenced by RepeatableRayStudyBase::generateRays(), ViewFactorRayStudy::generateRays(), and RayTracingStudyTest::RayTracingStudyTest().

1646 {
1647  if (!currentlyGenerating())
1648  mooseError("Can only reserve in Ray buffer during generateRays()");
1649 
1650  _parallel_ray_study->reserveBuffer(size);
1651 }
bool currentlyGenerating() const
Whether or not the study is generating.
const std::unique_ptr< ParallelRayStudy > _parallel_ray_study
The study that used is to actually execute (trace) the Rays.
void mooseError(Args &&... args) const

◆ resetReplicatedRayIDs()

void RayTracingStudy::resetReplicatedRayIDs ( )
protected

Resets the generation of unique replicated RayIDs accessed via generateReplicatedRayID().

Definition at line 1741 of file RayTracingStudy.C.

Referenced by RayTracingStudy().

1742 {
1743  libmesh_parallel_only(comm());
1744 
1745  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
1746 
1747  mooseAssert(!currentlyGenerating() && !currentlyPropagating(),
1748  "Cannot be reset during generation or propagation");
1749 
1751 }
bool currentlyGenerating() const
Whether or not the study is generating.
const Parallel::Communicator & comm() const
RayID _replicated_next_ray_id
Storage for the next available replicated RayID, obtained via generateReplicatedRayID() ...
Threads::spin_mutex _spin_mutex
Spin mutex object for locks.
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)

◆ resetUniqueRayIDs()

void RayTracingStudy::resetUniqueRayIDs ( )
protected

Resets the generation of unique RayIDs via generateUniqueRayID() to the beginning of the range.

Definition at line 1715 of file RayTracingStudy.C.

Referenced by RayTracingStudy().

1716 {
1717  libmesh_parallel_only(comm());
1718 
1719  Threads::spin_mutex::scoped_lock lock(_spin_mutex);
1720 
1721  mooseAssert(!currentlyGenerating() && !currentlyPropagating(),
1722  "Cannot be reset during generation or propagation");
1723 
1724  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
1726 }
bool currentlyGenerating() const
Whether or not the study is generating.
unsigned long int RayID
Type for a Ray&#39;s ID.
Definition: Ray.h:44
unsigned int n_threads()
const Parallel::Communicator & comm() const
std::vector< RayID > _threaded_next_ray_id
Storage for the next available unique RayID, obtained via generateUniqueRayID()
Threads::spin_mutex _spin_mutex
Spin mutex object for locks.
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)
const processor_id_type _pid
The rank of this processor (this actually takes time to lookup - so just do it once) ...
unsigned int THREAD_ID

◆ residualSetup()

void RayTracingStudy::residualSetup ( )
overridevirtual

Reimplemented from GeneralUserObject.

Definition at line 272 of file RayTracingStudy.C.

273 {
274  for (THREAD_ID tid = 0; tid < libMesh::n_threads(); ++tid)
275  mooseAssert(_num_cached[tid] == 0, "Cached residuals/Jacobians not empty");
276 
277  for (auto & rto : getRayTracingObjects())
278  rto->residualSetup();
279 }
unsigned int n_threads()
std::vector< RayTracingObject * > getRayTracingObjects()
Gets all of the currently active RayTracingObjects.
std::vector< std::size_t > _num_cached
Number of currently cached objects for Jacobian/residual for each thread.
unsigned int THREAD_ID

◆ sameLevelActiveElems()

bool RayTracingStudy::sameLevelActiveElems ( ) const
protected

Determine whether or not the mesh currently has active elements that are all the same level.

Definition at line 1677 of file RayTracingStudy.C.

Referenced by meshChanged().

1678 {
1679  unsigned int min_level = std::numeric_limits<unsigned int>::max();
1680  unsigned int max_level = std::numeric_limits<unsigned int>::min();
1681 
1682  for (const auto & elem : *_mesh.getActiveLocalElementRange())
1683  {
1684  const auto level = elem->level();
1685  min_level = std::min(level, min_level);
1686  max_level = std::max(level, max_level);
1687  }
1688 
1689  _communicator.min(min_level);
1690  _communicator.max(max_level);
1691 
1692  return min_level == max_level;
1693 }
libMesh::ConstElemRange * getActiveLocalElementRange()
MooseMesh & _mesh
The Mesh.
const Parallel::Communicator & _communicator
void min(const T &r, T &o, Request &req) const
void max(const T &r, T &o, Request &req) const

◆ segmentsOnCacheTraces()

bool RayTracingStudy::segmentsOnCacheTraces ( ) const
inline

Whether or not to cache individual element segments when _cache_traces = true.

Definition at line 577 of file RayTracingStudy.h.

Referenced by RayTracingMeshOutput::buildSegmentMesh(), RayTracingMeshOutput::neededNodes(), TraceRay::onContinueTrace(), TraceRay::onTrajectoryChanged(), RayTracingMeshOutput::RayTracingMeshOutput(), and RayTracingMeshOutput::startingIDs().

577 { return _segments_on_cache_traces; }
const bool _segments_on_cache_traces
Whether or not to cache individual element segments when caching.

◆ segmentSubdomainSetup()

void RayTracingStudy::segmentSubdomainSetup ( const SubdomainID  subdomain,
const THREAD_ID  tid,
const RayID  ray_id 
)
virtual

Setup for on subdomain change or subdomain AND ray change during ray tracing.

Parameters
subdomainThe subdomain changed to
tidThread id
ray_idID of the ray initiating the change

Definition at line 743 of file RayTracingStudy.C.

Referenced by TraceRay::onSubdomainChanged().

746 {
747  mooseAssert(currentlyPropagating(), "Should not call while not propagating");
748 
749  // Call subdomain setup on FE
750  _fe_problem.subdomainSetup(subdomain, tid);
751 
752  std::set<MooseVariableFEBase *> needed_moose_vars;
753  std::unordered_set<unsigned int> needed_mat_props;
754 
755  // Get RayKernels and their dependencies and call subdomain setup
756  getRayKernels(_threaded_current_ray_kernels[tid], subdomain, tid, ray_id);
757  for (auto & rkb : _threaded_current_ray_kernels[tid])
758  {
759  rkb->subdomainSetup();
760 
761  const auto & mv_deps = rkb->getMooseVariableDependencies();
762  needed_moose_vars.insert(mv_deps.begin(), mv_deps.end());
763 
764  const auto & mp_deps = rkb->getMatPropDependencies();
765  needed_mat_props.insert(mp_deps.begin(), mp_deps.end());
766  }
767 
768  // Prepare aux vars
769  for (auto & var : needed_moose_vars)
770  if (var->kind() == Moose::VarKindType::VAR_AUXILIARY)
771  var->prepareAux();
772 
773  _fe_problem.setActiveElementalMooseVariables(needed_moose_vars, tid);
774  _fe_problem.prepareMaterials(needed_mat_props, subdomain, tid);
775 }
void getRayKernels(std::vector< RayKernelBase *> &result, SubdomainID id, THREAD_ID tid)
Fills the active RayKernels associated with this study and a block into result.
std::vector< std::vector< RayKernelBase * > > _threaded_current_ray_kernels
The current RayKernel objects for each thread.
FEProblemBase & _fe_problem
virtual void setActiveElementalMooseVariables(const std::set< MooseVariableFEBase * > &moose_vars, const THREAD_ID tid) override
bool currentlyPropagating() const
Whether or not the study is propagating (tracing Rays)
virtual void subdomainSetup(SubdomainID subdomain, const THREAD_ID tid)
void prepareMaterials(const std::unordered_set< unsigned int > &consumer_needed_mat_props, const SubdomainID blk_id, const THREAD_ID tid)

◆ shouldCacheTrace()

virtual bool RayTracingStudy::shouldCacheTrace ( const std::shared_ptr< Ray > &  ) const
inlinevirtual

Virtual that allows for selection in if a Ray should be cached or not (only used when _cache_traces).

Definition at line 583 of file RayTracingStudy.h.

Referenced by initThreadedCachedTrace(), and TraceRay::trace().

584  {
585  return _always_cache_traces;
586  }
bool _always_cache_traces
Whether or not to cache traces on every trace execution.

◆ sideIsIncoming()

bool RayTracingStudy::sideIsIncoming ( const Elem *const  elem,
const unsigned short  side,
const Point direction,
const THREAD_ID  tid 
)

Whether or not side is incoming on element elem in direction direction.

Definition at line 1760 of file RayTracingStudy.C.

Referenced by TraceRay::findExternalBoundarySide(), TraceRay::onSegment(), ClaimRays::postClaimRay(), and TraceRay::trace().

1764 {
1765  const auto & normal = getSideNormal(elem, side, tid);
1766  const auto dot = normal * direction;
1767  return dot < TraceRayTools::TRACE_TOLERANCE;
1768 }
const Real TRACE_TOLERANCE
The standard tolerance to use in tracing.
Definition: TraceRayTools.h:48
virtual const Point & getSideNormal(const Elem *elem, const unsigned short side, const THREAD_ID tid)
Get the outward normal for a given element side.

◆ sideIsNonPlanar()

bool RayTracingStudy::sideIsNonPlanar ( const Elem *  elem,
const unsigned short  s 
) const
inline

Whether or not the side on elem elem is non-planar.

This is cached because checking whether or not a face is planar is costly

Definition at line 465 of file RayTracingStudy.h.

Referenced by TraceRay::exitsElem(), ViewFactorRayStudy::generateRays(), ReflectRayBC::onBoundary(), TraceRay::onSegment(), TraceRay::possiblyAddToBoundaryElems(), Ray::setStart(), and TraceRay::trace().

466  {
468  }
std::vector< std::vector< unsigned short > > _non_planar_sides
Non planar side data, which is for quick checking if an elem side is non-planar We use unsigned short...
libMesh::dof_id_type getIndex(const libMesh::Elem *elem) const
Get the index associated with the element elem.
ElemIndexHelper _elem_index_helper
Helper for defining a local contiguous index for each element.
bool _has_non_planar_sides
Whether or not the local mesh has elements with non-planar sides.

◆ subdomainHmax()

Real RayTracingStudy::subdomainHmax ( const SubdomainID  subdomain_id) const

Get the cached hmax for all elements in a subdomain.

Used for scaling tolerances in ray tracing.

Definition at line 1696 of file RayTracingStudy.C.

Referenced by RayTracingStudyTest::generateRays(), TraceRay::onSubdomainChanged(), TraceRay::subdomainHmax(), and subdomainHMaxSetup().

1697 {
1698  const auto find = _subdomain_hmax.find(subdomain_id);
1699  if (find == _subdomain_hmax.end())
1700  mooseError("Subdomain ", subdomain_id, " not found in subdomain hmax map");
1701  return find->second;
1702 }
KOKKOS_INLINE_FUNCTION const T * find(const T &target, const T *const begin, const T *const end)
std::unordered_map< SubdomainID, Real > _subdomain_hmax
The cached hmax for all elements in a subdomain.
void mooseError(Args &&... args) const

◆ subdomainHMaxSetup()

void RayTracingStudy::subdomainHMaxSetup ( )
private

Caches the hmax for all elements in each subdomain.

Definition at line 617 of file RayTracingStudy.C.

Referenced by initialSetup(), and meshChanged().

618 {
619  // Setup map with subdomain keys
620  _subdomain_hmax.clear();
621  for (const auto subdomain_id : _mesh.meshSubdomains())
622  _subdomain_hmax[subdomain_id] = std::numeric_limits<Real>::min();
623 
624  // Set local max for each subdomain
625  for (const auto & elem : *_mesh.getActiveLocalElementRange())
626  {
627  auto & entry = _subdomain_hmax.at(elem->subdomain_id());
628  entry = std::max(entry, elem->hmax());
629  }
630 
631  // Accumulate global max for each subdomain
633 
634  if (getParam<bool>("warn_subdomain_hmax"))
635  {
636  const auto warn_prefix = type() + " '" + name() + "': ";
637  const auto warn_suffix =
638  "\n\nRay tracing uses an approximate element size for each subdomain to scale the\n"
639  "tolerances used in computing ray intersections. This warning suggests that the\n"
640  "approximate element size is not a good approximation. This is likely due to poor\n"
641  "element aspect ratios.\n\n"
642  "This warning is only output for the first element affected.\n"
643  "To disable this warning, set warn_subdomain_hmax = false.\n";
644 
645  for (const auto & elem : *_mesh.getActiveLocalElementRange())
646  {
647  const auto hmin = elem->hmin();
648  const auto hmax = elem->hmax();
649  const auto max_hmax = subdomainHmax(elem->subdomain_id());
650 
651  const auto hmax_rel = hmax / max_hmax;
652  if (hmax_rel < 1.e-2 || hmax_rel > 1.e2)
653  mooseDoOnce(mooseWarning(warn_prefix,
654  "Element hmax varies significantly from subdomain hmax.\n",
655  warn_suffix,
656  "First element affected:\n",
657  Moose::stringify(*elem)););
658 
659  const auto h_rel = max_hmax / hmin;
660  if (h_rel > 1.e2)
661  mooseDoOnce(mooseWarning(warn_prefix,
662  "Element hmin varies significantly from subdomain hmax.\n",
663  warn_suffix,
664  "First element affected:\n",
665  Moose::stringify(*elem)););
666  }
667  }
668 }
libMesh::ConstElemRange * getActiveLocalElementRange()
MooseMesh & _mesh
The Mesh.
std::unordered_map< SubdomainID, Real > _subdomain_hmax
The cached hmax for all elements in a subdomain.
const Parallel::Communicator & _communicator
Real subdomainHmax(const SubdomainID subdomain_id) const
Get the cached hmax for all elements in a subdomain.
const std::string & name() const
const std::string & type() const
std::string stringify(const T &t)
void max(const T &r, T &o, Request &req) const
void mooseWarning(Args &&... args) const
const std::set< SubdomainID > & meshSubdomains() const

◆ timestepSetup()

void RayTracingStudy::timestepSetup ( )
overridevirtual

Reimplemented from GeneralUserObject.

Definition at line 292 of file RayTracingStudy.C.

293 {
294  for (auto & rto : getRayTracingObjects())
295  rto->timestepSetup();
296 }
std::vector< RayTracingObject * > getRayTracingObjects()
Gets all of the currently active RayTracingObjects.

◆ tolerateFailure()

bool RayTracingStudy::tolerateFailure ( ) const
inline

Whether or not to tolerate failure.

Definition at line 209 of file RayTracingStudy.h.

Referenced by TraceRay::exitsElem(), TraceRay::onBoundary(), and TraceRay::trace().

209 { return _tolerate_failure; }
const bool _tolerate_failure
Whether or not to tolerate a Ray Tracing failure.

◆ totalDistance()

Real RayTracingStudy::totalDistance ( ) const
inline

Total distance traveled by all Rays.

Definition at line 172 of file RayTracingStudy.h.

Referenced by RayTracingStudyResult::getValue().

172 { return _total_distance; }
Real _total_distance
Total distance traveled by all Rays.

◆ totalIntersections()

unsigned long long int RayTracingStudy::totalIntersections ( ) const
inline

Total number of Ray/element intersections.

Definition at line 155 of file RayTracingStudy.h.

155 { return _total_intersections; }
unsigned long long int _total_intersections
Total number of Ray/element intersections.

◆ totalProcessorCrossings()

unsigned long long int RayTracingStudy::totalProcessorCrossings ( ) const
inline

Total number of processor crossings.

Definition at line 138 of file RayTracingStudy.h.

Referenced by RayTracingStudyResult::getValue().

138 { return _total_processor_crossings; }
unsigned long long int _total_processor_crossings
Total number of processor crossings.

◆ totalVolume()

Real RayTracingStudy::totalVolume ( ) const
inline

Get the current total volume of the domain.

Definition at line 437 of file RayTracingStudy.h.

Referenced by isRectangularDomain().

437 { return _total_volume; }
Real _total_volume
The total volume of the domain.

◆ traceableMeshChecks()

void RayTracingStudy::traceableMeshChecks ( )
private

Check for if all of the element types in the mesh are supported by ray tracing.

Definition at line 382 of file RayTracingStudy.C.

Referenced by initialSetup().

383 {
384  for (const auto & elem : *_mesh.getActiveLocalElementRange())
385  {
386  if (_fe_problem.adaptivity().isOn())
387  {
389  mooseError("Element type ",
390  Utility::enum_to_string(elem->type()),
391  " is not supported in ray tracing with adaptivity");
392  }
393  else if (!TraceRayTools::isTraceableElem(elem))
394  mooseError("Element type ",
395  Utility::enum_to_string(elem->type()),
396  " is not supported in ray tracing");
397  }
398 }
libMesh::ConstElemRange * getActiveLocalElementRange()
MooseMesh & _mesh
The Mesh.
bool isTraceableElem(const Elem *elem)
bool isAdaptivityTraceableElem(const Elem *elem)
bool isOn()
FEProblemBase & _fe_problem
void mooseError(Args &&... args) const
Adaptivity & adaptivity()

◆ traceRay() [1/2]

TraceRay& RayTracingStudy::traceRay ( const THREAD_ID  tid)
inline

Gets the threaded TraceRay object for tid.

Definition at line 647 of file RayTracingStudy.h.

Referenced by BackfaceCullingStudyTest::BackfaceCullingStudyTest().

647 { return *_threaded_trace_ray[tid]; }
std::vector< std::shared_ptr< TraceRay > > _threaded_trace_ray
The TraceRay objects for each thread (they do the physical tracing)

◆ traceRay() [2/2]

const TraceRay& RayTracingStudy::traceRay ( const THREAD_ID  tid) const
inline

Definition at line 648 of file RayTracingStudy.h.

648 { return *_threaded_trace_ray[tid]; }
std::vector< std::shared_ptr< TraceRay > > _threaded_trace_ray
The TraceRay objects for each thread (they do the physical tracing)

◆ useRayRegistration()

bool RayTracingStudy::useRayRegistration ( ) const
inline

Whether or not ray registration is being used.

Definition at line 564 of file RayTracingStudy.h.

Referenced by Ray::getInfo(), RayIntegralValue::initialize(), and RayDataValue::RayDataValue().

564 { return _use_ray_registration; }
const bool _use_ray_registration
Whether or not to use Ray registration.

◆ validParams()

InputParameters RayTracingStudy::validParams ( )
static

Whether or not subdomain setup is dependent on the Ray

Definition at line 33 of file RayTracingStudy.C.

Referenced by StationaryRayStudyTest::validParams(), TestRay::validParams(), RayTracingStudyTest::validParams(), TestPICRayStudy::validParams(), ViewFactorRayStudy::validParams(), TestReuseRaysStudy::validParams(), TestTransientRaysStudy::validParams(), RepeatableRayStudyBase::validParams(), RayTracingStudyNoBankingTest::validParams(), and RayTracingStudyWithRegistrationTest::validParams().

34 {
35  auto params = GeneralUserObject::validParams();
36 
37  // Parameters for the execution of the rays
39 
40  params.addRangeCheckedParam<Real>("ray_distance",
41  std::numeric_limits<Real>::max(),
42  "ray_distance > 0",
43  "The maximum distance all Rays can travel");
44 
45  params.addParam<bool>(
46  "tolerate_failure", false, "Whether or not to tolerate a ray tracing failure");
47 
48  MooseEnum work_buffers("lifo circular", "circular");
49  params.addParam<MooseEnum>("work_buffer_type", work_buffers, "The work buffer type to use");
50 
51  params.addParam<bool>(
52  "ray_kernel_coverage_check", true, "Whether or not to perform coverage checks on RayKernels");
53  params.addParam<bool>("warn_non_planar",
54  true,
55  "Whether or not to produce a warning if any element faces are non-planar.");
56 
57  params.addParam<bool>(
58  "always_cache_traces",
59  false,
60  "Whether or not to cache the Ray traces on every execution, primarily for use in output. "
61  "Warning: this can get expensive very quick with a large number of rays!");
62  params.addParam<bool>("data_on_cache_traces",
63  false,
64  "Whether or not to also cache the Ray's data when caching its traces");
65  params.addParam<bool>("aux_data_on_cache_traces",
66  false,
67  "Whether or not to also cache the Ray's aux data when caching its traces");
68  params.addParam<bool>(
69  "segments_on_cache_traces",
70  true,
71  "Whether or not to cache individual segments when trace caching is enabled. If false, we "
72  "will instead cache a segment for each part of the trace where the direction is the same. "
73  "This minimizes the number of segments requied to represent the Ray's path, but removes the "
74  "ability to show Ray field data on each segment through an element.");
75 
76  params.addParam<bool>("use_internal_sidesets",
77  false,
78  "Whether or not to use internal sidesets for RayBCs in ray tracing");
79 
80  params.addParam<bool>("warn_subdomain_hmax",
81  true,
82  "Whether or not to warn if the approximated hmax (constant on subdomain) "
83  "varies significantly for an element");
84 
85  params.addParam<bool>(
86  "verify_rays",
87  true,
88  "Whether or not to verify the generated Rays. This includes checking their "
89  "starting information and the uniqueness of Rays before and after execution. This is also "
90  "used by derived studies for more specific verification.");
91  params.addParam<bool>("verify_trace_intersections",
92  true,
93  "Whether or not to verify the trace intersections in devel and dbg modes. "
94  "Trace intersections are not verified regardless of this parameter in "
95  "optimized modes (opt, oprof).");
96 
97  params.addParam<bool>("allow_other_flags_with_prekernels",
98  false,
99  "Whether or not to allow the list of execution flags to have PRE_KERNELS "
100  "mixed with other flags. If this parameter is not set then if PRE_KERNELS "
101  "is provided it must be the only execution flag.");
102 
103  ExecFlagEnum & exec_enum = params.set<ExecFlagEnum>("execute_on", true);
104  exec_enum.addAvailableFlags(EXEC_PRE_KERNELS);
105 
106  params.addParamNamesToGroup(
107  "always_cache_traces data_on_cache_traces aux_data_on_cache_traces segments_on_cache_traces",
108  "Trace cache");
109  params.addParamNamesToGroup("warn_non_planar warn_subdomain_hmax", "Tracing Warnings");
110  params.addParamNamesToGroup("ray_kernel_coverage_check verify_rays verify_trace_intersections",
111  "Checks and verifications");
112 
113  // Whether or not each Ray must be registered using the registerRay() API
114  params.addPrivateParam<bool>("_use_ray_registration", true);
115  // Whether or not to bank Rays on completion
116  params.addPrivateParam<bool>("_bank_rays_on_completion", true);
118  params.addPrivateParam<bool>("_ray_dependent_subdomain_setup", true);
119 
120  // Add a point neighbor relationship manager
121  params.addRelationshipManager("ElementPointNeighborLayers",
124  [](const InputParameters &, InputParameters & rm_params)
125  { rm_params.set<unsigned short>("layers") = 1; });
126 
127  return params;
128 }
static InputParameters validParams()
T & set(const std::string &name, bool quiet_mode=false)
void addAvailableFlags(const ExecFlagType &flag, Args... flags)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ verifyDependenciesExist()

void RayTracingStudy::verifyDependenciesExist ( const std::vector< RayTracingObject *> &  rtos)
private

Verifies that the dependencies exist for a set of RayTracingObjects.

Definition at line 363 of file RayTracingStudy.C.

Referenced by dependencyChecks().

364 {
365  for (const auto & rto : rtos)
366  for (const auto & dep_name : rto->getRequestedItems())
367  {
368  bool found = false;
369  for (const auto & rto_search : rtos)
370  if (rto_search->name() == dep_name)
371  {
372  found = true;
373  break;
374  }
375 
376  if (!found)
377  rto->paramError("depends_on", "The ", rto->getBase(), " '", dep_name, "' does not exist");
378  }
379 }

◆ verifyRays()

bool RayTracingStudy::verifyRays ( ) const
inline

Whether or not to verify if Rays have valid information before being traced.

Definition at line 654 of file RayTracingStudy.h.

Referenced by RayKernelBase::changeRayStartDirection(), ClaimRays::claim(), RepeatableRayStudyBase::defineRaysInternal(), executeStudy(), Ray::setStart(), Ray::setStartingEndPoint(), and TraceRay::trace().

654 { return _verify_rays; }
const bool _verify_rays
Whether or not to verify if Rays have valid information before being traced.

◆ verifyTraceIntersections()

bool RayTracingStudy::verifyTraceIntersections ( ) const
inline

Whether or not trace verification is enabled in devel/dbg modes.

Definition at line 659 of file RayTracingStudy.h.

Referenced by TraceRay::exitsElem(), TraceRay::onSegment(), TraceRay::onTrajectoryChanged(), and TraceRay::trace().

659 { return _verify_trace_intersections; }
const bool _verify_trace_intersections
Whether or not to verify the trace intersections in devel and dbg modes.

◆ verifyUniqueRayIDs()

void RayTracingStudy::verifyUniqueRayIDs ( const std::vector< std::shared_ptr< Ray >>::const_iterator  begin,
const std::vector< std::shared_ptr< Ray >>::const_iterator  end,
const bool  global,
const std::string &  error_suffix 
) const

Verifies that the Rays in the given range have unique Ray IDs.

Parameters
beginThe beginning of the range of Rays to check
endThe end of the range of Rays to check
globalIf true, this will complete the verification across all processors
error_suffixEntry point for additional information in the error message

Definition at line 1523 of file RayTracingStudy.C.

Referenced by executeStudy(), and RepeatableRayStudyBase::verifyReplicatedRays().

1527 {
1528  // Determine the unique set of Ray IDs on this processor,
1529  // and if not locally unique throw an error. Once we build this set,
1530  // we will send it to rank 0 to verify globally
1531  std::set<RayID> local_rays;
1532  for (const std::shared_ptr<Ray> & ray : as_range(begin, end))
1533  {
1534  mooseAssert(ray, "Null ray");
1535 
1536  // Try to insert into the set; the second entry in the pair
1537  // will be false if it was not inserted
1538  if (!local_rays.insert(ray->id()).second)
1539  {
1540  for (const std::shared_ptr<Ray> & other_ray : as_range(begin, end))
1541  if (ray.get() != other_ray.get() && ray->id() == other_ray->id())
1542  mooseError("Multiple Rays exist with ID ",
1543  ray->id(),
1544  " on processor ",
1545  _pid,
1546  " ",
1547  error_suffix,
1548  "\n\nOffending Ray information:\n\n",
1549  ray->getInfo(),
1550  "\n",
1551  other_ray->getInfo());
1552  }
1553  }
1554 
1555  // Send IDs from all procs to rank 0 and verify on rank 0
1556  if (global)
1557  {
1558  // Package our local IDs and send to rank 0
1559  std::map<processor_id_type, std::vector<RayID>> send_ids;
1560  if (local_rays.size())
1561  send_ids.emplace(std::piecewise_construct,
1562  std::forward_as_tuple(0),
1563  std::forward_as_tuple(local_rays.begin(), local_rays.end()));
1564  local_rays.clear();
1565 
1566  // Mapping on rank 0 from ID -> processor ID
1567  std::map<RayID, processor_id_type> global_map;
1568 
1569  // Verify another processor's IDs against the global map on rank 0
1570  const auto check_ids =
1571  [this, &global_map, &error_suffix](processor_id_type pid, const std::vector<RayID> & ids)
1572  {
1573  for (const RayID id : ids)
1574  {
1575  const auto emplace_pair = global_map.emplace(id, pid);
1576 
1577  // Means that this ID already exists in the map
1578  if (!emplace_pair.second)
1579  mooseError("Ray with ID ",
1580  id,
1581  " exists on ranks ",
1582  emplace_pair.first->second,
1583  " and ",
1584  pid,
1585  "\n",
1586  error_suffix);
1587  }
1588  };
1589 
1590  Parallel::push_parallel_vector_data(_communicator, send_ids, check_ids);
1591  }
1592 }
unsigned long int RayID
Type for a Ray&#39;s ID.
Definition: Ray.h:44
const Parallel::Communicator & _communicator
uint8_t processor_id_type
SimpleRange< IndexType > as_range(const std::pair< IndexType, IndexType > &p)
void mooseError(Args &&... args) const
const processor_id_type _pid
The rank of this processor (this actually takes time to lookup - so just do it once) ...

◆ verifyUniqueRays()

void RayTracingStudy::verifyUniqueRays ( const std::vector< std::shared_ptr< Ray >>::const_iterator  begin,
const std::vector< std::shared_ptr< Ray >>::const_iterator  end,
const std::string &  error_suffix 
)

Verifies that the Rays in the given range are unique.

That is, that there are not multiple shared_ptrs that point to the same Ray

Parameters
beginThe beginning of the range of Rays to check
endThe end of the range of Rays to check
error_suffixEntry point for additional information in the error message

Definition at line 1595 of file RayTracingStudy.C.

Referenced by executeStudy().

1598 {
1599  std::set<const Ray *> rays;
1600  for (const std::shared_ptr<Ray> & ray : as_range(begin, end))
1601  if (!rays.insert(ray.get()).second) // false if not inserted into rays
1602  mooseError("Multiple shared_ptrs were found that point to the same Ray ",
1603  error_suffix,
1604  "\n\nOffending Ray:\n",
1605  ray->getInfo());
1606 }
SimpleRange< IndexType > as_range(const std::pair< IndexType, IndexType > &p)
void mooseError(Args &&... args) const

◆ warnNonPlanar()

bool RayTracingStudy::warnNonPlanar ( ) const
inline

Whether or not to produce a warning when interacting with a non-planar mesh.

Definition at line 673 of file RayTracingStudy.h.

673 { return _warn_non_planar; }
const bool _warn_non_planar
Whether not to warn if non-planar faces are found.

◆ zeroAuxVariables()

void RayTracingStudy::zeroAuxVariables ( )
private

Zero the AuxVariables that the registered AuxRayKernels contribute to.

Definition at line 725 of file RayTracingStudy.C.

Referenced by executeStudy().

726 {
727  std::set<std::string> vars_to_be_zeroed;
728  std::vector<RayKernelBase *> ray_kernels;
729  getRayKernels(ray_kernels, 0);
730  for (auto & rk : ray_kernels)
731  {
732  AuxRayKernel * aux_rk = dynamic_cast<AuxRayKernel *>(rk);
733  if (aux_rk)
734  vars_to_be_zeroed.insert(aux_rk->variable().name());
735  }
736 
737  std::vector<std::string> vars_to_be_zeroed_vec(vars_to_be_zeroed.begin(),
738  vars_to_be_zeroed.end());
739  _fe_problem.getAuxiliarySystem().zeroVariables(vars_to_be_zeroed_vec);
740 }
void getRayKernels(std::vector< RayKernelBase *> &result, SubdomainID id, THREAD_ID tid)
Fills the active RayKernels associated with this study and a block into result.
FEProblemBase & _fe_problem
const std::string & name() const
MooseVariableFE< Real > & variable()
Gets the variable this AuxRayKernel contributes to.
Definition: AuxRayKernel.h:38
AuxiliarySystem & getAuxiliarySystem()
virtual void zeroVariables(std::vector< std::string > &vars_to_be_zeroed)

Member Data Documentation

◆ _always_cache_traces

bool RayTracingStudy::_always_cache_traces
protected

Whether or not to cache traces on every trace execution.

Definition at line 883 of file RayTracingStudy.h.

Referenced by shouldCacheTrace().

◆ _aux_data_on_cache_traces

const bool RayTracingStudy::_aux_data_on_cache_traces
protected

Whether or not to store the Ray aux data on the cache traces.

Definition at line 887 of file RayTracingStudy.h.

Referenced by auxDataOnCacheTraces().

◆ _b_box

libMesh::BoundingBox RayTracingStudy::_b_box
private

Nodal bounding box for the domain.

Definition at line 1067 of file RayTracingStudy.h.

Referenced by boundingBox(), isRectangularDomain(), and RayTracingStudy().

◆ _bank_rays_on_completion

const bool RayTracingStudy::_bank_rays_on_completion
protected

Whether or not to bank rays on completion.

Definition at line 878 of file RayTracingStudy.h.

Referenced by bankRaysOnCompletion(), onCompleteRay(), and rayBank().

◆ _cached_traces

std::vector<TraceData> RayTracingStudy::_cached_traces
private

Storage for the cached traces.

Definition at line 1047 of file RayTracingStudy.h.

Referenced by executeStudy(), and getCachedTraces().

◆ _called_initial_setup

bool RayTracingStudy::_called_initial_setup
private

Whether or not we've called initial setup - used to stop from late registration.

Definition at line 1139 of file RayTracingStudy.h.

Referenced by executeStudy(), getRayBCs(), getRayKernels(), initialSetup(), and registerRayDataInternal().

◆ _comm

const Parallel::Communicator& RayTracingStudy::_comm
protected

The Communicator.

Definition at line 863 of file RayTracingStudy.h.

Referenced by LotsOfRaysRayStudy::defineRays().

◆ _data_on_cache_traces

const bool RayTracingStudy::_data_on_cache_traces
protected

Whether or not to store the Ray data on the cache traces.

Definition at line 885 of file RayTracingStudy.h.

Referenced by dataOnCacheTraces().

◆ _domain_max_length

Real RayTracingStudy::_domain_max_length
private

An inflated max distance for the domain.

Definition at line 1071 of file RayTracingStudy.h.

Referenced by domainMaxLength().

◆ _elem_index_helper

ElemIndexHelper RayTracingStudy::_elem_index_helper
private

Helper for defining a local contiguous index for each element.

Definition at line 1142 of file RayTracingStudy.h.

Referenced by internalSidesetSetup(), localElemIndexSetup(), nonPlanarSideSetup(), and sideIsNonPlanar().

◆ _ending_distance

Real RayTracingStudy::_ending_distance
private

Total distance traveled by Rays that end on this processor.

Definition at line 1128 of file RayTracingStudy.h.

Referenced by endingDistance(), executeStudy(), and onCompleteRay().

◆ _ending_intersections

unsigned long long int RayTracingStudy::_ending_intersections
private

Total number of Ray/element intersections for Rays that finished on this processor.

Definition at line 1115 of file RayTracingStudy.h.

Referenced by endingIntersections(), executeStudy(), and onCompleteRay().

◆ _ending_max_intersections

unsigned int RayTracingStudy::_ending_max_intersections
private

Max number of intersections for Rays that finished on this processor.

Definition at line 1117 of file RayTracingStudy.h.

Referenced by endingMaxIntersections(), executeStudy(), and onCompleteRay().

◆ _ending_max_processor_crossings

unsigned int RayTracingStudy::_ending_max_processor_crossings
private

Max number of total processor crossings for Rays that finished on this processor.

Definition at line 1108 of file RayTracingStudy.h.

Referenced by endingMaxProcessorCrossings(), executeStudy(), and onCompleteRay().

◆ _ending_max_trajectory_changes

unsigned int RayTracingStudy::_ending_max_trajectory_changes
private

Max number of trajectory changes for Rays that finished on this processor.

Definition at line 1119 of file RayTracingStudy.h.

Referenced by executeStudy(), and onCompleteRay().

◆ _ending_processor_crossings

unsigned long long int RayTracingStudy::_ending_processor_crossings
private

Total number of processor crossings for Rays that finished on this processor.

Definition at line 1106 of file RayTracingStudy.h.

Referenced by endingProcessorCrossings(), executeStudy(), and onCompleteRay().

◆ _execution_start_time

std::chrono::steady_clock::time_point RayTracingStudy::_execution_start_time
private

Timing.

Definition at line 1025 of file RayTracingStudy.h.

Referenced by executeStudy().

◆ _execution_time

std::chrono::steady_clock::duration RayTracingStudy::_execution_time
private

Definition at line 1026 of file RayTracingStudy.h.

Referenced by executeStudy(), executionTime(), and executionTimeNano().

◆ _generation_time

std::chrono::steady_clock::duration RayTracingStudy::_generation_time
private

Definition at line 1027 of file RayTracingStudy.h.

Referenced by executeStudy(), and generationTime().

◆ _has_non_planar_sides

bool RayTracingStudy::_has_non_planar_sides
private

Whether or not the local mesh has elements with non-planar sides.

Definition at line 1059 of file RayTracingStudy.h.

Referenced by nonPlanarSideSetup(), and sideIsNonPlanar().

◆ _has_same_level_active_elems

bool RayTracingStudy::_has_same_level_active_elems
private

Whether or not the mesh has active elements of the same level.

Definition at line 1064 of file RayTracingStudy.h.

Referenced by hasSameLevelActiveElems(), and meshChanged().

◆ _internal_sidesets

std::set<BoundaryID> RayTracingStudy::_internal_sidesets
private

The BoundaryIDs on the local mesh that have internal RayBCs.

Definition at line 1055 of file RayTracingStudy.h.

Referenced by getInternalSidesets(), hasInternalSidesets(), and internalSidesetSetup().

◆ _internal_sidesets_map

std::vector<std::vector<std::vector<BoundaryID> > > RayTracingStudy::_internal_sidesets_map
private

Internal sideset data, if internal sidesets exist (indexed with getLocalElemIndex())

Definition at line 1057 of file RayTracingStudy.h.

Referenced by internalSidesetSetup().

◆ _local_trace_ray_results

std::vector<unsigned long long int> RayTracingStudy::_local_trace_ray_results
private

Cumulative results on this processor from the threaded TraceRay objects.

Definition at line 1133 of file RayTracingStudy.h.

Referenced by executeStudy(), and localTraceRayResult().

◆ _loose_b_box

libMesh::BoundingBox RayTracingStudy::_loose_b_box
private

Loose nodal bounding box for the domain.

Definition at line 1069 of file RayTracingStudy.h.

Referenced by looseBoundingBox(), and RayTracingStudy().

◆ _max_intersections

unsigned int RayTracingStudy::_max_intersections
private

Max number of intersections for a single Ray.

Definition at line 1123 of file RayTracingStudy.h.

Referenced by executeStudy(), and maxIntersections().

◆ _max_processor_crossings

unsigned int RayTracingStudy::_max_processor_crossings
private

Max number of processor crossings for all Rays.

Definition at line 1112 of file RayTracingStudy.h.

Referenced by executeStudy(), and maxProcessorCrossings().

◆ _max_trajectory_changes

unsigned int RayTracingStudy::_max_trajectory_changes
private

Max number of trajectory changes for a single Ray.

Definition at line 1125 of file RayTracingStudy.h.

Referenced by executeStudy(), and maxTrajectoryChanges().

◆ _mesh

MooseMesh& RayTracingStudy::_mesh
protected

◆ _non_planar_sides

std::vector<std::vector<unsigned short> > RayTracingStudy::_non_planar_sides
private

Non planar side data, which is for quick checking if an elem side is non-planar We use unsigned short here to avoid a std::vector<bool>; 0 = false, otherwise true.

Definition at line 1062 of file RayTracingStudy.h.

Referenced by nonPlanarSideSetup(), and sideIsNonPlanar().

◆ _num_cached

std::vector<std::size_t> RayTracingStudy::_num_cached
private

Number of currently cached objects for Jacobian/residual for each thread.

Definition at line 1052 of file RayTracingStudy.h.

Referenced by executeStudy(), jacobianSetup(), postOnSegment(), and residualSetup().

◆ _parallel_ray_study

const std::unique_ptr<ParallelRayStudy> RayTracingStudy::_parallel_ray_study
private

◆ _pid

const processor_id_type RayTracingStudy::_pid
protected

◆ _propagation_time

std::chrono::steady_clock::duration RayTracingStudy::_propagation_time
private

Definition at line 1028 of file RayTracingStudy.h.

Referenced by executeStudy(), and propagationTime().

◆ _ray_aux_data_map

std::unordered_map<std::string, RayDataIndex> RayTracingStudy::_ray_aux_data_map
private

The map from Ray aux data names to index.

Definition at line 1034 of file RayTracingStudy.h.

Referenced by getRayDataIndexInternal(), and registerRayDataInternal().

◆ _ray_aux_data_names

std::vector<std::string> RayTracingStudy::_ray_aux_data_names
private

The names for each Ray aux data entry.

Definition at line 1039 of file RayTracingStudy.h.

Referenced by getRayDataNameInternal(), hasRayAuxData(), rayAuxDataNames(), rayAuxDataSize(), and registerRayDataInternal().

◆ _ray_bank

std::vector<std::shared_ptr<Ray> > RayTracingStudy::_ray_bank
private

Cumulative Ray bank - stored only when _bank_rays_on_completion.

Definition at line 1093 of file RayTracingStudy.h.

Referenced by executeStudy(), onCompleteRay(), and rayBank().

◆ _ray_data_map

std::unordered_map<std::string, RayDataIndex> RayTracingStudy::_ray_data_map
private

The map from Ray data names to index.

Definition at line 1032 of file RayTracingStudy.h.

Referenced by getRayDataIndexInternal(), and registerRayDataInternal().

◆ _ray_data_names

std::vector<std::string> RayTracingStudy::_ray_data_names
private

The names for each Ray data entry.

Definition at line 1037 of file RayTracingStudy.h.

Referenced by getRayDataNameInternal(), hasRayData(), rayDataNames(), rayDataSize(), and registerRayDataInternal().

◆ _ray_dependent_subdomain_setup

const bool RayTracingStudy::_ray_dependent_subdomain_setup
protected

Whether or not subdomain setup is dependent on the Ray.

Definition at line 880 of file RayTracingStudy.h.

Referenced by rayDependentSubdomainSetup().

◆ _ray_kernel_coverage_check

const bool RayTracingStudy::_ray_kernel_coverage_check
protected

Whether or not to perform coverage checks on RayKernels.

Definition at line 868 of file RayTracingStudy.h.

Referenced by coverageChecks().

◆ _ray_max_distance

const Real RayTracingStudy::_ray_max_distance
protected

Max distance a Ray can travel before being killed (can change)

Definition at line 891 of file RayTracingStudy.h.

Referenced by rayMaxDistance().

◆ _registered_ray_map

std::unordered_map<std::string, RayID>& RayTracingStudy::_registered_ray_map
private

Map from registered Ray name to ID.

Definition at line 1042 of file RayTracingStudy.h.

Referenced by registeredRayID(), registeredRaySetup(), and registerRay().

◆ _replicated_next_ray_id

RayID RayTracingStudy::_replicated_next_ray_id
private

Storage for the next available replicated RayID, obtained via generateReplicatedRayID()

Definition at line 1097 of file RayTracingStudy.h.

Referenced by generateReplicatedRayID(), and resetReplicatedRayIDs().

◆ _reverse_registered_ray_map

std::vector<std::string>& RayTracingStudy::_reverse_registered_ray_map
private

Map from registered Ray ID to name.

Definition at line 1044 of file RayTracingStudy.h.

Referenced by registeredRayName(), and registerRay().

◆ _segment_qrule

std::unique_ptr<libMesh::QBase> RayTracingStudy::_segment_qrule
private

Quadrature rule for laying points across a 1D ray segment.

Definition at line 1103 of file RayTracingStudy.h.

Referenced by buildSegmentQuadrature(), and initialSetup().

◆ _segments_on_cache_traces

const bool RayTracingStudy::_segments_on_cache_traces
protected

Whether or not to cache individual element segments when caching.

Definition at line 889 of file RayTracingStudy.h.

Referenced by segmentsOnCacheTraces().

◆ _spin_mutex

Threads::spin_mutex RayTracingStudy::_spin_mutex
mutableprivate

◆ _subdomain_hmax

std::unordered_map<SubdomainID, Real> RayTracingStudy::_subdomain_hmax
private

The cached hmax for all elements in a subdomain.

Definition at line 1136 of file RayTracingStudy.h.

Referenced by subdomainHmax(), and subdomainHMaxSetup().

◆ _threaded_cache_ray_bc

std::vector<TheWarehouse::QueryCache<AttribBoundaries> > RayTracingStudy::_threaded_cache_ray_bc
private

Threaded cached boundary query for RayBC objects pertaining to this study.

Definition at line 1078 of file RayTracingStudy.h.

Referenced by getRayBCs().

◆ _threaded_cache_ray_kernel

std::vector<TheWarehouse::QueryCache<AttribSubdomains> > RayTracingStudy::_threaded_cache_ray_kernel
private

Threaded cached subdomain query for RayKernelBase objects pertaining to this study.

Definition at line 1076 of file RayTracingStudy.h.

Referenced by getRayKernels().

◆ _threaded_cached_normals

std::vector<std::unordered_map<std::pair<const Elem *, unsigned short>, Point> > RayTracingStudy::_threaded_cached_normals
private

Threaded cache for side normals that have been computed already during tracing.

Definition at line 1091 of file RayTracingStudy.h.

Referenced by executeStudy(), and getSideNormal().

◆ _threaded_cached_traces

std::vector<std::vector<TraceData> > RayTracingStudy::_threaded_cached_traces
private

The threaded storage for cached traces.

Definition at line 1049 of file RayTracingStudy.h.

Referenced by executeStudy(), and initThreadedCachedTrace().

◆ _threaded_current_ray_kernels

std::vector<std::vector<RayKernelBase *> > RayTracingStudy::_threaded_current_ray_kernels
private

The current RayKernel objects for each thread.

Definition at line 1082 of file RayTracingStudy.h.

Referenced by currentRayKernels(), currentRayKernelsWrite(), executeStudy(), and segmentSubdomainSetup().

◆ _threaded_elem_side_builders

std::vector<libMesh::ElemSideBuilder> RayTracingStudy::_threaded_elem_side_builders
private

Threaded helpers for building element sides without extraneous allocation.

Definition at line 1021 of file RayTracingStudy.h.

Referenced by elemSide().

◆ _threaded_fe_face

std::vector<std::unique_ptr<libMesh::FEBase> > RayTracingStudy::_threaded_fe_face
private

Face FE used for computing face normals for each thread.

Definition at line 1086 of file RayTracingStudy.h.

Referenced by getSideNormal(), and RayTracingStudy().

◆ _threaded_next_ray_id

std::vector<RayID> RayTracingStudy::_threaded_next_ray_id
private

Storage for the next available unique RayID, obtained via generateUniqueRayID()

Definition at line 1095 of file RayTracingStudy.h.

Referenced by generateUniqueRayID(), and resetUniqueRayIDs().

◆ _threaded_q_face

std::vector<std::unique_ptr<libMesh::QBase> > RayTracingStudy::_threaded_q_face
private

Face quadrature used for computing face normals for each thread.

Definition at line 1088 of file RayTracingStudy.h.

Referenced by RayTracingStudy().

◆ _threaded_ray_object_registration

std::vector<std::vector<std::set<const RayTracingObject *> > > RayTracingStudy::_threaded_ray_object_registration
private

Threaded storage for all of the RayTracingObjects associated with a single Ray.

Definition at line 1080 of file RayTracingStudy.h.

Referenced by getRayBCs(), getRayKernels(), and registeredRaySetup().

◆ _threaded_trace_ray

std::vector<std::shared_ptr<TraceRay> > RayTracingStudy::_threaded_trace_ray
private

The TraceRay objects for each thread (they do the physical tracing)

Definition at line 1084 of file RayTracingStudy.h.

Referenced by executeStudy(), meshChanged(), RayTracingStudy(), and traceRay().

◆ _tolerate_failure

const bool RayTracingStudy::_tolerate_failure
protected

Whether or not to tolerate a Ray Tracing failure.

Definition at line 876 of file RayTracingStudy.h.

Referenced by executeStudy(), and tolerateFailure().

◆ _total_distance

Real RayTracingStudy::_total_distance
private

Total distance traveled by all Rays.

Definition at line 1130 of file RayTracingStudy.h.

Referenced by executeStudy(), and totalDistance().

◆ _total_intersections

unsigned long long int RayTracingStudy::_total_intersections
private

Total number of Ray/element intersections.

Definition at line 1121 of file RayTracingStudy.h.

Referenced by executeStudy(), and totalIntersections().

◆ _total_processor_crossings

unsigned long long int RayTracingStudy::_total_processor_crossings
private

Total number of processor crossings.

Definition at line 1110 of file RayTracingStudy.h.

Referenced by executeStudy(), and totalProcessorCrossings().

◆ _total_volume

Real RayTracingStudy::_total_volume
private

The total volume of the domain.

Definition at line 1073 of file RayTracingStudy.h.

Referenced by totalVolume().

◆ _use_internal_sidesets

const bool RayTracingStudy::_use_internal_sidesets
protected

Whether or not to use the internal sidesets in ray tracing.

Definition at line 874 of file RayTracingStudy.h.

Referenced by internalSidesetSetup().

◆ _use_ray_registration

const bool RayTracingStudy::_use_ray_registration
protected

Whether or not to use Ray registration.

Definition at line 872 of file RayTracingStudy.h.

Referenced by getRayBCs(), getRayKernels(), registeredRayID(), registeredRayName(), registeredRaySetup(), registerRay(), and useRayRegistration().

◆ _verify_rays

const bool RayTracingStudy::_verify_rays
protected

Whether or not to verify if Rays have valid information before being traced.

Definition at line 893 of file RayTracingStudy.h.

Referenced by verifyRays().

◆ _verify_trace_intersections

const bool RayTracingStudy::_verify_trace_intersections
protected

Whether or not to verify the trace intersections in devel and dbg modes.

Definition at line 896 of file RayTracingStudy.h.

Referenced by verifyTraceIntersections().

◆ _warn_non_planar

const bool RayTracingStudy::_warn_non_planar
protected

Whether not to warn if non-planar faces are found.

Definition at line 870 of file RayTracingStudy.h.

Referenced by nonPlanarSideSetup(), and warnNonPlanar().


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