Python GDAL: find out spikes and pits in DEM

I’m working with on a DEM (Geotiff) in Python GDAL, and I’d like to find out the spikes or pits in my terrain.
I have no idea how to proceed, does anyone have good algorithm to help me out?
Thanks,
eo