The operators are written in a generic way appliable for a wide range of PDEs. Hence they should go into a custom package (i.e. FerriteOperators.jl).