subdomain_id fails with various dtypes

Issue #174 new
Jan Blechta created an issue

When using various dtypes for subdomain_id, tabulate_tensor() is still correctly generated but create_cell_integral() omits the case. MWE:

import numpy
c = Constant(interval)
#L = c*dx(42)             # works
#L = c*dx(numpy.int(42))  # works
L = c*dx(numpy.intc(42))  # fails; as well as many other dtypes

Might not be relevant in FFCX but it might turn out that the bug boils down to UFL.

Comments (0)

  1. Log in to comment