Found w(theta) bug; left older version present in wtheta.c