AccessorUnstructured< SolverType > Class Template Reference

#include <fvstg.h>

Inheritance diagram for AccessorUnstructured< SolverType >:

Collaboration diagram for AccessorUnstructured< SolverType >:

Public Member Functions
	AccessorUnstructured (MInt *sortedBndryCellIds, MInt size, SolverType *solver, const MPI_Comm &commStg, MBool cutOff)
MInt	getNghbrMapping (MInt stgId, MInt nghbr) const
MInt	getGhostIdFromStgId (MInt stgId) const
nDim_citerator	iterateSlopes ()
	Solver specific implementation of the interface defined in base class. More...
nDim_citerator	iterateSlopes_nDim_citerator_end ()
nDim_citerator	iterateB1 ()
nDim_citerator	iterateB1_nDim_citerator_end ()
nDim_citerator	iterateAll ()
nDim_citerator	iterateAll_nDim_citerator_end ()
MInt	getCellId (const nDim_citerator &it_) const
	Helper functions for the nDim_citerator as declared in base class. More...
MInt	getStgId (const nDim_citerator &it_) const
MInt	getNghbr (const nDim_citerator &it_, MInt dir) const
MFloat	a_coordinate (MInt cellId, MInt dim) const override
MInt	domainId () const
MInt	domainId (MInt, MInt) const override
MFloat &	a_pvariable (MInt cellId, MInt varId) override
MFloat &	UInfinity () const override
MFloat &	PInfinity () const override
MFloat &	rhoInfinity () const override
Public Member Functions inherited from Accessor< AccessorUnstructured< SolverType >, SolverType >
AccessorUnstructured< SolverType > *	d ()
virtual MFloat	a_coordinate (MInt, MInt) const=0
virtual MInt	domainId (MInt, MInt) const=0
virtual MFloat &	a_pvariable (MInt, MInt)=0
virtual MFloat &	UInfinity () const=0
virtual MFloat &	PInfinity () const=0
virtual MFloat &	rhoInfinity () const=0
MInt	sizeBC () const
MInt	sizeStg () const
nDim_citerator	iterateSlopes ()
nDim_citerator	iterateSlopes_nDim_citerator_end ()
nDim_citerator	iterateB1 ()
nDim_citerator	iterateB1_nDim_citerator_end ()
nDim_citerator	iterateAll ()
nDim_citerator	iterateAll_nDim_citerator_end ()
MInt	getCellId (const nDim_citerator &it_) const
MInt	getStgId (const nDim_citerator &it_) const
MInt	getNghbr (const nDim_citerator &it_, MInt dir) const
MInt	getNghbrStg (const nDim_citerator &it_, MInt dir) const

Public Attributes
std::vector< MInt >	m_stgBndry2stgGhost
std::vector< MInt >	m_stgGhost2stgBndry
Public Attributes inherited from Accessor< AccessorUnstructured< SolverType >, SolverType >
SolverType *const	m_solver
MPI_Comm	m_commStg
Storage	m_stgToCellId
Storage	m_bcStgId
Storage	m_stgId
MInt	m_noBcCells
MInt	m_stgSize

Private Attributes
MInt **	m_nghbrMapping = nullptr

Additional Inherited Members
Public Types inherited from Accessor< AccessorUnstructured< SolverType >, SolverType >
using	nDim_citerator = nDim_iterator_t< std::vector< MInt >::const_iterator, SolverType >
Static Public Attributes inherited from Accessor< AccessorUnstructured< SolverType >, SolverType >
static constexpr const MInt	m_nDim
Protected Types inherited from Accessor< AccessorUnstructured< SolverType >, SolverType >
using	Storage = std::vector< MInt >
using	iterator = std::vector< MInt >::iterator
using	const_iterator = std::vector< MInt >::const_iterator

Detailed Description

template<class SolverType>
class AccessorUnstructured< SolverType >

Definition at line 171 of file fvstg.h.

Constructor & Destructor Documentation

AccessorUnstructured()

template<class SolverType >

AccessorUnstructured< SolverType >::AccessorUnstructured ( MInt *  sortedBndryCellIds, MInt  size, SolverType *  solver, const MPI_Comm &  commStg, MBool  cutOff  )

inlineexplicit

Definition at line 182 of file fvstg.h.

    : Accessor<AccessorUnstructured<SolverType>, SolverType>(solver, commStg) {
    MInt myrank = 0;
    MPI_Comm_rank(commStg, &myrank);
 
    // First put the boundary cells into m_stgToCellId and then append the ghost cells and finally
    // the reconstruction neighbors, which haven't been appended already
    m_stgToCellId.assign(sortedBndryCellIds, sortedBndryCellIds + size);
    m_noBcCells = size;
    assert(m_stgToCellId.size() == static_cast<MUint>(m_noBcCells));
    m_stgSize = m_noBcCells;
 
    if(!cutOff) {
      // append ghost cells and create map from stgGhost to stgBndry
      m_stgBndry2stgGhost.resize(m_noBcCells);
      m_stgGhost2stgBndry.assign(2 * m_noBcCells, -1);
      m_stgToCellId.resize(2 * m_noBcCells);
      for(MInt id = 0; id < m_noBcCells; ++id) {
        const MInt cellId = m_stgToCellId[id];
        const MInt bndryId = solver->a_bndryId(cellId);
 
        if(bndryId < 0) {
          TERMM(1, "errror");
        }
 
        const MInt ghostCellId = m_solver->m_bndryCells->a[bndryId].m_srfcVariables[0]->m_ghostCellId;
        m_stgToCellId[m_noBcCells + id] = ghostCellId;
        m_stgBndry2stgGhost[id] = m_noBcCells + id;
        m_stgGhost2stgBndry[m_noBcCells + id] = id;
      }
    }
 
    m_stgSize = m_stgToCellId.size();
 
    mAlloc(m_nghbrMapping, m_noBcCells, m_solver->m_cells.noRecNghbrs(), "m_nghbrMapping_", AT_);
 
    if(!cutOff) {
      // Check if cellId has been already assigned a stgId
      using MyMap = std::map<MInt, MInt>;
      MyMap cellIdToStgId_map;
      auto cellIdToStgId = [this, &cellIdToStgId_map](MInt cellId) -> MInt {
        std::pair<MyMap::const_iterator, MBool> res = cellIdToStgId_map.insert(MyMap::value_type(cellId, m_stgSize));
        if(res.second) ++m_stgSize;
        return res.first->second;
      };
      for(MInt i = 0; i < m_stgSize; ++i) {
        cellIdToStgId_map.insert(MyMap::value_type(m_stgToCellId[i], i));
      }
 
      for(MInt i = 0; i < m_noBcCells; ++i) {
        for(MInt nghbr = 0; nghbr < m_solver->a_noReconstructionNeighbors(m_stgToCellId[i]); ++nghbr) {
          const MInt nghbrId = m_solver->a_reconstructionNeighborId(m_stgToCellId[i], nghbr);
          const MInt stgId = cellIdToStgId(nghbrId);
 
          m_nghbrMapping[i][nghbr] = stgId;
 
          if(m_stgToCellId.size() <= static_cast<MUint>(stgId)) {
            /* if(!cutOff){ */
            if(m_solver->a_isBndryGhostCell(nghbrId)) {
              std::cout << "NOO My rank " << myrank << "; " << nghbrId << "; " << stgId << "; "
                        << solver->a_coordinate(nghbrId, 0) << "|" << solver->a_coordinate(nghbrId, 1) << "|"
                        << solver->a_coordinate(nghbrId, 2) << std::endl;
              TERMM(1, "ERROR");
            }
            /* } */
            m_stgToCellId.push_back(nghbrId);
          }
        }
      }
 
      if(m_stgToCellId.size() != static_cast<MUint>(m_stgSize)) {
        TERMM(1, "ERROR");
      }
 
      /* if(!cutOff)  */ m_stgGhost2stgBndry.resize(m_stgSize, -1);
    }
 
    // Create list of stg Ids of only the boundary cells, which are stored at the beginning for unstructured case
    m_bcStgId.resize(m_noBcCells);
    std::iota(std::begin(m_bcStgId), std::end(m_bcStgId), 0); // Fill with 0, 1, ..., m_noBcCells-1
 
    // Create list of stg Ids of all cells
    m_stgId.resize(m_stgSize);
    std::iota(std::begin(m_stgId), std::end(m_stgId), 0); // Fill with 0, 1, ..., m_stgSize-1
 
 
    // Print summary
    /* #ifndef NDEBUG */
    /*     std::cout << "--- STG INFO ---" << std::endl */
    /*               << myrank << "# bcCells/stgCells: " << m_noBcCells << " / " << m_stgSize << std::endl; */
    /* #endif */
  } // Constructor ends

Member Function Documentation

a_coordinate()

template<class SolverType >

MFloat AccessorUnstructured< SolverType >::a_coordinate ( MInt  cellId, MInt  dim  ) const

inlineoverridevirtual

Implements Accessor< AccessorUnstructured< SolverType >, SolverType >.

Definition at line 309 of file fvstg.h.

{ return m_solver->a_coordinate(cellId, dim); }

a_pvariable()

template<class SolverType >

MFloat & AccessorUnstructured< SolverType >::a_pvariable ( MInt  cellId, MInt  varId  )

inlineoverridevirtual

Implements Accessor< AccessorUnstructured< SolverType >, SolverType >.

Definition at line 312 of file fvstg.h.

{ return m_solver->a_pvariable(cellId, varId); }

domainId() [1/2]

template<class SolverType >

MInt AccessorUnstructured< SolverType >::domainId ( ) const

inline

Definition at line 310 of file fvstg.h.

{ return m_solver->domainId(); }

domainId() [2/2]

template<class SolverType >

MInt AccessorUnstructured< SolverType >::domainId ( MInt  , MInt    ) const

inlineoverridevirtual

Implements Accessor< AccessorUnstructured< SolverType >, SolverType >.

Definition at line 311 of file fvstg.h.

{ TERMM(-1, "Invalid call!"); }

getCellId()

template<class SolverType >

MInt AccessorUnstructured< SolverType >::getCellId ( const nDim_citerator &  it_ ) const

inline

Definition at line 297 of file fvstg.h.

{ return m_stgToCellId[*it_]; }

getGhostIdFromStgId()

template<class SolverType >

MInt AccessorUnstructured< SolverType >::getGhostIdFromStgId ( MInt  stgId ) const

inline

Definition at line 279 of file fvstg.h.

                                             {
    const MInt cellId = m_stgToCellId[stgId];
    const MInt bndryId = m_solver->a_bndryId(cellId);
    // TODO labels:FV,totest check if it is a boundary cell
    return m_solver->m_bndryCells->a[bndryId].m_srfcVariables[0]->m_ghostCellId;
  }

getNghbr()

template<class SolverType >

MInt AccessorUnstructured< SolverType >::getNghbr ( const nDim_citerator &  it_, MInt  dir  ) const

inline

Definition at line 299 of file fvstg.h.

                                                           {
    const MInt cellId = m_stgToCellId[*it_];
    if(m_solver->a_hasNeighbor(cellId, dir) > 0) {
      return m_solver->c_neighborId(cellId, dir);
    } else {
      return -1;
    }
  }

getNghbrMapping()

template<class SolverType >

MInt AccessorUnstructured< SolverType >::getNghbrMapping ( MInt  stgId, MInt  nghbr  ) const

inline

Definition at line 277 of file fvstg.h.

{ return m_nghbrMapping[stgId][nghbr]; }

getStgId()

template<class SolverType >

MInt AccessorUnstructured< SolverType >::getStgId ( const nDim_citerator &  it_ ) const

inline

Definition at line 298 of file fvstg.h.

{ return *it_; }

iterateAll()

template<class SolverType >

nDim_citerator AccessorUnstructured< SolverType >::iterateAll ( )

inline

Definition at line 293 of file fvstg.h.

{ return {m_stgId.begin(), this}; }

iterateAll_nDim_citerator_end()

template<class SolverType >

nDim_citerator AccessorUnstructured< SolverType >::iterateAll_nDim_citerator_end ( )

inline

Definition at line 294 of file fvstg.h.

{ return {m_stgId.end(), this}; }

iterateB1()

template<class SolverType >

nDim_citerator AccessorUnstructured< SolverType >::iterateB1 ( )

inline

Definition at line 290 of file fvstg.h.

{ return {m_bcStgId.begin(), this}; }

iterateB1_nDim_citerator_end()

template<class SolverType >

nDim_citerator AccessorUnstructured< SolverType >::iterateB1_nDim_citerator_end ( )

inline

Definition at line 291 of file fvstg.h.

{ return {m_bcStgId.end(), this}; }

iterateSlopes()

template<class SolverType >

nDim_citerator AccessorUnstructured< SolverType >::iterateSlopes ( )

inline

Definition at line 287 of file fvstg.h.

{ return {m_bcStgId.begin(), this}; }

iterateSlopes_nDim_citerator_end()

template<class SolverType >

nDim_citerator AccessorUnstructured< SolverType >::iterateSlopes_nDim_citerator_end ( )

inline

Definition at line 288 of file fvstg.h.

{ return {m_bcStgId.end(), this}; }

PInfinity()

template<class SolverType >

MFloat & AccessorUnstructured< SolverType >::PInfinity ( ) const

inlineoverridevirtual

Implements Accessor< AccessorUnstructured< SolverType >, SolverType >.

Definition at line 314 of file fvstg.h.

{ return m_solver->m_PInfinity; }

rhoInfinity()

template<class SolverType >

MFloat & AccessorUnstructured< SolverType >::rhoInfinity ( ) const

inlineoverridevirtual

Implements Accessor< AccessorUnstructured< SolverType >, SolverType >.

Definition at line 315 of file fvstg.h.

{ return m_solver->m_rhoInfinity; }

UInfinity()

template<class SolverType >

MFloat & AccessorUnstructured< SolverType >::UInfinity ( ) const

inlineoverridevirtual

Implements Accessor< AccessorUnstructured< SolverType >, SolverType >.

Definition at line 313 of file fvstg.h.

{ return m_solver->m_UInfinity; }

Member Data Documentation

m_nghbrMapping

template<class SolverType >

MInt** AccessorUnstructured< SolverType >::m_nghbrMapping = nullptr

private

Definition at line 318 of file fvstg.h.

m_stgBndry2stgGhost

template<class SolverType >

std::vector<MInt> AccessorUnstructured< SolverType >::m_stgBndry2stgGhost

Definition at line 322 of file fvstg.h.

m_stgGhost2stgBndry

template<class SolverType >

std::vector<MInt> AccessorUnstructured< SolverType >::m_stgGhost2stgBndry

Definition at line 323 of file fvstg.h.

---

The documentation for this class was generated from the following file:

• /home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/FV/fvstg.h