Skip to content

Commit

Permalink
Merge with CPC and removing conflicts
Browse files Browse the repository at this point in the history
  • Loading branch information
mikaem committed Feb 18, 2014
1 parent 3be4705 commit 7a77f9e
Show file tree
Hide file tree
Showing 2 changed files with 0 additions and 23 deletions.
14 changes: 0 additions & 14 deletions problems/DrivenCavity.py
Original file line number Diff line number Diff line change
Expand Up @@ -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"
Expand Down Expand Up @@ -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:
Expand Down
9 changes: 0 additions & 9 deletions problems/TaylorGreen2D.py
Original file line number Diff line number Diff line change
Expand Up @@ -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,
Expand Down Expand Up @@ -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))
Expand Down

0 comments on commit 7a77f9e

Please sign in to comment.