The $q-$Onsager algebra and multivariable $q-$special functions

Two sets of mutually commuting $q-$difference operators $x_i$ and $y_j$, $i,j=1, ...,N$ such that $x_i$ and $y_i$ generate a homomorphic image of the $q-$Onsager algebra for each $i$ are introduced. The common polynomial eigenfunctions of each set are found to be entangled product of elementary Pochhammer functions in $N$ variables and $N+3$ parameters. Under certain conditions on the parameters, they form two `dual' bases of polynomials in $N$ variables. The action of each operator with respect to its dual basis is block tridiagonal. The overlap coefficients between the two dual bases are expressed as entangled products of $q-$Racah polynomials and satisfy an orthogonality relation. The overlap coefficients between either one of these bases and the multivariable monomial basis are also considered. One obtains in this case entangled products of dual $q-$Krawtchouk polynomials. Finally, the `split' basis in which the two families of operators act as block bidiagonal matrices is also provided.