def test_Refine_1D(self):
grid1d = RectangularGrid(3,1,1,1.0,1.0,1.0)
grid1dFine = RectangularGrid()
children = grid1dFine.refine(grid1d,2)
parent=0
child=0
for pN,cL in children.iteritems():
ok_(parent == pN)
for c in cL:
ok_(child == c.N)
child +=1
parent +=1
if GNUPLOT:
grid1dFine.writeEdgesGnuplot('grid1dFine')
grid1dFine.viewMeshGnuplotPipe('grid1dFine')
def test_RectangularGrid_2D():
grid2d = RectangularGrid(3, 3, 1, 1.0, 1.0, 1.0)
if GNUPLOT:
grid2d.writeEdgesGnuplot('grid2d')
grid2d.viewMeshGnuplotPipe('grid2d')
def test_RectangularGrid_3D(self):
grid3d = RectangularGrid(3,3,3,1.0,1.0,1.0)
if GNUPLOT:
grid3d.writeEdgesGnuplot('grid3d')
grid3d.viewMeshGnuplotPipe('grid3d')
def test_RectangularGrid_2D():
grid2d = RectangularGrid(3,3,1,1.0,1.0,1.0)
if GNUPLOT:
grid2d.writeEdgesGnuplot('grid2d')
grid2d.viewMeshGnuplotPipe('grid2d')