diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index a2b5b42bf8aabb46d96363319fa36da019e9d8f8..2067be0e554bc686fadf179163be42aa3be8fed7 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -269,6 +269,10 @@ int main(int argc, char *argv[]) {
               double const L = parset.get<double>("boundary.friction.ruina.L");
               double const unorm = u4_diff[i].two_norm();
 
+              // the (logarithmic) steady state corresponding to the current
+              // velocity
+              // std::cout << std::log(L/unorm * h) << std::endl;
+
               (*s4_new)[i] = compute_state_update(h, unorm, L, s4_old[i]);
             }
           }