How to use Bubble Functions

Yes! Try this:

from dolfin import *

mesh = UnitSquare(1, 1)
dim = mesh.topology().dim()

# Define bubble function space (polynomials of degree dim+1 that are
# zero on cell boundary):
B = FunctionSpace(mesh, "B", dim+1)
print "B = ", B

# Set all degrees of freedom to 1 and there you have your bubble
b_K = Function(B)
b_K.vector()[:] = 1.0

--
Marie

On 01/31/2012 11:41 AM, Thomas Fraunholz wrote:
> New question #186388 on DOLFIN:
> https://answers.launchpad.net/dolfin/+question/186388
>
> Hallo!
>
> Is there an easy way to define and use a Bubble function b_K on every triangle, eg. u*b_K*v*dx, where u and v are trial and test functions and b_K can be expressed in barycentric coordinates like \lambda_1*\lambda_2*\lambda_3 ?
>
> Thanks in advance,
> Thomas
>