The Toda equation and its variants are studied in the filed of integrable systems. One particularly generalized time discretisation of the Toda equation is known as the discrete hungry Toda (dhToda) equation, which has two main variants referred to as the dhTodaI equation and dhTodaII equation. The dhToda equations have both been shown to be applicable to the computation of eigenvalues of totally nonnegative (TN) matrices, which are matrices without negative minors. The dhTodaI equation has been investigated with respect to the properties of integrable systems, but the dhTodaII equation has not. Explicit solutions using determinants and matrix representations called Lax pairs are often considered as symbolic properties of discrete integrable systems. In this paper, we clarify the determinant solution and Lax pair of the dhTodaII equation by focusing on an infinite sequence. We show that the resulting determinant solution firmly covers the general solution to the dhTodaII equation, and provide an asymptotic analysis of the general solution as discrete-time variable goes to infinity.

We consider the unique solvability of an inverse-source problem with integral transmitting condition for a time-fractional mixed type equation in rectangular domain where the unknown source term depends only on the space variable. The solution is based on a series expansion using a bi-orthogonal basis in space, corresponding to a non-self-adjoint boundary value problem. Under certain regularity conditions on the given data, we prove the uniqueness and existence of the solution for the given problem. The influence of the transmitting condition on the solvability of the problem is also demonstrated. Two different transmitting conditions are considered — viz. a full integral form and a special case. In order to simplify the bulky expressions appearing in the proof of our main result, we establish a new property of the recently introduced Mittag-Leffler type function in two variables.

If the closed wire frame of a soap film having the shape of a Möbius strip is pulled apart and gradually deformed into a planar circle, the soap film transforms into a two-sided orientable surface. In the presence of a finite-time twist singularity, which changes the linking number of the film's Plateau border and the centreline of the wire, the topological transformation involves the collapse of the film toward the wire. In contrast to experimental studies of this process reported elsewhere, we use a numerical approach based on the immersed boundary method, which treats the soap film as a massless membrane in a Navier-Stokes fluid. In addition to known effects, we discover vibrating motions of the film arising after the topological change is completed, similar to the vibration of a circular membrane.

Preconditioned modified Hermitian and skew-Hermitian splitting method (PMHSS) is an unconditionally convergent iteration method for solving large sparse complex symmetric systems of linear equations, and uses one parameter α. Adding another parameter β, the generalized PMHSS method (GPMHSS) is essentially a twoparameter iteration method. In order to accelerate the GPMHSS method, using an unexpected way, we propose an accelerated GPMHSS method (AGPMHSS) for large complex symmetric linear systems. Numerical experiments show the numerical behavior of our new method.

The hybrid variational model for restoration of texture images corrupted by blur and Gaussian noise we consider combines total variation regularisation and a fractional-order regularisation, and is solved by an alternating minimisation direction algorithm. Numerical experiments demonstrate the advantage of this model over the adaptive fractional-order variational model in image quality and computational time.

Retinex theory explains how the human visual system perceives colors. The goal of retinex is to decompose the reflectance and the illumination from the given images and thereby compensating for non-uniform lighting. The existing methods for retinex usually use a single image with a fixed exposure to restore the reflectance of the image. In this paper, we propose a variational model for retinex problem by utilizing multi-exposure images of a given scene. The existence and uniqueness of the solutions of the proposed model have been elaborated. An alternating minimization method is constructed to solve the proposed model and its convergence is also demonstrated. The experimental results show that the proposed method is effective for reflectance recovery in retinex problem.

We propose a robust numerical algorithm for solving the nonlinear eigenvalue problem A(ƛ)x = 0. Our algorithm is based on the idea of finding the value of ƛ for which A(ƛ) is singular by computing the smallest eigenvalue or singular value of A(ƛ) viewed as a constant matrix. To further enhance computational efficiency, we introduce and use the concept of signed singular value. Our method is applicable when A(ƛ) is large and nonsymmetric and has strong nonlinearity. Numerical experiments on a nonlinear eigenvalue problem arising in the computation of scaling exponent in turbulent flow show robustness and effectiveness of our method.

Based on the shift-splitting technique and the idea of Hermitian and skew-Hermitian splitting, a fast shift-splitting iteration method is proposed for solving nonsingular and singular nonsymmetric saddle point problems in this paper. Convergence and semi-convergence of the proposed iteration method for nonsingular and singular cases are carefully studied, respectively. Numerical experiments are implemented to demonstrate the feasibility and effectiveness of the proposed method.

It is well known that numerical derivative contains two types of errors. One is truncation error and the other is rounding error. By evaluating variables with rounding error, together with step size and the unknown coefficient of the truncation error, the total error can be determined. We also know that the step size affects the truncation error very much, especially when the step size is large. On the other hand, rounding error will dominate numerical error when the step size is too small. Thus, to choose a suitable step size is an important task in computing the numerical differentiation. If we want to reach an accuracy result of the numerical difference, we had better estimate the best step size. We can use Taylor Expression to analyze the order of truncation error, which is usually expressed by the big O notation, that is, E(h) = Chk. Since the leading coefficient C contains the factor f(k)(ζ) for high order k and unknown ζ, the truncation error is often estimated by a roughly upper bound. If we try to estimate the high order difference f(k)(ζ), this term usually contains larger error. Hence, the uncertainty of ζ and the rounding errors hinder a possible accurate numerical derivative.

We will introduce the statistical process into the traditional numerical difference. The new method estimates truncation error and rounding error at the same time for a given step size. When we estimate these two types of error successfully, we can reach much better modified results. We also propose a genetic approach to reach a confident numerical derivative.

Based on the relaxed factorization techniques studied recently and the idea of the simple-like preconditioner, a modified relaxed positive-semidefinite and skew-Hermitian splitting (MRPSS) preconditioner is proposed for generalized saddle point problems. Some properties, including the eigenvalue distribution, the eigenvector distribution and the minimal polynomial of the preconditioned matrix are studied. Numerical examples arising from the mixed finite element discretization of the Oseen equation are illustrated to show the efficiency of the new preconditioner.

Some convergence bounds of the minimal residual (MINRES) method are studied when the method is applied for solving Hermitian indefinite linear systems. The matrices of these linear systems are supposed to have some properties so that their spectra are all clustered around ±1. New convergence bounds depending on the spectrum of the coefficient matrix are presented. Some numerical experiments are shown to demonstrate our theoretical results.

Based on a preconditioned shift-splitting of the (1,1)-block of non-Hermitian saddle-point matrix and the Uzawa iteration method, we establish a new Uzawa-type iteration method. The convergence properties of this iteration method are analyzed. In addition, based on this iteration method, a preconditioner is proposed. The spectral properties of the preconditioned saddle-point matrix are also analyzed. Numerical results are presented to verify the robustness and the efficiency of the new iteration method and the preconditioner.

Pinning effect on current-induced magnetic transverse domain wall dynamics in nanostrip is studied for its potential application to new magnetic memory devices. In this study, we carry out a series of calculations by solving generalized Landau-Lifshitz equation involving a current spin transfer torque in one and two dimensional models. The critical current for the transverse wall depinning in nanostrip depends on the size of artificial rectangular defects on the edges of nanostrip. We show that there is intrinsic pinning potential for a defect such that the transverse wall oscillates damply around the pinning site with an intrinsic frequency if the applied current is below critical value. The amplification of the transverse wall oscillation for both displacement and maximum value of m3 is significant by applying AC current and current pulses with appropriate frequency. We show that for given pinning potential, the oscillation amplitude as a function of the frequency of the AC current behaves like a Gaussian distribution in our numerical study, which is helpful to reduce strength of current to drive the transverse wall motion.

The numerical simulations for the blood flow in arteries by high order accurate schemes have a wide range of applications in medical engineering. The blood flow model admits the steady state solutions, in which the flux gradient is non-zero and is exactly balanced by the source term. In this paper, we present a high order finite volume weighted essentially non-oscillatory (WENO) scheme, which preserves the steady state solutions and maintains genuine high order accuracy for general solutions. The well-balanced property is obtained by a novel source term reformulation and discretisation, combined with well-balanced numerical fluxes. Extensive numerical experiments are carried out to verify well-balanced property, high order accuracy, as well as good resolution for smooth and discontinuous solutions.