Unify set and add dissipation + improve robustness flexibility
The user can
- Add one or more dissipation terms throughout the simulation
- Remove one or more dissipation terms throughout the simulation by setting the coefficient to None
- Remove all forms of dissipation by passing the formula None to add_dissipation
Automatic checks are performed to ensure that only one type (linear or quadratic) is possible at any time.