In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of convection-diffusion type of third order Ordinary Differential Equations (ODEs) in which the SPBVP is reduced into a weakly coupled system of two ODEs subject to suitable initial and boundary conditions. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. In order to get a numerical solution for the derivative of the solution, the domain is divided into two regions namely inner region and outer region. The shooting method is applied to the inner region while standard finite difference scheme (FD) is applied for the outer region. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing.

We present quantitative studies of transfer operators between finite element spaces associated with unrelated meshes. Several local approximations of the global L2-orthogonal projection are reviewed and evaluated computationally. The numerical studies in 3D provide the first estimates of the quantitative differences between a range of transfer operators between non-nested finite element spaces. We consider the standard finite element interpolation, Clément’s quasi-interpolation with different local polynomial degrees, the global L2-orthogonal projection, a local L2-quasi-projection via a discrete inner product, and a pseudo-L2-projection defined by a Petrov-Galerkin variational equation with a discontinuous test space. Understanding their qualitative and quantitative behaviors in this computational way is interesting per se; it could also be relevant in the context of discretization and solution techniques which make use of different non-nested meshes. It turns out that the pseudo-L2-projection approximates the actual L2-orthogonal projection best. The obtained results seem to be largely independent of the underlying computational domain; this is demonstrated by four examples (ball, cylinder, half torus and Stanford Bunny).

For the approximation in Lp-norm, we determine the weakly asymptotic orders for the simultaneous approximation errors of Sobolev classes by piecewise cubic Hermite interpolation with equidistant knots. For p = 1, ∞, we obtain its values. By these results we know that for the Sobolev classes, the approximation errors by piecewise cubic Hermite interpolation are weakly equivalent to the corresponding infinite-dimensional Kolmogorov widths. At the same time, the approximation errors of derivatives are weakly equivalent to the corresponding infinite-dimensional Kolmogorov widths.

In this paper, we study to use nonlocal bounded variation (NLBV) techniques to decompose an image intensity into the illumination and reflectance components. By considering spatial smoothness of the illumination component and nonlocal total variation (NLTV) of the reflectance component in the decomposition framework, an energy functional is constructed. We establish the theoretical results of the space of NLBV functions such as lower semicontinuity, approximation and compactness. These essential properties of NLBV functions are important tools to show the existence of solution of the proposed energy functional. Experimental results on both grey-level and color images are shown to illustrate the usefulness of the nonlocal total variation image decomposition model, and demonstrate the performance of the proposed method is better than the other testing methods.

In this paper we present a first supercloseness analysis for higher-order Galerkin FEM applied to a singularly perturbed convection-diffusion problem. Using a solution decomposition and a special representation of our finite element space, we are able to prove a supercloseness property of p + 1/4 in the energy norm where the polynomial order p ≥ 3 is odd.

In this paper, the G2 interpolation by Pythagorean-hodograph (PH) quintic curves in ℝd, d ≥ 2, is considered. The obtained results turn out as a useful tool in practical applications. Independently of the dimension d, they supply a G2 quintic PH spline that locally interpolates two points, two tangent directions and two curvature vectors at these points. The interpolation problem considered is reduced to a system of two polynomial equations involving only tangent lengths of the interpolating curve as unknowns. Although several solutions might exist, the way to obtain the most promising one is suggested based on a thorough asymptotic analysis of the smooth data case. The numerical algorithm traces this solution from a particular set of data to the general case by a homotopy continuation method. Numerical examples confirm the efficiency of the proposed method.

In this paper, we study shifted restated full orthogonalization method with deflation for simultaneously solving a number of shifted systems of linear equations. Theoretical analysis shows that with the deflation technique, the new residual of shifted restarted FOM is still collinear with each other. Hence, the new approach can solve the shifted systems simultaneously based on the same Krylov subspace. Numerical experiments show that the deflation technique can significantly improve the convergence performance of shifted restarted FOM.