StencilDataCollection.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 SRC_CORECOMPONENTS_PHYSICSSOLVERS_FLUIDFLOW_CELLTOCELLDATACOLLECTOR_HPP_
#define SRC_CORECOMPONENTS_PHYSICSSOLVERS_FLUIDFLOW_CELLTOCELLDATACOLLECTOR_HPP_
 
#include "events/tasks/TaskBase.hpp"
#include "mesh/ElementRegionManager.hpp"
namespace geos
{
 
class FlowSolverBase;
class MeshLevel;
class FluxApproximationBase;
 
class StencilDataCollection : public TaskBase
{
public:
 
  StencilDataCollection( const string & name,
                         dataRepository::Group * const parent );
 
  static string catalogName() { return "CellToCellDataCollection"; }
 
  virtual bool execute( real64 const time_n,
                        real64 const dt,
                        integer const cycleNumber,
                        integer const eventCounter,
                        real64 const eventProgress,
                        DomainPartition & domain ) override;
 
  struct viewKeyStruct
  {
    static constexpr char const * solverNameString() { return "flowSolverName"; }
    static constexpr char const * meshNameString() { return "meshBody"; }
    static constexpr char const * cellAGlobalIdString() { return "cellAGlobalId"; }
    static constexpr char const * cellBGlobalIdString() { return "cellBGlobalId"; }
    static constexpr char const * transmissibilityABString() { return "transmissibilityAB"; }
    static constexpr char const * transmissibilityBAString() { return "transmissibilityBA"; }
  };
 
  class Kernel;
 
  struct KernelConnectionData
  {
    real64 m_transmissibility[2];
    localIndex m_regionId[2];
    localIndex m_subRegionId[2];
    localIndex m_elementId[2];
  };
 
private:
 
  using Base = TaskBase;
  using LocalToGlobalMap = ElementRegionManager::ElementViewAccessor< arrayView1d< globalIndex const > >;
 
  struct ConnectionData
  {
    real64 m_transmissibility[2];
    globalIndex m_globalId[2];
 
 
    static ConnectionData fromKernel( KernelConnectionData const & kernelData,
                                      LocalToGlobalMap const & localToGlobalMap );
 
    bool operator<( ConnectionData const & other ) const
    {
      return m_globalId[0] != other.m_globalId[0] ?
             m_globalId[0] < other.m_globalId[0] :
             m_globalId[1] < other.m_globalId[1];
    }
  };
 
  array1d< globalIndex > m_cellAGlobalId;
  array1d< globalIndex > m_cellBGlobalId;
  array1d< real64 > m_transmissibilityAB;
  array1d< real64 > m_transmissibilityBA;
 
  string m_meshName;
  string m_solverName;
 
  FlowSolverBase const * m_solver = nullptr;
  MeshLevel const * m_meshLevel = nullptr;
  FluxApproximationBase const * m_discretization = nullptr;
 
 
 
  void postInputInitialization() override;
 
  void initializePostInitialConditionsPostSubGroups() override;
 
  template< typename STENCILWRAPPER_T >
  array1d< KernelConnectionData > gatherConnectionData( STENCILWRAPPER_T const & stencilWrapper ) const;
 
  void storeConnectionData( string_view stencilName,
                            arrayView1d< KernelConnectionData > const & kernelData );
 
  void logStoredConnections( string_view stencilName );
};
 
 
} /* namespace geos */
 
#endif /* SRC_CORECOMPONENTS_PHYSICSSOLVERS_FLUIDFLOW_CELLTOCELLDATACOLLECTOR_HPP_ */