BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Matthew M. Lin

Matthew M. Lin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NA

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2016arXiv

$AB$-algorithm and its application for solving matrix square roots

This work is to propose an iterative method of choice to compute a stable subspace of a regular matrix pencil. This approach is to define a sequence of matrix pencils via particular left null spaces. We show that this iteration preserves a discrete-type flow depending only on the initial matrix pencil. Via this recursion relationship, we propose an accelerated iterative method to compute the stable subspace and use it to provide a theoretical result to solve the principal square root of a given matrix, both nonsingular and singular. We show that this method can not only find out the matrix square root, but also construct an iterative approach which converges to the square root with any desired order.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Integer Low Rank Approximation of Integer matrices

Integer data sets frequently appear in many applications in sciences and technology. To analyze these, integer low rank approximation has received much attention due to its capacity of representing the results in integers preserving the meaning of the original data sets. To our knowledge, none of previously proposed techniques developed for real numbers can be successfully applied, since integers are discrete in nature. In this work, we start with a thorough review of algorithms for solving integer least squares problems, {and} then develop a block coordinate descent method based on the integer least squares estimation to obtain the integer low rank approximation of integer matrices. The numerical application on association analysis and numerical experiments on random integer matrices are presented. Our computed results seem to suggest that our method can find a more accurate solution than other existing methods for continuous data sets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Preconditioned iterative methods for space-time fractional advection-diffusion equations

In this paper we want to propose practical numerical methods to solve a class of initial-boundary problem of space-time fractional advection-diffusion equations. To start with, an implicit method based on two-sided Grünwald formulae is proposed with a discussion of the stability and consistency. Then, the preconditioned generalized minimal residual (preconditioned GMRES) method and the preconditioned conjugate gradient normal residual ({preconditioned} CGNR) method, with an easily constructed preconditioner, are developed. Importantly, because the resulting systems are Topelitz-like, the fast Fourier transform can be applied to significantly reduce the computational cost. Numerical experiments are implemented to show the efficiency of our preconditioner, even with cases of variable coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A fast numerical algorithm for constructing nonnegative matrices with prescribed real eigenvalues

The study of solving the inverse eigenvalue problem for nonnegative matrices has been around for decades. It is clear that an inverse eigenvalue problem is trivial if the desirable matrix is not restricted to a certain structure. Provided with the real spectrum, this paper presents a numerical procedure, based on the induction principle, to solve two kinds of inverse eigenvalue problems, one for nonnegative matrices and another for symmetric nonnegative matrices. As an immediate application, our approach can offer not only the sufficient condition for solving inverse eigenvalue problems for nonnegative or symmetric nonnegative matrices, but also a quick numerical way to solve inverse eigenvalue problem for stochastic matrices. Numerical examples are presented for problems of relatively larger size.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A note on Sylvester-type equations

This work is to provide a comprehensive treatment of the relationship between the theory of the generalized (palindromic) eigenvalue problem and the theory of the Sylvester-type equations. Under a regularity assumption for a specific matrix pencil, we show that the solution of the $\star$-Sylvester matrix equation is uniquely determined and can be obtained by considering its corresponding deflating subspace. We also propose an iterative method with quadratic convergence to compute the stabilizing solution of the $\star$-Sylvester matrix equation via the well-developed palindromic doubling algorithm. We believe that our discussion is the first which implements the tactic of the deflating subspace for solving Sylvester equations and could give rise to the possibility of developing an advanced and effective solver for different types of matrix equations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Eigenvalue Shift Technique and Its Eigenstructure Analysis of a Matrix

The eigenvalue shift technique is the most well-known and fundamental tool for matrix computations. Applications include the search of eigeninformation, the acceleration of numerical algorithms, the study of Google's PageRank. The shift strategy arises from the concept investigated by Brauer [1] for changing the value of an eigenvalue of a matrix to the desired one, while keeping the remaining eigenvalues and the original eigenvectors unchanged. The idea of shifting distinct eigenvalues can easily be generalized by Brauer's idea. However, shifting an eigenvalue with multiple multiplicities is a challenge issue and worthy of our investigation. In this work, we propose a new way for updating an eigenvalue with multiple multiplicities and thoroughly analyze its corresponding Jordan canonical form after the update procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Shifting Technique for Computing the Extreme Solutions of $X + A^\top X^{-1} A = Q$

We propose a new way for speeding up the search of the maximal solution $X_+$ of $X + A^\top X^{-1} A = Q$. It is known that the speed of convergence of traditional approaches for solving this problem depends highly on the spectral radius $ρ(X_+^{-1}A)$. If $ρ(X_+^{-1}A)$ is close to one or equal to one, the iterations of traditional approaches converges very slowly or does not converge. Our goal is to come up with a shifting tactic to remove the singularities embedded in $ρ(X_+^{-1}A)$. Finally, an example is used to demonstrate the capacity of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Shifting Technique for Solving a Nonsymmetric Algebraic Riccati Equation

This paper analyzes a special instance of nonsymmetric algebraic matrix Riccati equations arising from transport theory. Traditional approaches for finding the minimal nonnegative solution of the matrix Riccati equations are based on the fixed point iteration and the speed of the convergence is linear. Relying on simultaneously matrix computation, a structure-preserving doubling algorithm (SDA) with quadratic convergence is designed for improving the speed of convergence. The difficulty is that the double algorithm with quadratic convergence cannot guarantee to work all the time. Our main trust in this work is to show that applied with a suitable shifted technique, the SDA is guaranteed to converge quadratically with no breakdown. Also, we modify the conventional simple iteration algorithm in the critical case to dramatically improve the speed of convergence. Numerical experiments strongly suggest that the total number of computational steps can be significantly reduced via the shifting procedure.

0 citations0 reviewsNo code linkedOpen access