Put back the redundant step_length_positive_definite() parameter.

[octave.git] / tests / steepest_descent_tests.m

diff --git a/tests/steepest_descent_tests.m b/tests/steepest_descent_tests.m
index 93b2782bbd2278c8aee11f4837adb246c97ea5a5..d6de1ccaf4b013539d45833a70c3e293bb0ad6b6 100644 (file)
--- a/tests/steepest_descent_tests.m
+++ b/tests/steepest_descent_tests.m
@@ -23,7 +23,7 @@ q = @(x) (1/2)*x'*Q*x - b'*x;
 g = @(x) Q*x - b; % The gradient of q at x.
 
 % The step size algorithm to use in the steepest descent method.
-step_size = @(x) step_length_positive_definite(g(x), Q);
+step_size = @(x) step_length_positive_definite(g(x), Q, -g(x));
 sd = steepest_descent(g, x0, step_size, tolerance, max_iterations);
 
 diff = norm(cgm - sd, 'inf');
@@ -46,7 +46,7 @@ for n = [ 5, 10, 25, 50, 100 ]
   g = @(x) Q*x - b; % The gradient of q at x.
 
   % The step size algorithm to use in the steepest descent method.
-  step_size = @(x) step_length_positive_definite(g(x), Q);
+  step_size = @(x) step_length_positive_definite(g(x), Q, -g(x));
 
   ## pcg() stops when the /relative/ norm falls below tolerance. To
   ## eliminate the relativity, we divide the tolerance by the