I have some difficulties with expressing the following series $\sum\limits_{b=-\infty}^{+\infty}\sum\limits_{a=-\infty}^{+\infty}q^{1 + 3 a + 3 a^2 - 3 b - 3 a b + 3 b^2}$ using standart theta functions. I've tried to get the degree for the sum of squares: $1 + 3 a + 3 a^2 - 3 b - 3 a b + 3 b^2=3t^2+(\frac{3b-1}{2})^2$, where $a=t+\frac{b-1}{2}$ so as that the initial sum splits in $\sum\limits_{b=-\infty}^{+\infty}\sum\limits_{a=-\infty}^{+\infty}q^{1 + 3 a + 3 a^2 - 3 b - 3 a b + 3 b^2}$ with $b$ even and $\sum\limits_{t=-\infty}^{+\infty}\sum\limits_{b=-\infty}^{+\infty}q^{3t^2+(\frac{3b-1}{2})^2}$ , where $b$ is odd and then I tried analyzing cases with $\frac{3b-1}{2}\equiv 1 (mod3)$ and $\frac{3b-1}{2}\equiv 2 (mod3)$, but I were unsucceded.
It is also clear how to express $\sum\limits_{b=-\infty}^{+\infty}\sum\limits_{a=-\infty}^{+\infty}q^{3 a^2 - 3 a b + 3 b^2}$, so that we may decide it as being standart.
I will be really grateful for any suggestions.
Some useful definitions:
$\theta_2(z)=\sum\limits_{a=-\infty}^{+\infty}q^{(a+1/2)^2}$, where $q=e^{\pi i z}$;
$\theta_3(z)=\sum\limits_{a=-\infty}^{+\infty}q^{a^2}$;
$\theta_4(z)=\sum\limits_{a=-\infty}^{+\infty}(-q)^{a^2}$.
Some extra comments:
Let $\phi(z)=\sum\limits_{b=-\infty}^{+\infty}\sum\limits_{a=-\infty}^{+\infty}q^{a^2 - a b + b^2}$, and assume that $a^2 - a b + b^2=(a-\frac{b}{2})^2+\frac{3b^2}{4}$, then
$\sum\limits_{b=-\infty}^{+\infty}\sum\limits_{a=-\infty}^{+\infty}q^{(a-\frac{b}{2})^2+\frac{3b^2}{4}}=\sum\limits_{a=-\infty}^{+\infty}\sum\limits_{\substack{b=-\infty \\ \text{b even}}}^{+\infty}q^{(a-\frac{b}{2})^2+\frac{3b^2}{4}}+ \sum\limits_{a=-\infty}^{+\infty}\sum\limits_{\substack{b=-\infty \\ \text{b odd}}}^{+\infty}q^{(a-\frac{b}{2})^2+\frac{3b^2}{4}}= \sum\limits_{a=-\infty}^{+\infty}\sum\limits_{k=-\infty}^{+\infty}q^{(a-k)^2+3k^2} +\sum\limits_{a=-\infty}^{+\infty}\sum\limits_{k=-\infty}^{+\infty}q^{(a-k-1/2)^2+3(k+1/2)^2}= \sum\limits_{t=-\infty}^{+\infty}\sum\limits_{t=-\infty}^{+\infty}q^{t^2+3k^2}+ \sum\limits_{a=-\infty}^{+\infty}\sum\limits_{t=-\infty}^{+\infty}q^{(t+1/2)^2+3(k+1/2)^2}=\theta_2(z)\theta_2(3z)+\theta_3(z)\theta_3(3z)$.
And then $\sum\limits_{b=-\infty}^{+\infty}\sum\limits_{a=-\infty}^{+\infty}q^{3(a^2 - a b + b^2)}=\phi(3z)$.