Please wait ... |

Refine search result

CiteExportLink to result list
http://kth.diva-portal.org/smash/resultList.jsf?query=&language=en&searchType=SIMPLE&noOfRows=50&sortOrder=dateIssued_sort_desc&sortOrder2=title_sort_asc&onlyFullText=false&sf=all&aq=%5B%5B%7B%22author%22%3A%5B%22Jarlebring%22%5D%7D%5D%5D&aqe=%5B%5D&aq2=%5B%5B%5D%5D&af=%5B%5D $(function(){PrimeFaces.cw("InputTextarea","widget_formSmash_upper_j_idt482_recordPermLink",{id:"formSmash:upper:j_idt482:recordPermLink",widgetVar:"widget_formSmash_upper_j_idt482_recordPermLink",autoResize:true});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_upper_j_idt482_j_idt484",{id:"formSmash:upper:j_idt482:j_idt484",widgetVar:"widget_formSmash_upper_j_idt482_j_idt484",target:"formSmash:upper:j_idt482:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});

Permanent link

Cite

Citation styleapa harvard1 ieee modern-language-association-8th-edition vancouver Other style $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_upper_j_idt500",{id:"formSmash:upper:j_idt500",widgetVar:"widget_formSmash_upper_j_idt500",behaviors:{change:function(ext) {PrimeFaces.ab({s:"formSmash:upper:j_idt500",e:"change",f:"formSmash",p:"formSmash:upper:j_idt500",u:"formSmash:upper:otherStyle"},ext);}}});});

- apa
- harvard1
- ieee
- modern-language-association-8th-edition
- vancouver
- Other style

Languagede-DE en-GB en-US fi-FI nn-NO nn-NB sv-SE Other locale $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_upper_j_idt511",{id:"formSmash:upper:j_idt511",widgetVar:"widget_formSmash_upper_j_idt511",behaviors:{change:function(ext) {PrimeFaces.ab({s:"formSmash:upper:j_idt511",e:"change",f:"formSmash",p:"formSmash:upper:j_idt511",u:"formSmash:upper:otherLanguage"},ext);}}});});

- de-DE
- en-GB
- en-US
- fi-FI
- nn-NO
- nn-NB
- sv-SE
- Other locale

Output formathtml text asciidoc rtf $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_upper_j_idt521",{id:"formSmash:upper:j_idt521",widgetVar:"widget_formSmash_upper_j_idt521"});});

- html
- text
- asciidoc
- rtf

Rows per page

- 5
- 10
- 20
- 50
- 100
- 250

Sort

- Standard (Relevance)
- Author A-Ö
- Author Ö-A
- Title A-Ö
- Title Ö-A
- Publication type A-Ö
- Publication type Ö-A
- Issued (Oldest first)
- Issued (Newest first)
- Created (Oldest first)
- Created (Newest first)
- Last updated (Oldest first)
- Last updated (Newest first)
- Disputation date (earliest first)
- Disputation date (latest first)

- Standard (Relevance)
- Author A-Ö
- Author Ö-A
- Title A-Ö
- Title Ö-A
- Publication type A-Ö
- Publication type Ö-A
- Issued (Oldest first)
- Issued (Newest first)
- Created (Oldest first)
- Created (Newest first)
- Last updated (Oldest first)
- Last updated (Newest first)
- Disputation date (earliest first)
- Disputation date (latest first)

Select

The maximal number of hits you can export is 250. When you want to export more records please use the Create feeds function.

1. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_0_j_idt584",{id:"formSmash:items:resultList:0:j_idt584",widgetVar:"widget_formSmash_items_resultList_0_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.). KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:0:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:0:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); BROYDEN'S METHOD FOR NONLINEAR EIGENPROBLEMS2019In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 41, no 2, p. A989-A1012Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_0_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:0:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_0_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Broyden's method is a general method commonly used for nonlinear systems of equations when very little information is available about the problem. We develop an approach based on Broyden's method for the structure appearing in nonlinear eigenvalue problems. Our approach is designed for problems where the evaluation of a matrix vector product is computationally expensive, essentially as expensive as solving the corresponding linear system of equations. We show how the structure of the Jacobian matrix can be incorporated into the algorithm to improve convergence. The algorithm exhibits local superlinear convergence for simple eigenvalues, and we characterize the convergence. We show how deflation can be integrated and combined such that the method can be used to compute several eigenvalues. A specific problem in machine tool milling, coupled with a PDE, is used to illustrate the approach. The simulations were carried out using the Julia programming language, and the simulation software is provided publicly for reproducibility.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:0:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 2. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_1_j_idt584",{id:"formSmash:items:resultList:1:j_idt584",widgetVar:"widget_formSmash_items_resultList_1_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_1_j_idt587",{id:"formSmash:items:resultList:1:j_idt587",widgetVar:"widget_formSmash_items_resultList_1_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.). KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:1:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Poloni, F.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:1:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Iterative methods for the delay Lyapunov equation with T-Sylvester preconditioning2019In: Applied Numerical Mathematics, ISSN 0168-9274, E-ISSN 1873-5460, Vol. 135, p. 173-185Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_1_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:1:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_1_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The delay Lyapunov equation is an important matrix boundary-value problem which arises as an analogue of the Lyapunov equation in the study of time-delay systems x˙(t)=A0x(t)+A1x(t−τ)+B0u(t). We propose a new algorithm for the solution of the delay Lyapunov equation. Our method is based on the fact that the delay Lyapunov equation can be expressed as a linear system of equations, whose unknown is the value U(τ/2)∈Rn×n, i.e., the delay Lyapunov matrix at time τ/2. This linear matrix equation with n2 unknowns is solved by adapting a preconditioned iterative method such as GMRES. The action of the n2×n2 matrix associated to this linear system can be computed by solving a coupled matrix initial-value problem. A preconditioner for the iterative method is proposed based on solving a T-Sylvester equation MX+XTN=C, for which there are methods available in the literature. We prove that the preconditioner is effective under certain assumptions. The efficiency of the approach is illustrated by applying it to a time-delay system stemming from the discretization of a partial differential equation with delay. Approximate solutions to this problem can be obtained for problems of size up to n≈1000, i.e., a linear system with n2≈106 unknowns, a dimension which is outside of the capabi© 2018 IMACS

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:1:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 3. Koskela, Antti PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_2_j_idt584",{id:"formSmash:items:resultList:2:j_idt584",widgetVar:"widget_formSmash_items_resultList_2_j_idt584",onLabel:"Koskela, Antti ",offLabel:"Koskela, Antti ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_2_j_idt587",{id:"formSmash:items:resultList:2:j_idt587",widgetVar:"widget_formSmash_items_resultList_2_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:2:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Jarlebring, EliasKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:2:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); On a generalization of neumann series of bessel functions using Hessenberg matrices and matrix exponentials2019In: European Conference on Numerical Mathematics and Advanced Applications, ENUMATH 2017, Springer, 2019, Vol. 126, p. 205-214Conference paper (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_2_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:2:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_2_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The Neumann expansion of Bessel functions (of integer order) of a function g: ℂ→ ℂ corresponds to representing g as a linear combination of basis functions φ0, φ1, …, i.e., g(s)=∑ℓ=0 ∞wℓφℓ(s), where φi(s) = Ji(s), i = 0, …, are the Bessel functions. In this work, we study an expansion for a more general class of basis functions. More precisely, we assume that the basis functions satisfy an infinite dimensional linear ordinary differential equation associated with a Hessenberg matrix, motivated by the fact that these basis functions occur in certain iterative methods. A procedure to compute the basis functions as well as the coefficients is proposed. Theoretical properties of the expansion are studied. We illustrate that non-standard basis functions can give faster convergence than the Bessel functions.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:2:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 4. Araujo-Cabarcas, J. C. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_3_j_idt587",{id:"formSmash:items:resultList:3:j_idt587",widgetVar:"widget_formSmash_items_resultList_3_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:3:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Engström, C.Jarlebring, EliasKTH, Centres, SeRC - Swedish e-Science Research Centre. KTH, School of Engineering Sciences (SCI), Mathematics (Dept.).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:3:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Efficient resonance computations for Helmholtz problems based on a Dirichlet-to-Neumann map2018In: Journal of Computational and Applied Mathematics, ISSN 0377-0427, E-ISSN 1879-1778, Vol. 330, p. 177-192Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_3_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:3:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_3_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We present an efficient procedure for computing resonances and resonant modes of Helmholtz problems posed in exterior domains. The problem is formulated as a nonlinear eigenvalue problem (NEP), where the nonlinearity arises from the use of a Dirichlet-to-Neumann map, which accounts for modeling unbounded domains. We consider a variational formulation and show that the spectrum consists of isolated eigenvalues of finite multiplicity that only can accumulate at infinity. The proposed method is based on a high order finite element discretization combined with a specialization of the Tensor Infinite Arnoldi method (TIAR). Using Toeplitz matrices, we show how to specialize this method to our specific structure. In particular we introduce a pole cancellation technique in order to increase the radius of convergence for computation of eigenvalues that lie close to the poles of the matrix-valued function. The solution scheme can be applied to multiple resonators with a varying refractive index that is not necessarily piecewise constant. We present two test cases to show stability, performance and numerical accuracy of the method. In particular the use of a high order finite element discretization together with TIAR results in an efficient and reliable method to compute resonances.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:3:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 5. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_4_j_idt584",{id:"formSmash:items:resultList:4:j_idt584",widgetVar:"widget_formSmash_items_resultList_4_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_4_j_idt587",{id:"formSmash:items:resultList:4:j_idt587",widgetVar:"widget_formSmash_items_resultList_4_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.). KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:4:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Mele, GiampaoloKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.Palitta, DavideUniv Bologna, Dipartimento Matemat, Piazza Porta S Donato,5, I-40127 Bologna, Italy..Ringh, EmilKTH, School of Engineering Sciences (SCI), Mathematics (Dept.). KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:4:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Krylov methods for low-rank commuting generalized Sylvester equations2018In: Numerical Linear Algebra with Applications, ISSN 1070-5325, E-ISSN 1099-1506, Vol. 25, no 6, article id e2176Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_4_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:4:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_4_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We consider generalizations of the Sylvester matrix equation, consisting of the sum of a Sylvester operator and a linear operator pi with a particular structure. More precisely, the commutators of the matrix coefficients of the operator pi and the Sylvester operator coefficients are assumed to be matrices with low rank. We show (under certain additional conditions) low-rank approximability of this problem, that is, the solution to this matrix equation can be approximated with a low-rank matrix. Projection methods have successfully been used to solve other matrix equations with low-rank approximability. We propose a new projection method for this class of matrix equations. The choice of the subspace is a crucial ingredient for any projection method for matrix equations. Our method is based on an adaption and extension of the extended Krylov subspace method for Sylvester equations. A constructive choice of the starting vector/block is derived from the low-rank commutators. We illustrate the effectiveness of our method by solving large-scale matrix equations arising from applications in control theory and the discretization of PDEs. The advantages of our approach in comparison to other methods are also illustrated.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:4:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 6. Mele, Giampaolo PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_5_j_idt584",{id:"formSmash:items:resultList:5:j_idt584",widgetVar:"widget_formSmash_items_resultList_5_j_idt584",onLabel:"Mele, Giampaolo ",offLabel:"Mele, Giampaolo ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_5_j_idt587",{id:"formSmash:items:resultList:5:j_idt587",widgetVar:"widget_formSmash_items_resultList_5_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:5:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Jarlebring, EliasKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:5:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); On restarting the tensor infinite Arnoldi method2018In: BIT Numerical Mathematics, ISSN 0006-3835, E-ISSN 1572-9125, Vol. 58, no 1, p. 133-162Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_5_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:5:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_5_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); An efficient and robust restart strategy is important for any Krylov-based method for eigenvalue problems. The tensor infinite Arnoldi method (TIAR) is a Krylov-based method for solving nonlinear eigenvalue problems (NEPs). This method can be interpreted as an Arnoldi method applied to a linear and infinite dimensional eigenvalue problem where the Krylov basis consists of polynomials. We propose new restart techniques for TIAR and analyze efficiency and robustness. More precisely, we consider an extension of TIAR which corresponds to generating the Krylov space using not only polynomials, but also structured functions, which are sums of exponentials and polynomials, while maintaining a memory efficient tensor representation. We propose two restarting strategies, both derived from the specific structure of the infinite dimensional Arnoldi factorization. One restarting strategy, which we call semi-explicit TIAR restart, provides the possibility to carry out locking in a compact way. The other strategy, which we call implicit TIAR restart, is based on the Krylov–Schur restart method for the linear eigenvalue problem and preserves its robustness. Both restarting strategies involve approximations of the tensor structured factorization in order to reduce the complexity and the required memory resources. We bound the error introduced by some of the approximations in the infinite dimensional Arnoldi factorization showing that those approximations do not substantially influence the robustness of the restart approach. We illustrate the effectiveness of the approaches by applying them to solve large scale NEPs that arise from a delay differential equation and a wave propagation problem. The advantages in comparison to other restart methods are also illustrated.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:5:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 7. Ringh, Emil PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_6_j_idt584",{id:"formSmash:items:resultList:6:j_idt584",widgetVar:"widget_formSmash_items_resultList_6_j_idt584",onLabel:"Ringh, Emil ",offLabel:"Ringh, Emil ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_6_j_idt587",{id:"formSmash:items:resultList:6:j_idt587",widgetVar:"widget_formSmash_items_resultList_6_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Optimization and Systems Theory.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:6:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Mele, GiampaoloKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.Karlsson, JohanKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Optimization and Systems Theory.Jarlebring, EliasKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:6:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Sylvester-based preconditioning for the waveguide eigenvalue problem2018In: Linear Algebra and its Applications, ISSN 0024-3795, E-ISSN 1873-1856, Vol. 542, no 1, p. 441-463Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_6_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:6:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_6_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We consider a nonlinear eigenvalue problem (NEP) arising from absorbing boundary conditions in the study of a partial differential equation (PDE) describing a waveguide. We propose a new computational approach for this large-scale NEP based on residual inverse iteration (Resinv) with preconditioned iterative solves. Similar to many preconditioned iterative methods for discretized PDEs, this approach requires the construction of an accurate and efficient preconditioner. For the waveguide eigenvalue problem, the associated linear system can be formulated as a generalized Sylvester equation AX+XB+A1XB1+A2XB2+K(ring operator)X=C, where (ring operator) denotes the Hadamard product. The equation is approximated by a low-rank correction of a Sylvester equation, which we use as a preconditioner. The action of the preconditioner is efficiently computed by using the matrix equation version of the Sherman-Morrison-Woodbury (SMW) formula. We show how the preconditioner can be integrated into Resinv. The results are illustrated by applying the method to large-scale problems.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:6:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); The full text will be freely available from 2020-05-03 09:14$(function(){PrimeFaces.cw("Tooltip","widget_formSmash_items_resultList_6_j_idt853_0_j_idt856",{id:"formSmash:items:resultList:6:j_idt853:0:j_idt856",widgetVar:"widget_formSmash_items_resultList_6_j_idt853_0_j_idt856",showEffect:"fade",hideEffect:"fade",target:"formSmash:items:resultList:6:j_idt853:0:fullTextSvg"});}); 8. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_7_j_idt584",{id:"formSmash:items:resultList:7:j_idt584",widgetVar:"widget_formSmash_items_resultList_7_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_7_j_idt587",{id:"formSmash:items:resultList:7:j_idt587",widgetVar:"widget_formSmash_items_resultList_7_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:7:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Koskela, AnttiKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.Mele, GiampaoloKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:7:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); *Disguised and new quasi-Newton methods for nonlinear eigenvalue problems*2017In: Numerical Algorithms, ISSN 1017-1398, E-ISSN 1572-9265Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_7_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:7:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_7_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this paper, we take a quasi-Newton approach to nonlinear eigenvalue problems (NEPs) of the type M(λ)v = 0, where (Formula presented.) is a holomorphic function. We investigate which types of approximations of the Jacobian matrix lead to competitive algorithms, and provide convergence theory. The convergence analysis is based on theory for quasi-Newton methods and Keldysh’s theorem for NEPs. We derive new algorithms and also show that several well-established methods for NEPs can be interpreted as quasi-Newton methods, and thereby, we provide insight to their convergence behavior. In particular, we establish quasi-Newton interpretations of Neumaier’s residual inverse iteration and Ruhe’s method of successive linear problems.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:7:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 9. Gaaf, Sarah W. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_8_j_idt584",{id:"formSmash:items:resultList:8:j_idt584",widgetVar:"widget_formSmash_items_resultList_8_j_idt584",onLabel:"Gaaf, Sarah W. ",offLabel:"Gaaf, Sarah W. ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_8_j_idt587",{id:"formSmash:items:resultList:8:j_idt587",widgetVar:"widget_formSmash_items_resultList_8_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); TU Eindhoven.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:8:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Jarlebring, EliasKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:8:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); The infinite Bi-Lanczos method for nonlinear eigenvalue problems2017In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 39, p. S898-S919Article in journal (Refereed)10. Beeumen, R. V. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_9_j_idt587",{id:"formSmash:items:resultList:9:j_idt587",widgetVar:"widget_formSmash_items_resultList_9_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:9:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Jarlebring, EliasKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.Michiels, W.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:9:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A rank-exploiting infinite Arnoldi algorithm for nonlinear eigenvalue problems2016In: Numerical Linear Algebra with Applications, ISSN 1070-5325, E-ISSN 1099-1506Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_9_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:9:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_9_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Summary: We consider the nonlinear eigenvalue problem M(λ)x = 0, where M(λ) is a large parameter-dependent matrix. In several applications, M(λ) has a structure where the higher-order terms of its Taylor expansion have a particular low-rank structure. We propose a new Arnoldi-based algorithm that can exploit this structure. More precisely, the proposed algorithm is equivalent to Arnoldi's method applied to an operator whose reciprocal eigenvalues are solutions to the nonlinear eigenvalue problem. The iterates in the algorithm are functions represented in a particular structured vector-valued polynomial basis similar to the construction in the infinite Arnoldi method [Jarlebring, Michiels, and Meerbergen, Numer. Math., 122 (2012), pp. 169-195]. In this paper, the low-rank structure is exploited by applying an additional operator and by using a more compact representation of the functions. This reduces the computational cost associated with orthogonalization, as well as the required memory resources. The structure exploitation also provides a natural way in carrying out implicit restarting and locking without the need to impose structure in every restart. The efficiency and properties of the algorithm are illustrated with two large-scale problems.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:9:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 11. Koskela, Antti PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_10_j_idt584",{id:"formSmash:items:resultList:10:j_idt584",widgetVar:"widget_formSmash_items_resultList_10_j_idt584",onLabel:"Koskela, Antti ",offLabel:"Koskela, Antti ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_10_j_idt587",{id:"formSmash:items:resultList:10:j_idt587",widgetVar:"widget_formSmash_items_resultList_10_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:10:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Jarlebring, EliasKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.Hochstenbach, M. E.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:10:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Krylov approximation of linear odes with polynomial parameterization2016In: SIAM Journal on Matrix Analysis and Applications, ISSN 0895-4798, E-ISSN 1095-7162, Vol. 37, no 2, p. 519-538Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_10_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:10:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_10_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We propose a new numerical method to solve linear ordinary differential equations of the type δu/δt (t, ϵ) = A(ϵ) u(t,ϵ), where A: C → Cn×n is a matrix polynomial with large and sparse matrix coefficients. The algorithm computes an explicit parameterization of approximations of u(t, ϵ) such that approximations for many different values of ϵ and t can be obtained with a very small additional computational effort. The derivation of the algorithm is based on a reformulation of the parameterization as a linear parameter-free ordinary differential equation and on approximating the product of the matrix exponential and a vector with a Krylov method. The Krylov approximation is generated with Arnoldi's method and the structure of the coefficient matrix turns out to be independent of the truncation parameter so that it can also be interpreted as Arnoldi's method applied to an infinite dimensional matrix. We prove the super linear convergence of the algorithm and provide a posteriori error estimates to be used as termination criteria. The behavior of the algorithm is illustrated with examples stemming from spatial discretizations of partial differential equations.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:10:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 12. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_11_j_idt584",{id:"formSmash:items:resultList:11:j_idt584",widgetVar:"widget_formSmash_items_resultList_11_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_11_j_idt587",{id:"formSmash:items:resultList:11:j_idt587",widgetVar:"widget_formSmash_items_resultList_11_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:11:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Guettel, StefanPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:11:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A SPATIALLY ADAPTIVE ITERATIVE METHOD FOR A CLASS OF NONLINEAR OPERATOR EIGENPROBLEMS2014In: Electronic Transactions on Numerical Analysis, ISSN 1068-9613, E-ISSN 1068-9613, Vol. 41, p. 21-41Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_11_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:11:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_11_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We present a new algorithm for the iterative solution of nonlinear operator eigenvalue problems arising from partial differential equations (PDEs). This algorithm combines automatic spatial resolution of linear operators with the infinite Arnoldi method for nonlinear matrix eigenproblems proposed by Jarlebring et al. [Numer. Math., 122 (2012), pp. 169-195]. The iterates in this infinite Arnoldi method are functions, and each iteration requires the solution of an inhomogeneous differential equation. This formulation is independent of the spatial representation of the functions, which allows us to employ a dynamic representation with an accuracy of about the level of machine precision at each iteration similar to what is done in the Chebfun system with its chebop functionality although our function representation is entirely based on coefficients instead of function values. Our approach also allows nonlinearity in the boundary conditions of the PDE. The algorithm is illustrated with several examples, e.g., the study of eigenvalues of a vibrating string with delayed boundary feedback control.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:11:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 13. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_12_j_idt584",{id:"formSmash:items:resultList:12:j_idt584",widgetVar:"widget_formSmash_items_resultList_12_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_12_j_idt587",{id:"formSmash:items:resultList:12:j_idt587",widgetVar:"widget_formSmash_items_resultList_12_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:12:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Güttel, StefanPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:12:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A spatially adaptive iterative method for a class of nonlinear operator eigenvalue problems2014In: Electronic Transactions on Numerical Analysis, ISSN 1068-9613, E-ISSN 1068-9613, Vol. 41, p. 21-42-Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_12_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:12:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_12_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We present a new algorithm for the iterative solution of nonlinear operator eigenvalue problems arising from partial differential equations (PDEs). This algorithm combines automatic spatial resolution of linear operators with the infinite Arnoldi method for nonlinear matrix eigenproblems proposed by Jarlebring et al. [Numer.Math., 122 (2012), pp. 169–195]. The iterates in this infinite Arnoldi method are functions, and each iteration requires the solution of an inhomogeneous differential equation. This formulation is independent of the spatial representation of the functions, which allows us to employ a dynamic representation with an accuracy of about the level of machine precision at each iteration similar to what is done in the Chebfun system with its chebop functionality although our function representation is entirely based on coefficients instead of function values. Our approach also allows nonlinearity in the boundary conditions of the PDE. The algorithm is illustrated with several examples, e.g., the study of eigenvalues of a vibrating string with delayed boundary feedback control.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:12:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 14. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_13_j_idt584",{id:"formSmash:items:resultList:13:j_idt584",widgetVar:"widget_formSmash_items_resultList_13_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_13_j_idt587",{id:"formSmash:items:resultList:13:j_idt587",widgetVar:"widget_formSmash_items_resultList_13_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:13:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Kvaal, SimenMichiels, WimPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:13:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); An Inverse Iteration Method for Eigenvalue Problems with Eigenvector Nonlinearities2014In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 36, no 4, p. A1978-A2001Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_13_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:13:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_13_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Consider a symmetric matrix A(v) is an element of R-nxn depending on a vector v is an element of R-n and satisfying the property A(alpha v) = A(v) for any alpha is an element of R\{0}. We will here study the problem of finding (lambda,v) is an element of R x R-n\{0} such that (lambda,v) is an eigenpair of the matrix A(v) and we propose a generalization of inverse iteration for eigenvalue problems with this type of eigenvector nonlinearity. The convergence of the proposed method is studied and several convergence properties are shown to be analogous to inverse iteration for standard eigenvalue problems, including local convergence properties. The algorithm is also shown to be equivalent to a particular discretization of an associated ordinary differential equation, if the shift is chosen in a particular way. The algorithm is adapted to a variant of the Schrodinger equation known as the Gross-Pitaevskii equation. We use numerical simulations to illustrate the convergence properties, as well as the efficiency of the algorithm and the adaption.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:13:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 15. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_14_j_idt584",{id:"formSmash:items:resultList:14:j_idt584",widgetVar:"widget_formSmash_items_resultList_14_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_14_j_idt587",{id:"formSmash:items:resultList:14:j_idt587",widgetVar:"widget_formSmash_items_resultList_14_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:14:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Meerbergen, KarlKU Leuven.Michiels, WimKU Leuven.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:14:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Computing a partial schur factorization of nonlinear eigenvalue problems using the infinite Arnoldi method2014In: SIAM Journal on Matrix Analysis and Applications, ISSN 0895-4798, E-ISSN 1095-7162, Vol. 35, no 2, p. 411-436Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_14_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:14:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_14_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The partial Schur factorization can be used to represent several eigenpairs of a matrix in a numerically robust way. Different adaptions of the Arnoldi method are often used to compute partial Schur factorizations. We propose here a technique to compute a partial Schur factorization of a nonlinear eigenvalue problem (NEP). The technique is an extension of our algorithm from [E. Jarlebring, W. Michiels, and K. Meerbergen, Numer. Math., 122 (2012), pp. 169-195], now called the infinite Arnoldi method. The infinite Arnoldi method is a method designed for NEPs, and can be interpreted as Arnoldi's method applied to a linear infinite-dimensional operator, whose reciprocal eigenvalues are the solutions to the NEP. As a first result we show that the invariant pairs of the operator are equivalent to invariant pairs of the NEP. We characterize the structure of the invariant pairs of the operator and show how one can carry out a modification of the infinite Arnoldi method by respecting the structure. This also allows us to naturally add the feature known as locking. We nest this algorithm with an outer iteration, where the infinite Arnoldi method for a particular type of structured functions is appropriately restarted. The restarting exploits the structure and is inspired by the well-known implicitly restarted Arnoldi method for standard eigenvalue problems. The final algorithm is applied to examples from a benchmark collection, showing that both processing time and memory consumption can be considerably reduced with the restarting technique.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:14:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 16. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_15_j_idt584",{id:"formSmash:items:resultList:15:j_idt584",widgetVar:"widget_formSmash_items_resultList_15_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_15_j_idt587",{id:"formSmash:items:resultList:15:j_idt587",widgetVar:"widget_formSmash_items_resultList_15_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:15:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Damm, TobiasMichiels, WimPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:15:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Model reduction of time-delay systems using position balancing and delay Lyapunov equations2013In: MCSS. Mathematics of Control, Signals and Systems, ISSN 0932-4194, E-ISSN 1435-568X, Vol. 25, no 2, p. 147-166Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_15_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:15:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_15_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Balanced truncation is a standard and very natural approach to approximate dynamical systems. We present a version of balanced truncation for model order reduction of linear time-delay systems. The procedure is based on a coordinate transformation of the position and preserves the delay structure of the system. We therefore call it (structure-preserving) position balancing. To every position, we associate quantities representing energies for the controllability and observability of the position. We show that these energies can be expressed explicitly in terms of the solutions to corresponding delay Lyapunov equations. Apart from characterizing the energies, we show that one block of the (operator) controllability and observability Gramians in the operator formulation of the time-delay system can also be characterized with the delay Lyapunov equation. The delay Lyapunov equation undergoes a contragredient transformation when we apply the position coordinate transformation and we propose to truncate it in a classical fashion, such that positions which are only weakly connected to the input and the output in the sense of the energy concepts are removed.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:15:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 17. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_16_j_idt584",{id:"formSmash:items:resultList:16:j_idt584",widgetVar:"widget_formSmash_items_resultList_16_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_16_j_idt587",{id:"formSmash:items:resultList:16:j_idt587",widgetVar:"widget_formSmash_items_resultList_16_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:16:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Michiels, WimMeerbergen, KarlPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:16:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A linear eigenvalue algorithm for the nonlinear eigenvalue problem2012In: Numerische Mathematik, ISSN 0029-599X, E-ISSN 0945-3245, Vol. 122, no 1, p. 169-195Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_16_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:16:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_16_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The Arnoldi method for standard eigenvalue problems possesses several attractive properties making it robust, reliable and efficient for many problems. The first result of this paper is a characterization of the solutions to an arbitrary (analytic) nonlinear eigenvalue problem (NEP) as the reciprocal eigenvalues of an infinite dimensional operator denoted . We consider the Arnoldi method for the operator and show that with a particular choice of starting function and a particular choice of scalar product, the structure of the operator can be exploited in a very effective way. The structure of the operator is such that when the Arnoldi method is started with a constant function, the iterates will be polynomials. For a large class of NEPs, we show that we can carry out the infinite dimensional Arnoldi algorithm for the operator in arithmetic based on standard linear algebra operations on vectors and matrices of finite size. This is achieved by representing the polynomials by vector coefficients. The resulting algorithm is by construction such that it is completely equivalent to the standard Arnoldi method and also inherits many of its attractive properties, which are illustrated with examples.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:16:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 18. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_17_j_idt584",{id:"formSmash:items:resultList:17:j_idt584",widgetVar:"widget_formSmash_items_resultList_17_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA (closed 2012-06-30).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:17:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:17:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Convergence factors of Newton methods for nonlinear eigenvalue problems2012In: Linear Algebra and its Applications, ISSN 0024-3795, E-ISSN 1873-1856, Vol. 436, no 10, p. 3943-3953Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_17_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:17:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_17_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Consider a complex sequence convergent to

*λ*_{∗}∈*C*with order*p*∈*N*. The*convergence factor*is typically defined as the fraction*c*_{k}:=(*λ*_{k+1}-*λ*_{∗})/(*λ*_{k}-*λ*_{∗})^{p}in the limit*k*→*∞*. In this paper, we prove formulas characterizing*c*_{k}in the limit*k*→*∞*for two different Newton-type methods for nonlinear eigenvalue problems. The formulas are expressed in terms of the left and right eigenvectors.The two treated methods are called the

*method of successive linear problems*(MSLP) and*augmented Newton*and are widely used in the literature. We prove several explicit formulas for*c*_{k}for both methods. Formulas for both methods are found for simple as well as double eigenvalues. In some cases, we observe in examples that the limit*c*_{k}as*k*→*∞*does not exist. For cases where this limit does not appear to exist, we prove other limiting expressions such that a characterization of*c*_{k}in the limit is still possible.PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:17:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 19. Saadvandi, M. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_18_j_idt587",{id:"formSmash:items:resultList:18:j_idt587",widgetVar:"widget_formSmash_items_resultList_18_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:18:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Meerbergen, K.Jarlebring, EliasPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:18:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); On dominant poles and model reduction of second order time-delay systems2012In: Applied Numerical Mathematics, ISSN 0168-9274, E-ISSN 1873-5460, Vol. 62, no 1, p. 21-34Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_18_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:18:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_18_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The method known as the dominant pole algorithm (DPA) has previously been successfully used in combination with model order reduction techniques to approximate standard linear time-invariant dynamical systems and second order dynamical systems. In this paper, we show how this approach can be adapted to a class of second order delay systems, which are large scale nonlinear problems whose transfer functions have an infinite number of simple poles. Deflation is a very important ingredient for this type of methods. Because of the nonlinearity, many deflation approaches for linear systems are not applicable. We therefore propose an alternative technique that essentially removes computed poles from the systemʼs input and output vectors. In general, this technique changes the residues, and hence, modifies the order of dominance of the poles, but we prove that, under certain conditions, the residues stay near the original residues. The new algorithm is illustrated by numerical examples.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:18:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 20. Michiels, W. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_19_j_idt587",{id:"formSmash:items:resultList:19:j_idt587",widgetVar:"widget_formSmash_items_resultList_19_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:19:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Jarlebring, EliasKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.Meerbergen, K.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:19:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A projection approach for model reduction of large-scale time-delay systems2011In: Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference, 2011, p. 1-8Conference paper (Refereed)21. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_20_j_idt584",{id:"formSmash:items:resultList:20:j_idt584",widgetVar:"widget_formSmash_items_resultList_20_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_20_j_idt587",{id:"formSmash:items:resultList:20:j_idt587",widgetVar:"widget_formSmash_items_resultList_20_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Katholieke Univ Leuven, B-3001 Heverlee, Belgium .PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:20:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Michiels, WimKatholieke Univ Leuven, B-3001 Heverlee, Belgium .PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:20:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Analyzing the convergence factor of residual inverse iteration2011In: BIT Numerical Mathematics, ISSN 0006-3835, E-ISSN 1572-9125, Vol. 51, no 4, p. 937-957Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_20_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:20:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_20_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We will establish here a formula for the convergence factor of the method called residual inverse iteration, which is a method for nonlinear eigenvalue problems and a generalization of the well-known inverse iteration. The formula for the convergence factor is explicit and involves quantities associated with the eigenvalue to which the iteration converges, in particular the eigenvalue and eigenvector. Residual inverse iteration allows for some freedom in the choice of a vector w (k) and we can use the formula for the convergence factor to analyze how it depends on the choice of w (k) . We also use the formula to illustrate the convergence when the shift is close to the eigenvalue. Finally, we explain the slow convergence for double eigenvalues by showing that under generic conditions, the convergence factor is one, unless the eigenvalue is semisimple. If the eigenvalue is semisimple, it turns out that we can expect convergence similar to the simple case.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:20:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 22. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_21_j_idt584",{id:"formSmash:items:resultList:21:j_idt584",widgetVar:"widget_formSmash_items_resultList_21_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_21_j_idt587",{id:"formSmash:items:resultList:21:j_idt587",widgetVar:"widget_formSmash_items_resultList_21_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Katholieke Univ Leuven, B-3001 Louvain, Belgium .PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:21:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Vanbiervliet, JKatholieke Univ Leuven, B-3001 Louvain, Belgium .Michiels, WimKatholieke Univ Leuven, B-3001 Louvain, Belgium .PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:21:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Characterizing and Computing the H(2) Norm of Time-Delay Systems by Solving the Delay Lyapunov Equation2011In: IEEE Transactions on Automatic Control, ISSN 0018-9286, E-ISSN 1558-2523, Vol. 56, no 4, p. 814-825Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_21_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:21:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_21_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); It is widely known that the solutions of Lyapunov equations can be used to compute the norm of linear time-invariant (LTI) dynamical systems. In this paper, we show how this theory extends to dynamical systems with delays. The first result is that the norm can be computed from the solution of a generalization of the Lyapunov equation, which is known as the delay Lyapunov equation. From the relation with the delay Lyapunov equation we can prove an explicit formula for the norm if the system has commensurate delays, here meaning that the delays are all integer multiples of a basic delay. The formula is explicit and contains only elementary linear algebra operations applied to matrices of finite dimension. The delay Lyapunov equations are matrix boundary value problems. We show how to apply a spectral discretization scheme to these equations for the general, not necessarily commensurate, case. The convergence of spectral methods typically depends on the smoothness of the solution. To this end we describe the smoothness of the solution to the delay Lyapunov equations, for the commensurate as well as for the non-commensurate case. The smoothness properties allow us to completely predict the convergence order of the spectral method

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:21:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 23. Kvaal, Simen et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_22_j_idt587",{id:"formSmash:items:resultList:22:j_idt587",widgetVar:"widget_formSmash_items_resultList_22_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:22:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Jarlebring, EliasDepartement Computerwetenschappen, Katholieke Universiteit Leuven.Michiels, WimPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:22:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Computing singularities of perturbation series2011In: Physical Review A. Atomic, Molecular, and Optical Physics, ISSN 1050-2947, E-ISSN 1094-1622, Vol. 83, no 3, p. 032505-Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_22_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:22:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_22_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Many properties of current ab initio approaches to the quantum many-body problem, both perturbational and otherwise, are related to the singularity structure of the Rayleigh-Schrodinger perturbation series. A numerical procedure is presented that in principle computes the complete set of singularities, including the dominant singularity which limits the radius of convergence. The method approximates the singularities as eigenvalues of a certain generalized eigenvalue equation which is solved using iterative techniques. It relies on computation of the action of the Hamiltonian matrix on a vector and does not rely on the terms in the perturbation series. The method can be useful for studying perturbation series of typical systems of moderate size, for fundamental development of resummation schemes, and for understanding the structure of singularities for typical systems. Some illustrative model problems are studied, including a helium-like model with delta-function interactions for which Moller-Plesset perturbation theory is considered and the radius of convergence found.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:22:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 24. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_23_j_idt584",{id:"formSmash:items:resultList:23:j_idt584",widgetVar:"widget_formSmash_items_resultList_23_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_23_j_idt587",{id:"formSmash:items:resultList:23:j_idt587",widgetVar:"widget_formSmash_items_resultList_23_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:23:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Kvaal, SMichiels, WPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:23:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Computing all pairs (lambda,mu) such that lambdais a double eigenvalue of A plus mu B2011In: SIAM Journal on Matrix Analysis and Applications, ISSN 0895-4798, E-ISSN 1095-7162, Vol. 32, no 3, p. 902-927Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_23_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:23:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_23_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Double eigenvalues are not generic for matrices without any particular structure. A matrix depending linearly on a scalar parameter, A + mu B, will, however, generically have double eigenvalues for some values of the parameter mu. In this paper, we consider the problem of finding those values. More precisely, we construct a method to accurately find all scalar pairs (lambda, mu) such that A + mu B has a double eigenvalue lambda, where A and B are given arbitrary complex matrices. The general idea of the globally convergent method is that if mu is close to a solution, then A + mu B has two eigenvalues which are close to each other. We fix the relative distance between these two eigenvalues and construct a method to solve and study it by observing that the resulting problem can be stated as a two-parameter eigenvalue problem, which is already studied in the literature. The method, which we call the method of fixed relative distance (MFRD), involves solving a two-parameter eigenvalue problem which returns approximations of all solutions. It is unfortunately not possible to get full accuracy with MFRD. In order to compute solutions with full accuracy, we present an iterative method which returns a very accurate solution, for a sufficiently good starting value. The approach is illustrated with one academic example and one application to a simple problem in computational quantum mechanics.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:23:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 25. Michiels, W. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_24_j_idt587",{id:"formSmash:items:resultList:24:j_idt587",widgetVar:"widget_formSmash_items_resultList_24_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:24:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Jarlebring, EliasKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.Meerbergen, K.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:24:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Krylov-Based Model Order Reduction of Time-delay Systems2011In: SIAM Journal on Matrix Analysis and Applications, ISSN 0895-4798, E-ISSN 1095-7162, Vol. 32, no 4, p. 1399-1421Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_24_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:24:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_24_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); We present a model order reduction method which allows the construction of a reduced, delay-free model of a given dimension for linear time-delay systems, whose characteristic matrix is nonlinear due to the presence of exponential functions. The method builds on the equivalent representation of the time-delay system as an infinite-dimensional linear problem. It combines ideas from a finite-dimensional approximation via a spectral discretization, on the one hand, and a Krylov–Padé model reduction approach, on the other hand. The method exhibits a good spectral approximation of the original model, in the sense that the smallest characteristic roots are well approximated and the nonconverged eigenvalues of the reduced model have a favorable location, and it preserves moments at zero and at infinity. The spectral approximation is due to an underlying Arnoldi process that relies on building an appropriate Krylov space for the linear infinite-dimensional problem. The preservation of moments is guaranteed, because the chosen finite-dimensional approximation preserves moments and, in addition, the space on which one projects is constructed in such a way that the preservation of moments carries over to the reduced model. The implementation of the method is dynamic, since the number of grid points in the spectral discretization does not need to be chosen beforehand and the accuracy of the reduced model can always be improved by doing more iterations. It relies on a reformulation of the problem involving a companion-like system matrix and a highly structured input matrix, whose structure are fully exploited.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:24:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 26. Vanbiervliet, J. et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_25_j_idt587",{id:"formSmash:items:resultList:25:j_idt587",widgetVar:"widget_formSmash_items_resultList_25_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:25:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Michiels, W.Jarlebring, EliasKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:25:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Using spectral discretization for the optimal H2 design of time-delay systems2011In: International Journal of Control, ISSN 0020-7179, E-ISSN 1366-5820, Vol. 84, no 2, p. 228-241Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_25_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:25:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_25_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The stabilization and robustification of a time-delay system is the topic of this paper. More precisely, we want to minimize the H2 norm of the transfer function corresponding to this class of linear time-invariant input-output systems with fixed time delays in the states. Due to the presence of the delays, the transfer function is a nonrational, nonlinear function, and the classical procedure which involves solving Lyapunov equations is no longer applicable. We therefore propose an approach based on a spectral discretization applied to a reformulation of the time-delay system as an infinite-dimensional standard linear system. In this way, we obtain a large delay-free system, which serves as an approximation to the original time-delay system, and which allows the application of standard H2 norm optimization techniques. We give an interpretation of this approach in the frequency domain and relate it to the approximation of the nonlinear terms in the time-delay transfer function by means of a rational function. Using this property, we can provide some insight in the convergence be- haviour of the approximation, justifying its use for the purpose of H2 norm computation. Along with this, the easy availibility of derivatives with respect to the original matrices allows for an efficient integration into any standard optimization framework. A numerical example finally illustrates how the presented method can be employed to perform optimal H2 norm design using smooth optimization techniques.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:25:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 27. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_26_j_idt584",{id:"formSmash:items:resultList:26:j_idt584",widgetVar:"widget_formSmash_items_resultList_26_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_26_j_idt587",{id:"formSmash:items:resultList:26:j_idt587",widgetVar:"widget_formSmash_items_resultList_26_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Katholieke Univ Leuven, Dept Comp Sci, B-3001 Heverlee, Belgium .PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:26:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Meerbergen, KarlKatholieke Univ Leuven, Dept Comp Sci, B-3001 Heverlee, Belgium .Michiels, WimKatholieke Univ Leuven, Dept Comp Sci, B-3001 Heverlee, Belgium .PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:26:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A Krylov method for the delay eigenvalue problem2010In: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197, Vol. 32, no 6, p. 3278-3300Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_26_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:26:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_26_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The Arnoldi method is currently a very popular algorithm to solve large-scale eigenvalue problems. The main goal of this paper is to generalize the Arnoldi method to the characteristic equation of a delay-differential equation (DDE), here called a delay eigenvalue problem (DEP). The DDE can equivalently be expressed with a linear infinite-dimensional operator whose eigenvalues are the solutions to the DEP. We derive a new method by applying the Arnoldi method to the generalized eigenvalue problem associated with a spectral discretization of the operator and by exploiting the structure. The result is a scheme where we expand a subspace not only in the traditional way done in the Arnoldi method. The subspace vectors are also expanded with one block of rows in each iteration. More important, the structure is such that if the Arnoldi method is started in an appropriate way, it has the (somewhat remarkable) property that it is in a sense independent of the number of discretization points. It is mathematically equivalent to an Arnoldi method with an infinite matrix, corresponding to the limit where we have an infinite number of discretization points. We also show an equivalence with the Arnoldi method in an operator setting. It turns out that with an appropriately defined operator over a space equipped with scalar product with respect to which Chebyshev polynomials are orthonormal, the vectors in the Arnoldi iteration can be interpreted as the coefficients in a Chebyshev expansion of a function. The presented method yields the same Hessenberg matrix as the Arnoldi method applied to the operator.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:26:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 28. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_27_j_idt584",{id:"formSmash:items:resultList:27:j_idt584",widgetVar:"widget_formSmash_items_resultList_27_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_27_j_idt587",{id:"formSmash:items:resultList:27:j_idt587",widgetVar:"widget_formSmash_items_resultList_27_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:27:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Meerbergen, KMichiels, WPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:27:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); An Arnoldi method with structured starting vectors for the delay eigenvalue problem2010In: Proceedings of the 9th IFAC Workshop on Time Delay Systems, 2010Conference paper (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_27_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:27:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_27_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); The method called Arnoldi is currently a very popular method to solve largescale eigenvalue problems. The general purpose of this paper is to generalize Arnoldi to the characteristic equation of a time-delay system, here called adelay eigenvalue problem . The presented generalization is mathematically equivalent to Arnoldi applied to the problem corresponding to a Taylor approximation of the exponential. Even though the derivation of the result is with a Taylor approximation, the constructed method can be implemented in such a way that it is independent of the Taylor truncation paramater N . This is achieved by exploiting properties of vectors with a special structure, the vectorization of a rank one matrix plus the vectorization of a matrix which right-most columns are zero. It turns out that this set of vectors is closed under the matrix vector product as well as orthogonalization. Moreover, both operations can be efficiently computed. Since Arnoldi only consists of these operations, if Arnoldi is started with the special vector structure, the method can be efficiently executed. The presented numerical experiments indicate that the method is very efficient in comparison to methods in the literature.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:27:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 29. Rott, Oliver et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_28_j_idt587",{id:"formSmash:items:resultList:28:j_idt587",widgetVar:"widget_formSmash_items_resultList_28_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:28:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Jarlebring, EliasKTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:28:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); An iterative method for the multipliers of periodic delaydifferential equations and the analysis of a PDE milling model2010In: Proceedings of the 9th IFAC Workshop on Time Delay Systems / [ed] Vyhlidal, Tomas; Zitek, Pavel, 2010, p. 1-6Conference paper (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_28_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:28:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_28_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Locally convergent iterative schemes have turned out to be very useful in the analysis of the characteristic roots of delay-differential equations (DDEs) with constant coefficients. In this work we present a locally convergent iterative scheme for the characteristic multipliers of periodic-coefficient DDEs. The method is an adaption of an iterative method called residual inverse iteration. The possibility to use this method stems from an observation that the characteristic matrix can be expressed with the fundamental solution of a differential equation. We apply the method to a coupled milling model containing a partial and an ordinary differential equation. The conclusion of the numerical results is that the stability diagram of the coupled model differs significantly from the combined stability diagrams for each subsystem.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:28:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 30. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_29_j_idt584",{id:"formSmash:items:resultList:29:j_idt584",widgetVar:"widget_formSmash_items_resultList_29_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_29_j_idt587",{id:"formSmash:items:resultList:29:j_idt587",widgetVar:"widget_formSmash_items_resultList_29_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:29:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Vanbiervliet, JMichiels, WimPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:29:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Explicit expressions for the H2 norm of time-delay systems based on the delay lyapunov equation2010In: Proceedings of the 49th IEEE Conference on Decision and Control, 2010Conference paper (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_29_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:29:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_29_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); It is widely known that the solutions of Lyapunov equations can be used to compute the H

_{2}norm of linear time-invariant (LTI) dynamical systems. In this paper, we show how this theory extends to dynamical systems with delays. We establish, in the main result, that the H_{2}norm can be computed from the solution of a generalization of the Lyapunov equation, which is known as the delay Lyapunov equation. From the relation with the delay Lyapunov equation we can prove an explicit formula for the H_{2}norm if the system has commensurate delays, here meaning that the delays are all integer multiples of a basic delay. The formula is explicit and contains only elementary linear algebra operations applied to matrices of finite dimension.PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:29:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 31. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_30_j_idt584",{id:"formSmash:items:resultList:30:j_idt584",widgetVar:"widget_formSmash_items_resultList_30_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_30_j_idt587",{id:"formSmash:items:resultList:30:j_idt587",widgetVar:"widget_formSmash_items_resultList_30_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Katholieke Univ Leuven, Dept Comp Sci, B-3001 Heverlee, Belgium .PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:30:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Michiels, WimKatholieke Univ Leuven, Dept Comp Sci, B-3001 Heverlee, Belgium .PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:30:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Invariance properties in the root sensitivity of time-delay systems with double imaginary roots2010In: Automatica, ISSN 0005-1098, E-ISSN 1873-2836, Vol. 46, no 6, p. 1112-1115Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_30_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:30:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_30_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); If 1 omega 1R is an eigenvalue of a time-delay system for the delay tau(0) then 1 omega o is also an eigenvalue for the delays tau(k) := tau(0) + k2 pi/omega, for any k is an element of Z We investigate the sensitivity. periodicity and invariance properties of the root 1 omega for the case that 1 omega is a double eigenvalue for some tau(k) It turns out that under natural conditions (the condition that the root exhibits the completely regular splunng property if the delay is perturbed), the presence of a double imaginary root ico for some delay tau(0) implies that 1 omega is a simple root for the other delays tau(k), k not equal 0. Moreover, we show how to characterize the root locus around 1 omega. The entire local root locus picture can be completely determined from the square root splitting of the double root. We separate the general picture into two cases depending on the sign of a single scalar constant; the imaginary part of the first coefficient in the square root expansion of the double eigenvalue. (C) 2010 Elsevier Ltd All rights reserved.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:30:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 32. Larsson, Pia et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_31_j_idt587",{id:"formSmash:items:resultList:31:j_idt587",widgetVar:"widget_formSmash_items_resultList_31_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:31:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Sochor, JanaKTH, School of Architecture and the Built Environment (ABE), Transport and Economics (closed 20110301), Traffic and Logistics (closed 20110301).Udin, ChristianSjöström, ThomasJarlebring, IsakPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:31:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); ITS and Telematic Services: Different Implementation Aspects2010Report (Other academic)33. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_32_j_idt584",{id:"formSmash:items:resultList:32:j_idt584",widgetVar:"widget_formSmash_items_resultList_32_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Tech Univ Carolo Wilhelmina Braunschweig, Inst Computat Math, D-38023 Braunschweig, Germany.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:32:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:32:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Critical delays and Polynomial Eigenvalue Problems2009In: Journal of Computational and Applied Mathematics, ISSN 0377-0427, E-ISSN 1879-1778, Vol. 224, no 1, p. 296-306Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_32_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:32:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_32_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this work we present a new method to compute the delays of delay-differential equations (DDEs), such that the DDE has a purely imaginary eigenvalue. For delay-differential equations with multiple delays, the critical curves or critical surfaces in delay space (that is, the set of delays where the DDE has a purely imaginary eigenvalue) are parameterized. We show how the method is related to other works in the field by treating the case where the delays are integer multiples of some delay value, i.e., commensurate delays. The parameterization is done by solving a quadratic eigenvalue problem which is constructed from the vectorization of a matrix equation and hence typically of large size. For commensurate delay-differential equations, the corresponding equation is a polynomial eigenvalue problem. As a special case of the proposed method, we find a closed form for a parameterization of the critical surface for the scalar case. We provide several examples with visualizations where the computation is done with some exploitation of the structure of eigenvalue problems. (c) 2008 Elsevier B.V. All rights reserved.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:32:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 34. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_33_j_idt584",{id:"formSmash:items:resultList:33:j_idt584",widgetVar:"widget_formSmash_items_resultList_33_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_33_j_idt587",{id:"formSmash:items:resultList:33:j_idt587",widgetVar:"widget_formSmash_items_resultList_33_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Katholieke Univ Leuven, Dept Comp Sci, B-3001 Heverlee, Belgium.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:33:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Hochstenbach, M.E.Eindhoven Univ Technol, Dept Math & Comp Sci, NL-5600 MB Eindhoven, Netherlands .PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:33:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Polynomial two-parameter eigenvalue problems and matrix pencil methods for stability of delay-differential equations2009In: Linear Algebra and its Applications, ISSN 0024-3795, E-ISSN 1873-1856, Vol. 431, no 3-4, p. 369-380Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_33_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:33:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_33_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Several recent methods used to analyze asymptotic stability of delay-differential equations (DDEs) involve determining the eigenvalues of a matrix, a matrix pencil or a matrix polynomial constructed by Kronecker products. Despite some similarities between the different types of these so-called matrix pencil methods, the general ideas used as well as the proofs differ considerably. Moreover, the available theory hardly reveals the relations between the different methods. In this work, a different derivation of various matrix pencil methods is presented using a unifying framework of a new type of eigenvalue problem: the polynomial two-parameter eigenvalue problem, of which the quadratic two-parameter eigenvalue problem is a special case. This framework makes it possible to establish relations between various seemingly different methods and provides further insight in the theory of matrix pencil methods. We also recognize a few new matrix pencil variants to determine DDE stability. Finally, the recognition of the new types of eigenvalue problem opens a door to efficient computation of DDE stability. (C) 2009 Elsevier Inc. All rights reserved.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:33:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 35. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_34_j_idt584",{id:"formSmash:items:resultList:34:j_idt584",widgetVar:"widget_formSmash_items_resultList_34_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:34:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:34:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); On Critical Delays for Linear Neutral Delay Systems2007In: Proceedings of the European Control Conference, 2007Conference paper (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_34_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:34:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_34_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this work we address the problem of finding the critical delays of a linear neutral delay system, i.e., the delays such that the system has a purely imaginary eigenvalue. Even though neutral delay systems exhibit some discontinuity properties with respect to changes in the delays an essential part in a non-conservative stability analysis with respect to changes in the delays, is the computation of the critical delays. We generalize previous results on critical delays and stability switches for retarded time-delay systems, under some minor assumptions on the delay system. The work starts with stating a general equivalence theorem between the spectrum and a matrix function condition.We show how this theorem can be applied to the commensurate timedelay system to compute the critical delays. It turns out that the resulting method is closely related to parts of the results of Fu, Niculescu and Chen[6]. For the incommensurate case we present a scheme which allows the computation of the critical curves, i.e., the points in delay-space for which the system has a purely imaginary eigenvalue. We apply the method to previously investigated examples, in order to provide a verification of the results, as well as to an example for which the stability picture is, to our knowledge, not yet known.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:34:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 36. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_35_j_idt584",{id:"formSmash:items:resultList:35:j_idt584",widgetVar:"widget_formSmash_items_resultList_35_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_35_j_idt587",{id:"formSmash:items:resultList:35:j_idt587",widgetVar:"widget_formSmash_items_resultList_35_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC), Numerical Analysis, NA.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:35:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Damm, TobiasPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:35:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); The Lambert W function and the spectrum of some multidimensional time-delay systems2007In: Automatica, ISSN 0005-1098, E-ISSN 1873-2836, Vol. 43, no 12, p. 2124-2128Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_35_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:35:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_35_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this note we find an explicit expression for the eigenvalues of a retarded time-delay system with one delay, for the special case that the system matrices are simultaneously triangularizable, which includes the case where they commute. Using matrix function definitions we define a matrix version of the Lambert W function, from which we form the expression. We prove by counter-example that some expressions in other publications on Lambert W for time-delay systems do not always hold. (c) 2007 Elsevier Ltd. All rights reserved.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:35:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 37. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_36_j_idt584",{id:"formSmash:items:resultList:36:j_idt584",widgetVar:"widget_formSmash_items_resultList_36_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); Institut Computational Mathematics, TU Braunschweig.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:36:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:36:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); A quadratic eigenproblem in the analysis of a time delay system2006In: Proceedings in Applied Mathematics and Mechanics: PAMM, ISSN 1617-7061, E-ISSN 1617-7061, Vol. 6, no 1, p. 63-66Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_36_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:36:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_36_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this work we solve a

*quadratic eigenvalue problem*occurring in a method to compute the set of delays of a*linear time delay system*(TDS) such that the system has an imaginary eigenvalue. The computationally dominating part of the method is to find all eigenvalues*z*of modulus one of the quadratic eigenvalue problemBecause of its origin in the vectorization of a

*Lyapunov type matrix equation*, the quadratic eigenvalue problem is, even for moderate size problems, of very large size. We show one way to treat this problem by exploiting the Lyapunov type structure of the quadratic eigenvalue problem when constructing an iterative solver. More precisely, we show that the*shift-invert operation*for the companion form of the quadratic eigenvalue problem can be efficiently computed by solving a*Sylvester equation*. The usefulness of this exploitation is demonstrated with an example.PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:36:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 38. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_37_j_idt584",{id:"formSmash:items:resultList:37:j_idt584",widgetVar:"widget_formSmash_items_resultList_37_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); TU Braunschweig, Germany.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:37:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:37:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Computing Critical delays for time delay systems with multiple delays2006Conference paper (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_37_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:37:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_37_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this work we present a method to analyze the robustness of stability of a time-delay system (TDS) with respect to the delays. This is done by computing the delays for which the system has a purely imaginary eigenvalue. These delays, called critical delays , generate potential points for a stability switch, i.e., the point where the system switches from a stable to unstable. To derive the method, we find a Lyapunov-type equation , equivalent to the characteristic equation of the TDS. Unlike the characteristic equation, the Lyapunov-type equation does not have any non-exponential terms if the eigenvalue is imaginary. This allows us to solve the Lyapunov-type equation by rewriting it to a quadratic eigenvalue problem for which there are efficient numerical methods. For the scalar case, we find a new explicit expression for the curves in the stability chart. The method is applied to previously solved examples as well as previously unsolved problems of larger dimension.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:37:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 39. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_38_j_idt584",{id:"formSmash:items:resultList:38:j_idt584",widgetVar:"widget_formSmash_items_resultList_38_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:38:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:38:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Computing the stability region in delay-space of a tds using polynomial eigenproblems2006In: Proceedings of the 6th IFAC Workshop on Time-Delay Systems, 2006Conference paper (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_38_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:38:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_38_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In this work we analyze stability properties of retarded linear time invariant multi-dimensional, multi-delay, time delay systems with respect to perturbations in the delay parameters. We analyze two methods which allow the computation of the critical delays, i.e., the points in delay-space which causes the system to have a purely imaginary eigenvalue. The critical delays are potential stability boundaries as the boundaries of the stability region is necessarily a subset of the critical delays.The two methods originates from a Lyapunov-type condition, which is completely self-contained in this work. The first method corresponds to the case of commensurate delays, for which the the Lyapunov-type condition reduces to a polynomial eigenvalue problem for which the first companion form is exactly the eigenvalue problem occurring in Chen et al. (1995). The second method is the result of a simple substitution which allows the computation of the critical delays of an incommensurate system by solving a quadratic eigenvalue problem. For the scalar multi-delay case we find a closed expression for the critical curves using this method. We confirm the methods by comparing it to previous work and published examples.

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:38:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 40. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_39_j_idt584",{id:"formSmash:items:resultList:39:j_idt584",widgetVar:"widget_formSmash_items_resultList_39_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_39_j_idt587",{id:"formSmash:items:resultList:39:j_idt587",widgetVar:"widget_formSmash_items_resultList_39_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Computer Science and Communication (CSC).PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:39:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Voss, HeinrichPrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:39:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Rational Krylov for Nonlinear Eigenproblems, an Iterative Projection Method2005In: Applications of Mathematics, ISSN 0862-7940, E-ISSN 1572-9109, Vol. 50, no 6, p. 543-554Article in journal (Refereed)Abstract [en] PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_39_j_idt622_0_j_idt623",{id:"formSmash:items:resultList:39:j_idt622:0:j_idt623",widgetVar:"widget_formSmash_items_resultList_39_j_idt622_0_j_idt623",onLabel:"Abstract [en]",offLabel:"Abstract [en]",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); In recent papers Ruhe [10], [12] suggested rational Krylov method for nonlinear eigenproblems knitting together a secant method for linearizing the nonlinear problem and the Krylov method for the linearized problem. In this note we point out that the method can be understood as an iterative projection method. Similar to the Arnoldi method presented in [13], [14] the search space is expanded by the direction from residual inverse iteration. Numercial methods demonstrate that the rational Krylov method can be accelerated considerably by replacing an inner iteration by an explicit solver of projected problems

PrimeFaces.cw("Panel","tryPanel",{id:"formSmash:items:resultList:39:j_idt622:0:abstractPanel",widgetVar:"tryPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); 41. Jarlebring, Elias PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_40_j_idt584",{id:"formSmash:items:resultList:40:j_idt584",widgetVar:"widget_formSmash_items_resultList_40_j_idt584",onLabel:"Jarlebring, Elias ",offLabel:"Jarlebring, Elias ",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); et al. PrimeFaces.cw("SelectBooleanButton","widget_formSmash_items_resultList_40_j_idt587",{id:"formSmash:items:resultList:40:j_idt587",widgetVar:"widget_formSmash_items_resultList_40_j_idt587",onLabel:"et al.",offLabel:"et al.",onIcon:"ui-icon-triangle-1-s",offIcon:"ui-icon-triangle-1-e"}); KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:40:orgPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Koskela, AnttiKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.Mele, GiampaoloKTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, NA. KTH, Centres, SeRC - Swedish e-Science Research Centre.PrimeFaces.cw("Panel","testPanel",{id:"formSmash:items:resultList:40:etAlPanel",widgetVar:"testPanel",toggleable:true,toggleSpeed:500,collapsed:false,toggleOrientation:"vertical",closable:true,closeSpeed:500}); Disguised and new quasi-Newton methods for nonlinear eigenvalue problemsIn: SIAM Journal on Scientific Computing, ISSN 1064-8275, E-ISSN 1095-7197Article in journal (Refereed)

CiteExportLink to result list
http://kth.diva-portal.org/smash/resultList.jsf?query=&language=en&searchType=SIMPLE&noOfRows=50&sortOrder=dateIssued_sort_desc&sortOrder2=title_sort_asc&onlyFullText=false&sf=all&aq=%5B%5B%7B%22author%22%3A%5B%22Jarlebring%22%5D%7D%5D%5D&aqe=%5B%5D&aq2=%5B%5B%5D%5D&af=%5B%5D $(function(){PrimeFaces.cw("InputTextarea","widget_formSmash_lower_j_idt902_recordPermLink",{id:"formSmash:lower:j_idt902:recordPermLink",widgetVar:"widget_formSmash_lower_j_idt902_recordPermLink",autoResize:true});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_lower_j_idt902_j_idt904",{id:"formSmash:lower:j_idt902:j_idt904",widgetVar:"widget_formSmash_lower_j_idt902_j_idt904",target:"formSmash:lower:j_idt902:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});

Permanent link

Cite

Citation styleapa harvard1 ieee modern-language-association-8th-edition vancouver Other style $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_lower_j_idt920",{id:"formSmash:lower:j_idt920",widgetVar:"widget_formSmash_lower_j_idt920",behaviors:{change:function(ext) {PrimeFaces.ab({s:"formSmash:lower:j_idt920",e:"change",f:"formSmash",p:"formSmash:lower:j_idt920",u:"formSmash:lower:otherStyle"},ext);}}});});

- apa
- harvard1
- ieee
- modern-language-association-8th-edition
- vancouver
- Other style

Languagede-DE en-GB en-US fi-FI nn-NO nn-NB sv-SE Other locale $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_lower_j_idt931",{id:"formSmash:lower:j_idt931",widgetVar:"widget_formSmash_lower_j_idt931",behaviors:{change:function(ext) {PrimeFaces.ab({s:"formSmash:lower:j_idt931",e:"change",f:"formSmash",p:"formSmash:lower:j_idt931",u:"formSmash:lower:otherLanguage"},ext);}}});});

- de-DE
- en-GB
- en-US
- fi-FI
- nn-NO
- nn-NB
- sv-SE
- Other locale

Output formathtml text asciidoc rtf $(function(){PrimeFaces.cw("SelectOneMenu","widget_formSmash_lower_j_idt941",{id:"formSmash:lower:j_idt941",widgetVar:"widget_formSmash_lower_j_idt941"});});

- html
- text
- asciidoc
- rtf