Efficient QR Decomposition Using Low Complexity Column-wise Givens Rotation (CGR)

QR decomposition (QRD) is a widely used Numerical Linear Algebra (NLA) kernel with applications ranging from SONAR beamforming to wireless MIMO receivers. In this paper, we propose a novel Givens Rotation (GR) based QRD (GR QRD) where we reduce the computational complexity of GR and exploit higher degree of parallelism. This low complexity Column-wise GR (CGR) can annihilate multiple elements of a column of a matrix simultaneously. The algorithm is first realized on a Two-Dimensional (2 D) systolic array and then implemented on REDEFINE which is a Coarse Grained run-time Reconfigurable Architecture (CGRA). We benchmark the proposed implementation against state-of-the-art implementations to report better throughput, convergence and scalability.

Markov chains, R-trivial monoids and representation theory

We develop a general theory of Markov chains realizable as random walks on R-trivial monoids. It provides explicit and simple formulas for the eigenvalues of the transition matrix, for multiplicities of the eigenvalues via Mobius inversion along a lattice, a condition for diagonalizability of the transition matrix and some techniques for bounding the mixing time. In addition, we discuss several examples, such as Toom-Tsetlin models, an exchange walk for finite Coxeter groups, as well as examples previously studied by the authors, such as nonabelian sandpile models and the promotion Markov chain on posets. Many of these examples can be viewed as random walks on quotients of free tree monoids, a new class of monoids whose combinatorics we develop.

Forgotten records of Chrysopelea taprobanica Smith, 1943 (Squamata: Colubridae) from India

The colubrid snake Chrysopelea taprobanica Smith, 1943 was described from a holotype from Kanthali (= Kantalai) and paratypes from Kurunegala, both localities in Sri Lanka (formerly Ceylon) (Smith 1943). Since its description, literature pertaining to Sri Lankan snake fauna considered this taxon to be endemic to the island (Taylor 1950, Deraniyagala 1955, de Silva 1980, de Silva 1990, Somaweera 2004, Somaweera 2006, de Silva 2009, Pyron et al. 2013). In addition, earlier efforts on the Indian peninsula (e.g. Das 1994, 1997, Das 2003, Whitaker & Captain 2004, Aengals et al. 2012) and global data compilations (e.g. Wallach et al. 2014, Uetz & Hošek 2015) did not identify any record from mainland India until Guptha et al. (2015) recorded a specimen (voucher BLT 076 housed at Bio-Lab of Seshachalam Hills, Tirupathi, India) in the dry deciduous forest of Chamala, Seshachalam Biosphere Reserve in Andhra Pradesh, India in November 2013. Guptha et al. (2015) further mentioned an individual previously photographed in 2000 at Rishi Valley, Andhra Pradesh, but with no voucher specimen collected. Guptha’s record, assumed to be the first confirmed record of C. taprobanica in India, is noteworthy as it results in a large range extension, from northern Sri Lanka to eastern India with an Euclidean distance of over 400 km, as well as a change of status, i.e., species not endemic to Sri Lanka. However, at least three little-known previous records of this species from India evaded most literature and were overlooked by the researchers including ourselves.

Spontaneous Lorentz violation: the case of infrared QED

It is by now clear that the infrared sector of quantum electrodynamics (QED) has an intriguingly complex structure. Based on earlier pioneering work on this subject, two of us recently proposed a simple modification of QED by constructing a generalization of the U(1) charge group of QED to the ``Sky'' group incorporating the well-known spontaneous Lorentz violation due to infrared photons, but still compatible in particular with locality (Balachandran and Vaidya, Eur Phys J Plus 128:118, 2013). It was shown that the ``Sky'' group is generated by the algebra of angle-dependent charges and a study of its superselection sectors has revealed a manifest description of spontaneous breaking of the Lorentz symmetry. We further elaborate this approach here and investigate in some detail the properties of charged particles dressed by the infrared photons. We find that Lorentz violation due to soft photons may be manifestly codified in an angle-dependent fermion mass, modifying therefore the fermion dispersion relations. The fact that the masses of the charged particles are not Lorentz invariant affects their spin content, and time dilation formulas for decays should also get corrections.

Non-Abelian gauge redundancy and entropic ambiguities

The von Neumann entropy of a generic quantum state is not unique unless the state can be uniquely decomposed as a sum of extremal or pure states. Therefore one reaches the remarkable possibility that there may be many entropies for a given state. We show that this happens if the GNS representation (of the algebra of observables in some quantum state) is reducible, and some representations in the decomposition occur with non-trivial degeneracy. This ambiguity in entropy, which can occur at zero temperature, can often be traced to a gauge symmetry emergent from the non-trivial topological character of the configuration space of the underlying system. We also establish the analogue of an H-theorem for this entropy by showing that its evolution is Markovian, determined by a stochastic matrix. After demonstrating this entropy ambiguity for the simple example of the algebra of 2 x 2 matrices, we argue that the degeneracies in the GNS representation can be interpreted as an emergent broken gauge symmetry, and play an important role in the analysis of emergent entropy due to non-Abelian anomalies. We work out the simplest situation with such non-Abelian symmetry, that of an ethylene molecule.

L-infinity norms of holomorphic modular forms in the case of compact quotient

We prove a sub-convex estimate for the sup-norm of L-2-normalized holomorphic modular forms of weight k on the upper half plane, with respect to the unit group of a quaternion division algebra over Q. More precisely we show that when the L-2 norm of an eigenfunction f is one, parallel to f parallel to(infinity) <<(epsilon) k(1/2-1/33+epsilon) for any epsilon > 0 and for all k sufficiently large.

Minimization Problems Based on Relative alpha-Entropy I: Forward Projection

Minimization problems with respect to a one-parameter family of generalized relative entropies are studied. These relative entropies, which we term relative alpha-entropies (denoted I-alpha), arise as redundancies under mismatched compression when cumulants of compressed lengths are considered instead of expected compressed lengths. These parametric relative entropies are a generalization of the usual relative entropy (Kullback-Leibler divergence). Just like relative entropy, these relative alpha-entropies behave like squared Euclidean distance and satisfy the Pythagorean property. Minimizers of these relative alpha-entropies on closed and convex sets are shown to exist. Such minimizations generalize the maximum Renyi or Tsallis entropy principle. The minimizing probability distribution (termed forward I-alpha-projection) for a linear family is shown to obey a power-law. Other results in connection with statistical inference, namely subspace transitivity and iterated projections, are also established. In a companion paper, a related minimization problem of interest in robust statistics that leads to a reverse I-alpha-projection is studied.

On Locally Gabriel Geometric Graphs

Let be a set of points in the plane. A geometric graph on is said to be locally Gabriel if for every edge in , the Euclidean disk with the segment joining and as diameter does not contain any points of that are neighbors of or in . A locally Gabriel graph(LGG) is a generalization of Gabriel graph and is motivated by applications in wireless networks. Unlike a Gabriel graph, there is no unique LGG on a given point set since no edge in a LGG is necessarily included or excluded. Thus the edge set of the graph can be customized to optimize certain network parameters depending on the application. The unit distance graph(UDG), introduced by Erdos, is also a LGG. In this paper, we show the following combinatorial bounds on edge complexity and independent sets of LGG: (i) For any , there exists LGG with edges. This improves upon the previous best bound of . (ii) For various subclasses of convex point sets, we show tight linear bounds on the maximum edge complexity of LGG. (iii) For any LGG on any point set, there exists an independent set of size .

Co-Exploration of NLA Kernels and Specification of Compute Elements in Distributed Memory CGRAs

Coarse Grained Reconfigurable Architectures (CGRA) are emerging as embedded application processing units in computing platforms for Exascale computing. Such CGRAs are distributed memory multi- core compute elements on a chip that communicate over a Network-on-chip (NoC). Numerical Linear Algebra (NLA) kernels are key to several high performance computing applications. In this paper we propose a systematic methodology to obtain the specification of Compute Elements (CE) for such CGRAs. We analyze block Matrix Multiplication and block LU Decomposition algorithms in the context of a CGRA, and obtain theoretical bounds on communication requirements, and memory sizes for a CE. Support for high performance custom computations common to NLA kernels are met through custom function units (CFUs) in the CEs. We present results to justify the merits of such CFUs.

Anomalous transport at weak coupling

We evaluate the contribution of chiral fermions in d = 2, 4, 6, chiral bosons, a chiral gravitino like theory in d = 2 and chiral gravitinos in d = 6 to all the leading parity odd transport coefficients at one loop. This is done by using finite temperature field theory to evaluate the relevant Kubo formulae. For chiral fermions and chiral bosons the relation between the parity odd transport coefficient and the microscopic anomalies including gravitational anomalies agree with that found by using the general methods of hydrodynamics and the argument involving the consistency of the Euclidean vacuum. For the gravitino like theory in d = 2 and chiral gravitinos in d = 6, we show that relation between the pure gravitational anomaly and parity odd transport breaks down. From the perturbative calculation we clearly identify the terms that contribute to the anomaly polynomial, but not to the transport coefficient for gravitinos. We also develop a simple method for evaluating the angular integrals in the one loop diagrams involved in the Kubo formulae. Finally we show that charge diffusion mode of an ideal 2 dimensional Weyl gas in the presence of a finite chemical potential acquires a speed, which is equal to half the speed of light.

Action of Hopf on the twisted modular Hecke operator

Let Gamma subset of SL2(Z) be a principal congruence subgroup. For each sigma is an element of SL2(Z), we introduce the collection A(sigma)(Gamma) of modular Hecke operators twisted by sigma. Then, A(sigma)(Gamma) is a right A(Gamma)-module, where A(Gamma) is the modular Hecke algebra introduced by Connes and Moscovici. Using the action of a Hopf algebra h(0) on A(sigma)(Gamma), we define reduced Rankin-Cohen brackets on A(sigma)(Gamma). Moreover A(sigma)(Gamma) carries an action of H 1, where H 1 is the Hopf algebra of foliations of codimension 1. Finally, we consider operators between the levels A(sigma)(Gamma), sigma is an element of SL2(Z). We show that the action of these operators can be expressed in terms of a Hopf algebra h(Z).

Quantum stochastic calculus associated with quadratic quantum noises

We first study a class of fundamental quantum stochastic processes induced by the generators of a six dimensional non-solvable Lie dagger-algebra consisting of all linear combinations of the generalized Gross Laplacian and its adjoint, annihilation operator, creation operator, conservation, and time, and then we study the quantum stochastic integrals associated with the class of fundamental quantum stochastic processes, and the quantum Ito formula is revisited. The existence and uniqueness of solution of a quantum stochastic differential equation is proved. The unitarity conditions of solutions of quantum stochastic differential equations associated with the fundamental processes are examined. The quantum stochastic calculus extends the Hudson-Parthasarathy quantum stochastic calculus. (C) 2016 AIP Publishing LLC.

Counting zero kernel pairs over a finite field

Helmke et al. have recently given a formula for the number of reachable pairs of matrices over a finite field. We give a new and elementary proof of the same formula by solving the equivalent problem of determining the number of so called zero kernel pairs over a finite field. We show that the problem is, equivalent to certain other enumeration problems and outline a connection with some recent results of Guo and Yang on the natural density of rectangular unimodular matrices over F-qx]. We also propose a new conjecture on the density of unimodular matrix polynomials. (C) 2016 Elsevier Inc. All rights reserved.

Análisis espacial y multitemporal de la cobertura y uso del suelo con base en imágenes de satélites en la subcuenca Río Dipilto, Nueva Segovia, Nicaragua (1993, 2000, 2011)

La cuantificación del cambio de uso del suelo presenta aún altos niveles de incertidumbre, lo que repercute por ejemplo en la estimación de las emisiones de CO 2 . En este estudio se desarrollaron métodos, basados en imágenes de satélite y trabajo de campo, para estimar la tasa de cambio de la cobertura y uso del suelo, y las emisiones de CO 2 en la subcuenca río Dipilto, Nueva Segovia. La superficie de los tipos de vegetación se determinó con imágenes Landsat. Se utilizaron datos de carbono de nueve parcelas de muestreo en bosque de pino que fueron correlacionadas, para establecer un modelo de regresión lineal con el objetivo de estimar el Stock de Carbono. La sobreposición y algebra de mapas se utilizó para el escenario de emisiones de CO 2 . El análisis con imágenes de los años 1993, 2000 y 2011 reveló que durante es tos 18 años la velocidad a la que se perdieron los bosques latifoliados cerrado fue variable. Durante los primeros 7 años (1993 a 2000) se registró un aumento de 99.95 ha , que corresponde a una tasa de deforestación de - 1.45 % anual. Durante los últimos once años (2000 a 2011) esta cantidad cambió totalmente, ya que se eliminaron 331.76 h a , que corresponde a una tasa de deforestación anual de 3.41 %. Finalmente considerando el periodo de análisis, se transformaron más de 232.01 h a por año, correspondiente a u na tasa de deforestación anual de 1.55 %. La imagen de 2011 demostró que las reservas o Stock de C oscila entre 40 - 150 t/ha. Este intervalo de valores fue estimado por un modelo de regresión con razonable ajuste (R2 = 0.73 ), cuyas variables independientes fueron la reflectancia de las distintas bandas como índices de vegetación e infrarrojo cercano. Las pérdidas de C se estimaron en intervalos 1 - 191 t/h a en 20.76% del área. El 32.85% del área se mantuvo estable y 46.39% ganancias de 1 - 210 t/ha. La combinación de imágenes de resolución espacial media como son las de la serie Landsat para definir trayectorias de cambio de la cobertura del suelo, es una opción viable para la solución de interrogantes relacionadas con el cambio climático, tales como la estimación de las emisiones de CO 2 derivadas del cambio de uso del suelo.