Specific "scientific" data structures, and their processing

Programming physicists use, as all programmers, arrays, lists, tuples, records, etc., and this requires some change in their thought patterns while converting their formulae into some code, since the "data structures" operated upon, while elaborating some theory and its consequences, are rather: power series and Pad\'e approximants, differential forms and other instances of differential algebras, functionals (for the variational calculus), trajectories (solutions of differential equations), Young diagrams and Feynman graphs, etc. Such data is often used in a [semi-]numerical setting, not necessarily "symbolic", appropriate for the computer algebra packages. Modules adapted to such data may be "just libraries", but often they become specific, embedded sub-languages, typically mapped into object-oriented frameworks, with overloaded mathematical operations. Here we present a functional approach to this philosophy. We show how the usage of Haskell datatypes and - fundamental for our tutorial - the application of lazy evaluation makes it possible to operate upon such data (in particular: the "infinite" sequences) in a natural and comfortable manner.Comment: In Proceedings DSL 2011, arXiv:1109.032

Row-switched states in two-dimensional underdamped Josephson junction arrays

When magnetic flux moves across layered or granular superconductor structures, the passage of vortices can take place along channels which develop finite voltage, while the rest of the material remains in the zero-voltage state. We present analytical studies of an example of such mixed dynamics: the row-switched (RS) states in underdamped two-dimensional Josephson arrays, driven by a uniform DC current under external magnetic field but neglecting self-fields. The governing equations are cast into a compact differential-algebraic system which describes the dynamics of an assembly of Josephson oscillators coupled through the mesh current. We carry out a formal perturbation expansion, and obtain the DC and AC spatial distributions of the junction phases and induced circulating currents. We also estimate the interval of the driving current in which a given RS state is stable. All these analytical predictions compare well with our numerics. We then combine these results to deduce the parameter region (in the damping coefficient versus magnetic field plane) where RS states can exist.Comment: latex, 48 pages, 15 figs using psfi

Reusable Components of Semantic Specifications

Semantic specifications of programming languages typically have poor modularity. This hinders reuse of parts of the semantics of one language when specifying a different language – even when the two languages have many constructs in common – and evolution of a language may require major reformulation of its semantics. Such drawbacks have discouraged language developers from using formal semantics to document their designs. In the PLanCompS project, we have developed a component-based approach to semantics. Here, we explain its modularity aspects, and present an illustrative case study: a component-based semantics for Caml Light. We have tested the correctness of the semantics by running programs on an interpreter generated from the semantics, comparing the output with that produced on the standard implementation of the language. Our approach provides good modularity, facilitates reuse, and should support co-evolution of languages and their formal semantics. It could be particularly useful in connection with domain-specific languages and language-driven software development

Disorder and interference: localization phenomena

The specific problem we address in these lectures is the problem of transport and localization in disordered systems, when interference is present, as characteristic for waves, with a focus on realizations with ultracold atoms.Comment: Notes of a lecture delivered at the Les Houches School of Physics on "Ultracold gases and quantum information" 2009 in Singapore. v3: corrected mistakes, improved script for numerics, Chapter 9 in "Les Houches 2009 - Session XCI: Ultracold Gases and Quantum Information" edited by C. Miniatura et al. (Oxford University Press, 2011

Numerical techniques in linear duct acoustics

Both finite difference and finite element analyses of small amplitude (linear) sound propagation in straight and variable area ducts with flow, as might be found in a typical turboject engine duct, muffler, or industrial ventilation system, are reviewed. Both steady state and transient theories are discussed. Emphasis is placed on the advantages and limitations associated with the various numerical techniques. Examples of practical problems are given for which the numerical techniques have been applied

General-relativistic resistive magnetohydrodynamics in three dimensions: Formulation and tests

We present a new numerical implementation of the general-relativistic resistive magnetohydrodynamics (MHD) equations within the Whisky code. The numerical method adopted exploits the properties of implicit-explicit Runge-Kutta numerical schemes to treat the stiff terms that appear in the equations for large electrical conductivities. Using tests in one, two, and three dimensions, we show that our implementation is robust and recovers the ideal-MHD limit in regimes of very high conductivity. Moreover, the results illustrate that the code is capable of describing scenarios in a very wide range of conductivities. In addition to tests in flat spacetime, we report simulations of magnetized nonrotating relativistic stars, both in the Cowling approximation and in dynamical spacetimes. Finally, because of its astrophysical relevance and because it provides a severe testbed for general-relativistic codes with dynamical electromagnetic fields, we study the collapse of a nonrotating star to a black hole. We show that also in this case our results on the quasinormal mode frequencies of the excited electromagnetic fields in the Schwarzschild background agree with the perturbative studies within 0.7% and 5.6% for the real and the imaginary part of the l=1 mode eigenfrequency, respectively. Finally we provide an estimate of the electromagnetic efficiency of this process.Comment: 22 pages, 19 figure

A variational principle for actions on symmetric symplectic spaces

We present a definition of generating functions of canonical relations, which are real functions on symmetric symplectic spaces, discussing some conditions for the presence of caustics. We show how the actions compose by a neat geometrical formula and are connected to the hamiltonians via a geometrically simple variational principle which determines the classical trajectories, discussing the temporal evolution of such ``extended hamiltonians'' in terms of Hamilton-Jacobi-type equations. Simplest spaces are treated explicitly.Comment: 28 pages. Edited english translation of first author's PhD thesis (2000