I stumbled upon this using Sage, but developers there have traced it to a problem in GiNaC by running the following in ginsh:
> atan2(-Pi,0);
The response is:
power::eval(): division by zero
atan2 appears to work for other values, including (positive) (Pi,0), but (-Pi,0) fails.