import sys
from os import sep
sys.path.append('..'+sep+'pyvolution')

#========================================================================
from config import g, epsilon
from Numeric import allclose, array, zeros, ones, Float
from parallel_advection import *
from Numeric import array
from parallel_meshes import *

import pypar

numprocs = pypar.size()
myid = pypar.rank()
processor_name = pypar.Get_processor_name()



#N = 50
N = 25
M = 25

points, vertices, boundary, full_send_dict, ghost_recv_dict = parallel_rectangular(N, M, len1=1.0)

#Create advection domain with direction (1,-1)
domain = Parallel_Domain(points, vertices, boundary,
                         full_send_dict, ghost_recv_dict, velocity=[1.0, 0.0])

# Initial condition is zero by default

#domain.initialise_visualiser()

#Boundaries
T = Transmissive_boundary(domain)
D = Dirichlet_boundary(array([1.0]))

domain.default_order = 2


domain.set_boundary( {'left': T, 'right': T, 'bottom': T, 'top': T} )
domain.check_integrity()

class Set_Stage:
    """Set an initial condition with constant water height, for x<x0
    """

    def __init__(self, x0=0.25, x1=0.5, h=1.0):
        self.x0 = x0
        self.x1 = x1
        self.h  = h

    def __call__(self, x, y):
        return self.h*((x>self.x0)&(x<self.x1))

if myid == 0:
    domain.set_quantity('stage', Set_Stage(0.2,0.4,1.0))
    import time
    t0 = time.time()

#Check that the boundary value gets propagated to all elements
for t in domain.evolve(yieldstep = 0.1, finaltime = 2.0):
    if myid == 0:
        domain.write_time()

if myid == 0:
    print 'That took %.2f seconds' %(time.time()-t0)