fix(pzvec_extras): complex norm had a compile error

labmec · Jun 4, 2024 · 957e69c · 957e69c
1 parent cc75e1e
commit 957e69c
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/Util/pzvec_extras.h b/Util/pzvec_extras.h
@@ -273,9 +273,11 @@ RType(T) Norm(const TPZVec< T > &one) {
     RType(T) res = 0.;
     int size = one.NElements();
     if constexpr (is_complex<T>::value){
+        T cres{0};
         for (int i = 0; i < size; i++) {
-            res += (RType(T)) one[i] * std::conj(one[i]);
+            cres += one[i] * std::conj(one[i]);
         }
+        res = cres.real();
     }else{
         for (int i = 0; i < size; i++) {
             res += one[i]*one[i];