diff --git a/lib/lbfgs.c b/lib/lbfgs.c
index 88af94b..6f5ee4d 100644
--- a/lib/lbfgs.c
+++ b/lib/lbfgs.c
@@ -466,6 +466,7 @@ int lbfgs(
         step = 1.0 / sqrt(vecdot(d, d, n))
      */
     vec2norminv(&step, d, n);
+    step = 1.0;
 
     k = 1;
     end = 0;