incorrect assemble in subdomain for seond order Function space for some indexes

[Andrey Semenov]

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

‌