DebugRay.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

// The following can be uncommented for detailed output for
// the trace of a Ray or multiple Rays
// #define DEBUG_RAY_IF (*_current_ray)->id() == 0
// #define DEBUG_RAY_MESH_IF (*_current_ray)->id() == 0
// #define DEBUG_RAY_INTERSECTIONS

#include "MooseError.h"

#ifndef DEBUG_RAY_IF
#define debugRay(...) ((void)0)
#else
#define debugRay(...)                                                                              \
  if (DEBUG_RAY_IF)                                                                                \
  {                                                                                                \
    std::ostringstream ss;                                                                         \
    moose::internal::mooseStreamAll(ss, __VA_ARGS__);                                              \
    libMesh::err << "[pid/tid/id/line:" << _pid << "/" << _tid << "/" << (*_current_ray)->id()     \
                 << "/" << std::setw(4) << std::left << __LINE__ << "]: " << ss.str() << std::endl \
                 << std::flush;                                                                    \
  }
#endif

#ifndef DEBUG_RAY_INTERSECTIONS
#define debugRaySimple(...) ((void)0)
#else
#define debugRaySimple(...)                                                                        \
  if (debug)                                                                                       \
  {                                                                                                \
    std::ostringstream ss;                                                                         \
    moose::internal::mooseStreamAll(ss, __VA_ARGS__);                                              \
    libMesh::err << ss.str() << std::endl << std::flush;                                           \
  }
#endif

#ifndef DEBUG_RAY_MESH_IF
#define possiblySaveDebugRayMesh() ((void)0)
#else
#define possiblySaveDebugRayMesh()                                                                 \
  if (DEBUG_RAY_MESH_IF && _debug_mesh)                                                            \
  {                                                                                                \
    _debug_mesh->prepare_for_use();                                                                \
    _debug_mesh->write("debug_ray" + std::to_string((*_current_ray)->id()) + "_pc" +               \
                       std::to_string((*_current_ray)->processorCrossings()) + "_pid" +            \
                       std::to_string(_pid) + ".e");                                               \
    delete _debug_mesh;                                                                            \
  }
#endif

#ifndef DEBUG_RAY_MESH_IF
#define possiblyAddDebugRayMeshPoint(start, end) ((void)0)
#else
#define possiblyAddDebugRayMeshPoint(start, end)                                                   \
  if (_debug_mesh && DEBUG_RAY_MESH_IF)                                                            \
  {                                                                                                \
    _debug_mesh->add_point(start);                                                                 \
    _debug_node_count++;                                                                           \
    _debug_mesh->add_point(end);                                                                   \
    _debug_node_count++;                                                                           \
                                                                                                   \
    Elem * elem = Elem::build(EDGE2).release();                                                    \
    elem->subdomain_id() = 0;                                                                      \
                                                                                                   \
    elem = _debug_mesh->add_elem(elem);                                                            \
                                                                                                   \
    elem->set_node(0, _debug_mesh->node_ptr((_debug_node_count - 2)));                             \
    elem->set_node(1, _debug_mesh->node_ptr((_debug_node_count - 2) + 1));                         \
  }
#endif