Put back the redundant step_length_positive_definite() parameter.
authorMichael Orlitzky <michael@orlitzky.com>
Tue, 26 Mar 2013 00:48:04 +0000 (20:48 -0400)
committerMichael Orlitzky <michael@orlitzky.com>
Tue, 26 Mar 2013 00:48:04 +0000 (20:48 -0400)
tests/steepest_descent_tests.m

index 93b2782bbd2278c8aee11f4837adb246c97ea5a5..d6de1ccaf4b013539d45833a70c3e293bb0ad6b6 100644 (file)
@@ -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