Suppose $K:[0,1]^2 \times [0,1]^2 \to \mathbb{R}$ is continuous and positive semi-definite, and define the corresponding Hilbert-Schmidt integral operator
$$ [Cu](x) = \int_{[0,1]^2} K(x,y) u(y) dy.$$
Suppose further that $K(x,y)=k_1(x_1,y_1) k_2(x_2,y_2)$, and that one has a complete set of eigenfunctions and eigenvalues $\{ (f_n^{i},\lambda_n^{i}) \}_{n=1}^\infty$for the corresponding one-dimensional operators
$$ [C_i v](x_i) = \int_0^1 k_i(x_i,y_i)v(y_i)dy_i$$
for $i=1,2$. Clearly $u_{ij}(x)=f_i^1(x_1)f_j^2(x_2)$ is an eigenfunction of $C$ with eigenvalue $\lambda_i^1 \lambda_j^2$ for any $i,j \in \mathbb{N}$, due to the structure of the kernel and the decoupling of the integral.
My question is, does this exhaust all possibilities? Namely, can every eigenfunction of $C$ be obtained in this way? Impose any conditions that you'd like... for me, $k_1$ and $k_2$ are continuous and piecewise smooth. A reference would be greatly appreciated, as well as any responses... also I am happy to clarify if anything is unclear. Thanks in advance!