I kindly ask someone to join as a reviewer to my PR
https://github.com/numpy/numpy/pull/21814 about handling np.float32 arrays in ma.mean and ma.std functions. The PR is aimed to fix different behaviour between np.core and ma variants of these functions.