-
An extension of the low-rank Lyapunov ADI to non-zero initial values and its applications
Authors:
Jonas Schulze,
Jens Saak
Abstract:
We derive the Alternating-Direction Implicit (ADI) method based on a commuting operator split and apply the results to the continuous time algebraic Lyapunov equation with low-rank constant term and approximate solution. Previously, it has been mandatory to start the low-rank ADI (LR-ADI) with an all-zero initial value. Our approach retains the known efficient iteration schemes of low-rank increme…
▽ More
We derive the Alternating-Direction Implicit (ADI) method based on a commuting operator split and apply the results to the continuous time algebraic Lyapunov equation with low-rank constant term and approximate solution. Previously, it has been mandatory to start the low-rank ADI (LR-ADI) with an all-zero initial value. Our approach retains the known efficient iteration schemes of low-rank increments and residual to arbitrary low-rank initial values for the LR-ADI method. We further generalize some of the known properties of the LR-ADI for Lyapunov equations to larger classes of algorithms or problems.
We investigate the performance of arbitrary initial values using two outer iterations in which LR-ADI is typically called. First, we solve an algebraic Riccati equation with the Newton method. Second, we solve a differential Riccati equation with a first-order Rosenbrock method. Numerical experiments confirm that the proposed new initial value of the alternating-directions implicit (ADI) can lead to a significant reduction in the total number of ADI steps, while also showing a 17% and 8x speed-up over the zero initial value for the two equation types, respectively.
△ Less
Submitted 29 June, 2024; v1 submitted 19 June, 2024;
originally announced June 2024.
-
Numerical methods for closed-loop systems with non-autonomous data
Authors:
B. Baran,
P. Benner,
J. Saak,
T. Stillfjord
Abstract:
By computing a feedback control via the linear quadratic regulator (LQR) approach and simulating a non-linear non-autonomous closed-loop system using this feedback, we combine two numerically challenging tasks. For the first task, the computation of the feedback control, we use the non-autonomous generalized differential Riccati equation (DRE), whose solution determines the time-varying feedback g…
▽ More
By computing a feedback control via the linear quadratic regulator (LQR) approach and simulating a non-linear non-autonomous closed-loop system using this feedback, we combine two numerically challenging tasks. For the first task, the computation of the feedback control, we use the non-autonomous generalized differential Riccati equation (DRE), whose solution determines the time-varying feedback gain matrix. Regarding the second task, we want to be able to simulate non-linear closed-loop systems for which it is known that the regulator is only valid for sufficiently small perturbations. Thus, one easily runs into numerical issues in the integrators when the closed-loop control varies greatly. For these systems, e.g., the A-stable implicit Euler methods fails.\newline On the one hand, we implement non-autonomous versions of splitting schemes and BDF methods for the solution of our non-autonomous DREs. These are well-established DRE solvers in the autonomous case. On the other hand, to tackle the numerical issues in the simulation of the non-linear closed-loop system, we apply a fractional-step-theta scheme with time-adaptivity tuned specifically to this kind of challenge. That is, we additionally base the time-adaptivity on the activity of the control. We compare this approach to the more classical error-based time-adaptivity.\newline We describe techniques to make these two tasks computable in a reasonable amount of time and are able to simulate closed-loop systems with strongly varying controls, while avoiding numerical issues. Our time-adaptivity approach requires fewer time steps than the error-based alternative and is more reliable.
△ Less
Submitted 21 February, 2024;
originally announced February 2024.
-
Using $LDL^{T}$ factorizations in Newton's method for solving general large-scale algebraic Riccati equations
Authors:
Jens Saak,
Steffen W. R. Werner
Abstract:
Continuous-time algebraic Riccati equations can be found in many disciplines in different forms. In the case of small-scale dense coefficient matrices, stabilizing solutions can be computed to all possible formulations of the Riccati equation. This is not the case when it comes to large-scale sparse coefficient matrices. In this paper, we provide a reformulation of the Newton-Kleinman iteration sc…
▽ More
Continuous-time algebraic Riccati equations can be found in many disciplines in different forms. In the case of small-scale dense coefficient matrices, stabilizing solutions can be computed to all possible formulations of the Riccati equation. This is not the case when it comes to large-scale sparse coefficient matrices. In this paper, we provide a reformulation of the Newton-Kleinman iteration scheme for continuous-time algebraic Riccati equations using indefinite symmetric low-rank factorizations. This allows the application of the method to the case of general large-scale sparse coefficient matrices. We provide convergence results for several prominent realizations of the equation and show in numerical examples the effectiveness of the approach.
△ Less
Submitted 9 February, 2024;
originally announced February 2024.
-
Efficient solution of sequences of parametrized Lyapunov equations with applications
Authors:
Davide Palitta,
Zoran Tomljanović,
Ivica Nakić,
Jens Saak
Abstract:
Sequences of parametrized Lyapunov equations can be encountered in many application settings. Moreover, solutions of such equations are often intermediate steps of an overall procedure whose main goal is the computation of $\text{trace}(EX)$ where $X$ denotes the solution of a Lyapunov equation and $E$ is a given matrix. We are interested in addressing problems where the parameter dependency of th…
▽ More
Sequences of parametrized Lyapunov equations can be encountered in many application settings. Moreover, solutions of such equations are often intermediate steps of an overall procedure whose main goal is the computation of $\text{trace}(EX)$ where $X$ denotes the solution of a Lyapunov equation and $E$ is a given matrix. We are interested in addressing problems where the parameter dependency of the coefficient matrix is encoded as a low-rank modification to a \emph{seed}, fixed matrix. We propose two novel numerical procedures that fully exploit such a common structure. The first one builds upon the Sherman-Morrison-Woodbury (SMW) formula and recycling Krylov techniques, and it is well-suited for small dimensional problems as it makes use of dense numerical linear algebra tools. The second algorithm can instead address large-scale problems by relying on state-of-the-art projection techniques based on the extended Krylov subspace. We test the new algorithms on several problems arising in the study of damped vibrational systems and the analyses of output synchronization problems for multi-agent systems. Our results show that the algorithms we propose are superior to state-of-the-art techniques as they are able to remarkably speed up the computation of accurate solutions.
△ Less
Submitted 28 May, 2024; v1 submitted 13 December, 2023;
originally announced December 2023.
-
Diagonally-Addressed Matrix Nicknack: How to improve SpMV performance
Authors:
Jens Saak,
Jonas Schulze
Abstract:
We suggest a technique to reduce the storage size of sparse matrices at no loss of information. We call this technique Diagonally-Adressed (DA) storage. It exploits the typically low matrix bandwidth of matrices arising in applications. For memory-bound algorithms, this traffic reduction has direct benefits for both uni-precision and multi-precision algorithms. In particular, we demonstrate how to…
▽ More
We suggest a technique to reduce the storage size of sparse matrices at no loss of information. We call this technique Diagonally-Adressed (DA) storage. It exploits the typically low matrix bandwidth of matrices arising in applications. For memory-bound algorithms, this traffic reduction has direct benefits for both uni-precision and multi-precision algorithms. In particular, we demonstrate how to apply DA storage to the Compressed Sparse Rows (CSR) format and compare the performance in computing the Sparse Matrix Vector (SpMV) product, which is a basic building block of many iterative algorithms. We investigate 1367 matrices from the SuiteSparse Matrix Collection fitting into the CSR format using signed 32 bit indices. More than 95% of these matrices fit into the DA-CSR format using 16 bit column indices, potentially after Reverse Cuthill-McKee (RCM) reordering. Using IEEE 754 double precision scalars, we observe a performance uplift of 11% (single-threaded) or 17.5% (multithreaded) on average when the traffic exceeds the size of the last-level CPU cache. The predicted uplift in this scenario is 20%. For traffic within the CPU's combined level 2 and level 3 caches, the multithreaded performance uplift is over 40% for a few test matrices.
△ Less
Submitted 12 July, 2023;
originally announced July 2023.
-
Towards a Benchmark Framework for Model Order Reduction in the Mathematical Research Data Initiative (MaRDI)
Authors:
Peter Benner,
Kathryn Lund,
Jens Saak
Abstract:
The race for the most efficient, accurate, and universal algorithm in scientific computing drives innovation. At the same time, this healthy competition is only beneficial if the research output is actually comparable to prior results. Fairly comparing algorithms can be a complex endeavor, as the implementation, configuration, compute environment, and test problems need to be well-defined. Due to…
▽ More
The race for the most efficient, accurate, and universal algorithm in scientific computing drives innovation. At the same time, this healthy competition is only beneficial if the research output is actually comparable to prior results. Fairly comparing algorithms can be a complex endeavor, as the implementation, configuration, compute environment, and test problems need to be well-defined. Due to the increase in computer-based experiments, new infrastructure for facilitating the exchange and comparison of new algorithms is also needed. To this end, we propose a benchmark framework, as a set of generic specifications for comparing implementations of algorithms using test cases native to a community. Its value lies in its ability to fairly compare and validate existing methods for new applications, as well as compare newly developed methods with existing ones. As a prototype for a more general framework, we have begun building a benchmark tool for the model order reduction (MOR) community. The data basis of the tool is the collection of the Model Order Reduction Wiki (MORWiki). The wiki features three main categories: benchmarks, methods, and software. An editorial board curates submissions and patrols edited entries. Data sets for linear and parametric-linear models are already well represented in the existing collection. Data sets for non-linear or procedural models, for which only evaluation data, or codes / algorithmic descriptions, rather than equations, are available, are being added and extended. Properties and interesting characteristics used for benchmark selection and later assessments are recorded in the model metadata. Our tool, the Model Order Reduction Benchmarker (MORB) is under active development for linear time-invariant systems and solvers.
△ Less
Submitted 30 June, 2023;
originally announced July 2023.
-
Research-Data Management Planning in the German Mathematical Community
Authors:
Tobias Boege,
René Fritze,
Christiane Görgen,
Jeroen Hanselman,
Dorothea Iglezakis,
Lars Kastner,
Thomas Koprucki,
Tabea Krause,
Christoph Lehrenfeld,
Silvia Polla,
Marco Reidelbach,
Christian Riedel,
Jens Saak,
Björn Schembera,
Karsten Tabelow,
Marcus Weber
Abstract:
In this paper we discuss the notion of research data for the field of mathematics and report on the status quo of research-data management and planning. A number of decentralized approaches are presented and compared to needs and challenges faced in three use cases from different mathematical subdisciplines. We highlight the importance of tailoring research-data management plans to mathematicians'…
▽ More
In this paper we discuss the notion of research data for the field of mathematics and report on the status quo of research-data management and planning. A number of decentralized approaches are presented and compared to needs and challenges faced in three use cases from different mathematical subdisciplines. We highlight the importance of tailoring research-data management plans to mathematicians' research processes and discuss their usage all along the data life cycle.
△ Less
Submitted 22 November, 2022;
originally announced November 2022.
-
Riccati-feedback Control of a Two-dimensional Two-phase Stefan Problem
Authors:
Björn Baran,
Peter Benner,
Jens Saak
Abstract:
We discuss the feedback control problem for a two-dimensional two-phase Stefan problem. In our approach, we use a sharp interface representation in combination with mesh-movement to track the interface position. To attain a feedback control, we apply the linear-quadratic regulator approach to a suitable linearization of the problem. We address details regarding the discretization and the interface…
▽ More
We discuss the feedback control problem for a two-dimensional two-phase Stefan problem. In our approach, we use a sharp interface representation in combination with mesh-movement to track the interface position. To attain a feedback control, we apply the linear-quadratic regulator approach to a suitable linearization of the problem. We address details regarding the discretization and the interface representation therein. Further, we document the matrix assembly to generate a non-autonomous generalized differential Riccati equation. To numerically solve the Riccati equation, we use low-rank factored and matrix-valued versions of the non-autonomous backward differentiation formulas, which incorporate implicit index reduction techniques. For the numerical simulation of the feedback controlled Stefan problem, we use a time-adaptive fractional-step-theta scheme.
We provide the implementations for the developed methods and test these in several numerical experiments. With these experiments we show that our feedback control approach is applicable to the Stefan control problem and makes this large-scale problem computable. Also, we discuss the influence of several controller design parameters, such as the choice of inputs and outputs.
△ Less
Submitted 21 December, 2022; v1 submitted 12 September, 2022;
originally announced September 2022.
-
On an integrated Krylov-ADI solver for large-scale Lyapunov equations
Authors:
Peter Benner,
Davide Palitta,
Jens Saak
Abstract:
One of the most computationally expensive steps of the low-rank ADI method for large-scale Lyapunov equations is the solution of a shifted linear system at each iteration. We propose the use of the extended Krylov subspace method for this task. In particular, we illustrate how a single approximation space can be constructed to solve all the shifted linear systems needed to achieve a prescribed acc…
▽ More
One of the most computationally expensive steps of the low-rank ADI method for large-scale Lyapunov equations is the solution of a shifted linear system at each iteration. We propose the use of the extended Krylov subspace method for this task. In particular, we illustrate how a single approximation space can be constructed to solve all the shifted linear systems needed to achieve a prescribed accuracy in terms of Lyapunov residual norm. Moreover, we show how to fully merge the two iterative procedures in order to obtain a novel, efficient implementation of the low-rank ADI method, for an important class of equations. Many state-of-the-art algorithms for the shift computation can be easily incorporated into our new scheme, as well. Several numerical results illustrate the potential of our novel procedure when compared to an implementation of the low-rank ADI method based on sparse direct solvers for the shifted linear systems.
△ Less
Submitted 8 August, 2022; v1 submitted 31 March, 2022;
originally announced March 2022.
-
Optimization-based parametric model order reduction via $\mathcal{H}_2\otimes\mathcal{L}_2$ first-order necessary conditions
Authors:
Manuela Hund,
Tim Mitchell,
Petar Mlinarić,
Jens Saak
Abstract:
In this paper, we generalize existing frameworks for $\mathcal{H}_2\otimes\mathcal{L}_2$-optimal model order reduction to a broad class of parametric linear time-invariant systems. To this end, we derive first-order necessary ptimality conditions for a class of structured reduced-order models, and then building on those, propose a stability-preserving optimization-based method for computing locall…
▽ More
In this paper, we generalize existing frameworks for $\mathcal{H}_2\otimes\mathcal{L}_2$-optimal model order reduction to a broad class of parametric linear time-invariant systems. To this end, we derive first-order necessary ptimality conditions for a class of structured reduced-order models, and then building on those, propose a stability-preserving optimization-based method for computing locally $\mathcal{H}_2\otimes\mathcal{L}_2$-optimal reduced-order models. We also make a theoretical comparison to existing approaches in the literature, and in numerical experiments, show how our new method, with reasonable computational effort, produces stable optimized reduced-order models with significantly lower approximation errors.
△ Less
Submitted 31 March, 2022; v1 submitted 4 March, 2021;
originally announced March 2021.
-
Parametric model order reduction using pyMOR
Authors:
Petar Mlinarić,
Stephan Rave,
Jens Saak
Abstract:
pyMOR is a free software library for model order reduction that includes both reduced basis and system-theoretic methods. All methods are implemented in terms of abstract vector and operator interfaces, which allows direct integration of pyMOR's algorithms with a wide array of external PDE solvers. In this contribution, we give a brief overview of the available methods and experimentally compare t…
▽ More
pyMOR is a free software library for model order reduction that includes both reduced basis and system-theoretic methods. All methods are implemented in terms of abstract vector and operator interfaces, which allows direct integration of pyMOR's algorithms with a wide array of external PDE solvers. In this contribution, we give a brief overview of the available methods and experimentally compare them for the parametric instationary thermal-block benchmark defined in arXiv:2003.00846.
△ Less
Submitted 2 July, 2020; v1 submitted 12 March, 2020;
originally announced March 2020.
-
Matrix Equations, Sparse Solvers: M-M.E.S.S.-2.0.1 -- Philosophy, Features and Application for (Parametric) Model
Authors:
Peter Benner,
Martin Köhler,
Jens Saak
Abstract:
Matrix equations are omnipresent in (numerical) linear algebra and systems theory. Especially in model order reduction (MOR) they play a key role in many balancing based reduction methods for linear dynamical systems. When these systems arise from spatial discretizations of evolutionary partial differential equations, their coefficient matrices are typically large and sparse. Moreover, the numbers…
▽ More
Matrix equations are omnipresent in (numerical) linear algebra and systems theory. Especially in model order reduction (MOR) they play a key role in many balancing based reduction methods for linear dynamical systems. When these systems arise from spatial discretizations of evolutionary partial differential equations, their coefficient matrices are typically large and sparse. Moreover, the numbers of inputs and outputs of these systems are typically far smaller than the number of spatial degrees of freedom. Then, in many situations the solutions of the corresponding large-scale matrix equations are observed to have low (numerical) rank. This feature is exploited by M-M.E.S.S. to find successively larger low-rank factorizations approximating the solutions. This contribution describes the basic philosophy behind the implementation and the features of the package, as well as its application in the model order reduction of large-scale linear time-invariant (LTI) systems and parametric LTI systems.
△ Less
Submitted 9 May, 2020; v1 submitted 4 March, 2020;
originally announced March 2020.
-
A Non-stationary Thermal-Block Benchmark Model for Parametric Model Order Reduction
Authors:
Stephan Rave,
Jens Saak
Abstract:
In this contribution we aim to satisfy the demand for a publicly available benchmark for parametric model order reduction that is scalable both in degrees of freedom as well as parameter dimension.
In this contribution we aim to satisfy the demand for a publicly available benchmark for parametric model order reduction that is scalable both in degrees of freedom as well as parameter dimension.
△ Less
Submitted 16 April, 2020; v1 submitted 2 March, 2020;
originally announced March 2020.
-
Sustainable Research Software Hand-Over
Authors:
Jörg Fehr,
Christian Himpe,
Stephan Rave,
Jens Saak
Abstract:
Scientific software projects evolve rapidly in their initial development phase, yet at the end of a funding period, the completion of a research project, thesis, or publication, further engagement in the project may slow down or cease completely. To retain the invested effort for the sciences, this software needs to be preserved or handed over to a succeeding developer or team, such as the next ge…
▽ More
Scientific software projects evolve rapidly in their initial development phase, yet at the end of a funding period, the completion of a research project, thesis, or publication, further engagement in the project may slow down or cease completely. To retain the invested effort for the sciences, this software needs to be preserved or handed over to a succeeding developer or team, such as the next generation of (PhD) students.
Comparable guides provide top-down recommendations for project leads. This paper intends to be a bottom-up approach for sustainable hand-over processes from a developer's perspective. An important characteristic in this regard is the project's size, by which this guideline is structured. Furthermore, checklists are provided, which can serve as a practical guide for implementing the proposed measures.
△ Less
Submitted 30 September, 2020; v1 submitted 19 September, 2019;
originally announced September 2019.
-
A numerical comparison of solvers for large-scale, continuous-time algebraic Riccati equations and LQR problems
Authors:
Peter Benner,
Zvonimir Bujanović,
Patrick Kürschner,
Jens Saak
Abstract:
In this paper, we discuss numerical methods for solving large-scale continuous-time algebraic Riccati equations. These methods have been the focus of intensive research in recent years, and significant progress has been made in both the theoretical understanding and efficient implementation of various competing algorithms. There are several goals of this manuscript: first, to gather in one place a…
▽ More
In this paper, we discuss numerical methods for solving large-scale continuous-time algebraic Riccati equations. These methods have been the focus of intensive research in recent years, and significant progress has been made in both the theoretical understanding and efficient implementation of various competing algorithms. There are several goals of this manuscript: first, to gather in one place an overview of different approaches for solving large-scale Riccati equations, and to point to the recent advances in each of them. Second, to analyze and compare the main computational ingredients of these algorithms, to detect their strong points and their potential bottlenecks. And finally, to compare the effective implementations of all methods on a set of relevant benchmark examples, giving an indication of their relative performance.
△ Less
Submitted 10 April, 2020; v1 submitted 2 November, 2018;
originally announced November 2018.
-
Efficient Solution of Large-Scale Algebraic Riccati Equations Associated with Index-2 DAEs via the Inexact Low-Rank Newton-ADI Method
Authors:
Peter Benner,
Matthias Heinkenschloss,
Jens Saak,
Heiko K. Weichelt
Abstract:
This paper extends the algorithm of Benner, Heinkenschloss, Saak, and Weichelt: An inexact low-rank Newton-ADI method for large-scale algebraic Riccati equations, Applied Numerical Mathematics Vol.~108 (2016), pp.~125--142, doi:10.1016/j.apnum.2016.05.006 to Riccati equations associated with Hessenberg index-2 Differential Algebratic Equation (DAE) systems. Such DAE systems arise, e.g., from semi-…
▽ More
This paper extends the algorithm of Benner, Heinkenschloss, Saak, and Weichelt: An inexact low-rank Newton-ADI method for large-scale algebraic Riccati equations, Applied Numerical Mathematics Vol.~108 (2016), pp.~125--142, doi:10.1016/j.apnum.2016.05.006 to Riccati equations associated with Hessenberg index-2 Differential Algebratic Equation (DAE) systems. Such DAE systems arise, e.g., from semi-discretized, linearized (around steady state) Navier-Stokes equations. The solution of the associated Riccati equation is important, e.g., to compute feedback laws that stabilize the Navier-Stokes equations. Challenges in the numerical solution of the Riccati equation arise from the large-scale of the underlying systems and the algebraic constraint in the DAE system. These challenges are met by a careful extension of the inexact low-rank Newton-ADI method to the case of DAE systems. A main ingredient in the extension to the DAE case is the projection onto the manifold described by the algebraic constraints. In the algorithm, the equations are never explicitly projected, but the projection is only applied as needed. Numerical experience indicates that the algorithmic choices for the control of inexactness and line-search can help avoid subproblems with matrices that are only marginally stable. The performance of the algorithm is illustrated on a large-scale Riccati equation associated with the stabilization of Navier-Stokes flow around a cylinder.
△ Less
Submitted 12 September, 2019; v1 submitted 4 April, 2018;
originally announced April 2018.
-
A Connection Between Time Domain Model Order Reduction and Moment Matching for LTI Systems
Authors:
Manuela Hund,
Jens Saak
Abstract:
We investigate the time domain model order reduction (MOR) framework using general orthogonal polynomials by Jiang and Chen 2012 and extend their idea by exploiting the structure of the corresponding linear system of equations. Identifying an equivalent Sylvester equation, we show a connection to a rational Krylov subspace, and thus to moment matching. This theoretical link between the MOR techniq…
▽ More
We investigate the time domain model order reduction (MOR) framework using general orthogonal polynomials by Jiang and Chen 2012 and extend their idea by exploiting the structure of the corresponding linear system of equations. Identifying an equivalent Sylvester equation, we show a connection to a rational Krylov subspace, and thus to moment matching. This theoretical link between the MOR techniques is illustrated by three numerical examples. For linear time-invariant systems, the link also motivates that the time domain approach can be at best as accurate as moment matching, since the expansion points are fixed by the choice of the polynomial basis, while in moment matching they can be adapted to the system.
△ Less
Submitted 6 July, 2018; v1 submitted 22 January, 2018;
originally announced January 2018.
-
Predictive Modeling: An Optimized and Dynamic Solution Framework for Systematic Value Investing
Authors:
R. J. Sak
Abstract:
This paper defines systematic value investing as an empirical optimization problem. Predictive modeling is introduced as a systematic value investing methodology with dynamic and optimization features. A predictive modeling process is demonstrated using financial metrics from Gray & Carlisle and Buffett & Clark. A 31-year portfolio backtest (1985 - 2016) compares performance between predictive mod…
▽ More
This paper defines systematic value investing as an empirical optimization problem. Predictive modeling is introduced as a systematic value investing methodology with dynamic and optimization features. A predictive modeling process is demonstrated using financial metrics from Gray & Carlisle and Buffett & Clark. A 31-year portfolio backtest (1985 - 2016) compares performance between predictive models and Gray & Carlisle's Quantitative Value strategy. A 26-year portfolio backtest (1990 - 2016) uses an expanded set of predictor variables to show financial performance improvements. This paper includes secondary novel contributions. Quantitative definitions are provided for Buffett & Clark's value investing metrics. The "Sak ratio" is proposed as an extension to the Benjamini-Hochberg procedure for the inferential identification of false positive observations.
△ Less
Submitted 10 September, 2017;
originally announced September 2017.
-
Best Practices for Replicability, Reproducibility and Reusability of Computer-Based Experiments Exemplified by Model Reduction Software
Authors:
Jörg Fehr,
Jan Heiland,
Christian Himpe,
Jens Saak
Abstract:
Over the recent years the importance of numerical experiments has gradually been more recognized. Nonetheless, sufficient documentation of how computational results have been obtained is often not available. Especially in the scientific computing and applied mathematics domain this is crucial, since numerical experiments are usually employed to verify the proposed hypothesis in a publication. This…
▽ More
Over the recent years the importance of numerical experiments has gradually been more recognized. Nonetheless, sufficient documentation of how computational results have been obtained is often not available. Especially in the scientific computing and applied mathematics domain this is crucial, since numerical experiments are usually employed to verify the proposed hypothesis in a publication. This work aims to propose standards and best practices for the setup and publication of numerical experiments. Naturally, this amounts to a guideline for development, maintenance, and publication of numerical research software. Such a primer will enable the replicability and reproducibility of computer-based experiments and published results and also promote the reusability of the associated software.
△ Less
Submitted 5 July, 2016;
originally announced July 2016.
-
RADI: A low-rank ADI-type algorithm for large scale algebraic Riccati equations
Authors:
Peter Benner,
Zvonimir Bujanović,
Patrick Kürschner,
Jens Saak
Abstract:
This paper introduces a new algorithm for solving large-scale continuous-time algebraic Riccati equations (CARE). The advantage of the new algorithm is in its immediate and efficient low-rank formulation, which is a generalization of the Cholesky-factored variant of the Lyapunov ADI method. We discuss important implementation aspects of the algorithm, such as reducing the use of complex arithmetic…
▽ More
This paper introduces a new algorithm for solving large-scale continuous-time algebraic Riccati equations (CARE). The advantage of the new algorithm is in its immediate and efficient low-rank formulation, which is a generalization of the Cholesky-factored variant of the Lyapunov ADI method. We discuss important implementation aspects of the algorithm, such as reducing the use of complex arithmetic and shift selection strategies. We show that there is a very tight relation between the new algorithm and three other algorithms for CARE previously known in the literature -- all of these seemingly different methods in fact produce exactly the same iterates when used with the same parameters: they are algorithmically different descriptions of the same approximation sequence to the Riccati solution.
△ Less
Submitted 9 October, 2016; v1 submitted 30 September, 2015;
originally announced October 2015.