For a Schwartz function $f$ on $\mathbb R$, let $\Gamma(f,s)=\int_{\mathbb R^\times} |x|^sf(x){dx\over |x|}$. The usual Gamma factor is obtained by taking a Gaussian. The local functional equation (proven by changing variables in the defining integrals, in the range $0<{\rm Re}(s)<1$), is
$$\Gamma(f,s)\cdot \Gamma(\hat{g},1-s)= \Gamma(\hat{f},1-s)\cdot \Gamma(g,s)$$
for any two Schwartz functions $f,g$. And Riemann's argument proves
$$
\Gamma(f,s)\cdot \zeta(s) = \Gamma(\hat{f},1-s)\cdot \zeta(1-s)