Support [`DiagonalQN`](https://github.com/JuliaSmoothOptimizers/LinearOperators.jl/blob/main/src/DiagonalHessianApproximation.jl) for 1-variable nonlinear element functions merged together