Fastest phot() routine achieved