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

Jörg Liesen

Jörg Liesen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.CVastro-ph.GANumerical Analysisastro-ph.IM

Catalog footprint

What is connected

11works
5topics
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

11 published item(s)

preprint2022arXiv

On non-Hermitian positive (semi)definite linear algebraic systems arising from dissipative Hamiltonian DAEs

We discuss different cases of dissipative Hamiltonian differential-algebraic equations and the linear algebraic systems that arise in their linearization or discretization. For each case we give examples from practical applications. An important feature of the linear algebraic systems is that the (non-Hermitian) system matrix has a positive definite or semidefinite Hermitian part. In the positive definite case we can solve the linear algebraic systems iteratively by Krylov subspace methods based on efficient three-term recurrences. We illustrate the performance of these iterative methods on several examples. The semidefinite case can be challenging and requires additional techniques to deal with "singular part", while the "positive definite part" can still be treated with the three-term recurrence methods.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On conformal maps from multiply connected domains onto lemniscatic domains

We study conformal maps from multiply connected domains in the extended complex plane onto lemniscatic domains. Walsh proved the existence of such maps in 1956 and thus obtained a direct generalization of the Riemann mapping theorem to multiply connected domains. For polynomial pre-images of simply connected sets we derive a construction principle for Walsh's conformal map in terms of the Riemann map for the simply connected set. Moreover, we explicitly construct examples of Walsh's conformal map for certain radial slit domains and circular domains.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast Recovery and Approximation of Hidden Cauchy Structure

We derive an algorithm of optimal complexity which determines whether a given matrix is a Cauchy matrix, and which exactly recovers the Cauchy points defining a Cauchy matrix from the matrix entries. Moreover, we study how to approximate a given matrix by a Cauchy matrix with a particular focus on the recovery of Cauchy points from noisy data. We derive an approximation algorithm of optimal complexity for this task, and prove approximation bounds. Numerical examples illustrate our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Note on the Maximum Number of Zeros of $r(z) - \bar{z}$

An important theorem of Khavinson & Neumann (Proc. Amer. Math. Soc. 134(4), 2006) states that the complex harmonic function $r(z) - \bar{z}$, where $r$ is a rational function of degree $n \geq 2$, has at most $5 (n - 1)$ zeros. In this note we resolve a slight inaccuracy in their proof and in addition we show that for certain functions of the form $r(z) - \bar{z}$ no more than $5 (n - 1) - 1$ zeros can occur. Moreover, we show that $r(z) - \bar{z}$ is regular, if it has the maximal number of zeros.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Perturbing rational harmonic functions by poles

We study how adding certain poles to rational harmonic functions of the form $R(z)-\bar{z}$, with $R(z)$ rational and of degree $d\geq 2$, affects the number of zeros of the resulting functions. Our results are motivated by and generalize a construction of Rhie derived in the context of gravitational microlensing (ArXiv e-print 2003). Of particular interest is the construction and the behavior of rational functions $R(z)$ that are {\em extremal} in the sense that $R(z)-\bar{z}$ has the maximal possible number of $5(d-1)$ zeros.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pták's nondiscrete induction and its application to matrix iterations

Vlastimil Pták's method of nondiscrete induction is based on the idea that in the analysis of iterative processes one should aim at rates of convergence as functions rather than just numbers, because functions may give convergence estimates that are tight throughout the iteration rather than just asymptotically. In this paper we motivate and prove a theorem on nondiscrete induction originally due to Potra and Pták, and we apply it to the Newton iterations for computing the matrix polar decomposition and the matrix square root. Our goal is to illustrate the application of the method of nondiscrete induction in the finite dimensional numerical linear algebra context. We show the sharpness of the resulting convergence estimate analytically for the polar decomposition iteration and for special cases of the square root iteration, as well as on some numerical examples for the square root iteration. We also discuss some of the method's limitations and possible extensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sharp parameter bounds for certain maximal point lenses

Starting from an $n$-point circular gravitational lens having $3n+1$ images, Rhie (2003) used a perturbation argument to construct an $(n+1)$-point lens producing $5n$ images. In this work we give a concise proof of Rhie's result, and we extend the range of parameters in Rhie's model for which maximal lensing occurs. We also study a slightly different construction given by Bayer and Dyer (2007) arising from the $(3n+1)$-point lens. In particular, we extend their results and give sharp parameter bounds for their lens model. By a substitution of variables and parameters we show that both models are equivalent in a certain sense.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A framework for deflated and augmented Krylov subspace methods

We consider deflation and augmentation techniques for accelerating the convergence of Krylov subspace methods for the solution of nonsingular linear algebraic systems. Despite some formal similarity, the two techniques are conceptually different from preconditioning. Deflation (in the sense the term is used here) "removes" certain parts from the operator making it singular, while augmentation adds a subspace to the Krylov subspace (often the one that is generated by the singular operator); in contrast, preconditioning changes the spectrum of the operator without making it singular. Deflation and augmentation have been used in a variety of methods and settings. Typically, deflation is combined with augmentation to compensate for the singularity of the operator, but both techniques can be applied separately. We introduce a framework of Krylov subspace methods that satisfy a Galerkin condition. It includes the families of orthogonal residual (OR) and minimal residual (MR) methods. We show that in this framework augmentation can be achieved either explicitly or, equivalently, implicitly by projecting the residuals appropriately and correcting the approximate solutions in a final step. We study conditions for a breakdown of the deflated methods, and we show several possibilities to avoid such breakdowns for the deflated MINRES method. Numerical experiments illustrate properties of different variants of deflated MINRES analyzed in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Characterization of worst-case GMRES

Given a matrix $A$ and iteration step $k$, we study a best possible attainable upper bound on the GMRES residual norm that does not depend on the initial vector $b$. This quantity is called the worst-case GMRES approximation. We show that the worst case behavior of GMRES for the matrices $A$ and $A^T$ is the same, and we analyze properties of initial vectors for which the worst-case residual norm is attained. In particular, we show that such vectors satisfy a certain "cross equality", and we characterize them as right singular vectors of the corresponding GMRES residual matrix. We show that the worst-case GMRES polynomial may not be uniquely determined, and we consider the relation between the worst-case and the ideal GMRES approximations, giving new examples in which the inequality between the two quantities is sharp at all iteration steps $k\geq 3$. Finally, we give a complete characterization of how the values of the approximation problems in the context of worst-case and ideal GMRES for a real matrix change, when one considers complex (rather than real) polynomials and initial vectors in these problems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Max-min and min-max approximation problems for normal matrices revisited

We give a new proof for an equality of certain max-min and min-max approximation problems involving normal matrices. The previously published proofs of this equality apply tools from matrix theory, (analytic) optimization theory and constrained convex optimization. Our proof uses a classical characterization theorem from approximation theory and thus exploits the link between the two approximation problems with normal matrices on the one hand and approximation problems on compact sets in the complex plane on the other.

0 citations0 reviewsNo code linkedOpen access