We should also have a test for the divide-by-zero case and specify what the behavior is in the doc (not just here --- we have tons of boundary conditions like this to specify in the doc).  For now, I think it's OK to follow the std lib's behavior for 1 / 0.