couplingutils.h File Reference

Go to the source code of this file.

Functions
template<class SolverA , class SolverB >
MInt	convertId (SolverA &solverA, SolverB &solverB, const MInt solverAId)
	Conversion from solverA id to the solverB id on the same-level only! More...
template<class SolverA , class SolverB >
MInt	convertIdParent (SolverA &solverA, SolverB &solverB, const MInt solverAId)
	Conversion from solverA id to the solverB id If no cell on the same level is found, a connection through the Cell-parents is used! More...

Function Documentation

convertId()

template<class SolverA , class SolverB >

MInt convertId	(	SolverA &	solverA,
		SolverB &	solverB,
		const MInt	solverAId
	)

Author

Tim Wegmann, Julian Vorspohl

Definition at line 21 of file couplingutils.h.

                                                                         {
  if(solverAId >= solverA.c_noCells() || solverAId < 0) return -1;
  solverA.assertValidGridCellId(solverAId);
 
  const MInt gridId = solverA.grid().tree().solver2grid(solverAId);
  ASSERT(solverA.grid().solverFlag(gridId, solverA.solverId()), "");
 
  if(!solverB.grid().solverFlag(gridId, solverB.solverId())) return -1;
 
  const MInt solverBId = solverB.grid().tree().grid2solver(gridId);
 
  if(solverBId > 0) {
    ASSERT(solverA.a_level(solverAId) == solverB.a_level(solverBId),
           std::to_string(solverAId) + " " + std::to_string(solverA.a_level(solverAId)) + " "
               + std::to_string(solverB.a_level(solverBId)));
  }
 
  return solverBId;
}

convertIdParent()

template<class SolverA , class SolverB >

MInt convertIdParent	(	SolverA &	solverA,
		SolverB &	solverB,
		const MInt	solverAId
	)

Author

Tim Wegmann, Julian Vorspohl

Definition at line 46 of file couplingutils.h.

                                                                               {
  if(solverAId >= solverA.c_noCells() || solverAId < 0) return -1; // not even a grid-cell, still return -1
  solverA.assertValidGridCellId(solverAId);
 
  MInt gridId = solverA.grid().tree().solver2grid(solverAId);
  ASSERT(solverA.grid().solverFlag(gridId, solverA.solverId()), "");
 
  if(solverB.grid().solverFlag(gridId, solverB.solverId())) { // connection on the same level!
    return solverB.grid().tree().grid2solver(gridId);
 
  } else { // going through all parents! to find a connection
 
    MInt solverAParentId = solverA.c_parentId(solverAId);
    MInt solverBParentId = convertId(solverA, solverB, solverAParentId);
 
    while(solverBParentId < 0 && solverAParentId > -1 && solverAParentId < solverA.c_noCells()) {
      solverAParentId = solverA.c_parentId(solverAParentId);
      solverBParentId = convertId(solverA, solverB, solverAParentId);
    }
    return solverBParentId;
  }
}