(seems) a more efficient way to use diffEqOperator.jl

[ChrisRackauckas]

@Fred_He

ctually A1 and A2 are Diagonal matrix, we can further mod

Yes, it needs more optimizations. It should use @turbo internally.

Try @turbo on the convolution definitions