From: Michael Orlitzky Date: Mon, 25 Mar 2013 22:48:14 +0000 (-0400) Subject: Tiny refactoring in the PCGM. X-Git-Url: https://gitweb.michael.orlitzky.com/?a=commitdiff_plain;h=92116b34e755b3ef5de14a1777676bc09180f007;p=octave.git Tiny refactoring in the PCGM. --- diff --git a/optimization/preconditioned_conjugate_gradient_method.m b/optimization/preconditioned_conjugate_gradient_method.m index dec2eee..e2f841a 100644 --- a/optimization/preconditioned_conjugate_gradient_method.m +++ b/optimization/preconditioned_conjugate_gradient_method.m @@ -66,7 +66,8 @@ function [x, k] = preconditioned_conjugate_gradient_method(Q, ... % % We use this in the inner loop. - sqrt_n = floor(sqrt(length(x0))); + n = length(x0); + sqrt_n = floor(sqrt(n)); % Set k=0 first, that way the references to xk,rk,zk,dk which % immediately follow correspond (semantically) to x0,r0,z0,d0.