“…In addition, the method reduces the computational work of evaluating and assembling the stiffness matrix, which often constitutes a major portion of the overall computational work, from O(N 2 ) to O(N ). In fact, the fast collocation method is nothing but the conventional one, and so naturally inherits the stability, convergence, and other numerical advantages of conventional collocation methods [12,22,30,33,36,42,52]. Furthermore, any numerical techniques developed for traditional collocation methods, such as the numerical quadra-tures for the accurate evaluation of singular integrals [2,18,26,27,29,36,51], can be applied to the fast collocation method.…”