GitLab

Fix major bug affecting analysis when using several plans that shared the legendre tables (same normalization, same lmax and mres, mmax smaller than earlier plan).

Fix array overflow with --enable-magic-layout and padding not allowed (bug was introduced with padding)

Bugfix for initial conditions of Galewsky et al. benchmark.
Approved-by: Nathanaël Schaeffer <nathanael.schaeffer@gmail.com>

Gaussian bump was not centered at [-pi;pi], but located at 0 with
x \in [0;2*pi] and calculating exp(-x*x).
This created a discontinuity in the initial conditions and hence a high
frequency (hard to see) oscillation in the geopotential field.

Change equation according to equation (4) in the original paper
Galewsky, J., Scott, R. K., & Polvani, L. M. (2004). An initial-value problem for testing numerical models of the global shallow-water equations. Tellus, Series A: Dynamic Meteorology and Oceanography, 56(5), 429–440. https://doi.org/10.1111/j.1600-0870.2004.00071.x
Otherwise, no southerly/northerly waves will show up

New: compiler of SHT kernels defaults to gcc and is independent of the compiler used for the rest of the library (set by CC=). Use --enable-kernel-compiler to override the default.

    Enable at run-time using the `SHT_ALLOW_PADDING` flag.

only two variants are needed "std" and "m". The former "ltr" is removed, because it was redundant with "std".

proof-of-concept: alternate openmp mode that does not rely on the fftw3 threaded implementation. Significantly faster for small transforms.