Source author record

Alicja Smoktunowicz

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

ResearcherUnclaimed source record

math.NA math.RA Numerical Analysis

Catalog footprint

What is connected

6works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

On computing the symplectic $LL^T$ factorization

We analyze two algorithms for computing the symplectic $LL^T$ factorization $A=LL^T$ of a given symmetric positive definite symplectic matrix $A$. The first algorithm $W_1$ is an implementation of the $HH^T$ factorization from [Dopico et al., 2009], see Theorem 5.2. The second one, algorithm $W_2$ uses both Cholesky and Reverse Cholesky decompositions of symmetric positive definite matrices. We presents a comparison of these algorithms and illustrate their properties by numerical experiments in MATLAB. A particular emphasis is given on simplecticity properties of the computed matrices in floating-point arithmetic.

preprint2020arXiv

Involutive solutions of the Yang-Baxter-like matrix equation -- theory and algorithms

We find all explicit involutive solutions $X \in \mathbb C^{n \times n}$ of the Yang-Baxter-like matrix equation $AXA=XAX$, where $A \in \mathbb C^{n \times n}$ is a given involutory matrix. The construction is algorithmic.

preprint2015arXiv

Numerical stability of iterative refinement with a relaxation for linear systems

Stability analysis of Wilkinson's iterative refinement with a relaxation IR(omega) for solving linear systems is given. It extends existing results for omega=1, i.e., for Wilkinson's iterative refinement. We assume that all computations are performed in fixed (working) precision arithmetic. Numerical tests were done in MATLAB to illustrate our theoretical results. A particular emphasis is given on convergence of iterative refinement with a relaxation. Our tests confirm that the choice omega=1 is the best choice from the point of numerical stability.

preprint2013arXiv

Numerical solution of saddle point problems by block {Gram--Schmidt} orthogonalization

Saddle point problems arise in many important practical applications. In this paper we propose and analyze some algorithms for solving symmetric saddle point problems which are based upon the block Gram-Schmidt method. In particular, we prove that the algorithm BCGS2 (Reorthogonalized Block Classical Gram-Schmidt) using Householder Q-R decomposition implemented in floating point arithmetic is backward stable, under a mild assumption on the matrix $M$. This means that the computed vector $\tilde z$ is the exact solution to a slightly perturbed linear system of equations $Mz = f$.

preprint2013arXiv

On the accuracy and stability of algorithms most commonly used in the evaluation of Chebyshev polynomials of the first kind

This paper provides error analyses of the algorithms most commonly used for the evaluation of the Chebyshev polynomial of the first kind $T_N(x)$. Some of these algorithms are shown to be backward stable. This means that the computed value of $T_N(x)$ in floating point arithmetic by these algorithms can be interpreted as a slightly perturbed value of polynomial $T_N$, for slightly perturbed value of $x$.

preprint2011arXiv

Reorthogonalized Block Classical Gram--Schmidt

A new reorthogonalized block classical Gram--Schmidt algorithm is proposed that factorizes a full column rank matrix $A$ into $A=QR$ where $Q$ is left orthogonal (has orthonormal columns) and $R$ is upper triangular and nonsingular. With appropriate assumptions on the diagonal blocks of $R$, the algorithm, when implemented in floating point arithmetic with machine unit $\macheps$, produces $Q$ and $R$ such that $\| I- Q^{T} Q \|_2 =O(\macheps)$ and $\| A-QR \|_2 =O(\macheps \| A \|_2)$. The resulting bounds also improve a previous bound by Giraud et al. [Num. Math., 101(1):87-100,\ 2005] on the CGS2 algorithm originally developed by Abdelmalek [BIT, 11(4):354--367,\ 1971]. \medskip Keywords: Block matrices, Q--R factorization, Gram-Schmidt process, Condition numbers, Rounding error analysis.