#include <fvmbzonal.h>

Inheritance diagram for CouplerFvMbZonal< nDim, SysEqn >:

Collaboration diagram for CouplerFvMbZonal< nDim, SysEqn >:

Public Types
using	FvMbSolver = FvMbCartesianSolverXD< nDim, SysEqn >

using	FVBndryCnd = FvBndryCndXD< nDim, SysEqn >

Public Types inherited from CouplingFvMb< nDim, SysEqn >
using	solverType = FvMbCartesianSolverXD< nDim, SysEqn >

Public Member Functions
	CouplerFvMbZonal (const MInt couplingId, FvMbSolver upStream, FvMbSolver downStream)

	~CouplerFvMbZonal ()

void	init () override

void	finalizeSubCoupleInit (MInt)

void	postCouple (MInt) override
	unify the time-Step More...

void	cleanUp ()

void	checkProperties () override

void	readProperties () override

void	finalizeCouplerInit ()
	create first zonal mapping for the zonal exchange More...

void	finalizeBalance (const MInt) override
	update zonal mapping after balance More...

void	finalizeAdaptation (const MInt) override
	update zonal mapping after adaptation More...

void	subCouple (MInt, MInt, std::vector< MBool > &)
	exchange zonal variables each runge-Kutta Step More...

void	preCouple (MInt) override
	exchange zonal variables before each runge-Kutta Step More...

Public Member Functions inherited from CouplingFvMb< nDim, SysEqn >
	CouplingFvMb (const MInt couplingId, solverType *b)

solverType &	fvMbSolver () const

MInt	a_noFvCells () const

MInt	a_noFvGridCells () const

MInt	a_noLevelSetsMb () const

Public Member Functions inherited from Coupling
	Coupling (const MInt couplingId)

virtual	~Coupling ()=default

	Coupling (const Coupling &)=delete

Coupling &	operator= (const Coupling &)=delete

MInt	couplerId () const

virtual void	init ()=0

virtual void	finalizeSubCoupleInit (MInt solverId)=0

virtual void	finalizeCouplerInit ()=0

virtual void	preCouple (MInt recipeStep)=0

virtual void	subCouple (MInt recipeStep, MInt solverId, std::vector< MBool > &solverCompleted)=0

virtual void	postCouple (MInt recipeStep)=0

virtual void	cleanUp ()=0

virtual void	balancePre ()
	Load balancing. More...

virtual void	balancePost ()

virtual void	reinitAfterBalance ()

virtual void	prepareAdaptation ()

virtual void	postAdaptation ()

virtual void	finalizeAdaptation (const MInt)

virtual void	writeRestartFile (const MInt)

virtual MInt	noCellDataDlb () const
	Methods to inquire coupler data during balancing. More...

virtual MInt	cellDataTypeDlb (const MInt NotUsed(dataId)) const

virtual MInt	cellDataSizeDlb (const MInt NotUsed(dataId), const MInt NotUsed(cellId))

virtual void	getCellDataDlb (const MInt NotUsed(dataId), const MInt NotUsed(oldNoCells), const MInt const NotUsed(bufferIdToCellId), MInt const NotUsed(data))

virtual void	getCellDataDlb (const MInt NotUsed(dataId), const MInt NotUsed(oldNoCells), const MInt const NotUsed(bufferIdToCellId), MLong const NotUsed(data))

virtual void	getCellDataDlb (const MInt NotUsed(dataId), const MInt NotUsed(oldNoCells), const MInt const NotUsed(bufferIdToCellId), MFloat const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MInt *const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MLong *const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MFloat *const NotUsed(data))

virtual void	finalizeBalance (const MInt)

virtual MInt	noCouplingTimers (const MBool NotUsed(allTimings)) const
	Number of coupling timers. More...

virtual void	getCouplingTimings (std::vector< std::pair< MString, MFloat > > &NotUsed(timings), const MBool NotUsed(allTimings))
	Return coupling timings. More...

virtual void	getDomainDecompositionInformation (std::vector< std::pair< MString, MInt > > &NotUsed(domainInfo))
	Return information on current domain decomposition (e.g. number of coupled cells/elements/...) More...

void	setDlbTimer (const MInt timerId)

void	startLoadTimer (const MString &name) const
	Start the load timer of the coupler. More...

void	stopLoadTimer (const MString &name) const
	Stop the load timer of the coupler. More...

Private Member Functions
FvMbSolver &	upStream () const

FvMbSolver &	downStream () const

MInt	couplerId () const

void	createZonalMapping ()
	create the cell mapping which allows for a faster exchange during each RK-Step. This requires little additional memory. More...

void	updateZonalMapping ()
	update the zonal mapping before the FVMb-timeStep as the bndryCell order might have changed! NOTE: currently unused, as its not necessary, the bndry-ghostCell values are updated inside the solutionStep. More...

void	exchangeZonalValues (const MInt)
	exchange zonal variables around the zonal coordinate from upstream to downstream solver and the other way around! More...

void	unifyTimeStep ()
	set the same timeStep both FvMb-solvers More...

MInt	up2downId (MInt, const MBool)
	conversion from upStream solverId to the downStream solverId NOTE: also handles bndry-ghost cell conversion! More...

MInt	down2upId (MInt, const MBool)
	conversion from downStream to upStream cellId NOTE: also handles bndry-ghost cells More...

Private Attributes
FvMbSolver *	m_solverUp

FvMbSolver *	m_solverDown

MInt	m_upStreamId {}

MInt	m_downStreamId {}

MInt	m_couplingId {}

MString	m_zonalMethod {}

MInt	m_zonalDir {}

MFloat	m_zonalCoordinate {}

MBool	m_zonalDualTimeStepping = false

std::vector< std::pair< MInt, MInt > >	m_upDown

std::vector< std::pair< MInt, MInt > >	m_downUp

MInt	m_upStreamOffset = 0

MInt	m_downStreamOffset = 0

MInt	m_RKStep = 0

MInt	m_noRKSteps = -1

Additional Inherited Members
Protected Member Functions inherited from Coupling
MFloat	returnLoadRecord () const

MFloat	returnIdleRecord () const

Detailed Description

template<MInt nDim, class SysEqn>
class CouplerFvMbZonal< nDim, SysEqn >

Definition at line 31 of file fvmbzonal.h.

Member Typedef Documentation

◆ FVBndryCnd

template<MInt nDim, class SysEqn >

using CouplerFvMbZonal< nDim, SysEqn >::FVBndryCnd = FvBndryCndXD<nDim, SysEqn>

Definition at line 34 of file fvmbzonal.h.

◆ FvMbSolver

template<MInt nDim, class SysEqn >

using CouplerFvMbZonal< nDim, SysEqn >::FvMbSolver = FvMbCartesianSolverXD<nDim, SysEqn>

Definition at line 33 of file fvmbzonal.h.

Constructor & Destructor Documentation

◆ CouplerFvMbZonal()

template<MInt nDim, class SysEqn >

CouplerFvMbZonal< nDim, SysEqn >::CouplerFvMbZonal	(	const MInt	couplingId,
		FvMbSolver *	upStream,
		FvMbSolver *	downStream
	)

Definition at line 23 of file fvmbzonal.cpp.

  : Coupling(couplingId), CouplingFvMb<nDim, SysEqn>(couplingId, upStream_) {
  TRACE();
 
  // set solver pointer and ids
  m_solverUp = upStream_;
  m_solverDown = downStream_;
  m_couplingId = couplingId;
 
  m_upStreamId = upStream().m_solverId;
  m_downStreamId = downStream().m_solverId;
 
  readProperties();
}

◆ ~CouplerFvMbZonal()

template<MInt nDim, class SysEqn >

CouplerFvMbZonal< nDim, SysEqn >::~CouplerFvMbZonal ( )

inline

Definition at line 38 of file fvmbzonal.h.

38{};

Member Function Documentation

◆ checkProperties()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::checkProperties ( )

inlineoverridevirtual

Implements Coupling.

Definition at line 45 of file fvmbzonal.h.

45{};

◆ cleanUp()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::cleanUp ( )

inlinevirtual

Reimplemented from CouplingFvMb< nDim, SysEqn >.

Definition at line 44 of file fvmbzonal.h.

44{};

◆ couplerId()

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::couplerId ( ) const

inlineprivate

Definition at line 60 of file fvmbzonal.h.

60{ return m_couplingId; }

◆ createZonalMapping()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::createZonalMapping

private

Author: Tim Wegmann

Definition at line 201 of file fvmbzonal.cpp.

                                                        {
  TRACE();
 
  m_upDown.clear();
  m_downUp.clear();
 
  m_upStreamOffset = 0;
  m_downStreamOffset = 0;
 
  if(upStream().isActive()) {
    if(m_zonalMethod == "P") {
      // first create mapping for valid grid cells
      // this mapping will only change duing adaptation or balance
      for(MInt upId = 0; upId < upStream().c_noCells(); upId++) {
        if(upStream().a_coordinate(upId, m_zonalDir) > m_zonalCoordinate) {
          MBool lastLayer = false;
          if(!upStream().a_hasNeighbor(upId, m_zonalDir * 2 + 1, false) && !upStream().a_isHalo(upId)) {
            const MInt parentId = upStream().c_parentId(upId);
            if(parentId > -1 && !upStream().a_hasNeighbor(parentId, m_zonalDir * 2 + 1, false)) {
              lastLayer = true;
            }
          }
          const MInt downId = up2downId(upId, lastLayer);
          if(downId > -1) {
            m_downUp.push_back(make_pair(downId, upId));
          }
        }
      }
    } else if(m_zonalMethod == "B") {
    }
    m_downStreamOffset = m_downUp.size();
  }
 
  if(downStream().isActive()) {
    if(m_zonalMethod == "P") {
      // first create mapping for valid grid cells
      // this mapping will only change duing adaptation or balance
      for(MInt downId = 0; downId < downStream().c_noCells(); downId++) {
        if(downStream().a_coordinate(downId, m_zonalDir) < m_zonalCoordinate) {
          MBool lastLayer = false;
          if(!downStream().a_hasNeighbor(downId, m_zonalDir * 2, false) && !downStream().a_isHalo(downId)) {
            const MInt parentId = downStream().c_parentId(downId);
            if(parentId > -1 && !downStream().a_hasNeighbor(parentId, m_zonalDir * 2, false)) {
              lastLayer = true;
            }
          }
          const MInt upId = down2upId(downId, lastLayer);
          if(upId > -1) {
            m_upDown.push_back(make_pair(upId, downId));
          }
        }
      }
    } else if(m_zonalMethod == "B") {
    }
    m_upStreamOffset = m_upDown.size();
  }
}

◆ down2upId()

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::down2upId	(	MInt	downCell,
		const	MBool
	)

private

Author: Tim Wegmann

Definition at line 438 of file fvmbzonal.cpp.

                                                                             {
  TRACE();
 
  if(!upStream().isActive()) return -1;
 
  // to hanlde bndry ghostCells and others as well
  //   1: regular grid cell
  // < 1: bndry-ghost cell (index represents matching surface)
  MInt cellType = 1;
  if(downCell >= downStream().c_noCells()) {
    const MInt gridCell = downStream().getAssociatedInternalCell(downCell);
    if(downStream().a_isBndryGhostCell(downCell)) {
      // find the matching surface of the grid cell to which the upCell is the bndry-ghostcell
      const MInt bndryId = downStream().a_bndryId(gridCell);
      MInt srfc = -1;
      for(MInt surf = 0; surf < downStream().m_bndryCells->a[bndryId].m_noSrfcs; surf++) {
        if(downStream().m_bndryCells->a[bndryId].m_srfcVariables[surf]->m_ghostCellId == downCell) {
          srfc = surf;
          break;
        }
      }
      ASSERT(srfc >= 0, "");
      cellType = -srfc;
    } else if(downStream().a_hasProperty(downCell, FvCell::IsSplitChild)
              || downStream().a_hasProperty(downCell, FvCell::IsSplitClone)) {
      cellType = 2;
    }
    downCell = gridCell;
  }
 
  downStream().assertValidGridCellId(downCell);
 
  const MInt gridId = downStream().grid().tree().solver2grid(downCell);
  ASSERT(downStream().grid().solverFlag(gridId, m_downStreamId), "");
 
  if(!upStream().grid().solverFlag(gridId, m_upStreamId)) return -1;
 
  MInt upCell = upStream().grid().tree().grid2solver(gridId);
 
  if(upCell > 0) {
    ASSERT(downStream().a_level(downCell) == upStream().a_level(upCell), "");
  }
 
  if(cellType == 1 || upCell < 0) return upCell;
 
  if(!lastLayer) {
    ASSERT(upStream().a_isBndryCell(upCell), "");
  } else {
    return -1;
  }
 
  if(cellType < 1) {
    // bndry-ghost cell
    const MInt srfc = -cellType;
    const MInt bndryId = upStream().a_bndryId(upCell);
    upCell = upStream().m_bndryCells->a[bndryId].m_srfcVariables[srfc]->m_ghostCellId;
    return upCell;
 
  } else {
    // split child
    mTerm(1, AT_, "Not yet implemented for splitchilds");
  }
 
  return -1;
}

◆ downStream()

template<MInt nDim, class SysEqn >

FvMbSolver & CouplerFvMbZonal< nDim, SysEqn >::downStream ( ) const

inlineprivate

Definition at line 59 of file fvmbzonal.h.

59{ return *m_solverDown; }

◆ exchangeZonalValues()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::exchangeZonalValues ( const MInt solverId )

private

Author: Tim Wegmann

Definition at line 326 of file fvmbzonal.cpp.

                                                                            {
  TRACE();
 
  if(m_upDown.empty()) {
    ASSERT(m_downUp.empty(), "");
    return;
  }
 
  // transfer data from up-stream to down-stream
  if(solverId == m_upStreamId || solverId == -1) {
    for(auto it = m_upDown.begin(); it != m_upDown.end(); it++) {
      const MInt upId = it->first;
      const MInt downId = it->second;
 
      for(MInt var = 0; var < upStream().m_sysEqn.CV->noVariables; var++) {
        downStream().a_variable(downId, var) = upStream().a_variable(upId, var);
        downStream().a_pvariable(downId, var) = upStream().a_pvariable(upId, var);
      }
    }
  }
 
  // transfer data from down-stream to up-stream
  if(solverId == m_downStreamId || solverId == -1) {
    for(auto it = m_downUp.begin(); it != m_downUp.end(); it++) {
      const MInt downId = it->first;
      const MInt upId = it->second;
 
      for(MInt var = 0; var < downStream().m_sysEqn.CV->noVariables; var++) {
        upStream().a_variable(upId, var) = downStream().a_variable(downId, var);
        upStream().a_pvariable(upId, var) = downStream().a_pvariable(downId, var);
      }
    }
  }
}

◆ finalizeAdaptation()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::finalizeAdaptation ( const MInt )

overridevirtual

Author: Tim Wegmann

Reimplemented from Coupling.

Definition at line 123 of file fvmbzonal.cpp.

                                                                         {
  if(solverLoopId >= m_upStreamId && solverLoopId >= m_downStreamId) {
    createZonalMapping();
    unifyTimeStep();
  }
}

◆ finalizeBalance()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::finalizeBalance ( const MInt )

overridevirtual

Author: Tim Wegmann

Reimplemented from Coupling.

Definition at line 134 of file fvmbzonal.cpp.

                                                                      {
  if(solverLoopId >= m_upStreamId && solverLoopId >= m_downStreamId) {
    createZonalMapping();
  }
}

◆ finalizeCouplerInit()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::finalizeCouplerInit

virtual

Author: Tim Wegmann

Implements Coupling.

Definition at line 113 of file fvmbzonal.cpp.

                                                         {
  createZonalMapping();
 
  ASSERT(upStream().a_timeStepComputationInterval() == downStream().a_timeStepComputationInterval(), "");
}

◆ finalizeSubCoupleInit()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::finalizeSubCoupleInit ( MInt )

inlinevirtual

Implements Coupling.

Definition at line 42 of file fvmbzonal.h.

42{};

◆ init()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::init ( )

inlineoverridevirtual

Implements Coupling.

Definition at line 41 of file fvmbzonal.h.

41{};

◆ postCouple()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::postCouple ( MInt solverStep )

overridevirtual

Author: Tim Wegmann

Implements Coupling.

Definition at line 185 of file fvmbzonal.cpp.

                                                               {
  if(solverStep == 1) {
    if(upStream().isActive() && downStream().isActive()
       && fabs(upStream().timeStep() - downStream().timeStep()) > upStream().m_eps) {
      cerr << setprecision(16) << "UpStream-TimeStep " << upStream().timeStep() << " DownStream-TimeStep "
           << upStream().timeStep() << endl;
      mTerm(1, AT_, "Time in Fv-Mb-Solvers differs!.");
    }
  }
}

◆ preCouple()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::preCouple ( MInt solverStep )

overridevirtual

Author: Tim Wegmann

Implements Coupling.

Definition at line 171 of file fvmbzonal.cpp.

                                                              {
  // past the LS-solver
  // before the FvMb-solver
  if(solverStep == 1) {
    // updateZonalMapping();
    exchangeZonalValues(-1);
    m_RKStep = 0;
  }
}

◆ readProperties()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::readProperties ( )

overridevirtual

Implements Coupling.

◆ subCouple()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::subCouple	(	MInt	recipeStep,
		MInt	solverId,
		std::vector< MBool > &
	)

virtual

Author: Tim Wegmann

Implements Coupling.

Definition at line 144 of file fvmbzonal.cpp.

                                                                    {
  // only during FVMb computation
  if(recipeStep == 1 && (solverId == m_upStreamId || solverId == m_downStreamId)) {
    m_RKStep++;
 
    if(m_RKStep == 2 * m_noRKSteps) {
      // complete exchange during the last runge-Kutta Step
      exchangeZonalValues(-1);
    } else {
      // only update from the solver who just made its timeStep!
      exchangeZonalValues(solverId);
    }
 
    // check the timeStep contraint at the last-RK-step of both solvers
    // if the timeStep might have changed!
    if(m_RKStep == 2 * m_noRKSteps && upStream().a_timeStepComputationInterval() > 0
       && globalTimeStep % upStream().a_timeStepComputationInterval() == 0) {
      unifyTimeStep();
    }
  }
}

◆ unifyTimeStep()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::unifyTimeStep

private

Author: Tim Wegmann

Definition at line 509 of file fvmbzonal.cpp.

                                                   {
  const MFloat invalidTimeStep = std::numeric_limits<MFloat>::max();
  MFloat timeUp = invalidTimeStep;
  MFloat timeDown = invalidTimeStep;
  MInt maxLevelUp = -1;
  MInt maxLevelDown = -1;
  MFloat newTimeUp = invalidTimeStep;
  MFloat newTimeDown = invalidTimeStep;
 
  // determine the new timeStep on ranks which have both up-Stream and downStream solvers!
  if(upStream().isActive() && downStream().isActive()) {
    timeUp = upStream().timeStep(true);
    maxLevelUp = upStream().maxLevel();
    timeDown = downStream().timeStep(true);
    maxLevelDown = downStream().maxLevel();
 
    MInt levelDiff = maxLevelDown - maxLevelUp;
    if(levelDiff == 0) {
      newTimeUp = mMin(timeUp, timeDown);
      newTimeDown = newTimeUp;
    } else {
      ASSERT(m_zonalDualTimeStepping, "");
      if(levelDiff == 1) {
        // downStream has higher level => lower timeStep!
        newTimeDown = mMin(newTimeDown, newTimeUp / 2.0);
        newTimeUp = 2.0 * newTimeDown;
      } else if(levelDiff == -1) {
        // upStream has the higher level => lower timeStep!
        newTimeUp = mMin(newTimeUp, newTimeDown / 2.0);
        newTimeDown = 2.0 * newTimeUp;
      } else {
        mTerm(1, AT_, "Zonal dual-timeStepping currently only implemented for 1 level difference!");
      }
    }
  }
 
  if(downStream().isActive()) {
    MPI_Allreduce(MPI_IN_PLACE, &newTimeDown, 1, MPI_DOUBLE, MPI_MIN, downStream().mpiComm(), AT_, "MPI_IN_PLACE",
                  "newTimeDown");
  }
  if(upStream().isActive()) {
    MPI_Allreduce(MPI_IN_PLACE, &newTimeUp, 1, MPI_DOUBLE, MPI_MIN, upStream().mpiComm(), AT_, "MPI_IN_PLACE",
                  "newTimeUp");
  }
 
  if(upStream().isActive()) {
    upStream().forceTimeStep(newTimeUp);
  }
  if(downStream().isActive()) {
    downStream().forceTimeStep(newTimeDown);
  }
}

◆ up2downId()

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::up2downId	(	MInt	upCell,
		const	MBool
	)

private

Author: Tim Wegmann

Definition at line 367 of file fvmbzonal.cpp.

                                                                           {
  TRACE();
 
  if(!downStream().isActive()) return -1;
 
  // to hanlde bndry ghostCells and others as well
  //   1: regular grid cell
  // < 1: bndry-ghost cell (index represents matching surface)
  MInt cellType = 1;
  if(upCell >= upStream().c_noCells()) {
    const MInt gridCell = upStream().getAssociatedInternalCell(upCell);
    if(upStream().a_isBndryGhostCell(upCell)) {
      // find the matching surface of the grid cell to which the upCell is the bndry-ghostcell
      const MInt bndryId = upStream().a_bndryId(gridCell);
      MInt srfc = -1;
      for(MInt surf = 0; surf < upStream().m_bndryCells->a[bndryId].m_noSrfcs; surf++) {
        if(upStream().m_bndryCells->a[bndryId].m_srfcVariables[surf]->m_ghostCellId == upCell) {
          srfc = surf;
          break;
        }
      }
      ASSERT(srfc >= 0, "");
      cellType = -srfc;
    } else if(upStream().a_hasProperty(upCell, FvCell::IsSplitChild)
              || upStream().a_hasProperty(upCell, FvCell::IsSplitClone)) {
      cellType = 2;
    }
    upCell = gridCell;
  }
 
  upStream().assertValidGridCellId(upCell);
 
  const MInt gridId = upStream().grid().tree().solver2grid(upCell);
  ASSERT(upStream().grid().solverFlag(gridId, m_upStreamId), "");
 
  if(!downStream().grid().solverFlag(gridId, m_downStreamId)) return -1;
 
  MInt downId = downStream().grid().tree().grid2solver(gridId);
 
  if(downId > 0) {
    ASSERT(downStream().a_level(downId) == upStream().a_level(upCell), "");
  }
 
  if(cellType == 1 || downId < 0) return downId;
 
  if(!lastLayer) {
    ASSERT(downStream().a_isBndryCell(upCell), "");
  } else {
    return -1;
  }
 
  if(cellType < 1) {
    // bndry-ghost cell
    const MInt srfc = -cellType;
    const MInt bndryId = downStream().a_bndryId(downId);
    downId = downStream().m_bndryCells->a[bndryId].m_srfcVariables[srfc]->m_ghostCellId;
    return downId;
 
  } else {
    // split child
    mTerm(1, AT_, "Not yet implemented for splitchilds");
  }
 
  return -1;
}

◆ updateZonalMapping()

template<MInt nDim, class SysEqn >

void CouplerFvMbZonal< nDim, SysEqn >::updateZonalMapping

private

Author: Tim Wegmann

Definition at line 266 of file fvmbzonal.cpp.

                                                        {
  TRACE();
 
  // resize to grid-cell offset
  m_upDown.resize(m_upStreamOffset);
  m_downUp.resize(m_downStreamOffset);
 
  // re-build bndryCell mapping
  if(upStream().isActive()) {
    if(m_zonalMethod == "P") {
      for(MInt upId = upStream().c_noCells(); upId < upStream().a_noCells(); upId++) {
        if(upStream().a_coordinate(upId, m_zonalDir) > m_zonalCoordinate) {
          MBool lastLayer = false;
          const MInt associatedId = upStream().getAssociatedInternalCell(upId);
          if(!upStream().a_hasNeighbor(associatedId, m_zonalDir * 2 + 1, false) && !upStream().a_isHalo(associatedId)) {
            const MInt parentId = upStream().c_parentId(associatedId);
            if(parentId > -1 && !upStream().a_hasNeighbor(parentId, m_zonalDir * 2 + 1, false)) {
              lastLayer = true;
            }
          }
          const MInt downId = up2downId(upId, lastLayer);
          if(downId > -1) {
            m_upDown.push_back(make_pair(upId, downId));
          }
        }
      }
    } else if(m_zonalMethod == "B") {
    }
  }
 
 
  if(downStream().isActive()) {
    if(m_zonalMethod == "P") {
      for(MInt downId = downStream().c_noCells(); downId < downStream().a_noCells(); downId++) {
        if(downStream().a_coordinate(downId, m_zonalDir) < m_zonalCoordinate) {
          MBool lastLayer = false;
          const MInt associatedId = downStream().getAssociatedInternalCell(downId);
          if(!downStream().a_hasNeighbor(associatedId, m_zonalDir * 2, false) && !downStream().a_isHalo(associatedId)) {
            const MInt parentId = downStream().c_parentId(associatedId);
            if(parentId > -1 && !downStream().a_hasNeighbor(parentId, m_zonalDir * 2, false)) {
              lastLayer = true;
            }
          }
          const MInt upId = down2upId(downId, lastLayer);
          if(upId > -1) {
            m_upDown.push_back(make_pair(downId, upId));
          }
        }
      }
    } else if(m_zonalMethod == "B") {
    }
  }
}

◆ upStream()

template<MInt nDim, class SysEqn >

FvMbSolver & CouplerFvMbZonal< nDim, SysEqn >::upStream ( ) const

inlineprivate

Definition at line 58 of file fvmbzonal.h.

58{ return *m_solverUp; }

Member Data Documentation

◆ m_couplingId

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::m_couplingId {}

private

Definition at line 67 of file fvmbzonal.h.

◆ m_downStreamId

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::m_downStreamId {}

private

Definition at line 66 of file fvmbzonal.h.

◆ m_downStreamOffset

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::m_downStreamOffset = 0

private

Definition at line 79 of file fvmbzonal.h.

◆ m_downUp

template<MInt nDim, class SysEqn >

std::vector<std::pair<MInt, MInt> > CouplerFvMbZonal< nDim, SysEqn >::m_downUp

private

Definition at line 76 of file fvmbzonal.h.

◆ m_noRKSteps

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::m_noRKSteps = -1

private

Definition at line 81 of file fvmbzonal.h.

◆ m_RKStep

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::m_RKStep = 0

private

Definition at line 80 of file fvmbzonal.h.

◆ m_solverDown

template<MInt nDim, class SysEqn >

FvMbSolver* CouplerFvMbZonal< nDim, SysEqn >::m_solverDown

private

Definition at line 63 of file fvmbzonal.h.

◆ m_solverUp

template<MInt nDim, class SysEqn >

FvMbSolver* CouplerFvMbZonal< nDim, SysEqn >::m_solverUp

private

Definition at line 62 of file fvmbzonal.h.

◆ m_upDown

template<MInt nDim, class SysEqn >

std::vector<std::pair<MInt, MInt> > CouplerFvMbZonal< nDim, SysEqn >::m_upDown

private

Definition at line 75 of file fvmbzonal.h.

◆ m_upStreamId

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::m_upStreamId {}

private

Definition at line 65 of file fvmbzonal.h.

◆ m_upStreamOffset

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::m_upStreamOffset = 0

private

Definition at line 78 of file fvmbzonal.h.

◆ m_zonalCoordinate

template<MInt nDim, class SysEqn >

MFloat CouplerFvMbZonal< nDim, SysEqn >::m_zonalCoordinate {}

private

Definition at line 71 of file fvmbzonal.h.

◆ m_zonalDir

template<MInt nDim, class SysEqn >

MInt CouplerFvMbZonal< nDim, SysEqn >::m_zonalDir {}

private

Definition at line 70 of file fvmbzonal.h.

◆ m_zonalDualTimeStepping

template<MInt nDim, class SysEqn >

MBool CouplerFvMbZonal< nDim, SysEqn >::m_zonalDualTimeStepping = false

private

Definition at line 73 of file fvmbzonal.h.

◆ m_zonalMethod

template<MInt nDim, class SysEqn >

MString CouplerFvMbZonal< nDim, SysEqn >::m_zonalMethod {}

private

Definition at line 69 of file fvmbzonal.h.

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

/home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/COUPLER/fvmbzonal.h
/home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/COUPLER/fvmbzonal.cpp

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ FVBndryCnd

◆ FvMbSolver

Constructor & Destructor Documentation

◆ CouplerFvMbZonal()

◆ ~CouplerFvMbZonal()

Member Function Documentation

◆ checkProperties()

◆ cleanUp()

◆ couplerId()

◆ createZonalMapping()

◆ down2upId()

◆ downStream()

◆ exchangeZonalValues()

◆ finalizeAdaptation()

◆ finalizeBalance()

◆ finalizeCouplerInit()

◆ finalizeSubCoupleInit()

◆ init()

◆ postCouple()

◆ preCouple()

◆ readProperties()

◆ subCouple()

◆ unifyTimeStep()

◆ up2downId()

◆ updateZonalMapping()

◆ upStream()

Member Data Documentation

◆ m_couplingId

◆ m_downStreamId

◆ m_downStreamOffset

◆ m_downUp

◆ m_noRKSteps

◆ m_RKStep

◆ m_solverDown

◆ m_solverUp

◆ m_upDown

◆ m_upStreamId

◆ m_upStreamOffset

◆ m_zonalCoordinate

◆ m_zonalDir

◆ m_zonalDualTimeStepping

◆ m_zonalMethod