Optimising code generation with haggies

This article describes haggies, a program for the generation of optimised programs for the efficient numerical evaluation of mathematical expressions. It uses a multivariate Horner-scheme and Common Subexpression Elimination to reduce the overall number of operations.The package can serve as a back-end for virtually any general purpose computer algebra program. Built-in type inference that allows to deal with non-standard data types in strongly typed languages and a very flexible, pattern-based output specification ensure that haggies can produce code for a large variety of programming languages.We currently use haggies as part of an automated package for the calculation of one-loop scattering amplitudes in quantum field theories. The examples in this articles, however, demonstrate that its use is not restricted to the field of high energy physics.

Program summary

Program title: haggiesCatalogue identifier: AEGF_v1_0Program summary: URL: http://cpc.cs.qub.ac.uk/summaries/AEGF_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU GPL v3No. of lines in distributed program, including test data, etc.: 56 220No. of bytes in distributed program, including test data, etc.: 579 010Distribution format: tar.gzProgramming language: Java, JavaCCComputer: Any system that runs the Java Virtual MachineOperating system: Any system that runs the Java Virtual MachineRAM: Determined by the size of the problemClassification: 4.14, 5, 6.2, 6.5, 11.1Nature of problem: Generation of optimised programs for the evaluation of possibly large algebraic expressionsSolution method: Java implementationRunning time: Determined by the size of the problem     

A web-deployed interface for performing ab initio molecular dynamics, optimization, and electronic structure in Fireball

Fireball is an ab initio technique for fast local orbital simulations of nanotechnological, solid state, and biological systems. We have implemented a convenient interface for new users and software architects in the platform-independent Java language to access Fireball's unique and powerful capabilities. The graphical user interface can be run directly from a web server or from within a larger framework such as the Computational Science and Engineering Online (CSE-Online) environment or the Distributed Analysis of Neutron Scattering Experiments (DANSE) framework. We demonstrate its use for high-throughput electronic structure calculations and a multi-100 atom quantum molecular dynamics (MD) simulation.

Program summary

Program title: FireballUICatalogue identifier: AECF_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECF_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 279 784No. of bytes in distributed program, including test data, etc.: 12 836 145Distribution format: tar.gzProgramming language: JavaComputer: PC and workstationOperating system: The GUI will run under Windows, Mac and Linux. Executables for Mac and Linux are included in the package.RAM: 512 MBWord size: 32 or 64 bitsClassification: 4.14Nature of problem: The set up and running of many simulations (all of the same type), from the command line, is a slow process. But most research quality codes, including the ab initio tight-binding code FIREBALL, are designed to run from the command line. The desire is to have a method for quickly and efficiently setting up and running a host of simulations.Solution method: We have created a graphical user interface for use with the FIREBALL code. Once the user has created the files containing the atomic coordinates for each system that they are going to run a simulation on, the user can set up and start the computations of up to hundreds of simulations.Running time: 3 to 5 minutes on a 2 GHz Pentium IV processor.     

A general spectral method for the numerical simulation of one-dimensional interacting fermions

This work introduces a general framework for the direct numerical simulation of systems of interacting fermions in one spatial dimension. The approach is based on a specially adapted nodal spectral Galerkin method, where the basis functions are constructed to obey the antisymmetry relations of fermionic wave functions. An efficient Matlab program for the assembly of the stiffness and potential matrices is presented, which exploits the combinatorial structure of the sparsity pattern arising from this discretization to achieve optimal run-time complexity. This program allows the accurate discretization of systems with multiple fermions subject to arbitrary potentials, e.g., for verifying the accuracy of multi-particle approximations such as Hartree–Fock in the few-particle limit. It can be used for eigenvalue computations or numerical solutions of the time-dependent Schrödinger equation.Program summaryProgram title: assembleFermiMatrixCatalogue identifier: AEKO_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEKO_v1_0.htmlProgram obtainable from: CPC Program Library, Queen?s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 102No. of bytes in distributed program, including test data, etc.: 2294Distribution format: tar.gzProgramming language: MATLABComputer: Any architecture supported by MATLABOperating system: Any supported by MATLAB; tested under Linux (x86-64) and Mac OS X (10.6)RAM: Depends on the dataClassification: 4.3, 2.2Nature of problem: The direct numerical solution of the multi-particle one-dimensional Schrödinger equation in a quantum well is challenging due to the exponential growth in the number of degrees of freedom with increasing particles.Solution method: A nodal spectral Galerkin scheme is used where the basis functions are constructed to obey the antisymmetry relations of the fermionic wave function. The assembly of these matrices is performed efficiently by exploiting the combinatorial structure of the sparsity patterns.Restrictions: Only one-dimensional computational domains with homogeneous Dirichlet or periodic boundary conditions are supported.Running time: Seconds to minutes     

Real-time Java simulations of multiple interference dielectric filters

An interactive Java applet for real-time simulation and visualization of the transmittance properties of multiple interference dielectric filters is presented. The most commonly used interference filters as well as the state-of-the-art ones are embedded in this platform-independent applet which can serve research and education purposes. The Transmittance applet can be freely downloaded from the site http://cpc.cs.qub.ac.uk.

Program summary

Program title: TransmittanceCatalogue identifier: AEBQ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBQ_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 5778No. of bytes in distributed program, including test data, etc.: 90 474Distribution format: tar.gzProgramming language: JavaComputer: Developed on PC-Pentium platformOperating system: Any Java-enabled OS. Applet was tested on Windows ME, XP, Sun Solaris, Mac OSRAM: VariableClassification: 18Nature of problem: Sophisticated wavelength selective multiple interference filters can include some tens or even hundreds of dielectric layers. The spectral response of such a stack is not obvious. On the other hand, there is a strong demand from application designers and students to get a quick insight into the properties of a given filter.Solution method: A Java applet was developed for the computation and the visualization of the transmittance of multilayer interference filters. It is simple to use and the embedded filter library can serve educational purposes. Also, its ability to handle complex structures will be appreciated as a useful research and development tool.Running time: Real-time simulations     

ISICSoo: A class for the calculation of ionization cross sections from ECPSSR and PWBA theory

ISICS, originally a C language program for calculating K-, L- and M-shell ionization and X-ray production cross sections from ECPSSR and PWBA theory, has been reengineered into a C++ language class, named ISICSoo. The new software design enables the use of ISICS functionality in other software systems. The code, originally developed for Microsoft Windows operating systems, has been ported to Linux and Mac OS platforms to facilitate its use in a wider scientific environment. The reengineered software also includes some fixes to the original implementation, which ensure more robust computational results and a review of some physics parameters used in the computation. The paper describes the software design and the modifications to the implementation with respect to the previous version; it also documents the test process and provides some indications about the software performance.Program summaryProgram title: ISICSooCatalogue identifier: AEKN_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEKN_v1_0.htmlProgram obtainable from: CPC Program Library, Queen?s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 26 420No. of bytes in distributed program, including test data, etc.: 459 548Distribution format: tar.gzProgramming language: C++Computer: 80486 or higher-level PC or MacOperating system: Any OS with gcc compiler version 4.1 (or newer); tested on Scientific Linux 5 (gcc 4.1.2), Mac OS X 10.6.5 (gcc 4.2.1) and Windows XP (MS Visual C++ 2010 Express)Classification: 16.7Nature of problem: Ionization and X-ray production cross section calculations for ion–atom collisions.Solution method: Numerical integration of form factor using a logarithmic transform and Gaussian quadrature, plus exact integration limits.Additional comments: This program is a portable version of the program ADDS_v4_0.Reasons for the new version: Capability of using ISICS physics functionality in other software systems; porting the software to other platforms than Microsoft Windows; improved computational robustness and performance.Summary of revisions: Reengineering into a C++ class; several internal modifications to improve correctness and robustness; updated binding energies tabulations; performance improvements.Running time: The running time depends on the selected atomic shell and the number of polynomials used in the Gaussian quadrature integration. The examples provided only take seconds to run.     

JChainsAnalyser: an ImageJ-based stand-alone application for the study of magneto-rheological fluids

JChainsAnalyser is a Java-based program for the analysis of two-dimensional images of magneto-rheological fluids (MRF) at low concentration of particles obtained using the video-microscopy technique. MRF are colloidal dispersions of micron-sized polarizable particles in a carrier fluid with medium to low viscosity. When a magnetic field is applied to the suspension, the particles aggregate forming chains or clusters. Aggregation dynamics [P. Domínguez-García, S. Melle, J.M. Pastor, M.A. Rubio, Phys. Rev. E 76 (2007) 051403] and morphology of the aggregates [P. Domínguez-García, S. Melle, M.A. Rubio, J. Colloid Interface Sci. 333 (2009) 221-229] have been studied capturing images of the fluid and analyzing them by using this software. The program allows to analyze automatically the MRF images by means of an adequate combination of different imaging methods, while magnitudes and statistics are calculated and saved in data files. It is possible to run the program on a desktop computer, using the GUI (graphical user interface), or in a cluster of processors or remote computer by means of command-line instructions.

Program summary

Program title: JChainsAnalyserCatalogue identifier: AEDT_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDT_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC license, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 79 071No. of bytes in distributed program, including test data, etc.: 4 367 909Distribution format: tar.gzProgramming language: Java 2Computer: Any computer with Java Runtime Environment (JRE) installedOperating system: Any OS with Java Runtime Environment (JRE) installedRAM: Typically, 3.3 MBClassification: 23External routines: ImageJ, ij-imageIO, jdom, L2FProdNature of problem: The video-microscopy technique usually produces quite a big quantity of images to analyze. Although ImageJ gives the required filters and methods for image analysis, it fails when a large number of images is used. Moreover, an adequate combination of filters is needed for the segmentation and binarization of this kind of images.Solution method: JChainsAnalyser filters and analyses any quantity of MRF images automatically, so the application can be run on a desktop computer or using a cluster of processors. It can be run in a desktop computer using the GUI (graphical user interface) or by a command-line interface. JChainsAnalyser uses XML files to define input/output data and Java to ensure portability between operating systems. It also utilizes an image algorithm based on the application of different and adaptative ImageJ's filters.Running time: The test run provided takes only a few seconds.     

Fringe—A Java-based finite fringe analysis package

A package for analysing two-dimensional finite fringe interferograms is described. Through a combination of automatic and interactive routines, an interferogram can be processed to extract the phase shift imparted on the recording light by a transparent object. The package consists of routines to condition and pad the original image for Fourier transform analysis, to filter the image and obtain the phase, to unwrap the phase, and to remove the background phase ramp. A sample image recorded using holographic interferometry is successfully analysed.Program summaryProgram title: FRINGECatalogue identifier: AEMM_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEMM_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 134006No. of bytes in distributed program, including test data, etc.: 4029801Distribution format: tar.gzProgramming language: Java.Computer: Personal Computers.Operating system: Mac OS X, Windows XP, Linux and any other system that can run Java Jar files.RAM: 1GB recommendedClassification: 18.Nature of problem: A standalone multi-platform program to perform analysis of finite fringe interferograms.Solution method: Fourier filtering approach with phase unwrapping and background subtraction.Restrictions: Designed to analyse square images.Running time: Interactive processing takes several minutes. Minimal cpu time.     

SDECAY: a Fortran code for the decays of the supersymmetric particles in the MSSM

We present the Fortran code SDECAY, which calculates the decay widths and branching ratios of all the supersymmetric particles in the Minimal Supersymmetric Standard Model, including higher order effects. Besides the usual two-body decays of sfermions and gauginos and the three-body decays of charginos, neutralinos and gluinos, we have also implemented the three-body decays of stops and sbottoms, and even the four-body decays of the stop; the important loop-induced decay modes are also included. The QCD corrections to the two-body decays involving strongly interacting particles and the dominant components of the electroweak corrections to all decay modes are implemented.

Program summary

Title of program: SDECAY Version 1.1a (March 2005)Catalogue identifier: ADVJProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVJProgram obtainable: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputer for which the program is designed: Any with a Fortran77 systemOperating systems under which the program has been tested: Linux, UnixTypical running time: A few seconds on modern personal computers and workstationsProgramming language used: Fortran77No. of lines in distributed program, including test data, etc.: 59 621No. of bytes in distributed program, including test data, etc.: 338 478Distribution format: tar.gzMemory required to execute (with test data): 7.3 MBDistribution format: ASCIINature of physical problem: Numerical calculation of the decay widths and branching ratios of supersymmetric particles in the Minimal Supersymmetric Standard Model (MSSM). The program calculates two-, three- and four-body decays and loop decays. It includes the SUSY-QCD corrections to two-body decays involving strongly interacting particles. The top-quark decays within the MSSM are evaluated as well.Method of solution: Two-dimensional numerical integration of the analytic formulae for the double differential decay widths of the three-body decays. The other decay widths are calculated analytically.Restrictions on the complexity of the problem: In the higher order decay modes the total decay widths of the virtually exchanged (s)particles are not included in their respective propagators. The higher order decays are calculated when the two-body decays are kinematically closed.     

FeynEdit—a tool for drawing Feynman diagrams

We describe the FeynEdit tool for drawing Feynman diagrams. Input and output is done using the macros of FeynArts, which also implies that diagrams drawn by FeynArts can be edited with FeynEdit. The code can be conveniently transferred using copy-and-paste.

Program summary

Program title: FeynEditCatalogue identifier: AEBX_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBX_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 31 729No. of bytes in distributed program, including test data, etc.: 500 240Distribution format: tar.gzProgramming language: JavaComputer: All Java-capable platformsOperating system: Linux, Mac OS, WindowsRAM: 1-2 MBytesClassification: 4.4Nature of problem: Graphical editing of Feynman diagrams.Solution method: The user copy-and-pastes the LaTeX code of the Feynman diagram into the editor, clicks a button to visualize the diagram, modifies it using the mouse, and finally copy-and-pastes it back into the text.Restrictions: Propagators are presently drawn only as straight lines. This is largely for performance reasons and may be added in a future version. It is not a serious deficit because that information can easily be added in the LaTeX code.Unusual features: Uses FeynArts' LaTeX representation for input and outputRunning time: User-dependent     

Calculating the renormalisation group equations of a SUSY model with Susyno

Susyno is a Mathematica package dedicated to the computation of the 2-loop renormalisation group equations of a supersymmetric model based on any gauge group (the only exception being multiple U(1) groups) and for any field content.Program summaryProgram title: SusynoCatalogue identifier: AEMX_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEMX_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 30829No. of bytes in distributed program, including test data, etc.: 650170Distribution format: tar.gzProgramming language: Mathematica 7 or higher.Computer: All systems that Mathematica 7+ is available for (PC, Mac).Operating system: Any platform supporting Mathematica 7+ (Windows, Linux, Mac OS).Classification: 4.2, 5, 11.1.Nature of problem:Calculating the renormalisation group equations of a supersymmetric model involves using long and complicated general formulae [1, 2]. In addition, to apply them it is necessary to know the Lagrangian in its full form. Building the complete Lagrangian of models with small representations of SU(2) and SU(3) might be easy but in the general case of arbitrary representations of an arbitrary gauge group, this task can be hard, lengthy and error prone.Solution method:The Susyno package uses group theoretical functions to calculate the super-potential and the soft-SUSY-breaking Lagrangian of a supersymmetric model, and calculates the two-loop RGEs of the model using the general equations of [1, 2]. Susyno works for models based on any representation(s) of any gauge group (the only exception being multiple U(1) groups).Restrictions:As the program is based on the formalism of [1, 2], it shares its limitations. Running time can also be a significant restriction, in particular for models with many fields.Unusual features:Susyno contains functions that (a) calculate the Lagrangian of supersymmetric models and (b) calculate some group theoretical quantities. Some of these functions are available to the user and can be freely used. A built-in help system provides detailed information.Running time:Tests were made using a computer with an Intel Core i5 760 CPU, running under Ubuntu 11.04 and with Mathematica 8.0.1 installed. Using the option to suppress printing, the one- and two-loop beta functions of the MSSM were obtained in 2.5 s (NMSSM: 5.4 s). Note that the running time scales up very quickly with the total number of fields in the model.References:[1] S.P. Martin and M.T. Vaughn, Phys. Rev. D 50 (1994) 2282. [Erratum-ibid D 78 (2008) 039903] [arXiv:hep-ph/9311340].[2] Y. Yamada, Phys. Rev. D 50 (1994) 3537 [arXiv:hep-ph/9401241].     

Strongdeco: Expansion of analytical,strongly correlated quantum states into a many-body basis

We provide a Mathematica code for decomposing strongly correlated quantum states described by a first-quantized, analytical wave function into many-body Fock states. Within them, the single-particle occupations refer to the subset of Fock–Darwin functions with no nodes. Such states, commonly appearing in two-dimensional systems subjected to gauge fields, were first discussed in the context of quantum Hall physics and are nowadays very relevant in the field of ultracold quantum gases. As important examples, we explicitly apply our decomposition scheme to the prominent Laughlin and Pfaffian states. This allows for easily calculating the overlap between arbitrary states with these highly correlated test states, and thus provides a useful tool to classify correlated quantum systems. Furthermore, we can directly read off the angular momentum distribution of a state from its decomposition. Finally we make use of our code to calculate the normalization factors for Laughlin?s famous quasi-particle/quasi-hole excitations, from which we gain insight into the intriguing fractional behavior of these excitations.Program summaryProgram title: StrongdecoCatalogue identifier: AELA_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AELA_v1_0.htmlProgram obtainable from: CPC Program Library, Queen?s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 5475No. of bytes in distributed program, including test data, etc.: 31 071Distribution format: tar.gzProgramming language: MathematicaComputer: Any computer on which Mathematica can be installedOperating system: Linux, Windows, MacClassification: 2.9Nature of problem: Analysis of strongly correlated quantum states.Solution method: The program makes use of the tools developed in Mathematica to deal with multivariate polynomials to decompose analytical strongly correlated states of bosons and fermions into a standard many-body basis. Operations with polynomials, determinants and permanents are the basic tools.Running time: The distributed notebook takes a couple of minutes to run.     

ALOHA: Automatic libraries of helicity amplitudes for Feynman diagram computations

We present an application that automatically writes the Helas (HELicity Amplitude Subroutines) library corresponding to the Feynman rules of any quantum field theory Lagrangian. The code is written in Python and takes the Universal FeynRules Output (Ufo) as an input. From this input it produces the complete set of routines, wave-functions and amplitudes, that are needed for the computation of Feynman diagrams at leading as well as at higher orders. The representation is language independent and currently it can output routines in Fortran, C++, and Python. A few sample applications implemented in the MadGraph 5 framework are presented.Program summaryProgram title: ALOHACatalogue identifier: AEMS_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEMS_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: http://www.opensource.org/licenses/UoI-NCSA.phpNo. of lines in distributed program, including test data, etc.: 6094320No. of bytes in distributed program, including test data, etc.: 7479819Distribution format: tar.gzProgramming language: Python2.6Computer: 32/64 bitOperating system: Linux/Mac/WindowsRAM: 512 MbytesClassification: 4.4, 11.6Nature of problem:An effcient numerical evaluation of a squared matrix element can be done with the help of the helicity routines implemented in the HELAS library [1]. This static library contains a limited number of helicity functions and is therefore not always able to provide the needed routine in the presence of an arbitrary interaction. This program provides a way to automatically create the corresponding routines for any given model.Solution method:ALOHA takes the Feynman rules associated to the vertex obtained from the model information (in the Ufo format [2]), and multiplies it by the different wavefunctions or propagators. As a result the analytical expression of the helicity routines is obtained. Subsequently, this expression is automatically written in the requested language (Python, Fortran or C++)Restrictions: The allowed fields are currently spin 0, 1/2, 1 and 2, and the propagators of these particles are canonical.Running time: A few seconds for the SM and the MSSM, and up to a few minutes for models with spin 2 particles.References:[1] Murayama, H. and Watanabe, I. and Hagiwara, K., Helas: HELicity Amplitude Subroutines for Feynman diagram evaluations, KEK-91-11, (1992) http://www-lib.kek.jp/cgi-bin/img_index?199124011[2] C. Degrande, C. Duhr, B. Fuks, D. Grellscheid, O. Mattelaer, et al., Ufo— The Universal FeynRules Output, Comput. Phys. Commun. 183 (2012) 1201-1214. arXiv:1108.2040, doi:10.1016/j.cpc.2012.01.022.     

TSIL: a program for the calculation of two-loop self-energy integrals

TSIL is a library of utilities for the numerical calculation of dimensionally regularized two-loop self-energy integrals. A convenient basis for these functions is given by the integrals obtained at the end of O.V. Tarasov's recurrence relation algorithm. The program computes the values of all of these basis functions, for arbitrary input masses and external momentum. When analytical expressions in terms of polylogarithms are available, they are used. Otherwise, the evaluation proceeds by a Runge-Kutta integration of the coupled first-order differential equations for the basis integrals, using the external momentum invariant as the independent variable. The starting point of the integration is provided by known analytic expressions at (or near) zero external momentum. The code is written in C, and may be linked from C/C++ or Fortran. A Fortran interface is provided. We describe the structure and usage of the program, and provide a simple example application. We also compute two new cases analytically, and compare all of our notations and conventions for the two-loop self-energy integrals to those used by several other groups.

Program summary

Title of program:TSILVersion number: 1.0Catalogue identifier: ADWSProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADWSProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandProgramming language: CPlatform: Any platform supporting the GNU Compiler Collection (gcc), the Intel C compiler (icc), or a similar C compiler with support for complex mathematicsNo. of lines in distributed program, including test data, etc.: 42 730No. of bytes in distributed program, including test data, etc.: 297 101Distribution format: tar.gzNature of physical problem: Numerical evaluation of dimensionally regulated Feynman integrals needed in two-loop self-energy calculations in relativistic quantum field theory in four dimensions.Method of solution: Analytical evaluation in terms of polylogarithms when possible, otherwise through Runge-Kutta solution of differential equations.Limitations: Loss of accuracy in some unnatural threshold cases that do not have vanishing masses.Typical running time: Less than a second.     

NearFar: A computer program for nearside-farside decomposition of heavy-ion elastic scattering amplitude

The NearFar program is a package for carrying out an interactive nearside-farside decomposition of heavy-ion elastic scattering amplitude. The program is implemented in Java to perform numerical operations on the nearside and farside angular distributions. It contains a graphical display interface for the numerical results. A test run has been applied to the elastic scattering at Elab=1503 MeV.

Program summary

Title of program: NearFarCatalogue identifier: ADYP_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYP_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputers: designed for any machine capable of running Java, developed on PC-Pentium-4Operating systems under which the program has been tested: Microsoft Windows XP (Home Edition)Program language used: JavaNumber of bits in a word: 64Memory required to execute with typical data: case dependentNo. of lines in distributed program, including test data, etc.: 3484Number of bytes distributed program, including test data, etc.: 142 051Distribution format: tar.gzOther software required: A Java runtime interpreter, or the Java Development Kit, version 5.0Nature of physical problem: Interactive nearside-farside decomposition of heavy-ion elastic scattering amplitude.Method of solution: The user must supply a external data file or PPSM parameters which calculates theoretical values of the quantities to be decomposed.Typical running time: Problem dependent. In a test run, it is about 35 s on a 2.40 GHz Intel P4-processor machine.     

WIRED — World-Wide Web Interactive Remote Event Display

WIRED (World-Wide Web Interactive Remote Event Display) is a framework, written in the Java^™ language, for building High Energy Physics event displays. An event display based on the WIRED framework enables users of a HEP collaboration to visualise and analyse events remotely using ordinary WWW browsers, on any type of machine. In addition, event displays using WIRED may provide the general public with access to the research of high energy physics.The recent introduction of the object-oriented Java^™ language enables the transfer of machine independent code across the Internet, to be safely executed by a Java enhanced WWW browser. We have employed this technology to create a remote event display in WWW. The combined Java—WWW technology hence assures a world wide availability of such an event display, an always up-to-date program and a platform independent implementation, which is easy to use and to install.     

Program package for multicanonical simulations of U(1) lattice gauge theory

We document our Fortran 77 code for multicanonical simulations of 4D U(1) lattice gauge theory in the neighborhood of its phase transition. This includes programs and routines for canonical simulations using biased Metropolis heatbath updating and overrelaxation, determination of multicanonical weights via a Wang-Landau recursion, and multicanonical simulations with fixed weights supplemented by overrelaxation sweeps. Measurements are performed for the action, Polyakov loops and some of their structure factors. Many features of the code transcend the particular application and are expected to be useful for other lattice gauge theory models as well as for systems in statistical physics.

Program summary

Program title: STMC_U1MUCACatalogue identifier: AEET_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEET_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 18 376No. of bytes in distributed program, including test data, etc.: 205 183Distribution format: tar.gzProgramming language: Fortran 77Computer: Any capable of compiling and executing Fortran codeOperating system: Any capable of compiling and executing Fortran codeClassification: 11.5Nature of problem: Efficient Markov chain Monte Carlo simulation of U(1) lattice gauge theory close to its phase transition. Measurements and analysis of the action per plaquette, the specific heat, Polyakov loops and their structure factors.Solution method: Multicanonical simulations with an initial Wang-Landau recursion to determine suitable weight factors. Reweighting to physical values using logarithmic coding and calculating jackknife error bars.Running time: The prepared tests runs took up to 74 minutes to execute on a 2 GHz PC.     

Monte Carlo event generator for black hole production and decay in proton-proton collisions - QBH version 1.02

We describe the Monte Carlo event generator for black hole production and decay in proton-proton collisions - QBH version 1.02. The generator implements a model for quantum black hole production and decay based on the conservation of local gauge symmetries and democratic decays. The code in written entirely in C++ and interfaces to the PYTHIA 8 Monte Carlo code for fragmentation and decays.

Program summary

Program title: QBHCatalogue identifier: AEGU_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGU_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 10 048No. of bytes in distributed program, including test data, etc.: 118 420Distribution format: tar.gzProgramming language: C++Computer: x86Operating system: Scientific Linux, Mac OS XRAM: 1 GBClassification: 11.6External routines: PYTHIA 8130 (http://home.thep.lu.se/~torbjorn/pythiaaux/present.html) and LHAPDF (http://projects.hepforge.org/lhapdf/)Nature of problem: Simulate black hole production and decay in proton-proton collision.Solution method: Monte Carlo simulation using importance sampling.Running time: Eight events per second.