Get the FMM working with Bessel expansions on the 2D Helmholtz kernel.

- Add a H2DMultipoleExpansion class.
- Implement L2L, M2L, M2M operators.
- Teach loopy about casting to complex numbers (to handle a corner
  case where some multipole coefficients may look real valued).
- Update the tests.