search

CRPropa / CRPropa3

CRPropa is a public astrophysical simulation framework for propagating extraterrestrial ultra-high energy particles. https://crpropa.github.io/CRPropa3/
https://crpropa.desy.de
GNU General Public License v3.0
72 stars 69 forks source link

Turbulent Grids on Mac M1 not working #495

Open lukasmerten opened 3 months ago

lukasmerten commented 3 months ago

I have issues creating a turbulent grid on my Mac with M1Max chip, when using the following code from python.

lMin = 1 * parsec
lMax = 272 * parsec
Brms = 1
spacing = 0.5 * parsec
grid_n = 256
spectrum = SimpleTurbulenceSpectrum(Brms, lMin, lMax)
gp = GridProperties(Vector3d(0.), grid_n, spacing)
tf = SimpleGridTurbulence(spectrum, gp, seed)

The error message Trace/BPT trap: 5 is not super helpful but hints on some linking issues.

Using something like JF12Field().randomTurbulent(), which should call more or less the same code, works fine.

swig 4.2.1 numpy 1.26.0 clang-1500.3.9.4