Exception: Cell type prism not supported for integral type exterior_facet

General
1

Hello,

I am dealing with a heat conduction problem whose r.h.s is a surface heat source. Therefore I need to calculate the integral confined to a exterior surface. I encountered the following error

Traceback (most recent call last):
  File "/opt/project/field_mat.py", line 145, in <module>
    b_t = form(source_expr * v_t * ufl.ds)
          ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/usr/local/dolfinx-real/lib/python3.12/dist-packages/dolfinx/fem/forms.py", line 337, in form
    return _create_form(form)
           ^^^^^^^^^^^^^^^^^^
  File "/usr/local/dolfinx-real/lib/python3.12/dist-packages/dolfinx/fem/forms.py", line 331, in _create_form
    return _form(form)
           ^^^^^^^^^^^
  File "/usr/local/dolfinx-real/lib/python3.12/dist-packages/dolfinx/fem/forms.py", line 254, in _form
    ufcx_form, module, code = jit.ffcx_jit(
                              ^^^^^^^^^^^^^
  File "/usr/local/dolfinx-real/lib/python3.12/dist-packages/dolfinx/jit.py", line 62, in mpi_jit
    return local_jit(*args, **kwargs)
           ^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/usr/local/dolfinx-real/lib/python3.12/dist-packages/dolfinx/jit.py", line 212, in ffcx_jit
    r = ffcx.codegeneration.jit.compile_forms([ufl_object], options=p_ffcx, **p_jit)
        ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/dolfinx-env/lib/python3.12/site-packages/ffcx/codegeneration/jit.py", line 225, in compile_forms
    raise e
  File "/dolfinx-env/lib/python3.12/site-packages/ffcx/codegeneration/jit.py", line 205, in compile_forms
    impl = _compile_objects(
           ^^^^^^^^^^^^^^^^^
  File "/dolfinx-env/lib/python3.12/site-packages/ffcx/codegeneration/jit.py", line 330, in _compile_objects
    _, code_body = ffcx.compiler.compile_ufl_objects(
                   ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/dolfinx-env/lib/python3.12/site-packages/ffcx/compiler.py", line 113, in compile_ufl_objects
    ir = compute_ir(analysis, _object_names, _prefix, options, visualise)
         ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/dolfinx-env/lib/python3.12/site-packages/ffcx/ir/representation.py", line 135, in compute_ir
    _compute_integral_ir(
  File "/dolfinx-env/lib/python3.12/site-packages/ffcx/ir/representation.py", line 330, in _compute_integral_ir
    points, weights, tensor_factors = create_quadrature_points_and_weights(
                                      ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/dolfinx-env/lib/python3.12/site-packages/ffcx/ir/representationutils.py", line 80, in create_quadrature_points_and_weights
    raise Exception(f"Cell type {cell} not supported for integral type {integral_type}.")
Exception: Cell type prism not supported for integral type exterior_facet.

It looks as if the current version of dolfinx (== 0.9.0) does not yet support the integration of external facet on prism cells. Are there any possible solutions?

2

This is being worked on atm: Make ds integrals on prism/pyramids generate kernels for each facet type by mscroggs · Pull Request #739 · FEniCS/ffcx · GitHub

3

Glad to know. Thanks a lot.