diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 1aac9192059fa5a05faf1faf7a5e3125a0975710..9d2f7541cae1d9c8d3fa7d0d06c5c76aa51bad1c 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -243,13 +243,10 @@ int main(int argc, char *argv[]) {
     auto const &neumannFunction = functions.get("neumannCondition");
 
     // Find a (somewhat random) frictional node
-    size_t first_frictional_node;
-    for (size_t i = 0; i < frictionalNodes.size(); ++i) {
-      if (frictionalNodes[i][0]) {
-        first_frictional_node = i;
-        break;
-      }
-    }
+    size_t first_frictional_node = 0;
+    while (!frictionalNodes[first_frictional_node][0] &&
+           first_frictional_node < frictionalNodes.size())
+      ++first_frictional_node;
 
     for (size_t run = 1; run <= timesteps; ++run) {
       double const time = h * run;