RayTracingExodus.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

// TODO: remove ignore warnings once std::tuple<> is instantiated as a StandardType
// Using push_parallel_vector_data with std::tuple<> leads to a -Wextra error
// https://github.com/libMesh/TIMPI/issues/52
#include "libmesh/ignore_warnings.h"
#include "RayTracingExodus.h"
#include "libmesh/restore_warnings.h"

// Local includes
#include "RayTracingStudy.h"

#include "libmesh/exodusII_io.h"

registerMooseObject("RayTracingApp", RayTracingExodus);

InputParameters
RayTracingExodus::validParams()
{
  auto params = RayTracingMeshOutput::validParams();
  params.addClassDescription("Outputs ray segments and data as segments using the Exodus format.");
  return params;
}

RayTracingExodus::RayTracingExodus(const InputParameters & params) : RayTracingMeshOutput(params) {}

void
RayTracingExodus::outputMesh()
{
  TIME_SECTION("outputMesh", 3, "Writing Ray Mesh");

  libMesh::ExodusII_IO eio(*_segment_mesh);
  eio.set_hdf5_writing(false);

  // With nodal data, we need to output these variables in write_timestep
  if (_output_data_nodal)
    eio.set_output_variables(_study.rayDataNames());
  // Otherwise, there's no variables to write in write_timestep
  else
    eio.set_output_variables(std::vector<std::string>());
  // Write the timestep, which is the mesh + nodal vars (if any)
  eio.write_timestep(filename(), *_es, 1, time() + _app.getGlobalTimeOffset());

  // This will default to write_element_data getting all available elemental vars
  eio.set_output_variables(std::vector<std::string>());
  // Write the elemental variables, which are the variables with the constant Ray field data
  eio.write_element_data(*_es);

  // We write a new file every time as we don't keep track of whether or not the rays change
  ++_file_num;
}