AdaptiveInnerProduct.hpp

Go to the documentation of this file.

/*
 * ------------------------------------------------------------------------------------------------------------
 * SPDX-License-Identifier: LGPL-2.1-only
 *
 * Copyright (c) 2016-2024 Lawrence Livermore National Security LLC
 * Copyright (c) 2018-2024 TotalEnergies
 * Copyright (c) 2018-2024 The Board of Trustees of the Leland Stanford Junior University
 * Copyright (c) 2023-2024 Chevron
 * Copyright (c) 2019-     GEOS/GEOSX Contributors
 * All rights reserved
 *
 * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
 * ------------------------------------------------------------------------------------------------------------
 */
 
#ifndef GEOS_FINITEVOLUME_MIMETICINNERPRODUCTS_ADAPTIVEINNERPRODUCT_HPP_
#define GEOS_FINITEVOLUME_MIMETICINNERPRODUCTS_ADAPTIVEINNERPRODUCT_HPP_
 
#include "finiteVolume/mimeticInnerProducts/TPFAInnerProduct.hpp"
 
namespace geos
{
namespace mimeticInnerProduct
{
 
template< typename MFD_IP >
struct AdaptiveInnerProduct
{
  template< localIndex NF >
  GEOS_HOST_DEVICE
  static void
  computeM( arrayView2d< real64 const, nodes::REFERENCE_POSITION_USD > const & nodePosition,
            ArrayOfArraysView< localIndex const > const & faceToNodes,
            arraySlice1d< localIndex const > const & elemToFaces,
            arraySlice1d< real64 const > const & elemCenter,
            real64 const & elemVolume,
            real64 const (&elemPerm)[ 3 ],
            real64 const & lengthTolerance,
            real64 const & weight,
            arraySlice2d< real64 > const & M )
  {
    stackArray2d< real64, NF *NF > mfdMassMatrix( NF, NF );
    stackArray2d< real64, NF *NF > tpfaMassMatrix( NF, NF );
 
    MFD_IP::template computeM< NF >( nodePosition,
                                     faceToNodes,
                                     elemToFaces,
                                     elemCenter,
                                     elemVolume,
                                     elemPerm,
                                     lengthTolerance,
                                     mfdMassMatrix );
 
    TPFAInnerProduct::computeM< NF >( nodePosition,
                                      faceToNodes,
                                      elemToFaces,
                                      elemCenter,
                                      elemVolume,
                                      elemPerm,
                                      lengthTolerance,
                                      tpfaMassMatrix );
 
    for( localIndex i = 0; i < NF; ++i )
    {
      for( localIndex j = 0; j < NF; ++j )
      {
        M[i][j] = weight * mfdMassMatrix( i, j ) + ( 1.0 - weight ) * tpfaMassMatrix( i, j );
      }
    }
  }
};
 
} // end namespace mimeticInnerProduct
 
} // end namespace geos
 
#endif //GEOS_FINITEVOLUME_MIMETICINNERPRODUCTS_ADAPTIVEINNERPRODUCT_HPP_