Segmentation fault happened when use axisymmetric shtns functions on KNL nodes

magic-sph / magic

MagIC is a high-performance code that solves the magneto-hydrodynamics equations in rotating spherical shells

https://magic-sph.github.io/

80 stars 25 forks source link

Segmentation fault happened when use axisymmetric shtns functions on KNL nodes #17

Closed alexisreboulsalze closed 5 years ago

alexisreboulsalze commented 5 years ago

Some outputs (power, perpPar, TO, etc.) call shtns.axi_to_spat, shtns.spat_to_SH_axi and lead to a segmentation fault on KNL nodes. I'm using the latest version of the code and the SHTNS 3.2.2 version compiled with the configuration command : ./configure --enable-mkl --enable-knl --enable-magic-layout --prefix=$HOME/ Without these outputs, the code works well otherwise.

tgastine commented 5 years ago

This is an SHTns issue with the MagIC layout. Please clone the latest version on SHTns bitbucket repository it should work fine.

alexisreboulsalze commented 5 years ago

Thanks for the fast answer, it works now.