Fixed a bug: TAU was initialized to zero in the cyth version of the code (this should be fixed later in a more general way).