From 7a77f9eae84a8bf6d0b601ed4a86f4c6eb57aa07 Mon Sep 17 00:00:00 2001
From: Mikael Mortensen <mikaem@math.uio.no>
Date: Tue, 18 Feb 2014 07:10:58 -0800
Subject: [PATCH] Merge with CPC and removing conflicts

---
 problems/DrivenCavity.py  | 14 --------------
 problems/TaylorGreen2D.py |  9 ---------
 2 files changed, 23 deletions(-)

diff --git a/problems/DrivenCavity.py b/problems/DrivenCavity.py
index a02af293..5a3918bd 100644
--- a/problems/DrivenCavity.py
+++ b/problems/DrivenCavity.py
@@ -25,16 +25,6 @@ def mesh(Nx, Ny, **params):
     print_intermediate_info = 100,
     use_krylov_solvers = True)
 
-<<<<<<< HEAD
-def pre_solve_hook(mesh, velocity_degree, constrained_domain, **NS_namespace):    
-    # Declare a Function used for plotting in temporal_hook
-    Vv = VectorFunctionSpace(mesh, 'CG', velocity_degree, constrained_domain=constrained_domain)
-    return dict(uv=Function(Vv), Vv=Vv)
-
-noslip = "std::abs(x[0]*x[1]*(1-x[0]))<1e-8"
-lid    = "std::abs(x[1]-1) < 1e-8"
-=======
->>>>>>> CPC
 # Specify boundary conditions
 noslip = "std::abs(x[0]*x[1]*(1-x[0]))<1e-8"
 top    = "std::abs(x[1]-1) < 1e-8"
@@ -67,11 +57,7 @@ def theend_hook(u_, p_, uv, Vv, **NS_namespace):
     plot(p_, title='Pressure')
 
     try:
-<<<<<<< HEAD
         from fenicstools import StreamFunction
-=======
-        from fenicstools.Streamfunctions import StreamFunction
->>>>>>> CPC
         psi = StreamFunction(u_, [], use_strong_bc=True)
         plot(psi, title='Streamfunction', interactive=True)
     except:
diff --git a/problems/TaylorGreen2D.py b/problems/TaylorGreen2D.py
index cab855d3..b893567b 100644
--- a/problems/TaylorGreen2D.py
+++ b/problems/TaylorGreen2D.py
@@ -12,11 +12,6 @@
     dt = 0.001,
     Nx = 20, Ny = 20,
     folder = "taylorgreen2D_results",
-<<<<<<< HEAD
-=======
-    max_iter = 1,
-    iters_on_first_timestep = 2,
->>>>>>> CPC
     plot_interval = 1000,
     save_step = 10000,
     checkpoint = 10000,
@@ -106,12 +101,8 @@ def temporal_hook(q_, t, nu, VV, dt, plot_interval, initial_fields, tstep, sys_c
             error = norm(ue.vector())/uen
             err[ui] = "{0:2.6e}".format(norm(ue.vector()))
             total_error[i] += error*dt     
-<<<<<<< HEAD
         if MPI.process_number() == 0:
             print "Error is ", err, " at time = ", t 
-=======
-        print "Error is ", err, " at time = ", t 
->>>>>>> CPC
         
 def theend_hook(mesh, q_, t, dt, nu, VV, sys_comp, initial_fields, **NS_namespace):
     final_error = zeros(len(sys_comp))