From ce671e990e9a957ea1ecb727624ac8eca338bd6d Mon Sep 17 00:00:00 2001
From: Joan Sola <jsola@iri.upc.edu>
Date: Tue, 14 Jun 2016 17:48:44 +0200
Subject: [PATCH] better control of intial conditions for q and R rotations

---
 src/examples/test_matrix_prod.cpp | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/examples/test_matrix_prod.cpp b/src/examples/test_matrix_prod.cpp
index 777483084..36a10867e 100644
--- a/src/examples/test_matrix_prod.cpp
+++ b/src/examples/test_matrix_prod.cpp
@@ -140,7 +140,9 @@ int main()
     }
     t1 = clock();
     std::cout << "Time w = R * v: " << (double)(t1 - t0) * 1e9 / CLOCKS_PER_SEC / N << "ns" << std::endl;
+    std::cout << "v norm change: " << 10*log((long double)v.norm()/(long double)vn) << " dB" << std::endl;
 
+    v << 1,2,3; vn = v.norm();
     t0 = clock();
     for (int i = 0; i < N; i++)
     {
-- 
GitLab