From 59b9eccb0271f2ecf22b1a9d7ecf0b8b2744c82d Mon Sep 17 00:00:00 2001
From: simone-romiti <simone.romiti.1994@gmail.com>
Date: Tue, 7 Mar 2023 09:58:17 +0100
Subject: [PATCH] solved https://github.com/urbach/su2/issues/42

---
 include/su2.hh | 11 ++++++-----
 1 file changed, 6 insertions(+), 5 deletions(-)

diff --git a/include/su2.hh b/include/su2.hh
index cbcfdc5..5987751 100644
--- a/include/su2.hh
+++ b/include/su2.hh
@@ -4,9 +4,9 @@
  * @brief class representing an SU(2) matrix in the fundamental representation
  * @version 0.1
  * @date 2023-02-24
- * 
+ *
  * @copyright Copyright (c) 2023
- * 
+ *
  */
 
 #pragma once
@@ -79,7 +79,7 @@ public:
     b /= r;
   }
 
-  void print(){
+  void print() {
     std::cout << "--------------------\n";
     std::cout << a << " " << b << "\n";
     std::cout << -std::conj(b) << " " << std::conj(a) << "\n";
@@ -108,7 +108,9 @@ template <> inline _su2 dagger(const _su2 &u) {
 }
 
 template <> inline _su2 traceless_antiherm(const _su2 &x) {
-  return (_su2(0.5 * (x.geta() - std::conj(x.geta())), x.getb()));
+  const std::complex<double> a1 = 0.5 * (x.geta() - std::conj(x.geta()));
+  const std::complex<double> b1 = 0.5 * (x.getb() - std::conj(x.getb()));
+  return (_su2(a1, b1));
 }
 
 inline _su2 operator*(const _su2 &U1, const _su2 &U2) {
@@ -146,7 +148,6 @@ inline _su2 operator*(const _su2 &U1, const Complex &U2) {
   return (res);
 }
 
-
 using su2_accum = _su2;
 
 inline _su2 accum_to_Group(const su2_accum &x) {