#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.

309{ 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.

312{ 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.

310{ 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.

311{ 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.

297{ 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.

277{ 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.

298{ return *it_; }

◆ iterateAll()

template<class SolverType >

nDim_citerator AccessorUnstructured< SolverType >::iterateAll ( )

inline

Definition at line 293 of file fvstg.h.

293{ 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.

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

◆ iterateB1()

template<class SolverType >

nDim_citerator AccessorUnstructured< SolverType >::iterateB1 ( )

inline

Definition at line 290 of file fvstg.h.

290{ 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.

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

◆ iterateSlopes()

template<class SolverType >

nDim_citerator AccessorUnstructured< SolverType >::iterateSlopes ( )

inline

Definition at line 287 of file fvstg.h.

287{ 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.

288{ 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.

314{ 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.

315{ 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.

313{ 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

Public Member Functions

Public Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ AccessorUnstructured()

Member Function Documentation

◆ a_coordinate()

◆ a_pvariable()

◆ domainId() [1/2]

◆ domainId() [2/2]

◆ getCellId()

◆ getGhostIdFromStgId()

◆ getNghbr()

◆ getNghbrMapping()

◆ getStgId()

◆ iterateAll()

◆ iterateAll_nDim_citerator_end()

◆ iterateB1()

◆ iterateB1_nDim_citerator_end()

◆ iterateSlopes()

◆ iterateSlopes_nDim_citerator_end()

◆ PInfinity()

◆ rhoInfinity()

◆ UInfinity()

Member Data Documentation

◆ m_nghbrMapping

◆ m_stgBndry2stgGhost

◆ m_stgGhost2stgBndry