add diagop.jl for diagonal operating functions