1. A computer implemented method of performing projection based polynomial fitting, the method comprising: generating a plurality of sampling points as a function of variables; forming a polynomial model template representative of the plurality of sampling points, wherein the polynomial model template comprises at least one polynomial coefficient; and forming a low-rank matrix to approximate the polynomial coefficient.