Peano

Searching...

No Matches

#include <GridTraversalEventGeneratorTest.h>

Inheritance diagram for peano4::grid::tests::GridTraversalEventGeneratorTest:

[legend]

Collaboration diagram for peano4::grid::tests::GridTraversalEventGeneratorTest:

[legend]

Public Member Functions
	GridTraversalEventGeneratorTest ()

virtual void	run () override
	This routine is triggered by the TestCaseCollection.

Public Member Functions inherited from tarch::tests::TestCase
	TestCase (const std::string &testCaseName)
	Constructor.

virtual	~TestCase ()
	Destructor.

int	getNumberOfErrors () const

std::string	getTestCaseName () const

Private Member Functions
void	testAreFacesLocal1 ()

void	testAreFacesLocal2 ()
	A simple setup that I extracted from a bug in the code:

void	testGetFaceType ()
	Test face state identification.

void	testCreateLeaveCellTraversalEvent1 ()

void	testCreateEnterCellTraversalEvent1 ()

void	testCreateGenericCellTraversalEvent1 ()

Static Private Attributes
static tarch::logging::Log	_log
	Logging device.

Additional Inherited Members
Protected Member Functions inherited from tarch::tests::TestCase
	TestCase ()=delete

Protected Attributes inherited from tarch::tests::TestCase
const std::string	_testCaseName
	Name of the test case.

int	_errors
	Error counter.

bool	_error
	Indicate within macro if one specific validation failed.

Detailed Description

Definition at line 21 of file GridTraversalEventGeneratorTest.h.

Constructor & Destructor Documentation

◆ GridTraversalEventGeneratorTest()

peano4::grid::tests::GridTraversalEventGeneratorTest::GridTraversalEventGeneratorTest ( )

Definition at line 18 of file GridTraversalEventGeneratorTest.cpp.

Member Function Documentation

◆ run()

void peano4::grid::tests::GridTraversalEventGeneratorTest::run ( )

overridevirtual

This routine is triggered by the TestCaseCollection.

Implements tarch::tests::TestCase.

Definition at line 396 of file GridTraversalEventGeneratorTest.cpp.

References testMethod.

◆ testAreFacesLocal1()

void peano4::grid::tests::GridTraversalEventGeneratorTest::testAreFacesLocal1 ( )

private

Definition at line 76 of file GridTraversalEventGeneratorTest.cpp.

References peano4::grid::InvalidRank(), peano4::grid::GridVertex::setAdjacentRanks(), peano4::grid::GridVertex::setBackupOfAdjacentRanks(), peano4::grid::GridVertex::setState(), TwoPowerD, peano4::grid::GridVertex::Unrefined, and validateWithParams1.

Here is the call graph for this function:

◆ testAreFacesLocal2()

void peano4::grid::tests::GridTraversalEventGeneratorTest::testAreFacesLocal2 ( )

private

A simple setup that I extracted from a bug in the code:

In this example, the middle cell of the 3x3 coarsest pattern is deployed to rank 1. On rank 1, the marked face thus should not be local on rank 0.

Definition at line 24 of file GridTraversalEventGeneratorTest.cpp.

References peano4::grid::GridTraversalEventGenerator::areFacesLocal(), peano4::grid::GridVertex::HangingVertex, peano4::grid::GridVertex::setAdjacentRanks(), peano4::grid::GridVertex::setBackupOfAdjacentRanks(), peano4::grid::GridVertex::setState(), TwoPowerD, peano4::grid::GridVertex::Unrefined, and validateWithParams1.

Here is the call graph for this function:

◆ testCreateEnterCellTraversalEvent1()

void peano4::grid::tests::GridTraversalEventGeneratorTest::testCreateEnterCellTraversalEvent1 ( )

private

   Test to find out where faces go to and come from


   @image html GridTraversalEventGeneratorTest_testCreateEnterCellTraversalEvent.png


   This examples stems from the following (buggy) output:

   <pre>

18465645121 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::isVertexAdjacentToLocalSpacetree(...) in:vertex.toString():(state=Unrefined,adjacentRanks=[3,3,3,3],backupOfAdjacentRanks=[3,3,3,3],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.62963,0.407407],level=3),splittingIsConsideredLocal:1,joiningIsConsideredLocal:1 (file:peano4/grid/Spacetree.cpp,line:95) 18465660468 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::isVertexAdjacentToLocalSpacetree(...) out:result:0 (file:peano4/grid/Spacetree.cpp,line:108) 18465677543 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::isVertexAdjacentToLocalSpacetree(...) in:vertex.toString():(state=Unrefined,adjacentRanks=[3,2,3,2],backupOfAdjacentRanks=[3,2,3,2],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.666667,0.407407],level=3),splittingIsConsideredLocal:1,joiningIsConsideredLocal:1 (file:peano4/grid/Spacetree.cpp,line:95) 18465692370 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::isVertexAdjacentToLocalSpacetree(...) out:result:1 (file:peano4/grid/Spacetree.cpp,line:108) 18465709288 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::isVertexAdjacentToLocalSpacetree(...) in:vertex.toString():(state=New,adjacentRanks=[3,3,3,3],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=1684632135,x=[0.62963,0.444444],level=3),splittingIsConsideredLocal:1,joiningIsConsideredLocal:1 (file:peano4/grid/Spacetree.cpp,line:95) 18465723996 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::isVertexAdjacentToLocalSpacetree(...) out:result:0 (file:peano4/grid/Spacetree.cpp,line:108) 18465740694 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::isVertexAdjacentToLocalSpacetree(...) in:vertex.toString():(state=New,adjacentRanks=[3,2,3,2],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=1684632135,x=[0.666667,0.444444],level=3),splittingIsConsideredLocal:1,joiningIsConsideredLocal:1 (file:peano4/grid/Spacetree.cpp,line:95) 18466102574 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::isVertexAdjacentToLocalSpacetree(...) out:result:1 (file:peano4/grid/Spacetree.cpp,line:108) 18466130181 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::createGenericCellTraversalEvent(...) out:event.toString():(x=[0.648148,0.425926],h=[0.037037,0.037037],isRefined=0000,isVertexLocal=1010,isFaceLocal=0100,isCellLocal=0,vertexDataFrom=[0,0,0,0],vertexDataTo=[0,0,0,0],faceDataFrom=[0,0,0,0],faceDataTo=[0,0,0,0],cellData=2,relativePositionToFather=[2,2]) (file:peano4/grid/Spacetree.cpp,line:1735) 18466167134 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::createEnterCellTraversalEvent(...) out:state.toString():(level=3,x=[0.62963,0.407407],h=[0.037037,0.037037],inverted=1,evenFlags=11,accessNumber=[-1,-2,2,1]),event.toString():(x=[0.648148,0.425926],h=[0.037037,0.037037],isRefined=0000,isVertexLocal=1010,isFaceLocal=0100,isCellLocal=0,vertexDataFrom=[5,6,5,-2],vertexDataTo=[0,1,2,3],faceDataFrom=[4,3,-2,2],faceDataTo=[1,0,3,2],cellData=2,relativePositionToFather=[2,2]),_id:2 (file:peano4/grid/Spacetree.cpp,line:1826) 18466200951 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::receiveAndMergeUserDataAtHorizontalBoundary(...) in:state.toString():(level=3,x=[0.62963,0.407407],h=[0.037037,0.037037],inverted=1,evenFlags=11,accessNumber=[-1,-2,2,1]),enterCellTraversalEvent.toString():(x=[0.648148,0.425926],h=[0.037037,0.037037],isRefined=0000,isVertexLocal=1010,isFaceLocal=0100,isCellLocal=0,vertexDataFrom=[5,6,5,-2],vertexDataTo=[0,1,2,3],faceDataFrom=[4,3,-2,2],faceDataTo=[1,0,3,2],cellData=2,relativePositionToFather=[2,2]),_id:2 (file:peano4/grid/Spacetree.cpp,line:1832) 18466219530 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::getAdjacentDomainIds(...) in:faceNumber:2,calledByReceivingProcess:1,_id:2 (file:peano4/grid/Spacetree.cpp,line:1194) 18466237899 00:00:18 rank:0 core:0 trace peano4::grid::Spacetree::getAdjacentDomainIds(...) out:neighbour:0 (file:peano4/grid/Spacetree.cpp,line:1261) 18466259586 00:00:18 rank:0 core:0 info peano4::grid::Spacetree::receiveAndMergeUserDataAtHorizontalBoundary(...) receive and merge 3th face on stack 2 of tree 2 (relative position=0) with neighbour 0. Local position in cell=2, state=(level=3,x=[0.62963,0.407407],h=[0.037037,0.037037],inverted=1,evenFlags=11,accessNumber=[-1,-2,2,1]), inFaceStack=2 18466347886 00:00:18 rank:0 core:0 trace examples::exahype2::euler::observers::CreateGrid::receiveAndMergeFaceHorizontally(...) in:event.toString():(x=[0.648148,0.425926],h=[0.037037,0.037037],isRefined=0000,isVertexLocal=1010,isFaceLocal=0100,isCellLocal=0,vertexDataFrom=[5,6,5,-2],vertexDataTo=[0,1,2,3],faceDataFrom=[4,3,-2,2],faceDataTo=[1,0,3,2],cellData=2,relativePositionToFather=[2,2]),positionWithinCell:2,inOutStack:2,relativePositionOnInOutStack:0,fromTree:0,marker.toString():(x=[0.648148,0.425926],h=[0.037037,0.037037],select=2,is-cell-local=0),_spacetreeId:2 (file:observers/CreateGrid.cpp,line:1302) assertion in file observers/CreateGrid.cpp, line 1313 failed: tarch::la::equals( data.getDebugX(),incomingData.getDebugX() ) parameter data.getDebugX(): [0.666667,0.425926] parameter incomingData.getDebugX(): [0.833333,0.333333] parameter data.getDebugX(): [0.666667,0.425926] parameter incomingData.getDebugX(): [0.833333,0.333333] parameter fromTree: 0 parameter fromStack: 34 parameter inOutStack: 2 parameter relativePositionOnInOutStack: 0 parameter marker.toString(): (x=[0.648148,0.425926],h=[0.037037,0.037037],select=2,is-cell-local=0) parameter _spacetreeId: 2 peano4: observers/CreateGrid.cpp:1313: virtual void examples::exahype2::euler::observers::CreateGrid::receiveAndMergeFaceHorizontally(const peano4::grid::GridTraversalEvent&, int, int, int, int): Assertion ‘false’ failed. Aborted (core dumped)

 We run through this cell from left bottom to top right and this cell is
 created while we span this whole third level from scratch. That is, all
 grid elements are new. Indeed, only the bottom two vertices are unrefined,
 as they have been created by another cell before (the one below), while
 the top two vertices are new. The left top one has been created before
 (through the left neighbour), so the right top vertex is the one which is
 actually to be created by this event.

 The vertex indices in the output above are correct. Only one vertex (number
 three) is to be created. The flags for the faces are however wrong: The
 left and bottom one indeed are to be taken from temporary stacks, as we've
 visited their other adjacent neighbour/cell before. But the top one is to
 be created (correct) as well as the right one (wrong in dump above).

Definition at line 258 of file GridTraversalEventGeneratorTest.cpp.

◆ testCreateGenericCellTraversalEvent1()

void peano4::grid::tests::GridTraversalEventGeneratorTest::testCreateGenericCellTraversalEvent1 ( )

private

   A bug report where I observed that local is not valid.

   That's the input:

 57612233468  00:00:57     rank:0       core:1       trace        peano4::grid::GridTraversalEventGenerator::createGenericCellTraversalEvent(...)                 in:state.toString():(level=4,x=[0.37037,0.358025],h=[0.0123457,0.0123457],inverted=1,evenFlags=10,accessNumber=[-1,2,1,3]),relativePositionToFather:[0,2],fineGridVertices[0].toString():(state=HangingVertex,adjacentRanks=[0,0,0,0],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.37037,0.358025],level=4),fineGridVertices[1].toString():(state=Unrefined,adjacentRanks=[0,0,0,0],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.382716,0.358025],level=4),fineGridVertices[2].toString():(state=HangingVertex,adjacentRanks=[0,0,0,0],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.37037,0.37037],level=4),fineGridVertices[3].toString():(state=Unrefined,adjacentRanks=[0,0,0,0],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.382716,0.37037],level=4),_id:0 (file:peano4/grid/GridTraversalEventGenerator.cpp,line:229)

   which in the code gave us

 57612539282  00:00:57     rank:0       core:1       trace        peano4::grid::GridTraversalEventGenerator::createGenericCellTraversalEvent(...)                 out:event.toString():(x=[0.376543,0.364198],h=[0.0123457,0.0123457],hasBeenRefined=0000,willBeRefined=0000,i
sVertexLocal=1010,isFaceLocal=1111,isCellLocal=1,isParentCellLocal=1,isVertexAdjacentToParallelDomainBoundary=0000,isFaceAdjacentToParallelDomainBoundary=0000,vertexDataFrom=[30144400,0,30144496,0],vertexDataTo=[30144592,0,30144688,0],faceDataFrom=[0,0,30096912,0],faceD
ataTo=[30097328,0,30097744,0],cellData=30098160,relativePositionToFather=[0,2],invokingSpacetree=0,invokingSpacetreeIsNotInvolvedInAnyDynamicLoadBalancing=1),_id:0 (file:peano4/grid/GridTraversalEventGenerator.cpp,line:255)

   Here is the input in a nicer formatted way:

   state.toString():(level=4,x=[0.37037,0.358025],h=[0.0123457,0.0123457],inverted=1,evenFlags=10,accessNumber=[-1,2,1,3])
   relativePositionToFather:[0,2]
   fineGridVertices[0].toString():(state=HangingVertex,adjacentRanks=[0,0,0,0],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.37037,0.358025],level=4),
   fineGridVertices[1].toString():(state=Unrefined,adjacentRanks=[0,0,0,0],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.382716,0.358025],level=4),
   fineGridVertices[2].toString():(state=HangingVertex,adjacentRanks=[0,0,0,0],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.37037,0.37037],level=4),
   fineGridVertices[3].toString():(state=Unrefined,adjacentRanks=[0,0,0,0],backupOfAdjacentRanks=[0,0,0,0],hasBeenAntecessorOfRefinedVertexInPreviousTreeSweep=0,isAntecessorOfRefinedVertexInCurrentTreeSweep=0,numberOfAdjacentRefinedLocalCells=0,x=[0.382716,0.37037],level=4),
   _id:0 (file:peano4/grid/GridTraversalEventGenerator.cpp,line:229)

Definition at line 325 of file GridTraversalEventGeneratorTest.cpp.

References peano4::grid::GridVertex::HangingVertex, peano4::grid::GridVertex::Refined, peano4::grid::GridVertex::setAdjacentRanks(), peano4::grid::GridVertex::setBackupOfAdjacentRanks(), peano4::grid::GridVertex::setState(), state, TwoPowerD, peano4::grid::GridVertex::Unrefined, and validateEqualsWithParams1.

Here is the call graph for this function: