search

bempp / bempp-cl

A fast Python based just-in-time compiling boundary element library
http://www.bempp.com
MIT License
149 stars 38 forks source link

Bempp crashed when solving an electromagnetic scattering problem #190

Open xiamengqi2012 opened 1 year ago

xiamengqi2012 commented 1 year ago

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 commented 1 year ago

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 commented 1 year ago

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...