BoundaryStencil.hpp

Go to the documentation of this file.

/*
 * ------------------------------------------------------------------------------------------------------------
 * SPDX-License-Identifier: LGPL-2.1-only
 *
 * Copyright (c) 2018-2020 Lawrence Livermore National Security LLC
 * Copyright (c) 2018-2020 The Board of Trustees of the Leland Stanford Junior University
 * Copyright (c) 2018-2020 Total, S.A
 * Copyright (c) 2019-     GEOSX Contributors
 * All rights reserved
 *
 * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
 * ------------------------------------------------------------------------------------------------------------
 */

#ifndef GEOSX_FINITEVOLUME_BOUNDARYSTENCIL_HPP_
#define GEOSX_FINITEVOLUME_BOUNDARYSTENCIL_HPP_

#include "StencilBase.hpp"
#include "codingUtilities/traits.hpp"

namespace geosx
{

struct BoundaryStencil_Traits
{
  using IndexContainerType = array2d< localIndex >;

  using IndexContainerViewType = traits::ViewType< IndexContainerType >;

  using IndexContainerViewConstType = traits::ViewTypeConst< IndexContainerType >;

  using WeightContainerType = array2d< real64 >;

  using WeightContainerViewType = traits::ViewType< WeightContainerType >;

  using WeightContainerViewConstType = traits::ViewTypeConst< WeightContainerType >;

  static constexpr localIndex NUM_POINT_IN_FLUX = 2;

  static constexpr localIndex MAX_STENCIL_SIZE = 2;
};

class BoundaryStencil : public StencilBase< BoundaryStencil_Traits, BoundaryStencil >,
  public BoundaryStencil_Traits
{
public:

  struct Order
  {
    static constexpr localIndex ELEM = 0; 
    static constexpr localIndex FACE = 1; 
  };

  BoundaryStencil();

  virtual void add( localIndex const numPts,
                    localIndex const * const elementRegionIndices,
                    localIndex const * const elementSubRegionIndices,
                    localIndex const * const elementIndices,
                    real64 const * const weights,
                    localIndex const connectorIndex ) override final;

  virtual localIndex size() const override final
  { return m_elementRegionIndices.size( 0 ); }

  constexpr localIndex stencilSize( localIndex const index ) const
  {
    GEOSX_UNUSED_VAR( index )
    return MAX_STENCIL_SIZE;
  }

};

}

#endif //GEOSX_FINITEVOLUME_BOUNDARYSTENCIL_HPP_