MeshTally.h

Go to the documentation of this file.

/********************************************************************/
/*                  SOFTWARE COPYRIGHT NOTIFICATION                 */
/*                             Cardinal                             */
/*                                                                  */
/*                  (c) 2021 UChicago Argonne, LLC                  */
/*                        ALL RIGHTS RESERVED                       */
/*                                                                  */
/*                 Prepared by UChicago Argonne, LLC                */
/*               Under Contract No. DE-AC02-06CH11357               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*             Prepared by Battelle Energy Alliance, LLC            */
/*               Under Contract No. DE-AC07-05ID14517               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*                 See LICENSE for full restrictions                */
/********************************************************************/
 
#pragma once
 
#include "TallyBase.h"
#include "OpenMCCellAverageProblem.h"
 
#include "openmc/tallies/filter_mesh.h"
 
namespace libMesh
{
class ReplicatedMesh;
}
 
class MeshTally : public TallyBase
{
public:
  static InputParameters validParams();
 
  MeshTally(const InputParameters & parameters);
 
  virtual std::pair<unsigned int, openmc::Filter *> spatialFilter() override;
 
  virtual void resetTally() override;
 
  virtual void gatherLinkedSum() override;
 
  virtual bool addingGlobalTally() const override { return _needs_global_tally && _instance == 0; }
 
protected:
  virtual Real storeResultsInner(const std::vector<unsigned int> & var_numbers,
                                 unsigned int local_score,
                                 const std::vector<OMCTensor> & tally_vals,
                                 bool norm_by_src_rate = true) override;
  void checkMeshTemplateAndTranslations();
 
  const std::string * _mesh_template_filename = nullptr;
 
  Point _mesh_translation;
 
  const unsigned int _instance;
 
  unsigned int _mesh_index;
 
  openmc::MeshFilter * _mesh_filter;
 
  openmc::UnstructuredMesh * _mesh_template;
 
  const bool _use_dof_map;
 
  std::unique_ptr<libMesh::ReplicatedMesh> _libmesh_mesh_copy;
  std::vector<unsigned int> _bin_to_element_mapping;
};