Use a less scale-dependent assertion

src/compute_state.cc

@@ -57,6 +57,8 @@ double state_update_dieterich(double h, double uol, double old_state) {
 
      log((ret - old_state + uol)/h) = -ret
   */
-  assert(std::abs(std::log((ret - old_state + uol) / h) + ret) < 1e-13);
+  assert(std::min(std::abs(ret - old_state + uol - h * std::exp(-ret)),
+                  std::abs(std::log((ret - old_state + uol) / h) + ret)) <
+         1e-12);
   return ret;
 }