FORTRAN: Bessel functions & combining operators
Should we move the Bessel function generator to fortran? This would make it more in line with the way we deal with other Chebyshev coefficients (see the Fermi operator). Also, we should look at combining operators like the Hamiltonian and current operators into one generic operator.