52 
print "Time = %s seconds"%( time ) 
print "Time = %s seconds"%( time ) 
53 


54 
solution.initialize(fixed_u_mask=boundary_cond,eta=eta,f=Y) 
solution.initialize(fixed_u_mask=boundary_cond,eta=eta,f=Y) 
55 
velocity,pressure=solution.solve(velocity,pressure,max_iter=max_iter,verbose=verbose,useUzawa=useUzawa) 
velocity,pressure=solution.solve(velocity,pressure,max_iter=max_iter,verbose=verbose,usePCG=True) 
56 


57 
print "Max velocity =", Lsup(velocity), "m/s" 
print "Max velocity =", Lsup(velocity), "m/s" 
58 

