Smart Geodesic Tracing in GRworkbench

We have developed a new tool for numerical work in General Relativity: GRworkbench. We discuss how GRworkbench's implementation of a numerically-amenable analogue to Differential Geometry facilitates the development of robust and chart-independent numerical algorithms. We consider, as an example, geodesic tracing on two charts covering the exterior Schwarzschild space-time.     

Intended Treatments of Arithmetic Average in U.S. and Asian School Mathematics Textbooks

This study examined how selected U.S. and Asian mathematics curricula are designed to facilitate students' understanding of the arithmetic average. There is a consistency regarding the learning goals among these curriculum series, but the focuses are different between the Asian series and the U.S. reform series. The Asian series and the U.S. commercial series focus the arithmetic average more on conceptual and procedural understanding of the concept as a computational algorithm than on understanding the concept as a representative of a data set; however, the two U.S. reform series focus the concept more on the latter. Because of the different focuses, the Asian and the U.S. curriculum series treat the concept differently. In the Asian series, the concept is first introduced in the context of “equal‐sharing” or “per‐unit‐quantity,” and the averaging formula is formally introduced at a very early stage. In the U.S. reform series, the concept is discussed as a measure of central tendency, and after students have some intuitive ideas of the statistical aspect of the concept, the averaging algorithm is briefly introduced.     

Iterative methods for solvingAx=b,GMRES/FOM versus QMR/BiCG

We study the convergence of the GMRES/FOM and QMR/BiCG methods for solving nonsymmetric systems of equationsAx=b. We prove, in exact arithmetic, that any type of residual norm convergence obtained using BiCG can also be obtained using FOM but on a different system of equations. We consider practical comparisons of these procedures when they are applied to the same matrices. We use a unitary invariance shared by both methods, to construct test matrices where we can vary the nonnormality of the test matrix by variations in simplified eigenvector matrices. We used these test problems in two sets of numerical experiments. The first set of experiments was designed to study effects of increasing nonnormality on the convergence of GMRES and QMR. The second set of experiments was designed to track effects of the eigenvalue distribution on the convergence of QMR. In these tests the GMRES residual norms decreased significantly more rapidly than the QMR residual norms but without corresponding decreases in the error norms. Furthermore, as the nonnormality ofA was increased, the GMRES residual norms decreased more rapidly. This led to premature termination of the GMRES procedure on highly nonnormal problems. On the nonnormal test problems the QMR residual norms exhibited less sensitivity to changes in the nonnormality. The convergence of either type of procedure, as measured by the error norms, was delayed by the presence of large or small outliers and affected by the type of eigenvalues, real or complex, in the eigenvalue distribution ofA. For GMRES this effect can be seen only in the error norm plots.In honor of the 70th birthday of Ted RivlinThis work was supported by NSF grant GER-9450081.     

Latin Squares without Orthogonal Mates

In 1779 Euler proved that for every even n there exists a latin square of order n that has no orthogonal mate, and in 1944 Mann proved that for every n of the form 4k + 1, k ≥ 1, there exists a latin square of order n that has no orthogonal mate. Except for the two smallest cases, n = 3 and n = 7, it is not known whether a latin square of order n = 4k + 3 with no orthogonal mate exists or not. We complete the determination of all n for which there exists a mate-less latin square of order n by proving that, with the exception of n = 3, for all n = 4k + 3 there exists a latin square of order n with no orthogonal mate. We will also show how the methods used in this paper can be applied more generally by deriving several earlier non-orthogonality results.     

Restructuring ordered binary trees

We consider the problem of restructuring an ordered binary tree T, preserving the in-order sequence of its nodes, so as to reduce its height to some target value h. Such a restructuring necessarily involves the downward displacement of some of the nodes of T. Our results, focusing both on the maximum displacement over all nodes and on the maximum displacement over leaves only, provide (i) an explicit tradeoff between the worst-case displacement and the height restriction (including a family of trees that exhibit the worst-case displacements) and (ii) efficient algorithms to achieve height-restricted restructuring while minimizing the maximum node displacement.     

Block transitive Steiner systems with more than one point orbit

For all ‘reasonable’ finite t, k, and s, we construct a t‐(?₀, k, 1) design and a group of automorphisms which is transitive on blocks and has s orbits on points. In particular, there is a 2‐(?0, 4, 1) design with a block‐transitive group of automorphisms having two point orbits. This answers a question of P. J. Cameron and C. E. Praeger. The construction is presented in a purely combinatorial way, but is a by‐product of a new way of looking at a model‐theoretic construction of E. Hrushovski. © 2004 Wiley Periodicals, Inc.     

A Generally Covariant Wave Equation for Grand Unified Field Theory

A generally covariant wave equation is derived geometrically for grand unified field theory. The equation states most generally that the covariant d'Alembertian acting on the vielbein vanishes for the four fields which are thought to exist in nature: gravitation, electromagnetism, weak field and strong field. The various known field equations are derived from the wave equation when the vielbein is the eigenfunction. When the wave equation is applied to gravitation the wave equation is the eigenequation of wave mechanics corresponding to Einstein's field equation in classical mechanics, the vielbein eigenfunction playing the role of the quantized gravitational field. The three Newton laws, Newton's law of universal gravitation, and the Poisson equation are recovered in the classical and nonrelativistic, weak-field limits of the quantized gravitational field. The single particle wave-equation and Klein-Gordon equations are recovered in the relativistic, weak-field limit of the wave equation when scalar components are considered of the vielbein eigenfunction of the quantized gravitational field. The Schrödinger equation is recovered in the non-relativistec, weak-field limit of the Klein-Gordon equation). The Dirac equation is recovered in this weak-field limit of the quantized gravitational field (the nonrelativistic limit of the relativistic, quantezed gravitational field when the vielbein plays the role of the spinor. The wave and field equations of O(3) electrodynamics are recovered when the vielbein becomes the relativistic dreibein (triad) eigenfunction whose three orthonormal space indices become identified with the three complex circular indices (1), (2), (3), and whose four spacetime indices are the indices of non-Euclidean spacetime (the base manifold). This dreibein is the potential dreibein of the O(3) electromagnetic field (an electromagnetic potential four-vector for each index (1), (2), (3)). The wave equation of the parity violating weak field is recovered when the orthonormal space indices of the relativistic dreibein eigenfunction are identified with the indices of the three massive weak field bosons. The wave equation of the strong field is recovered when the orthonormal space indices of the relativistic vielbein eigenfunction become the eight indices defined by the group generators of the SU (3) group.