From 3ec90803acf4d3468a33874a07f0605a7ec0aa6c Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Fri, 14 Aug 2026 18:21:13 -0400
Subject: [PATCH] no compil patches

---
 tests/debug/Test_schur_dense_coarse.cc |  4 +++-
 tests/debug/Test_schur_inverse.cc      | 14 ++++++++++----
 2 files changed, 13 insertions(+), 5 deletions(-)

diff --git a/tests/debug/Test_schur_dense_coarse.cc b/tests/debug/Test_schur_dense_coarse.cc
index a2168cab6..912ee90a2 100644
--- a/tests/debug/Test_schur_dense_coarse.cc
+++ b/tests/debug/Test_schur_dense_coarse.cc
@@ -223,7 +223,9 @@ int main (int argc, char ** argv)
           Lexicographic::CoorFromIndex(gci, gsi, Coarse5d->GlobalDimensions());
           csobj si;
           peekSite(si, Me, gci);
-          eA(i,j) = ((ComplexD *)&si)[bi];
+          // Explicit re/im at the thrust/std boundary (HIP builds)
+          ComplexD zz = ((ComplexD *)&si)[bi];
+          eA(i,j) = std::complex<double>(zz.real(), zz.imag());
         }
       }
       Eigen::JacobiSVD<Eigen::MatrixXcd> svd(eA);
diff --git a/tests/debug/Test_schur_inverse.cc b/tests/debug/Test_schur_inverse.cc
index 3dd986c88..e0a2651e8 100644
--- a/tests/debug/Test_schur_inverse.cc
+++ b/tests/debug/Test_schur_inverse.cc
@@ -640,13 +640,19 @@ int main (int argc, char ** argv)
           for(int64_t d=0; d<b; d++) Aglob[(uint64_t)(d + d*N)] += ComplexD(bshift,0.0);
         }
 
-        // Eigen fp64 oracle
+        // Eigen fp64 oracle.  Explicit re/im conversion at the boundary:
+        // on HIP builds ComplexD is thrust::complex, which has no
+        // operators against Eigen's std::complex.
+        auto toStd = [](const ComplexD &z) -> std::complex<double>
+        {
+          return std::complex<double>(z.real(), z.imag());
+        };
         Eigen::MatrixXcd eA(N,N);
         for(int64_t j=0; j<N; j++)
         {
           for(int64_t i=0; i<N; i++)
           {
-            eA(i,j) = Aglob[(uint64_t)(i + j*N)];
+            eA(i,j) = toStd(Aglob[(uint64_t)(i + j*N)]);
           }
         }
         Eigen::MatrixXcd Xref = eA.inverse();
@@ -709,8 +715,8 @@ int main (int argc, char ** argv)
         {
           for(int64_t i=0; i<N; i++)
           {
-            maxref = std::max(maxref, abs(Xref(i,j)));
-            maxdif = std::max(maxdif, abs(Xfull[(uint64_t)(i + j*N)] - Xref(i,j)));
+            maxref = std::max(maxref, std::abs(Xref(i,j)));
+            maxdif = std::max(maxdif, std::abs(toStd(Xfull[(uint64_t)(i + j*N)]) - Xref(i,j)));
           }
         }
         double cert2 = maxdif / maxref;