Speed it up

[lamtung16]

I utilize numpy vectorization, resulting in significantly faster processing speeds.
My Test Code

[deepcharles]

Thanks for this comparison. The function from ruptures that you used is not really optimized for the piecewise constant model (it is pure python). You should compare to rpt.KernelCPD(kernel="linear") instead, which calls a C implementation.