Skip to content

Commit 48d0616

Browse files
geoffroylecontedpo
authored andcommitted
Adjoint opEye
1 parent 22cf605 commit 48d0616

File tree

2 files changed

+8
-0
lines changed

2 files changed

+8
-0
lines changed

src/special-operators.jl

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -25,6 +25,10 @@ struct opEye <: AbstractLinearOperator{Any} end
2525
*(T::AbstractLinearOperator, ::opEye) = T
2626
*(::opEye, T::opEye) = T
2727

28+
adjoint(A::opEye) = A
29+
transpose(A::opEye) = A
30+
conj(A::opEye) = A
31+
2832
function show(io::IO, op::opEye)
2933
println(io, "Identity operator")
3034
end

test/test_linop.jl

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -269,6 +269,10 @@ function test_linop()
269269
op2 = opEye()
270270
@test op === op2
271271
@test op === op * op2 === op2 * op
272+
273+
@test norm(transpose(op) * v - v) <= ϵ * norm(v)
274+
@test norm(adjoint(op) * v - v) <= ϵ * norm(v)
275+
@test norm(conj(op) * v - v) <= ϵ * norm(v)
272276
end
273277

274278
@testset "Ones" begin

0 commit comments

Comments
 (0)