In [35]:

    

from dolfin import *
mesh = UnitSquareMesh(int(2),int(2))
order = 1
Magnetic = FunctionSpace(mesh, "N1curl", order)
Lagrange = FunctionSpace(mesh, "CG", order)
W = Magnetic*Lagrange
def boundary(x, on_boundary):
    return on_boundary
bcW = DirichletBC(W.sub(0),Expression(("1.0","1.0")), boundary)
bcuW = DirichletBC(W.sub(1), Expression(("1.0")), boundary)


    

    




DEBUG:UFL:No integrals left after transformation, returning empty form.
DEBUG:FFC:Reusing form from cache.
DEBUG:UFL:No integrals left after transformation, returning empty form.
DEBUG:FFC:Reusing form from cache.
DEBUG:UFL:No integrals left after transformation, returning empty form.
DEBUG:FFC:Reusing form from cache.

In [8]:

    

print bcW.get_boundary_values()
print bcuW.markers()


    

    




{2L: 1.0, 4L: 0.5, 5L: 0.5, 7L: 1.0, 11L: 1.0, 12L: 0.5, 13L: 0.5, 15L: 1.0}
[ 2  4  5  7 11 12 13 15]

In [36]:

    

a = bcuW.get_boundary_values()


    

In [37]:

    

aa = a.keys()


    

In [43]:

    

aa
L = [1,2,3]       
" ".join(str(x).strip('L') for x in L)


    

    
Out[43]:





'1 2 3'

In [45]:

    

L


    

    
Out[45]:





[1, 2, 3]

In [ ]:

    

AAA.