Open jorgensd opened 9 months ago
I am trying this out on dolfinx=0.8.0 and running into the error
(8, np.asarray(integration_entities, dtype=np.int32))]) ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ ValueError: setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (20,) + inhomogeneous part.
I believe that the line integration_entities.append(local_index)
should instead be integration_entities.append(local_index[0])
@molel-gt I've updated the code to take this into account. Thanks for the feedback.
Source: https://fenicsproject.discourse.group/t/wrong-facetnormal-vector-on-internal-boundaries/12887/2?u=dokken