Initial code for handling tau=0 as a special case.
Add test for tau=0 and improve test for small tau.
Update docs; add proper docstrings to functions exposed to Python.