SetupDebugAction.C

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

#include "SetupDebugAction.h"
#include "FEProblem.h"
#include "ActionWarehouse.h"
#include "Factory.h"
#include "Output.h"
#include "MooseApp.h"
#include "MooseObjectAction.h"
#include "ActionFactory.h"
#include "AddAuxVariableAction.h"
#include "MooseUtils.h"
#include "BlockRestrictionDebugOutput.h"

using namespace libMesh;

registerMooseAction("MooseApp", SetupDebugAction, "add_output");

InputParameters
SetupDebugAction::validParams()
{
  InputParameters params = Action::validParams();
  params.addParam<unsigned int>(
      "show_top_residuals", 0, "The number of top residuals to print out (0 = no output)");
  params.addParam<bool>(
      "show_var_residual_norms",
      false,
      "Print the residual norms of the individual solution variables at each nonlinear iteration");
  params.addParam<bool>("show_action_dependencies", false, "Print out the action dependencies");
  params.addParam<bool>("show_actions", false, "Print out the actions being executed");
  params.addParam<bool>(
      "show_parser", false, "Shows parser block extraction and debugging information");
  params.addParam<bool>(
      "show_material_props",
      false,
      "Print out the material properties supplied for each block, face, neighbor, and/or sideset");
  params.addParam<bool>("show_chain_control_data",
                        false,
                        "Print out the chain control data on every time step setup");
  params.addParam<bool>("show_controllable",
                        false,
                        "Print out the controllable parameters from all input parameters");
  params.addParam<bool>("show_mesh_meta_data", false, "Print out the available mesh meta data");
  params.addParam<bool>(
      "show_reporters", false, "Print out information about the declared and requested Reporters");
  params.addParam<bool>(
      "show_mesh_generators", false, "Print out the mesh generators being executed");

  ExecFlagEnum print_on = MooseUtils::getDefaultExecFlagEnum();
  print_on.addAvailableFlags(EXEC_TRANSFER);
  print_on.addAvailableFlags(EXEC_FAILED);
  print_on.addAvailableFlags(EXEC_ALWAYS);
  params.addParam<ExecFlagEnum>(
      "show_execution_order",
      print_on,
      "Print more information about the order of execution during calculations");
  params.addDeprecatedParam<bool>(
      "pid_aux",
      "Add a AuxVariable named \"pid\" that shows the processors and partitioning",
      "pid_aux is deprecated, use output_process_domains");
  params.addParam<bool>(
      "output_process_domains",
      false,
      "Add a AuxVariable named \"pid\" that shows the partitioning for each process");
  params.addParam<bool>(
      "show_functors", false, "Whether to print information about the functors in the problem");
  params.addParam<MultiMooseEnum>(
      "show_block_restriction",
      BlockRestrictionDebugOutput::getScopes("none"),
      "Print out active objects like variables supplied for each block.");
  params.addParam<bool>("show_block_restriction_groups",
                        false,
                        "Print groups of objects with identical block restrictions.");
  params.addParam<bool>("show_boundary_restriction_groups",
                        false,
                        "Print groups of objects with identical boundary restrictions.");
  params.addParam<bool>(
      "error_on_residual_nan",
      false,
      "This option applies only to dbg and devel modes. If enabled, residual contributions are "
      "checked for NaN or Inf values; if found, an error is reported.");

  params.addClassDescription("Adds various debugging type output to the simulation system.");

  return params;
}

SetupDebugAction::SetupDebugAction(const InputParameters & parameters) : Action(parameters)
{
  _awh.showActionDependencies(getParam<bool>("show_action_dependencies"));
  _awh.showActions(getParam<bool>("show_actions"));
  _awh.showParser(getParam<bool>("show_parser"));
  _awh.mooseApp().getMeshGeneratorSystem().setVerbose(getParam<bool>("show_mesh_generators"));
}

void
SetupDebugAction::act()
{
  // Material properties
  if (_pars.get<bool>("show_material_props"))
  {
    const std::string type = "MaterialPropertyDebugOutput";
    auto params = _factory.getValidParams(type);
    _problem->addOutput(type, "_moose_material_property_debug_output", params);
  }

  // Variable residual norms
  if (_pars.get<bool>("show_var_residual_norms"))
  {
    const std::string type = "VariableResidualNormsDebugOutput";
    auto params = _factory.getValidParams(type);
    // Add one for every nonlinear system
    for (const auto & sys_name : _problem->getNonlinearSystemNames())
    {
      params.set<NonlinearSystemName>("nl_sys") = sys_name;
      _problem->addOutput(type, "_moose_variable_residual_norms_debug_output_" + sys_name, params);
    }
  }

  // Top residuals
  if (_pars.get<unsigned int>("show_top_residuals") > 0)
  {
    const std::string type = "TopResidualDebugOutput";
    auto params = _factory.getValidParams(type);
    params.set<unsigned int>("num_residuals") = _pars.get<unsigned int>("show_top_residuals");
    _problem->addOutput(type, "_moose_top_residual_debug_output", params);
  }

  // Print full names of mesh meta data
  if (getParam<bool>("show_mesh_meta_data"))
  {
    _console << "Mesh meta data:\n";
    for (auto it = _app.getRestartableDataMapBegin(); it != _app.getRestartableDataMapEnd(); ++it)
      if (it->first == MooseApp::MESH_META_DATA)
        for (auto & data : it->second.first)
          _console << " " << data.name() << std::endl;
  }

  // Print Reporter information
  if (getParam<bool>("show_reporters"))
  {
    const std::string type = "ReporterDebugOutput";
    auto params = _factory.getValidParams(type);
    _problem->addOutput(type, "_moose_reporter_debug_output", params);
  }

  // Print execution information in all loops
  if (parameters().isParamSetByUser("show_execution_order"))
    _problem->setExecutionPrinting(getParam<ExecFlagEnum>("show_execution_order"));

  // Add pid aux
  if (getParam<bool>("output_process_domains") ||
      (isParamValid("pid_aux") && getParam<bool>("pid_aux")))
  {
    if (_problem->hasVariable("pid"))
      paramError("output_process_domains", "Variable with the name \"pid\" already exists");

    auto fe_type = FEType(CONSTANT, MONOMIAL);
    auto type = AddAuxVariableAction::variableType(fe_type);
    auto var_params = _factory.getValidParams(type);
    _problem->addAuxVariable(type, "pid", var_params);

    InputParameters params = _factory.getValidParams("ProcessorIDAux");
    params.set<AuxVariableName>("variable") = "pid";
    params.set<ExecFlagEnum>("execute_on") = {EXEC_INITIAL, EXEC_TIMESTEP_BEGIN};
    _problem->addAuxKernel("ProcessorIDAux", "pid_aux", params);
  }

  // Add functor output
  if (getParam<bool>("show_functors"))
    _problem->setFunctorOutput(getParam<bool>("show_functors"));

  // Add chain control data output
  if (getParam<bool>("show_chain_control_data"))
    _problem->setChainControlDataOutput(true);

  // Block-restriction
  const MultiMooseEnum & block_restriction_scope =
      _pars.get<MultiMooseEnum>("show_block_restriction");
  const bool show_block_restriction_map =
      block_restriction_scope.isValid() && !block_restriction_scope.contains("none");
  const bool show_block_restriction_groups = getParam<bool>("show_block_restriction_groups");
  const bool show_boundary_restriction_groups = getParam<bool>("show_boundary_restriction_groups");
  if (show_block_restriction_map || show_block_restriction_groups ||
      show_boundary_restriction_groups)
  {
    const std::string type = "BlockRestrictionDebugOutput";
    auto params = _factory.getValidParams(type);
    params.set<MultiMooseEnum>("scope") = block_restriction_scope;
    params.set<bool>("show_block_restriction_map") = show_block_restriction_map;
    params.set<bool>("show_block_restriction_groups") = show_block_restriction_groups;
    params.set<bool>("show_boundary_restriction_groups") = show_boundary_restriction_groups;
    _problem->addOutput(type, "_moose_block_restriction_debug_output", params);
  }

  // Controllable output
  if (getParam<bool>("show_controllable"))
  {
    const std::string type = "ControlOutput";
    auto params = _factory.getValidParams(type);
    _problem->addOutput(type, "_moose_controllable_debug_output", params);
  }

  // Enable residual NaN/Inf-checking
  if (getParam<bool>("error_on_residual_nan"))
  {
#ifdef NDEBUG
    mooseError("The parameter 'error_on_residual_nan' may only be set to 'true' for 'dbg' and "
               "'devel' modes.");
#else
    _problem->setCheckResidualForNans(true);
#endif
  }
}