From d8379c13d5b0c0b983f4eab158a4ee59e1694c97 Mon Sep 17 00:00:00 2001
From: Tom Benson <30674819+benson31@users.noreply.github.com>
Date: Wed, 25 Jan 2023 21:01:34 -0500
Subject: [PATCH] Extend the HEEV threshold to 5000

This allows us to gracefully handle a particular application up to 4096 problem size.
---
 src/hydrogen/device/cuSOLVER_API.cpp | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/src/hydrogen/device/cuSOLVER_API.cpp b/src/hydrogen/device/cuSOLVER_API.cpp
index 642ec108c..4d1755e96 100644
--- a/src/hydrogen/device/cuSOLVER_API.cpp
+++ b/src/hydrogen/device/cuSOLVER_API.cpp
@@ -11,7 +11,10 @@
 // exactly performance should fall off. The asymptotics of the Jacobi
 // method aren't as favorable as QR, but for small matrices, it should
 // have some benefit.
-#define CUSOLVER_HEEV_JACOBI_THRESHOLD 1024
+#define CUSOLVER_HEEV_JACOBI_THRESHOLD 5000
+// This is set up to graciously handle matrix sizes up to 4096, but I was
+// noticing an application was doing (2^k + 1). And if I'm going to put
+// "4097" as a magic number, I may as well just round up.
 
 namespace hydrogen
 {