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

Zhonggang Zeng

Zhonggang Zeng contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysis

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
2topics
2close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

8 published item(s)

preprint2023arXiv

Computing multiple roots of inexact polynomials

We present a combination of two algorithms that accurately calculate multiple roots of general polynomials. Algorithm I transforms the singular root-finding into a regular nonlinear least squares problem on a pejorative manifold, and calculates multiple roots simultaneously from a given multiplicity structure and initial root approximations. To fulfill the input requirement of Algorithm I, we develop a numerical GCD-finder containing a successive singular value updating and an iterative GCD refinement as the main engine of Algorithm II that calculates the multiplicity structure and the initial root approximation. While limitations of our algorithm exist in identifying the multiplicity structure in certain situations, the combined method calculates multiple roots with high accuracy and consistency in practice without using multiprecision arithmetic even if the coefficients are inexact. This is perhaps the first blackbox-type root-finder with such capabilities. To measure the sensitivity of the multiple roots, a structure-preserving condition number is proposed and error bounds are established. According to our computational experiments and error analysis, a polynomial being ill-conditioned in the conventional sense can be well conditioned with the multiplicity structure being preserved, and its multiple roots can be computed with high accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A numerical method for computing the Jordan Canonical Form

The Jordan Canonical Form of a matrix is highly sensitive to perturbations, and its numerical computation remains a formidable challenge. This paper presents a regularization theory that establishes a well-posed least squares problem of finding the nearest staircase decomposition in the matrix bundle of the highest codimension. A two-staged algorithm is developed for computing the numerical Jordan Canonical Form. At the first stage, the method calculates the Jordan structure of the matrix and an initial approximation to the multiple eigenvalues. The staircase decomposition is then constructed by an iterative algorithm at the second stage. As a result, the numerical Jordan Canonical decomposition along with multiple eigenvalues can be computed with high accuracy even if the underlying matrix is perturbed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Geometric modeling and regularization of algebraic problems

Discontinuity with respect to data perturbations is common in algebraic computation where solutions are often highly sensitive. Such problems can be modeled as solving systems of equations at given data parameters. By appending auxiliary equations, the models can be formulated to satisfy four easily verifiable conditions so that the data form complex analytic manifolds on which the solutions maintain their structures and the Lipschitz continuity. When such a problem is given with empirical data, solving the system becomes a least squares problem whose solution uniquely exists and enjoys Lipschitz continuity as long as the data point is in a tubular neighborhood of the manifold. As a result, the singular problem is regularized as a well-posed computational problem.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the sensitivity of singular and ill-Conditioned linear systems

Solving a singular linear system for an individual vector solution is an ill-posed problem with a condition number infinity. From an alternative perspective, however, the general solution of a singular system is of a bounded sensitivity as a unique element in an affine Grassmannian. If a singular linear system is given through empirical data that are sufficiently accurate with a tight error bound, a properly formulated general numerical solution uniquely exists in the same affine Grassmannian, enjoys Lipschitz continuity and approximates the underlying exact solution with an accuracy in the same order as the data. Furthermore, any backward accurate numerical solution vector is an accurate approximation to one of the solutions of the underlying singular system.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Sensitivity and computation of a defective eigenvalue

A defective eigenvalue is well documented to be hypersensitive to data perturbations and round-off? errors, making it a formidable challenge in numerical computation particularly when the matrix is known through approximate data. This paper establishes a finitely bounded sensitivity of a defective eigenvalue with respect to perturbations that preserve the geometric multiplicity and the smallest Jordan block size. Based on this perturbation theory, numerical computation of a defective eigenvalue is regularized as a well-posed least squares problem so that it can be accurately carried out using floating point arithmetic even if the matrix is perturbed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Singular algebraic equations with empirical data

Singular equations with rank-deficient Jacobians arise frequently in algebraic computing applications. As shown in case studies in this paper, direct and intuitive modeling of algebraic problems often results in nonisolated singular solutions. The challenges become formidable when the problems need to be solved from empirical data of limited accuracy. A newly discovered low-rank Newton's iteration emerges as an effective regularization mechanism that enables solving singular equations accurately with an error bound in the same order as the data error. This paper elaborates applications of new methods on solving singular algebraic equations such as singular linear systems, polynomial GCD and factorizations as well as matrix defective eigenvalue problems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The numerical factorization of polynomials

Polynomial factorization in conventional sense is an ill-posed problem due to its discontinuity with respect to coefficient perturbations, making it a challenge for numerical computation using empirical data. As a regularization, this paper formulates the notion of numerical factorization based on the geometry of polynomial spaces and the stratification of factorization manifolds. Furthermore, this paper establishes the existence, uniqueness, Lipschitz continuity, condition number, and convergence of the numerical factorization to the underlying exact factorization, leading to a robust and efficient algorithm with a MATLAB implementation capable of accurate polynomial factorizations using floating point arithmetic even if the coefficients are perturbed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The numerical greatest common divisor of univariate polynomials

This paper presents a regularization theory for numerical computation of polynomial greatest common divisors and a convergence analysis, along with a detailed description of a blackbox-type algorithm. The root of the ill-posedness in conventional GCD computation is identified by its geometry where polynomials form differentiable manifolds entangled in a stratification structure. With a proper regularization, the numerical GCD is proved to be strongly well-posed. Most importantly, the numerical GCD solves the problem of finding the GCD accurately using floating point arithmetic even if the data are perturbed. A sensitivity measurement, error bounds at each computing stage, and the overall convergence are established rigorously. The computing results of selected test examples show that the algorithm and software appear to be robust and accurate.

0 citations0 reviewsNo code linkedOpen access