Dear Will, On 11/05/16 12:55, William Booker wrote:
Dear all,
I'm trying to assert that my projection in my incompressible case leaves my initial velocity to be discrete divergence free.
I've tried print an L2 norm but must be writing it wrong as it causes causes the compiler to break.
The relevant lines
solve(DIVU_perp == DIVU, u_perp, solver_parameters={'ksp_rtol': 1e-14, 'pc_type' : 'none', 'ksp_type': 'lsqr'} )
u0 += u_perp
# Print statement to verify initial velocity is divergence free
u_div = assemble( div_u( u0 , dFdPproject ) )
u_div_norm = norm (u_div, norm_type="L2")
print "l2 norm of discrete divergence on initial velocity is " % u_div_norm
How would I go about getting the script to print that the velocity is verified as divergence free?
We believe that this is a bug in COFFEE, but it was fixed yesterday afternoon. Can you please run firedrake-update followed by firedrake-clean and then try again. At least for me, your script runs. Cheers, Lawrence