IG-LSPIA: Least Squares Progressive Iterative Approximation for Isogeometric Collocation Method

The isogeometric collocation method (IGA-C), which is a promising branch of isogeometric analysis (IGA), can be considered fitting the load function with the combination of the numerical solution and its derivatives. In this study, we develop an iterative method, isogeometric least-squares progressive-iterative approximation (IG-LSPIA), to solve the fitting problem in the collocation method. IG-LSPIA starts with an initial blending function, where the control coefficients are combined with the B-spline basis functions and their derivatives. A new blending function is generated by constructing the differences for collocation points (DCP) and control coefficients (DCC), and then adding the DCC to the corresponding control coefficients. The procedure is performed iteratively until the stop criterion is reached. We prove the convergence of IG-LSPIA and show that the computation complexity in each iteration of IG-LSPIA is related only to the number of collocation points and unrelated to the number of control coefficients. Moreover, an incremental algorithm is designed; it alternates with knot refinement until the desired precision is achieved. After each knot refinement, the result of the last round of IG-LSPIA iterations is used to generate the initial blending function of the new round of iteration, thereby saving great computation. Experiments show that the proposed method is stable and efficient. In the three-dimensional case, the total computation time is saved twice compared to the traditional method.