Faster exact substrate thermorefractive noise calculation

For the "exact" calculation of substrate thermorefractive noise, pygwinc currently implements quadrature integration over wavevectors following P1400084 Eq. 17. However, BV04 Eq. E7 gives the solution in terms of a special function that is implemented in scipy.special. The special function expression takes about 2 ms to evaluate, compared to 600 ms for the current quadrature integration.

gwinc/noise/substratethermal.py

@@ -34,13 +34,11 @@ def substrate_thermorefractive(f, materials, wBeam, exact=False):
     omega = 2*pi*f
 
     if exact:
-        # arXiv:cond-mat/0402650, Eq. E7
         w = omega * r0**2 * rho * C / (2 * kappa)
         psd = np.abs(H * beta**2 * kBT * Temp / (2 * pi * kappa) * (exp(1j*w) * scipy.special.exp1(1j*w)
             + exp(-1j*w) * scipy.special.exp1(-1j*w)))
 
     else:
-        # arXiv:cond-mat/0402650, Eq. 5.3; and P1400084, Eq. 18
         psd = 4*H*beta**2*kappa*kBT*Temp/(pi*r0**4*omega**2*(rho*C)**2)
 
     return psd