QuasiRTInnerProduct.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 TotalEnergies
 * Copyright (c) 2019-     GEOSX Contributors
 * All rights reserved
 *
 * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
 * ------------------------------------------------------------------------------------------------------------
 */
 
#ifndef GEOS_FINITEVOLUME_MIMETICINNERPRODUCTS_QUASIRTINNERPRODUCT_HPP_
#define GEOS_FINITEVOLUME_MIMETICINNERPRODUCTS_QUASIRTINNERPRODUCT_HPP_
 
#include "finiteVolume/mimeticInnerProducts/MimeticInnerProductBase.hpp"
 
namespace geos
{
namespace mimeticInnerProduct
{
 
class QuasiRTInnerProduct : public MimeticInnerProductBase
{
public:
 
  template< localIndex NF >
  GEOS_HOST_DEVICE
  static void
  compute( arrayView2d< real64 const, nodes::REFERENCE_POSITION_USD > const & nodePosition,
           arrayView1d< real64 const > const & transMultiplier,
           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,
           arraySlice2d< real64 > const & transMatrix );
 
};
 
template< localIndex NF >
GEOS_HOST_DEVICE
void
QuasiRTInnerProduct::compute( arrayView2d< real64 const, nodes::REFERENCE_POSITION_USD > const & nodePosition,
                              arrayView1d< real64 const > const & transMultiplier,
                              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,
                              arraySlice2d< real64 > const & transMatrix )
{
  MimeticInnerProductBase::computeParametricInnerProduct< NF >( nodePosition,
                                                                transMultiplier,
                                                                faceToNodes,
                                                                elemToFaces,
                                                                elemCenter,
                                                                elemVolume,
                                                                elemPerm,
                                                                6.0,
                                                                lengthTolerance,
                                                                transMatrix );
}
 
} // end namespace mimeticInnerProduct
 
} // end namespace geos
 
#endif //GEOS_FINITEVOLUME_MIMETICINNERPRODUCTS_QUASIRTINNERPRODUCT_HPP_