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

#include "Moose.h"
#include "MooseTypes.h"
#include "MooseEnum.h"

#include <set>
#include <map>
#include <memory>

// Forward declarations
class MooseMesh;
class MooseBase;

namespace libMesh
{
class Elem;
class PointLocatorBase;
class Point;
}

class DiracKernelInfo
{
public:
  DiracKernelInfo();
  virtual ~DiracKernelInfo();

public:
  void addPoint(const Elem * elem, const Point & p, const Real & value = 1);

  void clearPoints();

  bool hasPoint(const Elem * elem, const Point & p);

  std::set<const Elem *> & getElements() { return _elements; }

  typedef std::map<const Elem *, std::pair<std::vector<Point>, std::vector<Real>>> MultiPointMap;

  MultiPointMap & getPoints() { return _points; }

  void updatePointLocator(const MooseMesh & mesh);

  CreateMooseEnumClass(PointNotFoundBehavior, ERROR, WARNING, IGNORE);

  const Elem * findPoint(const Point & p,
                         const MooseMesh & mesh,
                         const std::set<SubdomainID> & blocks,
                         const PointNotFoundBehavior point_not_found_behavior,
                         const MooseBase & consumer);

protected:
  bool pointsFuzzyEqual(const Point &, const Point &);

  std::set<const Elem *> _elements;

  MultiPointMap _points;

  std::unique_ptr<libMesh::PointLocatorBase> _point_locator;

  const Real _point_equal_distance_sq;
};