MIMS EPrints
Not a member yet
    2151 research outputs found

    Matrix Inverse Trigonometric and Inverse Hyperbolic Functions: Theory and Algorithms

    Get PDF
    Theoretical and computational aspects of matrix inverse trigonometric and inverse hyperbolic functions are studied. Conditions for existence are given, all possible values are characterized, and principal values \acos, \asin, \acosh, and \asinh are defined and shown to be unique primary matrix functions. Various functional identities are derived, some of which are new even in the scalar case, with care taken to specify precisely the choices of signs and branches. New results include a ``round trip'' formula that relates acos(cosA)\mathrm{acos}(\cos A) to AA and similar formulas for the other inverse functions. Key tools used in the derivations are the matrix unwinding function and the matrix sign function. A new inverse scaling and squaring type algorithm employing a Schur decomposition and variable-degree Pad\'e approximation is derived for computing acos, and it is shown how it can also be used to compute asin, acosh, and asinh. In numerical experiments the algorithm is found to behave in a forward stable fashion and to be superior to computing these functions via logarithmic formulas

    Performance Analysis of Asynchronous Parallel Jacobi

    No full text

    A Catalogue of Software for Matrix Functions. Version 2.0

    Get PDF
    A catalogue of software for computing matrix functions and their Fr\'echet derivatives is presented. For a wide variety of languages and for software ranging from commercial products to open source packages we describe what matrix function codes are available and which algorithms they implement

    Engel-type subgroups and length parameters of finite groups

    Get PDF
    Let gg be an element of a finite group GG. For a positive integer nn, let En(g)E_n(g) be the subgroup generated by all commutators [...[[x,g],g],,g][...[[x,g],g],\dots ,g] over xGx\in G, where gg is repeated nn times. By Baer's theorem, if En(g)=1E_n(g)=1, then gg belongs to the Fitting subgroup F(G)F(G). We generalize this theorem in terms of certain length parameters of En(g)E_n(g). For soluble GG we prove that if, for some nn, the Fitting height of En(g)E_n(g) is equal to kk, then gg belongs to the (k+1)(k+1)th Fitting subgroup Fk+1(G)F_{k+1}(G). For nonsoluble GG the results are in terms of nonsoluble length and generalized Fitting height. The generalized Fitting height h(H)h^*(H) of a finite group HH is the least number hh such that Fh(H)=HF^*_h(H)=H, where F0(H)=1F^*_0(H)=1, and Fi+1(H)F^*_{i+1}(H) is the inverse image of the generalized Fitting subgroup F(H/Fi(H))F^*(H/F^*_{i}(H)). Let mm be the number of prime factors of g|g| counting multiplicities. It is proved that if, for some nn, the generalized Fitting height of En(g)E_n(g) is equal to kk, then gg belongs to Ff(k,m)(G)F^*_{f(k,m)}(G), where f(k,m)f(k,m) depends only on kk and mm. The nonsoluble length~λ(H)\lambda (H) of a finite group~HH is defined as the minimum number of nonsoluble factors in a normal series each of whose factors either is soluble or is a direct product of nonabelian simple groups. It is proved that if λ(En(g))=k\lambda (E_n(g))=k, then gg belongs to a normal subgroup whose nonsoluble length is bounded in terms of kk and mm. We also state conjectures of stronger results independent of mm and show that these conjectures reduce to a certain question about automorphisms of direct products of finite simple groups

    Less is More II: an optimistic view of piecewise smooth bifurcation theory

    Get PDF
    The analysis of piecewise smooth bifurcations reveals an alarming proliferation of cases as the dimension of phase space increases. Rather than attempt the derivation of exhaustive lists of possibilities we describe ways of giving less detailed, but possibly more useful, results

    An inequality concerning the expected values of row-sum and column-sum products in Boolean matrices

    Get PDF
    We give a proof of a matrix inequality and indicate how it can be applied to establish a natural analogy principle in Pure Inductive Logic

    Odd Order Products of Conjugate Involutions in Linear Groups over GF(2^a)

    Get PDF
    Let GG be isomorphic to GLn(q)GL_n(q), SLn(q)SL_n(q), PGLn(q)PGL_n(q) or PSLn(q)PSL_n(q), where q=2aq=2^a. If tt is an involution lying in a GG-conjugacy class XX, then for arbitrary nn we show that as qq becomes large, the proportion of elements of XX which have odd-order product with tt tends to 11. Furthermore, for nn at most 44 we give formulae for the number of elements in XX which have odd-order product with tt, in terms of qq

    A Catalogue of Software for Matrix Functions. Version 2.0

    Get PDF
    A catalogue of software for computing matrix functions and their Fr\'echet derivatives is presented. For a wide variety of languages and for software ranging from commercial products to open source packages we describe what matrix function codes are available and which algorithms they implement

    A Catalogue of Software for Matrix Functions. Version 2.0

    Get PDF
    A catalogue of software for computing matrix functions and their Fr\'echet derivatives is presented. For a wide variety of languages and for software ranging from commercial products to open source packages we describe what matrix function codes are available and which algorithms they implement

    Computing the weighted geometric mean of two large-scale matrices and its inverse times a vector

    Get PDF
    We investigate different approaches for the computation of the action of the weighted geometric mean of two large-scale positive definite matrices on a vector. We derive several algorithms, based on numerical quadrature and the Krylov subspace, and compare them in terms of convergence speed and execution time. By exploiting an algebraic relation between the weighted geometric mean and its inverse, we show how these methods can be used for the solution of large linear system whose coefficient matrix is a weighted geometric mean. We derive two novel algorithms, based on Gauss�Jacobi quadrature, and tailor an existing technique based on contour integration. On the other hand, we adapt several existing Krylov subspace techniques to the computation of the weighted geometric mean. According to our experiments, both classes of algorithms perform well on some problems but there is no clear winner, while some problem-dependent recommendations are provided

    1,445

    full texts

    2,151

    metadata records
    Updated in last 30 days.
    MIMS EPrints
    Access Repository Dashboard
    Do you manage Open Research Online? Become a CORE Member to access insider analytics, issue reports and manage access to outputs from your repository in the CORE Repository Dashboard! 👇