Conjugation on Varieties of Nilpotent Matrices

Total Page:16

File Type:pdf, Size:1020Kb

Conjugation on Varieties of Nilpotent Matrices Conjugation on varieties of nilpotent matrices Dissertation zur Erlangung des Doktorgrades der Naturwissenschaften im Fachbereich C der Bergischen Universität Wuppertal vorgelegt von Magdalena Boos im März 2012 Die Dissertation kann wie folgt zitiert werden: urn:nbn:de:hbz:468-20120710-103808-3 [http://nbn-resolving.de/urn/resolver.pl?urn=urn%3Anbn%3Ade%3Ahbz%3A468-20120710-103808-3] For my mum 4 Contents Introduction7 1 Theoretical background 15 1.1 Methods from Algebraic geometry and Invariant theory 15 1.1.1 Algebraic group actions 16 1.1.2 Invariants and algebraic quotients 17 1.1.3 Semi-invariants and GIT-quotients 19 1.1.4 Toric varieties 21 1.2 Representation theory of finite-dimensional algebras 22 1.2.1 Covering theory of quiver algebras 24 1.2.2 Tame and wild algebras 26 1.2.3 Degenerations 27 2 The concrete setup 29 2.1 (Oriented) Link patterns 29 2.2 Known results 31 2.2.1 Results of M. Jordan and M. Gerstenhaber 31 2.2.2 Results of A. Melnikov 32 2.2.3 Results of L. Hille and G. Röhrle 33 2.3 Representation-theoretic approach 34 3 Nilpotency degree 2 37 3.1 Classification of the orbits 37 3.1.1 Parabolic orbits 37 3.1.2 Unipotent orbits 47 3.1.3 The examples n=3 and n=4 49 3.2 Homomorphisms, extensions and their combinatorial interpretation 52 3.3 Closures of Borel orbits 56 3.3.1 Minimal, disjoint pieces of degenerations 57 3.3.2 Minimal degenerations in general 67 3.3.3 Dimensions and the open orbit 78 3.3.4 Minimal singularities 80 3.4 Closures of parabolic orbits 84 3.4.1 Minimal, disjoint degenerations 84 3.4.2 Dimensions of orbits 86 6 Contents 4 Finite classifications in higher nilpotency degrees 89 4.1 Maximal parabolic action for x=3 89 4.1.1 Classification of the orbits 89 4.1.2 Orbit closures 98 4.1.3 Dimensions and the open orbit 100 4.1.4 The parabolic subgroup of block sizes (2,2) 104 4.2 A finiteness criterion 106 4.3 A wildness criterion 113 4.3.1 The representation type of the corresponding algebras 113 4.3.2 Concrete 2-parameter families via tree modules 114 5 Generic classification in the nilpotent cone 119 5.1 Generic normal forms 119 5.2 (Semi-) Invariants 127 5.3 Generation of semi-invariant rings 129 6 Towards an algebraic U-quotient of the nilpotent cone 135 6.1 A quotient criterion for unipotent actions 135 6.1.1 The examples n=2 and n=3 137 6.2 Toric invariants 140 6.2.1 Reductions 141 6.2.2 General description of toric invariants 148 6.3 Generic separation of the U-orbits 155 6.4 The toric variety X 157 6.4.1 Toric operation(s) on the U-invariant ring 157 6.4.2 Explicit description of toric invariants 158 6.4.3 Interrelation between N//U and X 159 6.5 The case n=4 161 7 Towards a GIT-quotient for the Borel action 167 7.1 The examples n=2 and n=3 167 7.2 Generic separation of the same weight 168 7.3 Translation to the language of quiver moduli 173 A Appendix 177 A.1 Singular computations 177 A.2 The parabolic subgroup of block sizes (4,3) 179 Bibliography 181 Introduction The study of algebraic group actions on affine varieties, especially the "vertical" study of orbits and their closures, and the "horizontal" study of parametric families of orbits and quotients, are a common topic in algebraic Lie theory. A well-known example is the study of the adjoint action of a reductive algebraic group on its Lie algebra and numerous variants thereof, in particular the conjugacy classes of com- plex (nilpotent) square matrices. In 1870, the classification of the orbits by so-called Jordan normal forms was described by M. Jordan [Jordan, 1989, Jordan, 1871]. Their closures were described by M. Gersten- haber [Gerstenhaber, 1959] and W. Hesselink [Hesselink, 1976] in the second half of the twentieth century in terms of partitions and visualized by combinatorial objects named Young Diagrams. Algebraic group actions of reductive groups have particularly been discussed elaborately in connection with orbit spaces and more generally algebraic quotients, even though their application to concrete examples is far from being trivial. In case of a non-reductive group, even most of these results fail to hold true immediately. For example, Hilbert’s theorem [Hilbert, 1890] yields that for reductive groups, the invari- ant ring is finitely generated; and a criterion for algebraic quotients is valid [Kraft, 1984]. In 1958, though, M. Nagata [Nagata, 1960, Nagata, 1959] constructed a counterexample of a not finitely generated invariant ring corresponding to a non-reductive algebraic group action, which answered Hilbert’s fourteenth problem in the negative. One exception are algebraic actions of unipotent subgroups that are induced by reductive groups, since the corresponding invariant ring is always finitely generated [Kraft, 1984]. Our main attention in this work is turned towards algebraic non-reductive group actions that are induced by the conjugation action of the general linear group GLn over C. For ex- ample, the standard parabolic subgroups P (and, therefore, the Borel subgroup B) and the unipotent subgroup U of GLn are not reductive. It suggests itself to consider their action (x) on the variety Nn of x-nilpotent matrices of square size n via conjugation. We discuss this setup in detail and, thereby, generalize certain known results. (x) The examination of a group action of P on Nn can be refined if we consider the P-action on a single nilpotent GLn-orbit O. This setup generalizes to arbitrary reductive groups G, where the classification of the orbits is equally interesting. Our considerations are, how- ever, restricted to the group GLn. 8 Introduction Recent studies A recent development in this field is A. Melnikov’s study of the B-action on the variety of upper-triangular 2-nilpotent matrices via conjugation [Melnikov, 2000, Melnikov, 2006, Melnikov, 2007] motivated by Springer Theory. The detailed description of the orbits and their closures is given in terms of so-called link patterns; these are combinatorial objects visualizing the set of involutions in the symmetric group S n. The Ph.D. thesis of B. D. Rothbach yields a description of the B-orbits and equations for their closures in the variety of 2-nilpotent matrices. There seems to be a desingularization of the B-orbit closures as well. However, this work was not available to the author of the present thesis, except for parts of it (see [Rothbach, 2009]). N(3) In her Bachelor thesis [Halbach, 2009], B. Halbach describes the B-action on 3 in all detail. The B-orbits as well as their closures, their minimal degenerations and their singularities are explicity given and a generic normal form is introduced. She general- izes the latter, obtaining a large set of pairwise non-B-conjugate matrices. These yield a (n) generic normal form in the nilpotent cone N B Nn for arbitrary n which is proven in [Boos and Reineke, 2011]. Another recent outcome is L. Hille’s and G. Röhrle’s study of an arbitrary parabolic ac- tion on the variety n of upper-triangular nilpotent matrices [Hille and Röhrle, 1999]. They obtain a criterion which determines whether the group action admits a finite or infinite number of orbits. Content of this work We translate the group action of P to a certain group action in the representation theory of finite dimensional algebras via an associated fibre product. In more detail, we “reduce” the classification of orbits of the action to the knowledge of certain isomorphism classes of representations of a quiver Q given by a linearly oriented quiver of Dynkin type An with a loop at the sink and an admissible ideal I given by just one relation. The quiver Q and the ideal I depend on the parabolic subgroup P and on the nilpotency degree x. Although in general, the classification of these isomorphism classes is far from well- known, there are several cases in which the algebra KQ=I is representation-finite and the representations can be classified using the decomposition theorem of W. Krull, R. Remak and O. Schmidt. Our first aim is to classify those cases in which a finite group action arises. The most obvi- ous case is an arbitrary parabolic action on the variety of 2-nilpotent matrices; the covering quiver of Q is of Dynkin type A2n, then. We classify the orbits in terms of “(enhanced) oriented link patterns”, a natural generaliza- tion of A. Melnikov’s link patterns, and provide the concrete structure of the orbit closures: By translating a description of orbit closures given by G. Zwara [Zwara, 1999], a natural generalization of the description A. Melnikov obtains in [Melnikov, 2006] is deduced. 9 By calculating the dimension of the stabilizer of a 2-nilpotent matrix (with respect to the chosen parabolic subgroup), we are able to compute the dimension of every orbit (closure); the description of the open orbit follows naturally. In case of the B-action, we give an explicit description of the minimal degenerations. Therefore, we describe all minimal, disjoint degenerations by using several results of K. Bongartz [Bongartz, 1994] and generalize them to arbitrary minimal degenerations after- wards. In this way, the precise degeneration diagram is obtained; this yields a concrete algorithm to derive a set-theoretic description of the orbit closure of a given matrix by turning around the arrows in the corresponding oriented link pattern.
Recommended publications
  • MATH 2370, Practice Problems
    MATH 2370, Practice Problems Kiumars Kaveh Problem: Prove that an n × n complex matrix A is diagonalizable if and only if there is a basis consisting of eigenvectors of A. Problem: Let A : V ! W be a one-to-one linear map between two finite dimensional vector spaces V and W . Show that the dual map A0 : W 0 ! V 0 is surjective. Problem: Determine if the curve 2 2 2 f(x; y) 2 R j x + y + xy = 10g is an ellipse or hyperbola or union of two lines. Problem: Show that if a nilpotent matrix is diagonalizable then it is the zero matrix. Problem: Let P be a permutation matrix. Show that P is diagonalizable. Show that if λ is an eigenvalue of P then for some integer m > 0 we have λm = 1 (i.e. λ is an m-th root of unity). Hint: Note that P m = I for some integer m > 0. Problem: Show that if λ is an eigenvector of an orthogonal matrix A then jλj = 1. n Problem: Take a vector v 2 R and let H be the hyperplane orthogonal n n to v. Let R : R ! R be the reflection with respect to a hyperplane H. Prove that R is a diagonalizable linear map. Problem: Prove that if λ1; λ2 are distinct eigenvalues of a complex matrix A then the intersection of the generalized eigenspaces Eλ1 and Eλ2 is zero (this is part of the Spectral Theorem). 1 Problem: Let H = (hij) be a 2 × 2 Hermitian matrix. Use the Min- imax Principle to show that if λ1 ≤ λ2 are the eigenvalues of H then λ1 ≤ h11 ≤ λ2.
    [Show full text]
  • Arxiv:1508.00183V2 [Math.RA] 6 Aug 2015 Atclrb Ed Yasemigroup a by field
    A THEOREM OF KAPLANSKY REVISITED HEYDAR RADJAVI AND BAMDAD R. YAHAGHI Abstract. We present a new and simple proof of a theorem due to Kaplansky which unifies theorems of Kolchin and Levitzki on triangularizability of semigroups of matrices. We also give two different extensions of the theorem. As a consequence, we prove the counterpart of Kolchin’s Theorem for finite groups of unipotent matrices over division rings. We also show that the counterpart of Kolchin’s Theorem over division rings of characteristic zero implies that of Kaplansky’s Theorem over such division rings. 1. Introduction The purpose of this short note is three-fold. We give a simple proof of Kaplansky’s Theorem on the (simultaneous) triangularizability of semi- groups whose members all have singleton spectra. Our proof, although not independent of Kaplansky’s, avoids the deeper group theoretical aspects present in it and in other existing proofs we are aware of. Also, this proof can be adjusted to give an affirmative answer to Kolchin’s Problem for finite groups of unipotent matrices over division rings. In particular, it follows from this proof that the counterpart of Kolchin’s Theorem over division rings of characteristic zero implies that of Ka- plansky’s Theorem over such division rings. We also present extensions of Kaplansky’s Theorem in two different directions. M arXiv:1508.00183v2 [math.RA] 6 Aug 2015 Let us fix some notation. Let ∆ be a division ring and n(∆) the algebra of all n × n matrices over ∆. The division ring ∆ could in particular be a field. By a semigroup S ⊆ Mn(∆), we mean a set of matrices closed under multiplication.
    [Show full text]
  • Chapter 6 Jordan Canonical Form
    Chapter 6 Jordan Canonical Form 175 6.1 Primary Decomposition Theorem For a general matrix, we shall now try to obtain the analogues of the ideas we have in the case of diagonalizable matrices. Diagonalizable Case General Case Minimal Polynomial: Minimal Polynomial: k k Y rj Y rj mA(λ) = (λ − λj) mA(λ) = (λ − λj) j=1 j=1 where rj = 1 for all j where 1 ≤ rj ≤ aj for all j Eigenspaces: Eigenspaces: Wj = Null Space of (A − λjI) Wj = Null Space of (A − λjI) dim(Wj) = gj and gj = aj for all j dim(Wj) = gj and gj ≤ aj for all j Generalized Eigenspaces: rj Vj = Null Space of (A − λjI) dim(Vj) = aj Decomposition: Decomposition: n n F = W1 ⊕ W2 ⊕ Wj ⊕ · · · ⊕ Wk W = W1⊕W2⊕Wj⊕· · ·⊕Wk is a subspace of F Primary Decomposition Theorem: n F = V1 ⊕ V2 ⊕ Vj ⊕ Vk n n x 2 F =) 9 unique x 2 F =) 9 unique xj 2 Wj, 1 ≤ j ≤ k such that vj 2 Vj, 1 ≤ j ≤ k such that x = x1 + x2 + ··· + xj + ··· + xk x = v1 + v2 + ··· + vj + ··· + vk 176 6.2 Analogues of Lagrange Polynomial Diagonalizable Case General Case Define Define k k Y Y ri fj(λ) = (λ − λi) Fj(λ) = (λ − λi) i=1 i=1 i6=j i6=j This is a set of coprime polynomials This is a set of coprime polynomials and and hence there exist polynomials hence there exist polynomials Gj(λ), 1 ≤ gj(λ), 1 ≤ j ≤ k such that j ≤ k such that k k X X fj(λ)gj(λ) = 1 Fj(λ)Gj(λ) = 1 j=1 j=1 The polynomials The polynomials `j(λ) = fj(λ)gj(λ); 1 ≤ j ≤ k Lj(λ) = Fj(λ)Gj(λ); 1 ≤ j ≤ k are the Lagrange Interpolation are the counterparts of the polynomials Lagrange Interpolation polynomials `1(λ) + `2(λ) + ··· + `k(λ) = 1 L1(λ) + L2(λ) + ··· + Lk(λ) = 1 λ1`1(λ)+λ2`2(λ)+···+λk`k(λ) = λ λ1L1(λ) + λ2L2(λ) + ··· + λkLk(λ) = d(λ) (a polynomial but it need not be = λ) For i 6= j the product `i(λ)`j(λ) has For i 6= j the product Li(λ)Lj(λ) has a a factor mA (λ) factor mA (λ) 6.3 Matrix Decomposition In the diagonalizable case, using the Lagrange polynomials we obtained a decomposition of the matrix A.
    [Show full text]
  • An Effective Lie–Kolchin Theorem for Quasi-Unipotent Matrices
    Linear Algebra and its Applications 581 (2019) 304–323 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa An effective Lie–Kolchin Theorem for quasi-unipotent matrices Thomas Koberda a, Feng Luo b,∗, Hongbin Sun b a Department of Mathematics, University of Virginia, Charlottesville, VA 22904-4137, USA b Department of Mathematics, Rutgers University, Hill Center – Busch Campus, 110 Frelinghuysen Road, Piscataway, NJ 08854, USA a r t i c l e i n f oa b s t r a c t Article history: We establish an effective version of the classical Lie–Kolchin Received 27 April 2019 Theorem. Namely, let A, B ∈ GLm(C)be quasi-unipotent Accepted 17 July 2019 matrices such that the Jordan Canonical Form of B consists Available online 23 July 2019 of a single block, and suppose that for all k 0the Submitted by V.V. Sergeichuk matrix ABk is also quasi-unipotent. Then A and B have a A, B < C MSC: common eigenvector. In particular, GLm( )is a primary 20H20, 20F38 solvable subgroup. We give applications of this result to the secondary 20F16, 15A15 representation theory of mapping class groups of orientable surfaces. Keywords: © 2019 Elsevier Inc. All rights reserved. Lie–Kolchin theorem Unipotent matrices Solvable groups Mapping class groups * Corresponding author. E-mail addresses: [email protected] (T. Koberda), fl[email protected] (F. Luo), [email protected] (H. Sun). URLs: http://faculty.virginia.edu/Koberda/ (T. Koberda), http://sites.math.rutgers.edu/~fluo/ (F. Luo), http://sites.math.rutgers.edu/~hs735/ (H.
    [Show full text]
  • Dimensions of Nilpotent Algebras
    DIMENSIONS OF NILPOTENT ALGEBRAS AMAJOR QUALIFYING PROJECT REPORT SUBMITTED TO THE FACULTY OF WORCESTER POLYTECHNIC INSTITUTE IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF BACHELOR OF SCIENCE WRITTEN BY KWABENA A. ADWETEWA-BADU APPROVED BY: MARCH 14, 2021 THIS REPORT REPRESENTS THE WORK OF WPI UNDERGRADUATE STUDENTS SUBMITTED TO THE FACULTY AS EVIDENCE OF COMPLETION OF A DEGREE REQUIREMENT. 1 Abstract In this thesis, we study algebras of nilpotent matrices. In the first chapter, we give complete proofs of some well-known results on the structure of a linear operator. In particular, we prove that every nilpotent linear transformation can be represented by a strictly upper triangular matrix. In Chapter 2, we generalise this result to algebras of nilpotent matrices. We prove the famous Lie-Kolchin theorem, which states that an algebra is nilpotent if and only if it is conjugate to an algebra of strictly upper triangular matrices. In Chapter 3 we introduce a family of algebras constructed from graphs and directed graphs. We characterise the directed graphs which give nilpotent algebras. To be precise, a digraph generates an algebra of nilpotency class d if and only if it is acyclic with no paths of length ≥ d. Finally, in Chapter 4, we give Jacobson’s proof of a Theorem of Schur on the maximal dimension of a subalgebra of Mn(k) of nilpotency class 2. We relate this result to problems in external graph theory, and give lower bounds on the dimension of subalgebras of nilpotency class d in Mn(k) for every pair of integers d and n.
    [Show full text]
  • Jordan Canonical Form of a Nilpotent Matrix Math 422 Schur's Triangularization Theorem Tells Us That Every Matrix a Is Unitari
    Jordan Canonical Form of a Nilpotent Matrix Math 422 Schur’s Triangularization Theorem tells us that every matrix A is unitarily similar to an upper triangular matrix T . However, the only thing certain at this point is that the the diagonal entries of T are the eigenvalues of A. The off-diagonal entries of T seem unpredictable and out of control. Recall that the Core-Nilpotent Decomposition of a singular matrix A of index k produces a block diagonal matrix C 0 0 L ∙ ¸ similar to A in which C is non-singular, rank (C)=rank Ak , and L is nilpotent of index k.Isitpossible to simplify C and L via similarity transformations and obtain triangular matrices whose off-diagonal entries are predictable? The goal of this lecture is to do exactly this¡ ¢ for nilpotent matrices. k 1 k Let L be an n n nilpotent matrix of index k. Then L − =0and L =0. Let’s compute the eigenvalues × k k 16 k 1 k 1 k 2 of L. Suppose x =0satisfies Lx = λx; then 0=L x = L − (Lx)=L − (λx)=λL − x = λL − (Lx)= 2 k 2 6 k λ L − x = = λ x; thus λ =0is the only eigenvalue of L. ··· 1 Now if L is diagonalizable, there is an invertible matrix P and a diagonal matrix D such that P − LP = D. Since the diagonal entries of D are the eigenvalues of L, and λ =0is the only eigenvalue of L,wehave 1 D =0. Solving P − LP =0for L gives L =0. Thus a diagonalizable nilpotent matrix is the zero matrix, or equivalently, a non-zero nilpotent matrix L is not diagonalizable.
    [Show full text]
  • A Polynomial in a of the Diagonalizable and Nilpotent Parts of A
    Western Oregon University Digital Commons@WOU Honors Senior Theses/Projects Student Scholarship 6-1-2017 A Polynomial in A of the Diagonalizable and Nilpotent Parts of A Kathryn Wilson Western Oregon University Follow this and additional works at: https://digitalcommons.wou.edu/honors_theses Recommended Citation Wilson, Kathryn, "A Polynomial in A of the Diagonalizable and Nilpotent Parts of A" (2017). Honors Senior Theses/Projects. 145. https://digitalcommons.wou.edu/honors_theses/145 This Undergraduate Honors Thesis/Project is brought to you for free and open access by the Student Scholarship at Digital Commons@WOU. It has been accepted for inclusion in Honors Senior Theses/Projects by an authorized administrator of Digital Commons@WOU. For more information, please contact [email protected], [email protected], [email protected]. A Polynomial in A of the Diagonalizable and Nilpotent Parts of A By Kathryn Wilson An Honors Thesis Submitted in Partial Fulfillment of the Requirements for Graduation from the Western Oregon University Honors Program Dr. Scott Beaver, Thesis Advisor Dr. Gavin Keulks, Honor Program Director June 2017 Acknowledgements I would first like to thank my advisor, Dr. Scott Beaver. He has been an amazing resources throughout this entire process. Your support was what made this project possible. I would also like to thank Dr. Gavin Keulks for his part in running the Honors Program and allowing me the opportunity to be part of it. Further I would like to thank the WOU Mathematics department for making my education here at Western immensely fulfilling, and instilling in me the quality of thinking that will create my success outside of school.
    [Show full text]
  • 6 Sep 2016 a Note on Products of Nilpotent Matrices
    A note on Products of Nilpotent Matrices Christiaan Hattingh∗† September 7, 2016 Abstract Let F be a field. A matrix A ∈ Mn(F ) is a product of two nilpotent matrices if and only if it is singular, except if A is a nonzero nilpotent matrix of order 2×2. This result was proved independently by Sourour [6] and Laffey [4]. While these results remain true and the general strategies and principles of the proofs correct, there are certain problematic details in the original proofs which are resolved in this article. A detailed and rigorous proof of the result based on Laffey’s original proof [4] is provided, and a problematic detail in the Proposition which is the main device underpinning the proof by Sourour is resolved. 1 Introduction First I will fix some notation. A standard basis vector, which is a column vector with a one in position i and zeros elsewhere, is indicated as ei. A matrix with a one in entry (i, j) and zeros elsewhere is indicated as E(i,j). Scripted letters such as F indicate a field, Mn(F ) indicate the set of all matrices of order n×n over the field F . A block diagonal matrix with diagonal blocks A and B (in that order from left top to bottom right) is indicated as Dg[A, B]. The notation [a,b,...,k] indicates a matrix with the column vectors a,b,...,k as columns. A simple Jordan block matrix of order k, and with eigenvalue λ is indicated as Jk(λ). In this text the convention arXiv:1608.04666v3 [math.RA] 6 Sep 2016 of ones on the subdiagonal is assumed for the Jordan canonical form, unless otherwise indicated.
    [Show full text]
  • Spring 2015 Complex and Linear Algebra
    Department of Mathematics OSU Qualifying Examination Spring 2015 PART II: COMPLEX ANALYSIS and LINEAR ALGEBRA • Do any of the two problems in Part CA, use the correspondingly marked blue book and indicate on the selection sheet with your identification number those problems that you wish graded. Similarly for Part LA. • Your solutions should contain all mathematical details. Please write them up as clearly as possible. • Explicitly state any standard theorems, including hypotheses, that are necessary to justify your reasoning. • On problems with multiple parts, individual parts may be weighted differently in grading. • You have three hours to complete Part II. • When you are done with the examination, place examination blue books and selec- tion sheets back into the envelope in which the test materials came. You will hand in all materials. If you use extra examination books, be sure to place your code number on them. 1 PART CA : COMPLEX ANALYSIS QUALIFYING EXAM 1. Let f : C → C be an entire function such that f(R) ⊆ R and f(iR) ⊆ iR, where iR = {it | t ∈ R}. Show that f(−z) = −f(z) for all z ∈ C. 2. Let D = {z ∈ C | |z| < 1}, and let f : D → D be a holomorphic map. Show that for any distinct z, w ∈ D, f(z) − f(w) z − w ≤ . f(w)f(z) − 1 wz − 1 3. Find the number of roots of z7 − 4z3 − 11 = 0 contained in the annulus 1 < |z| < 2. Exam continues on next page ... 2 PART LA: LINEAR ALGEBRA QUALIFYING EXAM ∗ 1. Let A = (aij) be a nonsingular n × n matrix with entries in C, and let A = (aji) 2 2 1/2 be its complex conjugate transpose (or adjoint).
    [Show full text]
  • High-Order Automatic Differentiation of Unmodified Linear Algebra
    High-Order Automatic Differentiation of Unmodified Linear Algebra Routines via Nilpotent Matrices by Benjamin Z Dunham B.A., Carroll College, 2009 M.S., University of Colorado, 2011 A thesis submitted to the Faculty of the Graduate School of the University of Colorado in partial fulfillment of the requirements for the degree of Doctor of Philosophy Department of Aerospace Engineering Sciences 2017 This thesis entitled: High-Order Automatic Differentiation of Unmodified Linear Algebra Routines via Nilpotent Matrices written by Benjamin Z Dunham has been approved for the Department of Aerospace Engineering Sciences Prof. Kurt K. Maute Prof. Alireza Doostan Date The final copy of this thesis has been examined by the signatories, and we find that both the content and the form meet acceptable presentation standards of scholarly work in the above mentioned discipline. iii Dunham, Benjamin Z (Ph.D., Aerospace Engineering Sciences) High-Order Automatic Differentiation of Unmodified Linear Algebra Routines via Nilpotent Matrices Thesis directed by Prof. Kurt K. Maute This work presents a new automatic differentiation method, Nilpotent Matrix Differentiation (NMD), capable of propagating any order of mixed or univariate derivative through common linear algebra functions – most notably third-party sparse solvers and decomposition routines, in addition to basic matrix arithmetic operations and power series – without changing data-type or modifying code line by line; this allows differentiation across sequences of ar- bitrarily many such functions with minimal implementation effort. NMD works by enlarging the matrices and vectors passed to the routines, replacing each original scalar with a matrix block augmented by derivative data; these blocks are constructed with special sparsity structures, termed “stencils,” each designed to be isomorphic to a particular mul- tidimensional hypercomplex algebra.
    [Show full text]
  • Construction of a System of Linear Differential Equations from a Scalar
    ISSN 0081-5438, Proceedings of the Steklov Institute of Mathematics, 2010, Vol. 271, pp. 322–338. c Pleiades Publishing, Ltd., 2010. Original Russian Text c I.V. V’yugin, R.R. Gontsov, 2010, published in Trudy Matematicheskogo Instituta imeni V.A. Steklova, 2010, Vol. 271, pp. 335–351. Construction of a System of Linear Differential Equations from a Scalar Equation I. V. V’yugin a andR.R.Gontsova Received January 2010 Abstract—As is well known, given a Fuchsian differential equation, one can construct a Fuch- sian system with the same singular points and monodromy. In the present paper, this fact is extended to the case of linear differential equations with irregular singularities. DOI: 10.1134/S008154381004022X 1. INTRODUCTION The classical Riemann–Hilbert problem, i.e., the question of existence of a system dy = B(z)y, y(z) ∈ Cp, (1) dz of p linear differential equations with given Fuchsian singular points a1,...,an ∈ C and monodromy χ: π1(C \{a1,...,an},z0) → GL(p, C), (2) has a negative solution in the general case. Recall that a singularity ai of system (1) is said to be Fuchsian if the matrix differential 1-form B(z) dz has simple poles at this point. The monodromy of a system is a representation of the fundamental group of a punctured Riemann sphere in the space of nonsingular complex matrices of dimension p. Under this representation, a loop γ is mapped to a matrix Gγ such that Y (z)=Y (z)Gγ ,whereY (z) is a fundamental matrix of the system in the neighborhood of the point z0 and Y (z) is its analytic continuation along γ.
    [Show full text]
  • Math 265, Midterm 2
    Math 265, Midterm 2 Name: This exam consists of 7 pages including this front page. Ground Rules 1. No calculator is allowed. 2. Show your work for every problem unless otherwise stated. Score 1 16 2 16 3 16 4 18 5 16 6 18 Total 100 1 Notations: R denotes the set of real number; C denotes the set of complex numbers. Rm×n denotes the set of m × n-matrices with entries in R; Rn = Rn×1 denotes the set of n-column vectors; Similar notations for matrices of complex numbers C; Pn[t] denotes the set of polynomials with coefficients in R and the most degree n, that is, n n−1 Pn[t] = ff(t) = ant + an−1t + ··· + a1t + a0; ai 2 R; 8ig: 1. The following are true/false questions. You don't have to justify your an- swers. Just write down either T or F in the table below. A; B; C; X; b are always matrices here. (2 points each) (a) If a matrix A has r linearly dependent rows then it has also r linearly independent columns. (b) Let A be 4×5-matrix with rank 4 then the linear system AX = b always has infinity many solutions for any b 2 R4. (c) If A is similar to B then A and B share the same eigenvectors. (d) Let v1 and v2 be eigenvectors with eigenvalues 0 and 1 respectively then v1 and v2 must be linearly independent. (e) Let A 2 Cn×n be a complex matrix. Then v 2 Cn is an eigenvector of A if and only if the complex conjugationv ¯ is.
    [Show full text]