Can there be compiler bugs in structure padding?  Simplify State
Use c.n_nu instead of N_NU in getdt() and rhd()
Fix bug: initialize I and tau to zero
Configure beta threshold at runtime
c.threshold is defined as 1 / beta_threshold
Enable extra warnings but disable the unused-function warning
Check return values of fread() and fscan() so glib doesn't complain
Follow the convention used in the paper to allow for td and tf