Adjoint opEye

JuliaSmoothOptimizers · Jan 8, 2024 · 48d0616 · 48d0616
1 parent 22cf605
commit 48d0616
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 0 deletions.
diff --git a/src/special-operators.jl b/src/special-operators.jl
@@ -25,6 +25,10 @@ struct opEye <: AbstractLinearOperator{Any} end
 *(T::AbstractLinearOperator, ::opEye) = T
 *(::opEye, T::opEye) = T
 
+adjoint(A::opEye) = A
+transpose(A::opEye) = A
+conj(A::opEye) = A
+
 function show(io::IO, op::opEye)
   println(io, "Identity operator")
 end

diff --git a/test/test_linop.jl b/test/test_linop.jl
@@ -269,6 +269,10 @@ function test_linop()
       op2 = opEye()
       @test op === op2
       @test op === op * op2 === op2 * op
+
+      @test norm(transpose(op) * v - v) <= ϵ * norm(v)
+      @test norm(adjoint(op) * v - v) <= ϵ * norm(v)
+      @test norm(conj(op) * v - v) <= ϵ * norm(v)
     end
 
     @testset "Ones" begin