We discuss an efficient preconditioner and iterative numerical method to solve large complex linear algebraic systems of the form (W + iT)u = c, where W and T are symmetric matrices, and at least one of them is nonsingular. When the real part W is dominantly stronger or weaker than the imaginary part T, we propose a block multiplicative (BM) preconditioner or its variant (VBM), respectively. The BM and VBM preconditioned iteration methods are shown to be parameter-free, in terms of eigenvalue distributions of the preconditioned matrix. Furthermore, when the relationship between W and T is obscure, we propose a new preconditioned BM method (PBM) to overcome this difficulty. Both convergent properties of these new iteration methods and spectral properties of the corresponding preconditioned matrices are discussed. The optimal value of iteration parameter for the PBM method is determined. Numerical experiments involving the Helmholtz equation and some other applications show the effectiveness and robustness of the proposed preconditioners and corresponding iterative methods.

Hexagonal grids are valuable in two-dimensional applications involving Laplacian. The methods and analysis are investigated in current work in both linear and nonlinear problems related to anisotropic Laplacian. Ordinary and compact hexagonal grid finite difference methods are developed by elementary arguments, and then analyzed by perturbation for standard Laplacian. In the anisotropic case, analysis is done through reduction to the standard one by using Fourier vectors of mixed types. These hexagonal seven-point methods, with established theoretic stabilities and accuracies, are numerically confirmed in linear and semi-linear anisotropic Poisson problems, and can be applied also in time-dependent problems and in many applications in two-dimensional irregular domains.

In this paper we generalize the notion of the C-numerical range of a matrix to operators in arbitrary tracial von Neumann algebras. For each self-adjoint operator C, the C-numerical range of such an operator is defined; it is a compact, convex subset of ℂ. We explicitly describe the C-numerical ranges of several operators and classes of operators.

The distribution for eigenvalues of Schur complement of matrices plays an important role in many mathematical problems. In this paper, we firstly present some criteria for H-matrix. Then as application, for two class matrices whose submatrices are γ-diagonally dominant and product γ-diagonally dominant, we show that the eigenvalues of the Schur complement are located in the Geršgorin discs and the Ostrowski discs of the original matrices under certain conditions.

In this paper, we consider two innovative structured matrices, CUPL-Toeplitz matrix and CUPL-Hankel matrix. The inverses of CUPL-Toeplitz and CUPL-Hankel matrices can be expressed by the Gohberg-Heinig type formulas, and the stability of the inverse matrices is verified in terms of 1-, ∞- and 2-norms, respectively. In addition, two algorithms for the inverses of CUPL-Toeplitz and CUPL-Hankel matrices are given and examples are provided to verify the feasibility of these algorithms.

There has been a lot of study on the SOR-like methods for solving the augmented system of linear equations since the outstanding work of Golub, Wu and Yuan (BIT 41(2001)71-85) was presented fifteen years ago. Based on the SOR-like methods, we establish a class of accelerated SOR-like methods for large sparse augmented linear systems by making use of optimization technique, which will find the optimal relaxation parameter ω by optimization models. We demonstrate the convergence theory of the new methods under suitable restrictions. The numerical examples show these methods are effective.

For symmetric eigenvalue problems, we constructed a three-term recurrence polynomial filter by means of Chebyshev polynomials. The new filtering technique does not need to solve linear systems and only needs matrix-vector products. It is a memory conserving filtering technique for its three-term recurrence relation. As an application, we use this filtering strategy to the Davidson method and propose the filtered-Davidson method. Through choosing suitable shifts, this method can gain cubic convergence rate locally. Theory and numerical experiments show the efficiency of the new filtering technique.

In this paper, we analyse the convergence rates of several different predictor-corrector iterations for computing the minimal positive solution of the nonsymmetric algebraic Riccati equation arising in transport theory. We have shown theoretically that the new predictor-corrector iteration given in [Numer. Linear Algebra Appl., 21 (2014), pp. 761–780] will converge no faster than the simple predictor-corrector iteration and the nonlinear block Jacobi predictor-corrector iteration. Moreover the last two have the same asymptotic convergence rate with the nonlinear block Gauss-Seidel iteration given in [SIAM J. Sci. Comput., 30 (2008), pp. 804–818]. Preliminary numerical experiments have been reported for the validation of the developed comparison theory.

We propose and investigate a novel solution strategy to efficiently and accurately compute approximate solutions to semilinear optimal control problems, focusing on the optimal control of phase field formulations of geometric evolution laws. The optimal control of geometric evolution laws arises in a number of applications in fields including material science, image processing, tumour growth and cell motility. Despite this, many open problems remain in the analysis and approximation of such problems. In the current work we focus on a phase field formulation of the optimal control problem, hence exploiting the well developed mathematical theory for the optimal control of semilinear parabolic partial differential equations. Approximation of the resulting optimal control problemis computationally challenging, requiring massive amounts of computational time and memory storage. The main focus of this work is to propose, derive, implement and test an efficient solution method for such problems. The solver for the discretised partial differential equations is based upon a geometric multigrid method incorporating advanced techniques to deal with the nonlinearities in the problem and utilising adaptive mesh refinement. An in-house two-grid solution strategy for the forward and adjoint problems, that significantly reduces memory requirements and CPU time, is proposed and investigated computationally. Furthermore, parallelisation as well as an adaptive-step gradient update for the control are employed to further improve efficiency. Along with a detailed description of our proposed solution method together with its implementation we present a number of computational results that demonstrate and evaluate our algorithms with respect to accuracy and efficiency. A highlight of the present work is simulation results on the optimal control of phase field formulations of geometric evolution laws in 3-D which would be computationally infeasible without the solution strategies proposed in the present work.

In this paper we establish operator quasilinearity properties of some functionals associated with Davis–Choi–Jensen’s inequality for positive maps and operator convex or concave functions. Applications for the power function and the logarithm are provided.

We present DASHMM, a general library implementing multipole methods (including both Barnes-Hut and the Fast Multipole Method). DASHMM relies on dynamic adaptive runtime techniques provided by the HPX-5 system to parallelize the resulting multipole moment computation. The result is a library that is easy-to-use, extensible, scalable, efficient, and portable. We present both the abstractions defined by DASHMM as well as the specific features of HPX-5 that allow the library to execute scalably and efficiently.

We introduce some new refinements of numerical radius inequalities for Hilbert space invertible operators. More precisely, we prove that if $T\in {\mathcal{B}}({\mathcal{H}})$ is an invertible operator, then $\Vert T\Vert \leq \sqrt{2}\unicode[STIX]{x1D714}(T)$ .

Assume that $A$ is a finite-dimensional algebra over some field, and assume that $A$ is weakly symmetric and indecomposable, with radical cube zero and radical square nonzero. We show that such an algebra of wild representation type does not have a nonprojective module $M$ whose ext-algebra is finite dimensional. This gives a complete classification of weakly symmetric indecomposable algebras which have a nonprojective module whose ext-algebra is finite dimensional. This shows in particular that existence of ext-finite nonprojective modules is not equivalent with the failure of the finite generation condition (Fg), which ensures that modules have support varieties.

This paper proposes improvements to the modified Fletcher–Reeves conjugate gradient method (FR-CGM) for computing $Z$ -eigenpairs of symmetric tensors. The FR-CGM does not need to compute the exact gradient and Jacobian. The global convergence of this method is established. We also test other conjugate gradient methods such as the modified Polak–Ribière–Polyak conjugate gradient method (PRP-CGM) and shifted power method (SS-HOPM). Numerical experiments of FR-CGM, PRP-CGM and SS-HOPM show the efficiency of the proposed method for finding $Z$ -eigenpairs of symmetric tensors.

We present RECFMM, a program representation and implementation of a recursive scheme for parallelizing the adaptive fast multipole method (FMM) on shared-memory computers. It achieves remarkable high performance while maintaining mathematical clarity and flexibility. The parallelization scheme signifies the recursion feature that is intrinsic to the FMM but was not well exploited. The program modules of RECFMM constitute a map between numerical computation components and advanced architecture mechanisms. The mathematical structure is preserved and exploited, not obscured nor compromised, by parallel rendition of the recursion scheme. Modern software system—CILK in particular, which provides graph-theoretic optimal scheduling in adaptation to the dynamics in parallel execution—is employed. RECFMM supports multiple algorithm variants that mark the major advances with low-frequency interaction kernels, and includes the asymmetrical version where the source particle ensemble is not necessarily the same as the target particle ensemble. We demonstrate parallel performance with Coulomb and screened Coulomb interactions.

The anisotropic perfectly matched layer (APML) defines a continuous vector field outside a rectangle domain and performs the complex coordinate stretching along the vector field. Inspired by [Z. Chen et al., Inverse Probl. Imag., 7, (2013):663–678] and based on the idea of the shortest distance, we propose a new approach to construct the vector field which still allows us to prove the exponential decay of the stretched Green function without the constraint on the thickness of the PML layer. Moreover, by using the reflection argument, we prove the stability of the PML problem in the PML layer and the convergence of the PML method. Numerical experiments are also included.

A lopsided alternating direction iteration (LADI) method and an induced block diagonal preconditioner for solving block two-by-two generalised saddle point problems are presented. The convergence of the LADI method is analysed, and the block diagonal preconditioner can accelerate the convergence rates of Krylov subspace iteration methods such as GMRES. Our new preconditioned method only requires a solver for two linear equation sub-systems with symmetric and positive definite coefficient matrices. Numerical experiments show that the GMRES with the new preconditioner is quite effective.

An efficient high order numerical method is presented to solve the mobile-immobile advection-dispersion model with the Coimbra time variable-order fractional derivative, which is used to simulate solute transport in watershed catchments and rivers. On establishing an efficient recursive algorithm based on the properties of Jacobi polynomials to approximate the Coimbra variable-order fractional derivative operator, we use spectral collocation method with both temporal and spatial discretisation to solve the time variable-order fractional mobile-immobile advection-dispersion model. Numerical examples then illustrate the effectiveness and high order convergence of our approach.

In this paper a determinant identity is established, from which a simple proof of the multivariate Lagrange–Good inversion formula follows directly. Further discussion on a discrete analogue of the Lagrange–Good inversion formula is also presented.

We investigate the central moments of (regular) hexagons and derive accordingly a discrete approximation to definite integrals on hexagons. The seven-point cubature rule makes use of interior and neighbor center nodes, and is of fourth order by construction. The result is expected to be useful in two-dimensional (open-field) applications of integral equations or image processing.