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

#pragma once

// MOOSE includes
#include "ThreadedElementLoop.h"

#include "libmesh/elem_range.h"

class FEProblemBase;
class NonlinearSystemBase;
class MaterialPropertyStorage;
class MaterialData;
class Assembly;

class ProjectMaterialProperties : public ThreadedElementLoop<ConstElemPointerRange>
{
public:
  ProjectMaterialProperties(bool refine,
                            FEProblemBase & fe_problem,
                            MaterialPropertyStorage & material_props,
                            MaterialPropertyStorage & bnd_material_props,
                            std::vector<std::vector<std::unique_ptr<Assembly>>> & assembly);

  // Splitting Constructor
  ProjectMaterialProperties(ProjectMaterialProperties & x, Threads::split split);

  virtual ~ProjectMaterialProperties();

  virtual void subdomainChanged() override;
  virtual void onElement(const Elem * elem) override;
  virtual void onBoundary(const Elem * elem,
                          unsigned int side,
                          BoundaryID bnd_id,
                          const Elem * lower_d_elem = nullptr) override;
  virtual void onInternalSide(const Elem * elem, unsigned int side) override;

  void join(const ProjectMaterialProperties & /*y*/);

protected:
  bool _refine;
  FEProblemBase & _fe_problem;
  MaterialPropertyStorage & _material_props;
  MaterialPropertyStorage & _bnd_material_props;
  std::vector<std::vector<std::unique_ptr<Assembly>>> & _assembly;
  bool _need_internal_side_material;

  const MaterialWarehouse & _materials;
  const MaterialWarehouse & _discrete_materials;

private:
  bool shouldComputeInternalSide(const Elem & /*elem*/, const Elem & /*neighbor*/) const override;
};

inline bool
ProjectMaterialProperties::shouldComputeInternalSide(const Elem & /*elem*/,
                                                     const Elem & /*neighbor*/) const
{
  return true;
}