commit 7006a1f2b1701b8d89b8d1525fc0101943802221
author Sameer Agarwal <sameeragarwal@google.com> Thu Aug 16 18:04:22 2012 -0700
committer Sameer Agarwal <sameeragarwal@google.com> Thu Aug 16 18:04:22 2012 -0700
tree 1b49e775eb31fc29bc77ecae09c565dae73ca0cb
parent 57a44b27bc6fc95b4e70fdc25c25c9925a2072a0

Correct example code in Powell's function example.

Thanks to Petter Strandmark for pointing this out.

Change-Id: I967632235dccdb481396e94904bb911c9a1efe1e

docs/powell.tex

diff --git a/docs/powell.tex b/docs/powell.tex
index 7fc94a3..5ff6ddc 100644
--- a/docs/powell.tex
+++ b/docs/powell.tex
@@ -20,14 +20,14 @@
                         double* residuals,
                         double** jacobians) const {
     double x1 = parameters[0][0];
-    double x4 = parameters[1][0];
+    double x4 = parameters[0][3];
     // $f_4 = \sqrt{10} * (x_1 - x_4)^2$
     residuals[0] = sqrt(10.0) * (x1 - x4) * (x1 - x4)
     if (jacobians != NULL) {
-      jacobians[0][0] = 2.0 * sqrt(10.0) * (x1 - x4);   // $\partial_{x_1}f_1(x)$
-      jacobians[0][1] = 0.0;                            // $\partial_{x_2}f_1(x)$
-      jacobians[0][2] = 0.0;                            // $\partial_{x_3}f_1(x)$
-      jacobians[0][3] = -2.0 * sqrt(10.0) * (x1 - x4);  // $\partial_{x_4}f_1(x)$
+      jacobians[0][0] = 2.0 * sqrt(10.0) * (x1 - x4);   // $\partial_{x_1}f_4(x)$
+      jacobians[0][1] = 0.0;                            // $\partial_{x_2}f_4(x)$
+      jacobians[0][2] = 0.0;                            // $\partial_{x_3}f_4(x)$
+      jacobians[0][3] = -2.0 * sqrt(10.0) * (x1 - x4);  // $\partial_{x_4}f_4(x)$
     }
     return true;
   }