BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kui Du

Kui Du appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysiscond-mat.mtrl-scicond-mat.mes-hall

Catalog footprint

What is connected

12works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2026arXiv

TriCG with deflated restarting for symmetric quasi-definite linear systems

TriCG is a short-recurrence iterative method recently introduced by Montoison and Orban [SIAM J. Sci. Comput., 43 (2021), pp. A2502--A2525] for solving symmetric quasi-definite (SQD) linear systems. TriCG takes advantage of the inherent block structure of SQD linear systems and performs substantially better than SYMMLQ. However, numerical experiments have revealed that the convergence of TriCG can be notably slow when the off-diagonal block contains a substantial number of large elliptic singular values. To address this limitation, we introduce a deflation strategy tailored for TriCG to improve its convergence behavior. Specifically, we develop a generalized Saunders--Simon--Yip process with deflated restarting to construct the deflation subspaces. Building upon this process, we propose a novel method termed TriCG with deflated restarting. The deflation subspaces can also be utilized to solve SQD linear systems with multiple right-hand sides. Numerical experiments are provided to illustrate the superior performance of the proposed methods.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

GPBiLQ and GPQMR: Two iterative methods for unsymmetric partitioned linear systems

We introduce two iterative methods, GPBiLQ and GPQMR, for solving unsymmetric partitioned linear systems. The basic mechanism underlying GPBiLQ and GPQMR is a novel simultaneous tridiagonalization via biorthogonality that allows for short-recurrence iterative schemes. Similar to the biconjugate gradient method, it is possible to develop another method, GPBiCG, whose iterate (if it exists) can be obtained inexpensively from the GPBiLQ iterate. Whereas the iterate of GPBiCG may not exist, the iterates of GPBiLQ and GPQMR are always well defined as long as the biorthogonal tridiagonal reduction process does not break down. We discuss connections between the proposed methods and some existing methods, and give numerical experiments to illustrate the performance of the proposed methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Preferred corrosion pathways for oxygen in Al2Ca-twin boundaries and dislocations

With an ongoing discussion on the oxygen diffusion along crystal defects remaining, it is difficult to study this phenomenon in Al containing intermetallic materials due to its rapid and passivating oxide formation. We report here the observation of enhanced oxygen diffusion along crystal defects, i.e. dislocations and twin boundaries, in the C15 Al 2 Ca Laves phase and how the presence of oxygen induces structural changes at these defects. Three main phases were identified and characterized structurally by aberration-corrected, atomic resolution scanning transmission electron microscopy, analytically by energy dispersive X-ray spectroscopy and electron energy loss spectroscopy. Unlike the C15 bulk phase, the twin boundary and dislocation transformed into a few nanometer wide amorphous phase, which depletes in Al and Ca but is highly enriched in oxygen. The dislocation even shows coexistence of the amorphous phase with a simple Al-rich A1 fcc phase. This A1 phase only depletes in Ca, not in Al (Al remains at bulk concentration), and is also enriched in oxygen. The Al-rich A1 phase is coherent with the C15 matrix. Electron energy loss spectroscopy revealed the amorphous phase to be Al 2 O 3 . We thereby show as one of the first studies that oxygen diffusion along crystal defects, especially also at the twin boundary can induce the formation of an amorphous oxide along themselves. The identification of oxygen-induced transformation at strained defects has to be considered when the material is exposed to air during plastic deformation at elevated temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A doubly stochastic block Gauss-Seidel algorithm for solving linear equations

We propose a simple doubly stochastic block Gauss--Seidel algorithm for solving linear systems of equations. By varying the row partition parameter and the column partition parameter of the coefficient matrix, we recover the Landweber algorithm, the randomized Kaczmarz algorithm, the randomized Gauss--Seidel algorithm, and the doubly stochastic Gauss--Seidel algorithm. For general (consistent or inconsistent) linear systems, we show the exponential convergence of the {\it norms of the expected iterates} via exact formulas. For consistent linear systems, we prove the exponential convergence of the {\it expected norms of the error and the residual}. Numerical experiments are given to illustrate the efficiency of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Randomized extended block Kaczmarz for solving least squares

Randomized iterative algorithms have recently been proposed to solve large-scale linear systems. In this paper, we present a simple randomized extended block Kaczmarz algorithm that exponentially converges in the mean square to the unique minimum $\ell_2$-norm least squares solution of a given linear system of equations. The proposed algorithm is pseudoinverse-free and therefore different from the projection-based randomized double block Kaczmarz algorithm of Needell, Zhao, and Zouzias. We emphasize that our method works for all types of linear systems (consistent or inconsistent, overdetermined or underdetermined, full-rank or rank-deficient). Moreover, our approach can utilize efficient implementations on distributed computing units, yielding remarkable improvements in computational time. Numerical examples are given to show the efficiency of the new algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stochastic gradient descent for linear least squares problems with partially observed data

We propose a novel stochastic gradient descent method for solving linear least squares problems with partially observed data. Our method uses submatrices indexed by a randomly selected pair of row and column index sets to update the iterate at each step. Theoretical convergence guarantees in the mean square sense are provided. Numerical experiments are reported to demonstrate the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On well-conditioned spectral collocation and spectral methods by the integral reformulation

Well-conditioned spectral collocation and spectral methods have recently been proposed to solve differential equations. In this paper, we revisit the well-conditioned spectral collocation methods proposed in [T.~A. Driscoll, {\it J. Comput. Phys.}, 229 (2010), pp.~5980-5998] and [L.-L. Wang, M.~D. Samson, and X.~Zhao, {\it SIAM J. Sci. Comput.}, 36 (2014), pp.~A907--A929], and the ultraspherical spectral method proposed in [S.~Olver and A.~Townsend, {\it SIAM Rev.}, 55 (2013), pp.~462--489] for an $m$th-order ordinary differential equation from the viewpoint of the integral reformulation. Moreover, we propose a Chebyshev spectral method for the integral reformulation. The well-conditioning of these methods is obvious by noting that the resulting linear operator is a compact perturbation of the identity. The adaptive QR approach for the ultraspherical spectral method still applies to the almost-banded infinite-dimensional system arising in the Chebyshev spectral method for the integral reformulation. Numerical examples are given to confirm the well-conditioning of the Chebyshev spectral method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Preconditioning fractional spectral collocation

Fractional spectral collocation (FSC) method based on fractional Lagrange interpolation has recently been proposed to solve fractional differential equations. Numerical experiments show that the linear systems in FSC become extremely ill-conditioned as the number of collocation points increases. By introducing suitable fractional Birkhoff interpolation problems, we present fractional integration preconditioning matrices for the ill-conditioned linear systems in FSC. The condition numbers of the resulting linear systems are independent of the number of collocation points. Numerical examples are given.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Preconditioning rectangular spectral collocation

Rectangular spectral collocation (RSC) methods have recently been proposed to solve linear and nonlinear differential equations with general boundary conditions and/or other constraints. The involved linear systems in RSC become extremely ill-conditioned as the number of collocation points increases. By introducing suitable Birkhoff-type interpolation problems, we present pseudospectral integration preconditioning matrices for the ill-conditioned linear systems in RSC. The condition numbers of the preconditioned linear systems are independent of the number of collocation points. Numerical examples are given.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two spectral methods for 2D quasi-periodic scattering problems

We consider the 2D quasi-periodic scattering problem in optics, which has been modelled by a boundary value problem governed by Helmholtz equation with transparent boundary conditions. A spectral collocation method and a tensor product spectral method are proposed to numerically solve the problem on rectangles. The discretization parameters can be adaptively chosen so that the numerical solution approximates the exact solution to a high accuracy. Our methods also apply to solve general partial differential equations in two space dimensions, one of which is periodic. Numerical examples are presented to illustrate the accuracy and efficiency of our methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interstitial-Boron Solution Strengthened WB$_{3+x}$

By means of variable-composition evolutionary algorithm coupled with density functional theory and in combination with aberration-corrected high-resolution transmission electron microscopy experiments, we have studied and characterized the composition, structure and hardness properties of WB$_{3+x}$ ($x$ $<$ 0.5). We provide robust evidence for the occurrence of stoichiometric WB$_3$ and non-stoichiometric WB$_{3+x}$ both crystallizing in the metastable $hP$16 ($P6_3/mmc$) structure. No signs for the formation of the highly debated WB$_4$ (both $hP$20 and $hP$10) phases were found. Our results rationalize the seemingly contradictory high-pressure experimental findings and suggest that the interstitial boron atom is located in the tungsten layer and vertically interconnect with four boron atoms, thus forming a typical three-center boron net with the upper and lower boron layers in a three-dimensional covalent network, which thereby strengthen the hardness.

0 citations0 reviewsNo code linkedOpen access