I’m making an attempt to resolve the Einstein-Klein-Gordon system (equations in file) for various piecewise potentials however the kernel retains crashing after repeated use of NDSolve.
Apparently, every little thing works advantageous at first, which means that I do not get any related warning and I receive options with the options I anticipate. Nonetheless, after a while (not all the time the identical time period) the kernel crashes with none earlier warning.
I attempted to maintain monitor of the reminiscence used at every step, however it by no means reaches a big worth, so it would not appear to me that it’s a reminiscence associated crash.

Let me present you a file with brief model of my code that reproduces the crash. I am sorry I’ve to supply a hyperlink however I am unable to identify the a part of the code that results in the crash…
As you’ll be able to see within the file, the perform bosonstar can be utilized to seek out options with none obvious drawback. However whether it is ran a number of instances, it will definitely results in a crash.

Btw, I am utilizing Mathematica 12.1 on a macOS Mojave model 10.14.6

I might vastly recognize any remark, recommendation or trace! Thanks!!