incorrect assemble in subdomain for seond order Function space for some indexes
Issue #1118
new
dolfin.version '2019.2.0.dev0' ubuntu 20.10
for subdomain indexes n in range 9:20 assemble equal to 0 but must be some (1. in this example).
for FunctionSpace(mesh, 'P', 1),
FunctionSpace(mesh, 'P', 3),FunctionSpace(mesh, 'P', 4) calculations is correct assemble == 1.0
from fenics import *
mesh = UnitSquareMesh(8, 8)
subv = MeshFunction('size_t', mesh, 2)
V = FunctionSpace(mesh, 'P', 2)
u = Function(V)
u.assign(Constant(1.))
def fun(n):
subv.set_all(n)
DX = dx(subdomain_data=subv)
return [assemble(u*dx),assemble(u*DX(n))]
print(fun(1),fun(9))
>>> [1.0, 1.0] [1.0, 0.0]
# 0.0 - is incorrect