From ec6b597a5aa0eeb361d4afcf6c18760cf397c8ff Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Mon, 30 Jul 2012 10:24:04 +0200
Subject: [PATCH] Calculate Ax - b numerically
---
src/test-gradient-trivial.cc | 11 ++---------
1 file changed, 2 insertions(+), 9 deletions(-)
diff --git a/src/test-gradient-trivial.cc b/src/test-gradient-trivial.cc
index 3ad7394f..83d8c30c 100644
--- a/src/test-gradient-trivial.cc
+++ b/src/test-gradient-trivial.cc
@@ -37,15 +37,8 @@ int main() {
Functional::SmallVector start = b;
start *= 17;
- /*
- j(x)
- = Ax - b
- = 17*(6, 9.5) - (1, 2)
- = (102 - 1, 161.5 - 2)
- */
- Functional::SmallVector analytic_descent;
- analytic_descent[0] = -101;
- analytic_descent[1] = -159.5;
+ Functional::SmallVector analytic_descent = b;
+ A.mmv(start, analytic_descent);
Functional::SmallVector numerical_descent;
J.descentDirection(start, numerical_descent);
assert(two_distance<dim>(numerical_descent, analytic_descent) < 1e-10);
--
GitLab