Hi Thomas.  Thanks so much, the second option seems to match scipy.signal.argrelmax() perfectly. It’s not superfast, but definitely faster than scipy and scales better with increased data length.  I’ll also link the stackoverflow answer here - nobody there seemed to know how to do this. Cheers, Cat