I am trying to interpolate a function defined over a three-dimensional real space: $$f: R^3\rightarrow R\\(x,y,z)\rightarrow f(x,y,z)$$ Let assume I have $N_1 N_2 N_3$ points in the space which form my grid for this interpolation, and the multivariate series $$F(x,y,z)=\sum_{a=0}^{N_1-1} \sum_{b=0}^{N_2-1} \sum_{c=0}^{N_3-1} C_{abc} x^a y^b z^c$$ is the chosen interpolator. In order to find the coefficients I should form the following sets of equations: $$\sum_{a=0}^{N_1-1} \sum_{b=0}^{N_2-1} \sum_{c=0}^{N_3-1} C_{abc} x_i^a y_i^b z_i^c=f_i=f(x_i,y_i,z_i)\,\mbox{, where, }i=1,\ldots,N_1 N_2 N_3\,,$$ and then solve for the coefficients. However, for this I first need to write the above set of linear equations in the standard form $$A\vec{x}=\vec{b}$$ wherein, $A$ is the matrix of coefficients, $\vec{x}$ is the vector of unknowns, and $b=\{f_i\}_1^{N_1 N_2 N_3}$ is the known vector. For this to be done I would require to expand the multivariate power series in the form of a single variable series, that is, $$\sum_{a=0}^{N_1-1} \sum_{b=0}^{N_2-1} \sum_{c=0}^{N_3-1} C_{abc} x_i^a y_i^b z_i^c=\sum_{m=0}^{(N_1-1)(N_2-1)(N_3-1)} C_m \eta_i^m$$ wherein $\eta_i=\eta_i(x_i,y_i,z_i)$ and probably $m=a b c$. OF course the expansion needs to be nontrivial and useful.
Is it possible at all? Any suggestion to find the coefficients more practically?
Regards, owari
UPDATE.
Maybe it appears that the most natural way for solving this problem is comprised of the following steps:
first solve for the coefficients of $$\sum_{a=0}^{N_1-1} E_a(y,z) x_i^a=f(x_i,y,z)\,\mbox{, where, }i=1,\ldots,N_1\,,$$
then solve for the coefficients of $$\sum_{b=0}^{N_2-1} D_{ab}(z) y_i^b=E_a(y_i,z)\,\mbox{, where, }i=1,\ldots,N_2\,\mbox{, and, }a=0,\ldots,N_1-1$$
and finally solve for the coefficients of $$\sum_{c=0}^{N_2-1} C_{abc} z_i^c=D_{ab}(z_i)\,\mbox{, where, }i=1,\ldots,N_3\,,\,a=0,\ldots,N_1-1\mbox{, and, }b=0,\ldots,N_2-1$$
which gives the coefficients we were looking for. However, this way, the first two steps will be working with function-type coefficients instead of numerical coefficients and that will prevent efficient usage of the available codes in numerical analysis. Solving for each set of coefficients at each step for different grid points would also inevitably increase the number of equations drastically, so any better suggestion to solve for this problem?