Cleanup

src/test-gradient-parabola.cc

@@ -42,6 +42,8 @@ int main() {
   M[0][0] += 2;
   M[1][1] += 2;
 
+  double ret1;
+  {
     Functional::SmallVector start = b;
     start *= 17;
 
@@ -52,15 +54,26 @@ int main() {
     J.descentDirection(start, numerical_descent);
     assert(two_distance<dim>(numerical_descent, analytic_descent) < 1e-10);
 
-  double const ret1 = functionTester(J, start, 6);
+    ret1 = functionTester(J, start, 6);
     assert(two_distance<dim>(start, solution) < 1e-6);
+  }
 
+  double ret2;
+  {
     // Something random
+    Functional::SmallVector start;
     start[0] = 279;
     start[1] = -96;
 
-  double const ret2 = functionTester(J, start, 15);
-  assert(two_distance<dim>(start, solution) < 1e-6);
+    Functional::SmallVector analytic_descent = b;
+    M.mmv(start, analytic_descent);
+
+    Functional::SmallVector numerical_descent;
+    J.descentDirection(start, numerical_descent);
+    assert(two_distance<dim>(numerical_descent, analytic_descent) < 1e-10);
 
+    ret2 = functionTester(J, start, 15);
+    assert(two_distance<dim>(start, solution) < 1e-6);
+  }
   assert(std::abs(ret1 - ret2) < 1e-11);
 }