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

New uncertainty relations for tomographic entropy: Application to squeezed states and solitons

Using the tomographic probability distribution (symplectic tomogram) describing the quantum state (instead of the wave function or density matrix) and properties of recently introduced tomographic entropy associated with the probability distribution, the new uncertainty relation for the tomographic entropy is obtained. Examples of the entropic uncertainty relation for squeezed states and solitons of the Bose--Einstein condensate are considered.Comment: 18 pages, 2 figures, to be published in European Physical Journal

Extension of the SIESTA MHD equilibrium code to free-plasma-boundary problems

is a recently developed MHD equilibrium code designed to perform fast and accurate calculations of ideal MHD equilibria for three-dimensional magnetic configurations. Since SIESTA does not assume closed magnetic surfaces, the solution can exhibit magnetic islands and stochastic regions. In its original implementation SIESTA addressed only fixed-boundary problems. That is, the shape of the plasma edge, assumed to be a magnetic surface, was kept fixed as the solution iteratively converges to equilibrium. This condition somewhat restricts the possible applications of SIESTA. In this paper, we discuss an extension that will enable SIESTA to address free-plasma-boundary problems, opening up the possibility of investigating problems in which the plasma boundary is perturbed either externally or internally. As an illustration, SIESTA is applied to a configuration of the W7-X stellarator.This research was funded in part by the Ministerio de Economía, Industria y Competitividad of Spain, Grant No. ENE2015-68265. This research was carried out in part at the Max-Planck-Institute for Plasma Physics in Greifswald (Germany), whose hospitality is gratefully acknowledged. This research was supported in part by the U.S. Department of Energy, Office of Fusion Energy Sciences under Award DE-AC05-00OR22725. SIESTA runs have been carred out in Uranus, a supercomputer cluster located at Universidad Carlos III de Madrid and funded jointly by the European Regional Development Funds (EU-FEDER) Project No. UNC313-4E-2361, and by the Ministerio de Economía, Industria y Competitividad via the National Project Nos. ENE2009-12213-C03-03, ENE2012-33219, and ENE2012-31753