_assert(dxij.shape == pdf.shape,
                'Shape of pdf does not match the arguments')