Circulant weighing matrices of weight 22t

A weighing matrix of weight k is a square matrix M with entries 0, ± 1 such that MM ^T = kI _n . We study the case that M is a circulant and k = 2^2t for some positive integer t. New structural results are obtained. Based on these results, we make a complete computer search for all circulant weighing matrices of order 16.      

Permutation Matrices, Wreath Products, and the Distribution of Eigenvalues

We consider a class of random matrix ensembles which can be constructed from the random permutation matrices by replacing the nonzero entries of the n×n permutation matrix matrix with M×M diagonal matrices whose entries are random Kth roots of unity or random points on the unit circle. Let X be the number of eigenvalues lying in a specified arc I of the unit circle, and consider the standardized random variable (X–E[X])/(Var(X))^1/2. We show that for a fixed set of arcs I ₁,...,I _N , the corresponding standardized random variables are jointly normal in the large n limit, and compare the covariance structures which arise with results for other random matrix ensembles.     

Locally Optimal Block Preconditioned Conjugate Gradient Method for Hierarchical Matrices

We present a method of almost linear complexity to approximate some (inner) eigenvalues of symmetric self-adjoint integral or differential operators. Using ℋ-arithmetic the discretisation of the operator leads to a large hierarchical (ℋ-) matrix M. We assume that M is symmetric, positive definite. Then we compute the smallest eigenvalues by the locally optimal block preconditioned conjugate gradient method (LOBPCG), which has been extensively investigated by Knyazev and Neymeyr. Hierarchical matrices were introduced by W. Hackbusch in 1998. They are data-sparse and require only O(nlog₂ n) storage. There is an approximative inverse, besides other matrix operations, within the set of ℋ-matrices, which can be computed in linear-polylogarithmic complexity. We will use the approximative inverse as preconditioner in the LOBPCG method. Further we combine the LOBPCG method with the folded spectrum method to compute inner eigenvalues of M. This is equivalent to the application of LOBPCG to the matrix M_μ = (M − μI)² . The matrix M_μ is symmetric, positive definite, too. Numerical experiments illustrate the behavior of the suggested approach. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hadamard Matrices Modulo p and Small Modular Hadamard Matrices

We use modular symmetric designs to study the existence of Hadamard matrices modulo certain primes. We solve the 7‐modular and 11‐modular versions of the Hadamard conjecture for all but a finite number of cases. In doing so, we state a conjectural sufficient condition for the existence of a p‐modular Hadamard matrix for all but finitely many cases. When 2 is a primitive root of a prime p, we conditionally solve this conjecture and therefore the p‐modular version of the Hadamard conjecture for all but finitely many cases when , and prove a weaker result for . Finally, we look at constraints on the existence of m‐modular Hadamard matrices when the size of the matrix is small compared to m.     

Some New Results on Circulant Weighing Matrices

We obtain a few structural theorems for circulant weighing matrices whose weight is the square of a prime number. Our results provide new schemes to search for these objects. We also establish the existence status of several previously open cases of circulant weighing matrices. More specifically we show their nonexistence for the parameter pairs (n, k) (here n is the order of the matrix and k its weight) = (147, 49), (125, 25), (200, 25), (55, 25), (95, 25), (133, 49), (195, 25), (11 w, 121) for w < 62.     

Tournament matrices and their generalizations,I.

If M is any complex matrix with rank (M + M ^* + I) = 1, we show that any eigenvalue of M that is not geometrically simple has 1/2 for its real part. This generalizes a recent finding of de Caen and Hoffman: the rank of any n × n tournament matrix is at least n ? 1. We extend several spectral properties of tournament matrices to this and related types of matrices. For example, we characterize the singular real matrices M with 0 diagonal for which rank (M + M^T + I) = 1 and we characterize the vectors that can be in the kernels of such matrices. We show that singular, irreducible n × n tournament matrices exist if and only n? {2,3,4,5} and exhibit many infinite families of such matrices. Connections with signed digraphs are explored and several open problems are presented.     

Non-skew symmetric orthogonal matrices with constant diagonals

A matrix C of order n is orthogonal if CC^T=dI. In this paper, we restrict the study to orthogonal matrices with a constant m > 1 on the diagonal and ±1's off the diagonal. It is observed that all skew symmetric orthogonal matrices of this type are constructed from skew symmetric Hadamard matrices and vice versa. Some simple necessary conditions for the existence of non-skew orthogonal matrices are derived. Two basic construction techniques for non-skew orthogonal matrices are given. Several families of non-skew orthogonal matrices are constructed by applying the basic techniques to well-known combinatorial objects like balanced incomplete block designs. It is also shown that if m is even and n=0 (mod 4), then an orthogonal matrix must be skew symmetric. The structure of a non-skew orthogonal matrix in the special case of m odd,n=2 (mod 4) and m?1/6n is also studied in detail. Finally, a list of cases with n?50 is given where the existence of non-skew orthogonal matrices are unknown.     

Inverse-preserving Linear Maps Between Spaces of Matrices over Fields

Suppose F is a field different from F2, the field with two elements. Let Mn（F） and Sn（F） be the space of n × n full matrices and the space of n ×n symmetric matrices over F, respectively. For any G1, G2 ∈ {Sn（F）, Mn（F）}, we say that a linear map f from G1 to G2 is inverse-preserving if f（X）^-1 = f（X^-1） for every invertible X ∈ G1. Let L （G1, G2） denote the set of all inverse-preserving linear maps from G1 to G2. In this paper the sets .L（Sn（F）,Mn（F））, L（Sn（F）,Sn（F））, L （Mn（F）,Mn（F）） and L（Mn （F）, Sn （F）） are characterized.     

On Skew E–W Matrices

An E–W matrix M is a ( ? 1, 1)‐matrix of order , where t is a positive integer, satisfying that the absolute value of its determinant attains Ehlich–Wojtas' bound. M is said to be of skew type (or simply skew) if is skew‐symmetric where I is the identity matrix. In this paper, we draw a parallel between skew E–W matrices and skew Hadamard matrices concerning a question about the maximal determinant. As a consequence, a problem posted on Cameron's website [7] has been partially solved. Finally, codes constructed from skew E–W matrices are presented. A necessary and sufficient condition for these codes to be self‐dual is given, and examples are provided for lengths up to 52.     

Transposing a matrix by similarity operations

A matrix T is said to co-transpose a square matrix A if T^?1AT=A′ and T^?1A′T=A. For every n?3 there exists a real n×n matrix which cannot be co-transposed by any matrix. However, it is shown that the following classes of real matrices can be co-transposed by a symmetric matrix of order two: 2×2 matrices, normal matrices, and matrices whose square is symmetric.     

Parallel Dixon Matrices by Bracket

It is known that the Dixon matrix can be constructed in parallel either by entry or by diagonal. This paper presents another parallel matrix construction, this time by bracket. The parallel by bracket algorithm is the fastest among the three, but not surprisingly it requires the highest number of processors. The method also shows analytically that the Dixon matrix has a total of m(m+1)²(m+2)n(n+1)²(n+2)/36 brackets but only mn(m+1)(n+1)(mn+2m+2n+1)/6 of them are distinct.     

Hadamard Matrices and Dihedral Groups

Let D _2p be a dihedral group of order 2p, where p is an odd integer. Let ZD _2p be the group ring of D _2p over the ring Z of integers. We identify elements of ZD _2p and their matrices of the regular representation of ZD _2p . Recently we characterized the Hadamard matrices of order 28 ([6] and [7]). There are exactly 487 Hadamard matrices of order 28, up to equivalence. In these matrices there exist matrices with some interesting properties. That is, these are constructed by elements of ZD ₆. We discuss relation of ZD _2p and Hadamard matrices of order n=8p+4, and give some examples of Hadamard matrices constructed by dihedral groups.     

The DJL Conjecture for CP Matrices over Special Inclines

Drew, Johnson, and Loewy conjectured that for n ≥ 4, the CP-rank of every n × n completely positive real matrix is at most [n²/4]. While this conjecture has recently been disproved for completely positive real matrices, we show that this conjecture is true for n × n completely positive matrices over certain special types of inclines. In addition, we prove an incline version of Markham's theorems which gives sufficient conditions for completely positive matrices over special inclines to have triangular factorizations.     

A Class of Sparse Invertible Matrices and Their Use for Nonlinear Prediction of Nearly Periodic Time Series with Fixed Period

We introduce a class of sparse matrices U _m (A _{p 1} ) of order m by m, where m is a composite natural number, p ₁ is a divisor of m, and A _{p 1} is a set of nonzero real numbers of length p ₁. The construction of U _m (A _{p 1} ) is achieved by iteration, involving repetitive dilation operations and block-matrix operations. We prove that the matrices U _m (A _{p 1} ) are invertible and we compute the inverse matrix (U _m (A _{p 1} ))^?1 explicitly. We prove that each row of the inverse matrix (U _m (A _{p 1} ))^?1 has only two nonzero entries with alternative signs, located at specific positions, related to the divisors of m. We use the structural properties of the matrix (U _m (A _{p 1} ))^?1 in order to build a nonlinear estimator for prediction of nearly periodic time series of length m with fixed period.     

On the matrix inequality (I + X)−1 ≤ I

In this note we consider the question under which conditions all entries of the matrix I???(I?+?X)^?1 are nonnegative in case matrix X is a real positive definite matrix. Sufficient conditions are presented as well as some necessary conditions. One sufficient condition is that matrix X ^?1 is an inverse M-matrix. A class of matrices for which the inequality holds is presented.     

Distribution of Subdominant Eigenvalues of Random Matrices

We mainly investigate the behavior of the subdominant eigenvalue of matrices B= (b _i,j)^n,n whose entries are independent random variables with an expectation Eb _i,j=1/n and with a variance n c/n ² for some constant c 0. For such matrices we show that for large n, the subdominant eigenvalue is, with great probability, in a small neighborhood of 0. We also show that for large n, the spectral radius of such matrices is, with great probability, in a small neighborhood of 1.     

A technique for constructing non‐embeddable quasi‐residual designs

We propose a technique for constructing two infinite families of non‐embeddable quasi‐residual designs as soon as one such design satisfying certain conditions exists. The main tools are generalized Hadamard matrices and balanced generalized weighing matrices. Starting with a specific non‐embeddable quasi‐residual 2‐(27,9,4) design, we construct for every positive integer m a non‐embeddable 2‐(3^m,3^m?1,(3^m?1?1)/2)‐design, and, if r_m=(3^m?1)/2 is a prime power, we construct for every positive integer n a non‐embeddable design. For each design in these families, a symmetric design with the corresponding parameters is known to exist. © 2002 Wiley Periodicals, Inc. J Combin Designs 10: 160–172, 2002; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/jcd.900     

Constructions for Circulant and Group‐Developed Generalized Weighing Matrices

An elementary construction yields a new class of circulant (so‐called “Butson‐type”) generalized weighing matrices, which have order and weight n², all of whose entries are nth roots of unity, for all positive integers , where . The idea is extended to a wider class of constructions giving various group‐developed generalized weighing matrices.     

A backward stability analysis of diagonal pivoting methods for solving unsymmetric tridiagonal systems without interchanges

This paper concerns the LBM ^T factorization of unsymmetric tridiagonal matrices, where L and M are unit lower triangular matrices and B is block diagonal with 1×1 and 2×2 blocks. In some applications, it is necessary to form this factorization without row or column interchanges while the tridiagonal matrix is formed. Bunch and Kaufman proposed a pivoting strategy without interchanges specifically for symmetric tridiagonal matrices, and more recently, Bunch and Marcia proposed pivoting strategies that are normwise backward stable for linear systems involving such matrices. In this paper, we extend these strategies to the unsymmetric tridiagonal case and demonstrate that the proposed methods both exhibit bounded growth factors and are normwise backward stable. Copyright © 2009 John Wiley & Sons, Ltd.     

线性流形上的广义反射矩阵反问题

设 R∈C^m×m 及 S∈C^n×n 是非平凡Hermitian酉矩阵, 即 R^H=R=R^-1≠±I_m ,S^H=S=S^-1≠±I_n.若矩阵 A∈C^m×n 满足 RAS=A, 则称矩阵 A 为广义反射矩阵.该文考虑线性流形上的广义反射矩阵反问题及相应的最佳逼近问题.给出了反问题解的一般表示, 得到了线性流形上矩阵方程AX₂=Z₂, Y₂^H A=W₂^H 具有广义反射矩阵解的充分必要条件, 导出了最佳逼近问题唯一解的显式表示.