In [1]:

    

%matplotlib inline
import numpy as np
from Tiling import Domain, Tile, Point


    

In [2]:

    

dom = Domain(lo=[0,0], hi=[10,10])


    

In [3]:

    

atile = Tile(lo=[0,0], hi=[5,5], virtual=False)
btile = Tile(lo=[5,0], hi=[10,10], virtual=False)
ctile = Tile(lo=[0,5], hi=[5,10], virtual=False)


    

In [4]:

    

dom.tiles.append(atile)
dom.tiles.append(btile)
dom.tiles.append(ctile)


    

In [5]:

    

dom.plot_domain_slice(show_plot=True)


    

In [ ]:

    

ta = dom.tiles[0]
tb = dom.tiles[1]
print(ta.whether_occludes_tile(tb, 0))
print(ta.whether_occludes_tile(tb, 1))


    

In [ ]:

    

sf, ct = ta.whether_osculates_tile(tb,0)
print('Dimension 0')
print('ta surface: {}'.format(sf))
print('ct surface: {}'.format(ct))
sf, ct = ta.whether_osculates_tile(tb,1)
print('Dimension 1')
print('ta surface: {}'.format(sf))
print('ct surface: {}'.format(ct))


    

In [ ]: