def METHOD_NAME(self):
    cube_small = self.cube_a[:, 0, 0]
    cube_small_masked = cube_small.copy()
    cube_small_masked.data = ma.array(
        cube_small.data, mask=np.array([0, 0, 0, 1, 1, 1], dtype=bool)
    )
    r = stats.pearsonr(cube_small, cube_small_masked, common_mask=True)
    self.assertArrayAlmostEqual(r.data, np.array([1.0]))