Fix

Author: kshyatt

test/testsuite/ad_utils.jl

@@ -165,6 +165,22 @@ end
 
 is_cpu(A) = typeof(parent(A)) <: Array
 
+"""
+    project_hermitian_inplace!(A, alg)
+
+Wrapper for `project_hermitian!(A, A, alg)`, invoked this way
+to avoid Enzyme's finite differences comparison getting confused.
+"""
+project_hermitian_inplace!(A, alg) = project_hermitian!(A, A, alg)
+
+"""
+    project_antihermitian_inplace!(A, alg)
+
+Wrapper for `project_hermitian!(A, A, alg)`, invoked this way
+to avoid Enzyme's finite differences comparison getting confused.
+"""
+project_antihermitian_inplace!(A, alg) = project_antihermitian!(A, A, alg)
+
 
 enzyme_fdm(T) = eltype(T) <: Union{Float32, ComplexF32} ? EnzymeTestUtils.FiniteDifferences.central_fdm(5, 1, max_range = 1.0e-2) : EnzymeTestUtils.FiniteDifferences.central_fdm(5, 1)
 

test/testsuite/enzyme/projections.jl

@@ -27,8 +27,8 @@ function test_enzyme_project_hermitian(
         B = instantiate_matrix(T, sz)
         alg = MatrixAlgebraKit.select_algorithm(project_hermitian, A)
         test_reverse(project_hermitian, RT, (A, TA), (alg, Const); atol, rtol, fdm)
-        #test_reverse(project_hermitian!, RT, (A, TA), (A, TA), (alg, Const); atol, rtol, fdm)
         test_reverse(project_hermitian!, RT, (A, TA), (B, TA), (alg, Const); atol, rtol, fdm)
+        test_reverse(project_hermitian_inplace!, RT, (A, TA), (alg, Const); atol, rtol, fdm)
     end
 end
 
@@ -47,7 +47,7 @@ function test_enzyme_project_antihermitian(
         B = instantiate_matrix(T, sz)
         alg = MatrixAlgebraKit.select_algorithm(project_hermitian, A)
         test_reverse(project_antihermitian, RT, (A, TA), (alg, Const); atol, rtol, fdm)
-        #test_reverse(project_antihermitian!, RT, (A, TA), (A, TA), (alg, Const); atol, rtol, fdm)
         test_reverse(project_antihermitian!, RT, (A, TA), (B, TA), (alg, Const); atol, rtol, fdm)
+        test_reverse(project_antihermitian_inplace!, RT, (A, TA), (alg, Const); atol, rtol, fdm)
     end
 end