We introduce a multiple interval Chebyshev-Gauss-Lobatto spectral collocation method for the initial value problems of the nonlinear ordinary differential equations (ODES). This method is easy to implement and possesses the high order accuracy. In addition, it is very stable and suitable for long time calculations. We also obtain the hp-version bound on the numerical error of the multiple inter...

A new adaptive quadrature algorithm that places a greater emphasis on cost reduction while still maintaining an acceptable accuracy is demonstrated. The different needs of science and engineering applications are highlighted as the existing algorithms are shown to be inadequate. The performance of the new algorithm is compared with the well known adaptive Simpson, Gauss-Lobatto and GaussKronrod...

Gauss-Lobatto quadrature formulae associated with symmetric weight functions are considered. The kernel of the remainder term for classes of analytic functions is investigated on elliptical contours. Sufficient conditions are found ensuring that the kernel attains its maximal absolute value at the intersection point of the contour with either the real or the imaginary axis. The results obtained...

The classical overlapping Schwarz algorithm is here extended to the triangular/tetrahedral spectral element (TSEM) discretization of elliptic problems. This discretization, based on Fekete nodes, is a generalization to nontensorial elements of the tensorial Gauss–Lobatto–Legendre quadrilateral spectral elements (QSEM). The overlapping Schwarz preconditioners are based on partitioning the domain...

This paper is concerned with the hp-version of the finite element method (hp-FEM) to treat a variational inequality that models frictional contact in linear elastostatics. Such an approximation of higher order leads to a nonconforming discretization scheme. We employ Gauss-Lobatto quadrature for the approximation of the nonsmooth frictiontype functional and take the resulting quadrature error i...