From dcbdd5bc547c7d0f631a44c8a5be76724e4f7324 Mon Sep 17 00:00:00 2001
From: Olivier Parcollet <olivier.parcollet@cea.fr>
Date: Thu, 4 Sep 2014 11:58:20 +0200
Subject: [PATCH] Clean the eigenelements tests.

- Clean the output of the test.
- Separate from previous commit to be able to retest previous test
  with same lib code.
---
 test/triqs/arrays/eigenelements.cpp    | 57 ++++++++------------------
 test/triqs/arrays/eigenelements.output | 31 --------------
 2 files changed, 16 insertions(+), 72 deletions(-)
 delete mode 100644 test/triqs/arrays/eigenelements.output

diff --git a/test/triqs/arrays/eigenelements.cpp b/test/triqs/arrays/eigenelements.cpp
index 14cdc03d..9507d6ec 100644
--- a/test/triqs/arrays/eigenelements.cpp
+++ b/test/triqs/arrays/eigenelements.cpp
@@ -29,19 +29,20 @@ using namespace triqs::arrays::linalg;
 using dcomplex = std::complex<double>;
 
 template <typename T> void check_eig(matrix<T> M, matrix<T> vectors, array<double, 1> values) {
-  auto _ = range();
- for (auto i : range(0,first_dim(M)))  {
-  std::cerr << "check "<< i << std::endl;
-  std::cerr  << (M  -values(i))* vectors(i, _)<<std::endl;
+ auto _ = range();
+ for (auto i : range(0, first_dim(M))) {
+  std::cerr << "check " << i << std::endl;
+  std::cerr << (M - values(i)) * vectors(i, _) << std::endl;
   assert_all_close(M * vectors(i, _), values(i) * vectors(i, _), 1.e-14);
- }}
+ }
+}
 
-template<typename M> void test(M A) { 
-  auto w = eigenelements(make_clone(A));
-  std::cerr << "A = " << A << std::endl;
-  std::cerr << " values = " <<w.first << std::endl;
-  std::cerr << " vectors = " << w.second << std::endl;
-  check_eig (A, w.second, w.first);
+template <typename M> void test(M A) {
+ auto w = eigenelements(make_clone(A));
+ std::cerr << "A = " << A << std::endl;
+ std::cerr << " values = " << w.first << std::endl;
+ std::cerr << " vectors = " << w.second << std::endl;
+ check_eig(A, w.second, w.first);
 }
 
 int main(int argc, char **argv) {
@@ -54,18 +55,8 @@ int main(int argc, char **argv) {
     A(i, j) = (i > j ? i + 2 * j : i - j);
     A(j, i) = A(i, j);
    }
-  std::cerr << "A = " << A << std::endl;
-  
-  auto B = A;
-  auto w = eigenelements(B);
-  std::cout << "A = " << B << std::endl;
-  std::cout << " vectors = " << w.first << std::endl;
-  std::cout << " values = " << w.second << std::endl;
-  check_eig (A, w.second, w.first);
- }
+  test(A);
 
- {
-  matrix<double> A(3, 3);
   A() = 0;
   A(0, 1) = 1;
   A(1, 0) = 1;
@@ -73,31 +64,15 @@ int main(int argc, char **argv) {
   A(0, 2) = 2;
   A(2, 0) = 2;
 
-  auto B = A;
-  std::cout << "A = " << A << std::endl;
-  auto w =  eigenelements(B);
-  std::cout << " values = " <<w.first << std::endl;
-  std::cout << " vectors = " << w.second << std::endl;
-  check_eig (A, w.second, w.first);
- }
-
- {
-  matrix<double> A(3, 3);
+  test(A);
 
   A() = 0;
   A(0, 1) = 1;
   A(1, 0) = 1;
   A(2, 2) = 8;
 
-  auto B = A;
-  std::cout << "A = " << A << std::endl;
-  auto w =  eigenelements(B);
-  std::cout << " vectors = " << w.second << std::endl;
-  std::cout << " values = " <<w.first << std::endl;
-  std::cout << "A = " << A << std::endl;
-  check_eig (A, w.second, w.first);
+  test(A);
  }
-
  { // the complex case
 
   matrix<dcomplex> M(2, 2);
@@ -110,7 +85,7 @@ int main(int argc, char **argv) {
   test(M);
  }
 
-{ // the complex case
+ { // the complex case
 
   matrix<dcomplex> M(2, 2, FORTRAN_LAYOUT);
 
diff --git a/test/triqs/arrays/eigenelements.output b/test/triqs/arrays/eigenelements.output
deleted file mode 100644
index 0cc9804e..00000000
--- a/test/triqs/arrays/eigenelements.output
+++ /dev/null
@@ -1,31 +0,0 @@
-A = 
-[[-0.193187,-0.65727,0.728474]
- [-0.896252,0.420344,0.141577]
- [-0.399264,-0.625545,-0.670284]]
- vectors = [-4.13941,-0.784934,4.92434]
- values = 
-[[-0.193187,-0.65727,0.728474]
- [-0.896252,0.420344,0.141577]
- [-0.399264,-0.625545,-0.670284]]
-A = 
-[[0,1,2]
- [1,0,0]
- [2,0,8]]
- values = [-1.23962,0.761178,8.47844]
- vectors = 
-[[-0.767503,0.619144,0.166133]
- [-0.597371,-0.784798,0.165046]
- [0.232568,0.0274306,0.972193]]
-A = 
-[[0,1,0]
- [1,0,0]
- [0,0,8]]
- vectors = 
-[[-0.707107,0.707107,0]
- [0.707107,0.707107,0]
- [0,0,1]]
- values = [-1,1,8]
-A = 
-[[0,1,0]
- [1,0,0]
- [0,0,8]]