fixing Florian's bug with overflow of norm for ints in test code

m-a-d-n-e-s-s · Sep 22, 2023 · 9050cb0 · 9050cb0
1 parent 704aaf9
commit 9050cb0
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 10 deletions.
diff --git a/src/madness/tensor/tensor.h b/src/madness/tensor/tensor.h
@@ -257,11 +257,11 @@ namespace madness {
     }
 
     namespace detail {
-        template <typename T> T mynorm(T t) {
-            return t*t;
+        template <typename T> double mynorm(T t) {
+            return double(t)*t;
         }
 
-        template <typename T> T mynorm(std::complex<T> t) {
+        template <typename T> double mynorm(std::complex<T> t) {
             return std::norm(t);
         }
     }

diff --git a/src/madness/tensor/test_tensor.cc b/src/madness/tensor/test_tensor.cc
@@ -53,13 +53,13 @@ namespace {
         return ((double) x)*x;
     }
 
-//    template <> double mynorm<float_complex>(float_complex x) {
-//        return (double) std::norm(x);
-//    }
-//
-//    template <> double mynorm<double_complex>(double_complex x) {
-//        return (double) std::norm(x);
-//    }
+    template <> double mynorm<float_complex>(float_complex x) {
+        return (double) std::norm(x);
+    }
+    
+    template <> double mynorm<double_complex>(double_complex x) {
+        return (double) std::norm(x);
+    }
 
     template <typename T, typename Q>
     inline