The operators returned by get_operator are usually wrapped as R*M*R' which is fine if R'R=I, but this is generally not the case. Instead, the expression should be R*M*inv(R), which is dependent on ContinuumArrays.jl implementing inv(::AbstractQuasiMatrix) and that inv(R)*R=I is a no-op.
The operators returned by
get_operatorare usually wrapped asR*M*R'which is fine ifR'R=I, but this is generally not the case. Instead, the expression should beR*M*inv(R), which is dependent on ContinuumArrays.jl implementinginv(::AbstractQuasiMatrix)and thatinv(R)*R=Iis a no-op.