One of the things I'm working on has required me to look into the literature of multidimensional theta functions, and I've gotten a bit confused on a few naming details.

A look at the DLMF says that "the" multidimensional theta function is the Riemann theta function,

$$\mathop{\Theta}\left(\mathbf{z}\mid\boldsymbol{{\Omega}}\right)=\sum_{{\mathbf{n}\in{\mathbb Z}^g}}e^{{2\pi i\left(\frac{1}{2}\mathbf{n}\cdot\boldsymbol{{\Omega}}\cdot\mathbf{n}+\mathbf{n}\cdot\mathbf{z}\right)}}$$

with $\mathbf{z}$ a complex vector and $\mathbf \Omega$ a complex symmetric matrix with symmetric positive definite imaginary part. I find that there is a nice implementation for numerically computing this function, but it is only available in Maple and Java. Since I work in *Mathematica*, I decided to see if there was already something in *Mathematica* before going through the trouble of translating code.

Peering at the docs for *Mathematica* netted the Siegel theta function. However, looking at how the Siegel theta function was defined, this and the Riemann theta function *seem* to be the same thing!

Not having Maple to check if the results from their implementation of Riemann theta and the results from *Mathematica*'s implementation of Siegel theta agree, and not being able to access the (older) references to these functions, I now wish to ask: is there no difference whatsoever between these two except in name (and if there are, how are these two different)? Did Riemann and Siegel independently study the exact same multidimensional theta function? How come it seems that there is no reference to these two being synonymous?