FECellDealii.h

Go to the documentation of this file.

/******************************************************************************
 * Copyright (c) 2021.                                                        *
 * The Regents of the University of Michigan and DFT-EFE developers.          *
 *                                                                            *
 * This file is part of the DFT-EFE code.                                     *
 *                                                                            *
 * DFT-EFE is free software: you can redistribute it and/or modify            *
 *   it under the terms of the Lesser GNU General Public License as           *
 *   published by the Free Software Foundation, either version 3 of           *
 *   the License, or (at your option) any later version.                      *
 *                                                                            *
 * DFT-EFE is distributed in the hope that it will be useful, but             *
 *   WITHOUT ANY WARRANTY; without even the implied warranty                  *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.                     *
 *   See the Lesser GNU General Public License for more details.              *
 *                                                                            *
 * You should have received a copy of the GNU Lesser General Public           *
 *   License at the top level of DFT-EFE distribution.  If not, see           *
 *   <https://www.gnu.org/licenses/>.                                         *
 ******************************************************************************/
 
/*
 * @author Bikash Kanungo, Vishal Subramanian
 */
 
#ifndef dftefeFECellDealii_h
#define dftefeFECellDealii_h
 
#include <utils/TypeConfig.h>
#include <deal.II/dofs/dof_accessor.h>
#include <utils/Point.h>
#include "FECellBase.h"
namespace dftefe
{
  namespace basis
  {
    template <size_type dim>
    class FECellDealii : public FECellBase
    {
      using DealiiFECellIterator =
        typename dealii::DoFHandler<dim>::active_cell_iterator;
 
    public:
      FECellDealii(DealiiFECellIterator dealiiFECellIter);
 
      void
      getVertices(std::vector<utils::Point> &points) const override;
 
      void
      getVertex(size_type i, utils::Point &point) const override;
 
      std::vector<std::shared_ptr<utils::Point>>
      getNodalPoints() const override;
 
      size_type
      getId() const override;
 
      bool
      isPointInside(const utils::Point &point) const override;
 
      bool
      isAtBoundary(const size_type i) const override;
 
      bool
      isAtBoundary() const override;
 
      virtual bool
      hasPeriodicNeighbor(const size_type i) const override;
 
      double
      diameter() const override;
 
      void
      center(dftefe::utils::Point &centerPoint) const override;
      void
      setRefineFlag() override;
 
      void
      clearRefineFlag() override;
 
      double
      minimumVertexDistance() const override;
 
      double
      distanceToUnitCell(dftefe::utils::Point &parametricPoint) const override;
 
      void
      setCoarsenFlag() override;
 
      void
      clearCoarsenFlag() override;
 
      bool
      isActive() const override;
 
      bool
      isLocallyOwned() const override;
 
      bool
      isGhost() const override;
 
      bool
      isArtificial() const override;
 
      size_type
      getDim() const override;
 
      void
      getParametricPoint(const utils::Point &   realPoint,
                         const CellMappingBase &cellMapping,
                         utils::Point &         parametricPoint) const override;
 
      void
      getRealPoint(const utils::Point &   parametricPoint,
                   const CellMappingBase &cellMapping,
                   utils::Point &         realPoint) const override;
 
      void
      cellNodeIdtoGlobalNodeId(
        std::vector<global_size_type> &vecId) const override;
 
      size_type
      getFaceBoundaryId(size_type faceId) const override;
 
      void
      getFaceDoFGlobalIndices(
        size_type                      faceId,
        std::vector<global_size_type> &vecNodeId) const override;
 
 
      size_type
      getFEOrder() const override;
 
      DealiiFECellIterator &
      getDealiiFECellIter();
 
    private:
      DealiiFECellIterator d_dealiiFECellIter;
 
 
    }; // end of class FECellDealii
  }    // namespace basis
} // namespace dftefe
#include "FECellDealii.t.cpp"
#endif // dftefeFECellDealii_h