@@ -0,0 +1,416 @@
+from __future__ import absolute_import, print_function, division%0A%0Afrom firedrake import *%0A%0Amesh = UnitSquareMesh(2, 2)%0A%0AV = FunctionSpace(mesh, %22RT%22, 1)%0Af1 = Function(V)%0Af2 = Function(V)%0A%0A# Apply a strong condition at the top of the square domain%0Abc1 = DirichletBC(V, Constant((0.0, 10.0)), 4)%0A%0Abc1.apply(f1)%0A%0Abc2 = DirichletBC(V, Expression((%220.0%22, %2210.0%22)), 4)%0Abc2.apply(f2)%0A%0Aprint(f1.dat.data)%0Aprint(f2.dat.data)%0A

Add Constant vs Expression testing