BCYCLIC: A parallel block tridiagonal matrix cyclic solver

13 pages, 6 figures.A block tridiagonal matrix is factored with minimal fill-in using a cyclic reduction algorithm that is easily parallelized. Storage of the factored blocks allows the application of the inverse to multiple right-hand sides which may not be known at factorization time. Scalability with the number of block rows is achieved with cyclic reduction, while scalability with the block size is achieved using multithreaded routines (OpenMP, GotoBLAS) for block matrix manipulation. This dual scalability is a noteworthy feature of this new solver, as well as its ability to efficiently handle arbitrary (non-powers-of-2) block row and processor numbers. Comparison with a state-of-the art parallel sparse solver is presented. It is expected that this new solver will allow many physical applications to optimally use the parallel resources on current supercomputers. Example usage of the solver in magneto-hydrodynamic (MHD), three-dimensional equilibrium solvers for high-temperature fusion plasmas is cited.This research has been sponsored by the US Department of Energy under Contract DE-AC05-00OR22725 with UT-Battelle, LLC. This research used resources of the National Center for Computational Sciences at Oak Ridge National Laboratory, which is supported by the Office of Science of the Department of Energy under Contract DE-AC05-00OR22725.Publicad

Local Production and Developing Core Regions: Ceramic Characterization in the Lake Pátzcuaro Basin, Western Mexico

A core region is the first place for expected shifts in archaeological materials before, during, and after political changes like state emergence and imperial consolidation. Yet, studies of ceramic production have shown that there are sometimes limited or more subtle changes in the ceramic economy throughout such political fluctuations. This article synthesizes recent efforts to address political economic changes via geochemical characterization (neutron activation analysis; NAA) in the Lake Pátzcuaro Basin in western Mexico. This region was home to the Purépecha state and then empire (Tarascan; ca. AD 1350-1530), one of the most powerful kingdoms in the Americas before European arrival. The combined ceramic dataset from four sites in the region result in eight geochemical groups. Our analysis indicates that the region experienced long-term and relatively stable ceramic production that was not substantially altered by the emergence of the state and empire. In addition, we find evidence for (1) dispersed, localized production; (2) long-lived compositional ceramic recipes; and (3) a complex ceramic economy with differential community participation. We discuss why documenting local ceramic production and craft production more generally is important for the study of past political economies

Plasmas and Controlled Nuclear Fusion

Contains reports on three research projects.U. S. Atomic Energy Commission (Contract AT(11-1)-3070

Strong "quantum" chaos in the global ballooning mode spectrum of three-dimensional plasmas

The spectrum of ideal magnetohydrodynamic (MHD) pressure-driven (ballooning) modes in strongly nonaxisymmetric toroidal systems is difficult to analyze numerically owing to the singular nature of ideal MHD caused by lack of an inherent scale length. In this paper, ideal MHD is regularized by using a $k$-space cutoff, making the ray tracing for the WKB ballooning formalism a chaotic Hamiltonian billiard problem. The minimum width of the toroidal Fourier spectrum needed for resolving toroidally localized ballooning modes with a global eigenvalue code is estimated from the Weyl formula. This phase-space-volume estimation method is applied to two stellarator cases.Comment: 4 pages typeset, including 2 figures. Paper accepted for publication in Phys. Rev. Letter

Perpendicular momentum injection by lower hybrid wave in a tokamak

The injection of lower hybrid waves for current drive into a tokamak affects the profile of intrinsic rotation. In this article, the momentum deposition by the lower hybrid wave on the electrons is studied. Due to the increase in the poloidal momentum of the wave as it propagates into the tokamak, the parallel momentum of the wave increases considerably. The change of the perpendicular momentum of the wave is such that the toroidal angular momentum of the wave is conserved. If the perpendicular momentum transfer via electron Landau damping is ignored, the transfer of the toroidal angular momentum to the plasma will be larger than the injected toroidal angular momentum. A proper quasilinear treatment proves that both perpendicular and parallel momentum are transferred to the electrons. The toroidal angular momentum of the electrons is then transferred to the ions via different mechanisms for the parallel and perpendicular momentum. The perpendicular momentum is transferred to ions through an outward radial electron pinch, while the parallel momentum is transferred through collisions.Comment: 22 pages, 4 figure

Plasmas and Controlled Nuclear Fusion

Contains research objectives, summary of research and reports on six research projects.U. S. Atomic Energy Commission (Contract AT(11-1)-3070

Designs Studies of Low-Aspect Ratio Quasi-Omnigenous Stellarators

Significant progress has been made in the development of new modest-size compact stellarator devices that could test optimization principles for the design of a more attractive reactor. These are 3 and 4 field period low-aspect-ratio quasi-omnigenous (QO) stellarators based on an optimization method that targets improved confinement, stability, ease of coil design, low-aspect-ratio, and low bootstrap current

Plasma Dynamics

Contains research objectives and summary of research on eighteen research projects split into seven sections and reports on four research projects.U.S. Atomic Energy Commission (Contract AT(l1-1)-3070)National Science Foundation (Grant GK-37979X1