astropy.utils
astropy.utils.masked
Like numpy.nancumprod, skipping masked values as well.
numpy.nancumprod
Masked/NaN values are replaced with 1. The output is not masked.