Hamburger Beiträge zur Angewandten Mathematik Fast Givens Transformation for Quaternion Valued Matrices Applied to Hessenberg Reductions
نویسندگان
چکیده
In a previous paper we investigated Givens transformations applied to quaternion valued matrices. Since arithmetic operations with quaternions are very costly it is desirable to reduce the number of arithmetic operations with quaternions. We show that the Fast Givens transformation, known for the real case, can also be defined for quaternion valued matrices, and we apply this technique to the reduction of an arbitrary quaternion valued matrix to upper Hessenberg form and also include a numerical example. We offer two algorithms. One is based on the classical real case using dynamically two transformation matrices, while the other is based on four transformation matrices where in each step that matrix is selected that has the smallest condition number. For the first algorithm we show that the essential information (namely the two numbers s and c which define the Givens transformation) can be stored in only one variable. This is apparently even new for the real case. We include, necessarily, some investigations on the determination of the relevant condition numbers. We show that in general the application of the Fast Givens transformation in the quaternion case is not as favorable as in the real case with respect to (relative) savings in arithmetic operations. We begin with some introduction into the field of quaternions. In the end in an appendix we present some results concerning the computation of roots of quaternions which in some cases are needed. AMS subject classifications. 12E15, 65F25, 65F30
منابع مشابه
Electronic Transactions on Numerical Analysis
In a previous paper we investigated Givens transformations applied to quaternion valued matrices. Since arithmetic operations with quaternions are very costly it is desirable to reduce the number of arithmetic operations with quaternions. We show that the Fast Givens transformation, known for the real case, can also be defined for quaternion valued matrices, and we apply this technique to the r...
متن کاملHamburger Beiträge zur Angewandten Mathematik
The well known conjugate gradient algorithm (cg-algorithm), introduced by Hestenes & Stiefel, [1952] intended for real, symmetric, positive definite matrices works as well for complex matrices and has the same typical convergence behavior. It will also work, not generally, but in many cases for hermitean, but not necessarily positive definite matrices. We shall show, that the same behavior is s...
متن کاملHamburger Beiträge zur Angewandten Mathematik
We consider a simplified accoustic model to describe nonlinear phenomena occuring in loudspeakers. The simplification is that we restrict to the one–dimensional isentropic Euler equations in a slab, where on the right end a membrane is moving periodically with frequency ω and maximal displacement ε 1. The asymptotic model based on the small parameter ε yields hyperbolic first order systems, w...
متن کاملHamburger Beiträge zur Angewandten Mathematik A Nonlinear Model Predictive Concept for Control of Two-Phase Flows Governed by the Cahn-Hilliard Navier-Stokes System
We present a nonlinear model predictive framework for closedloop control of two-phase flows governed by the Cahn-Hilliard NavierStokes system. We adapt the concept for instantaneous control from [6, 12, 16] to construct distributed closed-loop control strategies for twophase flows. It is well known that distributed instantaneous control is able to stabilize the Burger’s equation [16] and also t...
متن کاملHamburger Beiträge zur Angewandten Mathematik A Priori Error Estimates for a Finite Element Discretization of Parabolic Optimization Problems with Pointwise Constraints in Time on Mean Values of the Gradient of the State
This article is concerned with the discretization of parabolic optimization problems subject to pointwise in time constraints on mean values of the derivative of the state variable. Central component of the analysis are a priori error estimates for the dG(0)-cG(1) discretization of the parabolic partial differential equation (PDE) in the L∞(0, T ;H1 0 (Ω))-norm, together with corresponding esti...
متن کامل