Sitelet https://github.com/DedalusProject/dedalus/pull/337
Skip to content

Copy azimuthal-average m=0 modes by remaining group keys - #337

Open
shaneraphel wants to merge 1 commit into
DedalusProject:masterfrom
shaneraphel:azavg-m0-group-join-334
Open

shaneraphel wants to merge 1 commit into
DedalusProject:masterfrom
shaneraphel:azavg-m0-group-join-334

Conversation

@shaneraphel

@shaneraphel shaneraphel commented Sep 17, 2026 •

Copy link
Copy Markdown

Fixes #334

SphereAzimuthalAverage and SphericalAzimuthalAverage copy m = 0 with comp_out[m0_out] = comp_in[m0_in]. When M_max < L_max the input and output layouts have different local m=0 occupancy, so NumPy boolean assignment raises (96 into 64 in the report).

The current is the mask-length copy. This change copies m = 0 by the remaining group coordinates from local_group_arrays on both domains. PolarAzimuthalAverage.radial_matrix is not rewritten.

On unequal local masks the current raises; this change copies the overlapping groups. On equal masks this change matches current.

Sphere and spherical averages assigned m=0 with a boolean mask. When
M_max < L_max the local m=0 occupancy disagrees and NumPy raises.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

azimuthal average crashes in parallel when M_max<L_max

1 participant