From Nekrasov functions to theta functions

I'm going to demonstrate how to take an autonomous limit of the general solution of some (q-)isomonodromic system. Such solutions are usually given as Fourier transformations of Nekrasov functions. One can show using Seiberg-Witten equations that in the autonomous limit these Fourier transformations turn into Riemann theta functions, and thus satisfy Fay bilinear relations