def test_outflow_X1Left(): boundaryLeft = boundaryPy.OUTFLOW boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim) varsNumpy = prim.getVars() if (prim.iLocalStart == 0): for var in xrange(prim.numVars): for i in xrange(numGhost): np.testing.assert_equal(varsNumpy[var, :, :, i], varsNumpy[var, :, :, numGhost])
def test_mirror_X1Right(): boundaryRight = boundaryPy.MIRROR boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim) varsNumpy = prim.getVars() if (prim.iLocalEnd == N1): for var in xrange(prim.numVars): for i in xrange(numGhost): np.testing.assert_equal( varsNumpy[var, :, :, prim.N1Local + numGhost + i], varsNumpy[var, :, :, prim.N1Local + numGhost - 1 - i])
def test_mirror_X3Back(): boundaryBack = boundaryPy.MIRROR boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim) varsNumpy = prim.getVars() if (prim.kLocalStart == 0): for var in xrange(prim.numVars): for k in xrange(numGhost): np.testing.assert_equal( varsNumpy[var, k, :, :], varsNumpy[var, 2 * numGhost - k - 1, :, :])
def test_outflow_X3Front(): boundaryFront = boundaryPy.OUTFLOW boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim) varsNumpy = prim.getVars() if (prim.kLocalEnd == N3): for var in xrange(prim.numVars): for k in xrange(numGhost): np.testing.assert_equal( varsNumpy[var, prim.N3Local + numGhost + k, :, :], varsNumpy[var, prim.N3Local + numGhost - 1, :, :])
def test_mirror_X2Top(): boundaryTop = boundaryPy.MIRROR boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim) varsNumpy = prim.getVars() if (prim.jLocalEnd == N2): for var in xrange(prim.numVars): for j in xrange(numGhost): np.testing.assert_equal( varsNumpy[var, :, prim.N2Local + numGhost + j, :], varsNumpy[var, :, prim.N2Local + numGhost - 1 - j, :])
def test_mirror_X1Left(): boundaryLeft = boundaryPy.MIRROR boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim ) varsNumpy = prim.getVars() if (prim.iLocalStart==0): for var in xrange(prim.numVars): for i in xrange(numGhost): np.testing.assert_equal(varsNumpy[var, :, :, i], varsNumpy[var, :, :, 2*numGhost-i-1] )
def test_outflow_X1Right(): boundaryRight = boundaryPy.OUTFLOW boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim ) varsNumpy = prim.getVars() if (prim.iLocalEnd==N1): for var in xrange(prim.numVars): for i in xrange(numGhost): np.testing.assert_equal(varsNumpy[var, :, :, prim.N1Local+numGhost+i], varsNumpy[var, :, :, prim.N1Local+numGhost-1] )
def test_mirror_X3Front(): boundaryFront = boundaryPy.MIRROR boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim ) varsNumpy = prim.getVars() if (prim.kLocalEnd==N3): for var in xrange(prim.numVars): for k in xrange(numGhost): np.testing.assert_equal(varsNumpy[var, prim.N3Local+numGhost+k, :, :], varsNumpy[var, prim.N3Local+numGhost-1-k, :, :] )
def test_outflow_X3Back(): boundaryBack = boundaryPy.OUTFLOW boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim ) varsNumpy = prim.getVars() if (prim.kLocalStart==0): for var in xrange(prim.numVars): for k in xrange(numGhost): np.testing.assert_equal(varsNumpy[var, k, :, :], varsNumpy[var, numGhost, :, :] )
def test_mirror_X2Top(): boundaryTop = boundaryPy.MIRROR boundaryPy.applyBoundaryConditionsPy(boundaryLeft, boundaryRight, boundaryTop, boundaryBottom, boundaryFront, boundaryBack, prim ) varsNumpy = prim.getVars() if (prim.jLocalEnd==N2): for var in xrange(prim.numVars): for j in xrange(numGhost): np.testing.assert_equal(varsNumpy[var, :, prim.N2Local+numGhost+j, :], varsNumpy[var, :, prim.N2Local+numGhost-1-j, :] )