MWE :
from dolfin import *
mesh = UnitSquareMesh(3,3)
left = CompiledSubDomain("near(x[0], 0.) && on_boundary")
top = CompiledSubDomain("near(x[1], 1.) && on_boundary")
bottom = CompiledSubDomain("near(x[1], 0.) && on_boundary")
right = CompiledSubDomain("near(x[0], 1.) && on_boundary")
funcMarks = MeshFunction("size_t", mesh, 1)
bottom.mark(funcMarks, 1)
top.mark(funcMarks, 1)
left.mark(funcMarks, 2)
right.mark(funcMarks, 2)
mesh_LR = MeshView.create(funcMarks, 2)
mesh_TB = MeshView.create(funcMarks, 1)
subDomainLR = MeshFunction("size_t", mesh_LR, mesh_LR.topology().dim())
subDomainLR.set_all(0)
subDomainTB = MeshFunction("size_t", mesh_TB, mesh_TB.topology().dim())
subDomainTB.set_all(0)
left.mark(subDomainLR, 1)
right.mark(subDomainLR, 2)
bottom.mark(subDomainTB, 3)
top.mark(subDomainTB, 4)
ds_LR = Measure("dx", domain=mesh_LR, subdomain_data=subDomainLR)
ds_TB = Measure("dx", domain=mesh_TB, subdomain_data=subDomainTB)
print(assemble(1*ds_TB(3))) #gives incorrect result
Is marking of a submesh supported (it seems so at least looking at this post)
The above code gives incorrect result, however when I print
print(assemble(1*ds_TB))
I get the desired result (=2), which makes me wonder if there is an issue with the above code.