Example #1
0
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])
Example #2
0
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])
Example #3
0
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, :, :])
Example #4
0
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, :, :])
Example #5
0
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, :])
Example #6
0
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]
                               )
Example #7
0
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]
                               )
Example #8
0
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, :, :]
                               )
Example #9
0
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, :, :]
                               )
Example #10
0
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, :]
                               )