Many applications based on finite element and finite difference methods include the solution of large sparse linear systems using preconditioned iterative methods. Matrix vector multiplication is one of the key operations that has a significant impact on the performance of any iterative solver. In this paper, recent developments in sparse storage formats on vector machines are reviewed. Then, s...

After a brief survey on the parametric deformable models, we develop an iterative method based on the finite difference schemes in order to obtain energy-minimizing snakes. We estimate the approximation error, the residue, and the truncature error related to the corresponding algorithm, then we discuss its convergence, consistency, and stability. Some aspects regarding the prosthetic sugical me...

The finite volume method combines the advantage of the finite element method and the finite difference method, while overcoming their shortcomings. This paper describes the development of the finite volume method and the process of using finite volume method for fluid calculations. Finally, this paper introduces the process of using improved finite volume method for solving the flow equations o...

We analyze the discretization of initial and boundary value problems with a stationary interface in one space dimension for the heat equation, the Schrödinger equation, and the wave equation by finite difference methods. Extending the concept of the elliptic projection, well known from the analysis of Galerkin finite element methods, to our finite difference case, we prove second-order error es...

In this paper, we derive a highly accurate numerical method for the solution of one-dimensional wave equation with Neumann boundary conditions. This hyperbolic problem is solved by using semidiscrete approximations. The space direction is discretized by wavelet-Galerkin method and the time variable is discretized by using various classical finite difference schemes. The numerical results show t...

We study effects of finite height and surrounding material on photonic crystal slabs of oneand two-dimensional photonic crystals with a pseudo-spectral method and finite difference time domain simulation methods. The band gap is shown to be strongly modified by the boundary material. As an application we suggest reflection and guiding of light by patterning the material on top/below the slab.

A technique to incorporate directivity in the Pseudospectral time-domain method (PSTD) for auralization purposes is presented. The Pseudospectral time-domain method is a numerical method for solving the wave equation that requires only two spatial points per wavelength, making it more efficient compared to other wave based methods. The effect of source directivity as well as the head-related di...

by, for example, a finite difference method in the computational (μ, χ, φ) space. The above standard dipole coordinates is convenient and certainly appropriate for analytical studies in which the Earth’s dipolar field plays central roles. It also works as a base coordinates for the node and cell generation of the finite element method in the dipole geometry [2, 3]. However, when one tries to us...