Changeset 92 for trunk/test/test_wead.cc

Timestamp:

May 30, 2004, 12:44:39 AM (19 years ago)

Author:

Jari Häkkinen

Message:

Added test data for Wead, also wrote test for Wead.

File:

• trunk/test/test_wead.cc (modified) (2 diffs)

trunk/test/test_wead.cc

@@ -6 +6 @@
 
 #include "FileIO.h"
+#include "matrix.h"
 #include "Wead.h"
 
@@ -100 +101 @@
   theplu::cpptools::Wead wead(matrix_stream,flag_stream,option.neighbours());
   wead.estimate();
-  std::cout << wead.matrix() << std::endl;
-  return 0; // normal exit
+  std::ifstream control_stream("data/wead_result.data");
+  theplu::gslapi::matrix control(control_stream);
+  control-=wead.matrix();
+  for (unsigned int i=0; i<control.rows(); i++)
+    for (unsigned int j=0; j<control.columns(); j++)
+      // Jari, should we use GSL defined round off errors? Anyway, the
+      // hardcoded number below should be changed.
+      if (fabs(control(i,j))>5e-13)
+        exit(-1); // calculation result out of accepted error bounds
+  exit(0);        // normal exit
 }