Bempp crashed when solving an electromagnetic scattering problem

[xiamengqi2012]

Hi, I tried Bempp on electromagnetic scattering problems and the jupyter kernel crashed. The simulation was running for several hours, and I don't think I ran out of memory. It's kind of hard for me to locate the cause of such a crash.

I encountered this problem earlier at a smaller resolution and I was able to bypass it by forcing specific versions of libraries (see below). But now I have this issue again on a larger resolution. I attached the jupyter log file and the links to the notebook file and mesh file I used. I'm wondering if there is something suspicious to you in the script or the log file; or if there is a specific python environment setting I should try. Thank you!

I'm using
python=3.8.5
ipykernel=5.3.4
ipython=7.19.0
numba=0.51.2
nbformat=5.0.8
Bempp-cl=0.2.4

jupyter.log
The script I used is here: https://drive.google.com/file/d/1cRnAL52rX7Jh_L7hqdwyxdKelu5PW8Tf/view?usp=sharing
The mesha I used i here: https://drive.google.com/file/d/1U7Oe2UzfIphigjpFjTrNY90U2CFRn-nH/view?usp=sharing

[xiamengqi2012]

Though the memory usage at the beginning of matrix assembly seemed fine (20Gb), the memory usage increased over the matrix assembly process. I was able to run this problem on a more powerful machine with 256Gb memory.

[xiamengqi2012]

After further increasing the resolution, I ran into another issue that cause the program to stop: python: ../nptl/pthread_mutex_lock.c:81: __pthread_mutex_lock: Assertion `mutex->__data.__owner == 0' . I wan't able to test convergence on this problem yet...