The paper Fusion rules and modular transformations in 2D conformal field theory by Erik Verlinde mentions a simple case of rational conformal field theory, where the fusion algebra is just $(\mathbb{Z}_p,+)$:
$$ \phi_p \times \phi_{p'} = \phi_{p + p'}$$
Through some mysterious physics logic he obtains an identity relating various characters of the Virasoro algebra:
$$ S: \chi_p \mapsto \sum_{p'\in \mathbb{Z}_n} e^{\frac{2\pi i pp'}{N}}\chi_{p'}$$
In this particular case, the characters are just theta functions, but Verlinde doesn't specify which ones. The partition function is:
$$ Z = \mathrm{tr}[q^{L_0 - \frac{1}{24} }] = \frac{1}{\eta(q)}\sum q^{\frac{1}{2}\left(\tfrac{n}{R}+ \frac{1}{2}mR\right)^2} \tag{$\ast$}$$
The characters are given as partial traces over the different primary fields:
$$ \chi_i = \mathrm{tr}_{[\phi_i]}[q^{L_0 - \frac{1}{24} }]$$
The exercise for now is to re-write $(\ast)$ as an identity of theta functions and prove it for this very special case. How is it theta functions transform as the Discrete Fourier Transform
The trouble is I don't know how the traces for $\chi_i$ were computed. Also, the series Verlinde generates even more general than theta functions and can be defined for general Riemann surfaces, so I have been to understand them.
I was able to find work of Arnaud Beauville where he discusses the Verlinde formula which interprets this as the dimension of a line bundle over the moduli space of curves. Still trying to understand if this may be related, but I sense this particular case should be known.
- Arnaud Beauville, Vector bundles on Riemann surfaces and Conformal Field Theory