PointValueAtXFEMInterface.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "GeneralUserObject.h"
#include "ElementPairLocator.h"
#include "MooseVariableFE.h"
 
#include "libmesh/vector_value.h"
 
// Forward Declarations
class PointValueAtXFEMInterface;
class XFEM;
class LineSegmentCutSetUserObject;
 
template <>
InputParameters validParams<PointValueAtXFEMInterface>();
 
class PointValueAtXFEMInterface : public GeneralUserObject
{
public:
  PointValueAtXFEMInterface(const InputParameters & parameters);
 
  virtual ~PointValueAtXFEMInterface() {}
 
  virtual void initialize() override;
  virtual void execute() override;
  virtual void finalize() override;
 
  std::map<unsigned int, Real> getValueAtPositiveLevelSet() const
  {
    return _values_positive_level_set_side;
  };
 
  std::map<unsigned int, Real> getValueAtNegativeLevelSet() const
  {
    return _values_negative_level_set_side;
  };
 
  std::map<unsigned int, RealVectorValue> getGradientAtPositiveLevelSet() const
  {
    return _grad_values_positive_level_set_side;
  };
 
  std::map<unsigned int, RealVectorValue> getGradientAtNegativeLevelSet() const
  {
    return _grad_values_negative_level_set_side;
  };
 
  unsigned int numberPoints() const { return _points.size(); };
 
protected:
  const Elem * getElemContainingPoint(const Point & p, bool positive_level_set);
 
  MooseMesh & _mesh;
 
  std::vector<Point> _points;
 
  std::unique_ptr<PointLocatorBase> _pl;
 
  std::shared_ptr<XFEM> _xfem;
 
  const ElementPairLocator::ElementPairList * _elem_pairs;
 
  const LineSegmentCutSetUserObject * _geo_cut;
 
  MooseVariableFEBase * _var;
 
  const unsigned int _level_set_var_number;
 
  const System & _system;
 
  const NumericVector<Number> * _solution;
 
  std::map<unsigned int, Real> _values_positive_level_set_side;
 
  std::map<unsigned int, Real> _values_negative_level_set_side;
 
  std::map<unsigned int, RealVectorValue> _grad_values_positive_level_set_side;
 
  std::map<unsigned int, RealVectorValue> _grad_values_negative_level_set_side;
};