Source author record

Alicja Smoktunowicz

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.NA math.RA Numerical Analysis

Catalog footprint

Research overview

6works

3topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

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.

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.

Alicja Smoktunowicz

Research overview

Is this your work?

6 publications

On computing the symplectic $LL^T$ factorization

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

Numerical stability of iterative refinement with a relaxation for linear systems

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

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

Reorthogonalized Block Classical Gram--Schmidt

See the researcher in context

Building this map preview

Alicja Smoktunowicz

Research overview

Is this your work?

6 publications

On computing the symplectic LLTLL^TLLT factorization

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

Numerical stability of iterative refinement with a relaxation for linear systems

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

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

Reorthogonalized Block Classical Gram--Schmidt

See the researcher in context

Building this map preview

On computing the symplectic $LL^T$ factorization