Quasiseparable Approach to Evaluating Cubic Splines

The development of fast and efficient algorithms is crucial not only for computer scientists, but also for mathematicians and engineers as those algorithms lead to reduce complexity. Another common interest of these professionals is to construct models using existing data. This leads numerical analysts to explore interpolation techniques. One such technique is called cubic spline interpolation. In here, we will propose a cubic spline solver aiming to bridge the gap between numerical linear algebra, electrical engineering, systems engineering, sensor processing, and parallel processing. We will use quasiseparable structure to evaluate cubic splines by deriving a fast and stable algorithm. The derivation is carried through a specific factorization of the inverse of tridiagonal matrices. This factorization leads to an alternative method to solve the system of tridiagonal matrices as opposed to the existing methods. The proposed algorithm has the lowest computational complexity compared to existing algorithms

On post-Lie algebras, Lie--Butcher series and moving frames

Pre-Lie (or Vinberg) algebras arise from flat and torsion-free connections on differential manifolds. They have been studied extensively in recent years, both from algebraic operadic points of view and through numerous applications in numerical analysis, control theory, stochastic differential equations and renormalization. Butcher series are formal power series founded on pre-Lie algebras, used in numerical analysis to study geometric properties of flows on euclidean spaces. Motivated by the analysis of flows on manifolds and homogeneous spaces, we investigate algebras arising from flat connections with constant torsion, leading to the definition of post-Lie algebras, a generalization of pre-Lie algebras. Whereas pre-Lie algebras are intimately associated with euclidean geometry, post-Lie algebras occur naturally in the differential geometry of homogeneous spaces, and are also closely related to Cartan's method of moving frames. Lie--Butcher series combine Butcher series with Lie series and are used to analyze flows on manifolds. In this paper we show that Lie--Butcher series are founded on post-Lie algebras. The functorial relations between post-Lie algebras and their enveloping algebras, called D-algebras, are explored. Furthermore, we develop new formulas for computations in free post-Lie algebras and D-algebras, based on recursions in a magma, and we show that Lie--Butcher series are related to invariants of curves described by moving frames.Comment: added discussion of post-Lie algebroid

Reduction of Algebraic Parametric Systems by Rectification of their Affine Expanded Lie Symmetries

Lie group theory states that knowledge of a $m$-parameters solvable group of symmetries of a system of ordinary differential equations allows to reduce by $m$ the number of equations. We apply this principle by finding some \emph{affine derivations} that induces \emph{expanded} Lie point symmetries of considered system. By rewriting original problem in an invariant coordinates set for these symmetries, we \emph{reduce} the number of involved parameters. We present an algorithm based on this standpoint whose arithmetic complexity is \emph{quasi-polynomial} in input's size.Comment: Before analysing an algebraic system (differential or not), one can generally reduce the number of parameters defining the system behavior by studying the system's Lie symmetrie

Exact ZF Analysis and Computer-Algebra-Aided Evaluation in Rank-1 LoS Rician Fading

We study zero-forcing detection (ZF) for multiple-input/multiple-output (MIMO) spatial multiplexing under transmit-correlated Rician fading for an N_R X N_T channel matrix with rank-1 line-of-sight (LoS) component. By using matrix transformations and multivariate statistics, our exact analysis yields the signal-to-noise ratio moment generating function (m.g.f.) as an infinite series of gamma distribution m.g.f.'s and analogous series for ZF performance measures, e.g., outage probability and ergodic capacity. However, their numerical convergence is inherently problematic with increasing Rician K-factor, N_R , and N_T. We circumvent this limitation as follows. First, we derive differential equations satisfied by the performance measures with a novel automated approach employing a computer-algebra tool which implements Groebner basis computation and creative telescoping. These differential equations are then solved with the holonomic gradient method (HGM) from initial conditions computed with the infinite series. We demonstrate that HGM yields more reliable performance evaluation than by infinite series alone and more expeditious than by simulation, for realistic values of K , and even for N_R and N_T relevant to large MIMO systems. We envision extending the proposed approaches for exact analysis and reliable evaluation to more general Rician fading and other transceiver methods.Comment: Accepted for publication by the IEEE Transactions on Wireless Communications, on April 7th, 2016; this is the final revision before publicatio

Computer Algebra Applications for Numerical Relativity

We discuss the application of computer algebra to problems commonly arising in numerical relativity, such as the derivation of 3+1-splits, manipulation of evolution equations and automatic code generation. Particular emphasis is put on working with abstract index tensor quantities as much as possible.Comment: 5 pages, no figures. To appear in the Proceedings of the Spanish Relativity Meeting (ERE 2002), Mao, Menorca, Spain, 22-24 Sept 200

On the Lie enveloping algebra of a post-Lie algebra

We consider pairs of Lie algebras $g$ and $\bar{g}$, defined over a common vector space, where the Lie brackets of $g$ and $\bar{g}$ are related via a post-Lie algebra structure. The latter can be extended to the Lie enveloping algebra $U(g)$. This permits us to define another associative product on $U(g)$, which gives rise to a Hopf algebra isomorphism between $U(\bar{g})$ and a new Hopf algebra assembled from $U(g)$ with the new product. For the free post-Lie algebra these constructions provide a refined understanding of a fundamental Hopf algebra appearing in the theory of numerical integration methods for differential equations on manifolds. In the pre-Lie setting, the algebraic point of view developed here also provides a concise way to develop Butcher's order theory for Runge--Kutta methods.Comment: 25 page