Iteratively computing a limiting function and its consecutive derivatives, a Hermite subdivision scheme is of particular interest and importance in CAGD for generating smooth subdivision curves and in numerical PDEs for constructing Hermite wavelets. This paper characterizes the convergence and smoothness of a Hermite subdivision scheme, provides simple factorization of Hermite masks through the normal form of matrix-valued masks, and presents an algorithm for constructing all masks for Hermite subdivision schemes.

We introduce bivariate dual sqrt{2}-subdivision schemes using 1D stencils. Such subdivision schemes enjoy the quasi-interpolating property (i.e., they interpolate bivariate polynomials of high orders) and are intrinsically linked to both 1D primal interpolating subdivision schemes and 1D masks having linear-phase moments. Using only 1D stencils, such subdivision schemes can be straightforwardly implemented on any quadrilateral meshes and there is no need to design special subdivision rules near extraordinary or boundary vertices.

Using newly modified cubic B-splines with accuracy order 4, we proposed a numerical scheme for solving the 1D and 2D hyperbolic telegraph equation using a differential quadrature method. Our modified cubic B-splines retain the tridiagonal structure and achieve the fourth order convergence rate. The solution of the associated ODEs is advanced in the time domain by the SSPRK scheme.

Numerically solving the 2D Helmholtz equation is widely known to be \ very difficult largely due to its highly oscillatory solution. A very fine mesh\ size is necessary to deal with a large wavenumber leading to a severely ill-co\ nditioned huge coefficient matrix. To understand and tackle such challenges, it\ is crucial to analyze how the solution of the 2D Helmholtz equation depends on\ boundary and source data for large wavenumbers. In this paper, we analyze and\ derive several new sharp wavenumber-explicit stability bounds for the 2D Helmh\ oltz equation with inhomogeneous mixed boundary conditions: Dirichlet, Neumann,\ and impedance.

The elliptic interface problems with discontinuous and high-contrast coefficients have many applications but are very challenging to numerically solve, because their solutions often have low regularity or even discontinuity (which result in huge highly ill-conditioned linear systems). Thus, it is highly desired to construct high order schemes to solve the elliptic interface problems with discontinuous and high-contrast coefficients. To solve such elliptic interface problem with a smooth interface curve, we propose a fourth order compact finite difference method for computing the solution and its gradient by using uniform meshes. We showe that our proposed FDM numerically satisfies the discrete maximum principle when the mesh size is small and then we use it to prove their convergence of four order.