ArbitraryQuadrature.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 "Moose.h"

#include "libmesh/quadrature.h"

class ArbitraryQuadrature : public libMesh::QBase
{
public:
  ArbitraryQuadrature(const unsigned int _dim,
                      const libMesh::Order _order = libMesh::INVALID_ORDER);

  ArbitraryQuadrature(const ArbitraryQuadrature &) = default;
  ArbitraryQuadrature(ArbitraryQuadrature &&) = default;
  ArbitraryQuadrature & operator=(const ArbitraryQuadrature &) = default;
  ArbitraryQuadrature & operator=(ArbitraryQuadrature &&) = default;
  virtual ~ArbitraryQuadrature() = default;

  libMesh::QuadratureType type() const override;

  void setPoints(const std::vector<libMesh::Point> & points);

  void setWeights(const std::vector<libMesh::Real> & weights);

  virtual bool shapes_need_reinit() override { return true; }

  virtual std::unique_ptr<libMesh::QBase> clone() const override;

private:
#ifdef LIBMESH_QBASE_INIT_ARGUMENTS_REMOVED
  void init_1D() override;
  void init_2D() override;
  void init_3D() override;
#else
  void init_1D(const libMesh::ElemType _type = libMesh::INVALID_ELEM,
               unsigned int p_level = 0) override;
  void init_2D(const libMesh::ElemType _type = libMesh::INVALID_ELEM,
               unsigned int p_level = 0) override;
  void init_3D(const libMesh::ElemType _type = libMesh::INVALID_ELEM,
               unsigned int p_level = 0) override;
#endif // LIBMESH_QBASE_INIT_ARGUMENTS_REMOVED
};