For $z\in \mathbb{C}$ and $\tau \in i\mathbb{R}_{+}$ consider the function $$ \theta_3(z;\tau)=\sum_{n\in\mathbb{Z}}\exp\left(2\pi i n z +\pi i \tau n^2\right) $$ This function satisfies the well known quasiperiodicity properties $$\theta_3(z+k;\tau)=\theta_3(z;\tau)$$ and $$\theta_3(z+k\tau;\tau)=\exp\left(-\pi i \tau k^2-2\pi i k z \right)\theta_3(z;\tau) $$ Asymptotics in the parameter $\tau$ are well known ( and easily found through a search engine ). Namely
$$\theta_3(z;\tau) \to 1 \text{, } \Im(\tau)\to +\infty$$ and $$ \theta_3(z;\tau) \sim \sqrt{\frac{1}{-i\tau}}\exp\left( \frac{z^2}{i\tau}\right) \text{, } \Im(\tau)\to 0^{+} $$
Asymptotics in $z$ are a little less straightforward.
Since for any $k\in \mathbb{Z}$, $$\theta_3(z+k;\tau)=\theta_3(z;\tau)$$ holds, if $z=x$ is real there is no growth as $x\to \pm \infty$, just oscillations.
The question is how to apply the second identity to compute an asymptotic formula for $\theta_3(\pm ix;\tau)$ as $x\to \infty$.
Writing $ix=\epsilon_x \tau +k_x\tau$ with $k_x=\lfloor \frac{ix}{\tau} \rfloor$ and $\epsilon_x=\frac{ix}{\tau}-k$, we immediately have that $$\theta_3(\pm ix;\tau)=\exp\left( -\pi i k_x^2 \pm 2\pi k_x \epsilon_x \tau \right)\theta_3(\epsilon_x \tau;\tau)$$
The right most term is bounded and well defined. My calculations on the exponential term have been somewhat lengthy and did not give a nice formula.
I would like to, first, invite suggestions for references, this seems to be too simple to not be published somewhere!
If not, then, any suggestions on how to deal with the floor terms would be much appreciated.