I am looking to run the code from the solution to the question asked in this previous forum from 2017:
The issue is when I run this code in the latest fenics version, there is a segmentation fault in the function:
vb_v = np.array(bmesh.entity_map(0), dtype=int)
This code runs just fine when I use a docker container quay io/fenicsproject/stable:2017.1.0, but the segmentation fault occurs when using the latest version quay io/fenicsproject/stable:latest.
The segmentation fault occurs when trying to convert the entity_map object into a numpy array. Does anyone have any ideas on the cause of this and/or possible solutions for the latest version? Many thanks.