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     

TiReX: Replica-exchange molecular dynamics using Tinker

We present a driver program for performing replica-exchange molecular dynamics simulations with the Tinker package. Parallelization is based on the Message Passing Interface, with every replica assigned to a separate process. The algorithm is not communication intensive, which makes the program suitable for running even on loosely coupled cluster systems. Particular attention is paid to the practical aspects of analyzing the program output.

Program summary

Program title: TiReXCatalogue identifier: AEEK_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEEK_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.: 43 385No. of bytes in distributed program, including test data, etc.: 502 262Distribution format: tar.gzProgramming language: Fortran 90/95Computer: Most UNIX machinesOperating system: LinuxHas the code been vectorized or parallelized?: parallelized with MPIClassification: 16.13External routines: TINKER version 4.2 or 5.0, built as a libraryNature of problem: Replica-exchange molecular dynamics.Solution method: Each replica is assigned to a separate process; temperatures are swapped between replicas at regular time intervals.Running time: The sample run may take up to a few minutes.     

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     

Motion4D - A library for lightrays and timelike worldlines in the theory of relativity

The Motion4D-library solves the geodesic equation as well as the parallel- and Fermi-Walker-transport in four-dimensional Lorentzian spacetimes numerically. Initial conditions are given with respect to natural local tetrads which are adapted to the symmetries or the coordinates of the spacetime. Beside some already implemented metrics like the Schwarzschild and Kerr metric, the object oriented structure of the library permits to implement other metrics or integrators in a straight forward manner.

Program summary

Program title: Motion4D-libraryCatalogue identifier: AEEX_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEEX_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.: 150 425No. of bytes in distributed program, including test data, etc.: 5 139 407Distribution format: tar.gzProgramming language: C++Computer: All platforms with a C++ compilerOperating system: Linux, Unix, WindowsRAM: 39 MBytesClassification: 1.5External routines: Gnu Scientific Library (GSL) (http://www.gnu.org/software/gsl/)Nature of problem: Solve geodesic equation, parallel and Fermi-Walker transport in four-dimensional Lorentzian spacetimes.Solution method: Integration of ordinary differential equationsRunning time: The test runs provided with the distribution require only a few seconds to run.     

MDVRY: a polarizable classical molecular dynamics package for biomolecules

The MDVRY classical molecular dynamics package is presented for the study of biomolecules in the gas and liquid phase. Electrostatic polarization has been implemented in the formalism of point induced dipoles following the model of Thole. Two schemes have been implemented for the calculation of induced dipoles, i.e. resolution of the self-consistent equations and a ‘Car-Parrinello’ dynamical approach. In this latter, the induced dipoles are calculated at each time step of the dynamics through the dynamics of additional degrees of freedom associated with the dipoles. This method saves computer time and allows to study polarized solvated proteins at a very low CPU cost. The program is written in C-language and runs on LINUX machines. A detailed manual of the code is given. The main features of the package are illustrated taking on examples of proteins in the gas phase or immersed in liquid water.

Program summary

Program title: MDVRYCatalogue identifier: AEBY_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBY_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.: 39 156No. of bytes in distributed program, including test data, etc.: 277 197Distribution format: tar.bz2Programming language: CComputer: Linux machines with FFTW Fourier Transform package installedOperating system: Linux machines, SUSE & RedHat distributionsClassification: 3, 16.13, 23External routines: FFTW (http://www.fftw.org/)Nature of problem: Molecular Dynamics Software package.Solution method: Velocity Verlet algorithm. The implemented force field is composed of intra-molecular interactions and inter-molecular interactions (electrostatics, polarization, van der Waals). Polarization is accounted through induced point dipoles at each atomic site. Supplementary degrees of freedom are associated to the induced dipoles so that a modified Hamiltonian of the dynamics is written. This allows to calculate the induced dipoles with a very fast ‘Car-Parrinello’ type of dynamics.Running time: The test run provided takes approximately 6 minutes to run.     

Solving surface structures from normal incidence X-ray standing wave data

A program is provided to determine structural parameters of atoms in or adsorbed on surfaces by refinement of atomistic models towards experimentally determined data generated by the normal incidence X-ray standing wave (NIXSW) technique. The method employs a combination of Differential Evolution Genetic Algorithms and Steepest Descent Line Minimisations to provide a fast, reliable and user friendly tool for experimentalists to interpret complex multidimensional NIXSW data sets.

Program summary

Program title: NIXSW Planewave SolverCatalogue identifier: ADZE_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZE_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.: 16 874No. of bytes in distributed program, including test data, etc.: 1 631 874Distribution format: tar.gzProgramming language: Borland C++ Builder 5Computer: Any Windows CompatibleOperating system: Windows 2000 and XPRAM: <10 MBClassification: 7.4Nature of problem: Using NIXSW experimental data to calculate atomic positions of adsorbates.Restrictions: Data from substrates must have cubic, tetragonal or orthorhombic crystal structures i.e. with 90° between conventional cell axes.Running time: Seconds-minutes dependant on the number of plane waves and the number of atomic sites.     

A completely open source based computing system for computer generation of Fourier holograms

Computer generated holograms are usually generated using commercial software like MATLAB, MATHCAD, Mathematica, etc. This work is an approach in doing the same using freely distributed open source packages and Operating System. A Fourier hologram is generated using this method and tested for simulated and optical reconstruction. The reconstructed images are in good agreement with the objects chosen. The significance of using such a system is also discussed.

Program summary

Program title: FHOLOCatalogue identifier: AEDS_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDS_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.: 176 336No. of bytes in distributed program, including test data, etc.: 4 294 872Distribution format: tar.gzProgramming language: C++Computer: any X86 micro computerOperating system: Linux (Debian Etch)RAM: 512 MBClassification: 18Nature of problem: To generate a Fourier Hologram in micro computer only by using open source operating system and packages.Running time: Depends on the matrix size. 10 sec for a matrix of size 256×256.     

Object-oriented design patterns in Fortran 90/95: mazev1, mazev2 and mazev3

This paper discusses the concept, application, and usefulness of software design patterns for scientific programming in Fortran 90/95. An example from the discipline of object-oriented design patterns, that of a game based on navigation through a maze, is used to describe how some important patterns can be implemented in Fortran 90/95 and how the progressive introduction of design patterns can usefully restructure Fortran software as it evolves. This example is complemented by a discussion of how design patterns have been used in a real-life simulation of Particle-in-Cell plasma physics. The following patterns are mentioned in this paper: Factory, Strategy, Template, Abstract Factory and Facade.

Program summary

Program title: mazev1, mazev2, mazev3Catalogue identifier: AEAI_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAI_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.: 1958No. of bytes in distributed program, including test data, etc.: 17 100Distribution format: tar.gzProgramming language: Fortran 95Computer: PC/MacOperating system: Unix/Linux/Mac (FreeBSD)/Windows (Cygwin)RAM: These are interactive programs with small (KB) memory requirementsClassification: 6.5, 20Nature of problem: A sequence of programs which demonstrate the use of object oriented design patterns for the restructuring of Fortran 90/95 software. The programs implement a simple maze game similar to that described in [1].Solution method: Restructuring uses versions of the Template, Strategy and Factory design patterns.Running time: Interactive.References:

[1] 

E. Gamma, R. Helm, R. Johnson, J. Vlissides, Design Patterns: Elements of Reusable Object Oriented Software, Addison-Wesley, 1995, ISBN 0201633612.

     

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.     

An object-oriented C++ implementation of Davidson method for finding a few selected extreme eigenpairs of a large, sparse, real, symmetric matrix

A C++ class named Davidson is presented for determining a few eigenpairs with lowest or alternatively highest values of a large, real, symmetric matrix. The algorithm described by Stathopoulos and Fischer is used. The exception mechanism is involved to report the errors. The class is written in ANSI C++, so it is fully portable. In addition a console program as well as a program with graphical user interface for Microsoft Windows is attached, which allow one to calculate the lowest eigenstates of time-independent Schrödinger equation for a given binding potential in one, two or three spatial dimensions. The package contains the classes providing often used potential functions (model atom potential, Coulomb potential, square well potential and Kramers-Henneberger well potential) as well as a possibility to use any potential stored in a file (then any dimensionality of the problem is allowed).The described code is the subject of M.Sc. thesis of T.D. prepared under the supervision of J.M.

Program summary

Program title: DavidsonCatalogue identifier: ADZM_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZM_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.: 3 037 055No. of bytes in distributed program, including test data, etc.: 20 002 609Distribution format: tar.gzProgramming language: C++Computer: AllOperating system: AnyRAM: User's parameters dependentWord size: 32 and 64 bitsSupplementary material: Test results for the 2D and 3D cases is availableClassification: 4, 4.8Nature of problem: Finding a few extreme eigenpairs of a real, symmetric, sparse matrix. Examples in quantum optics (interaction of matter with a laser field).Solution method: Davidson algorithmRunning time: The test example included in the distribution package (1D matrix) takes approximately 30 minutes to run. 2D matrix calculations can take hours and 3D, days, to run.     

Gauss-Legendre and Chebyshev quadratures for singular integrals

Exact expressions are presented for efficient computation of the weights in Gauss-Legendre and Chebyshev quadratures for selected singular integrands. The singularities may be of Cauchy type, logarithmic type or algebraic-logarithmic end-point branching points. We provide Fortran 90 routines for computing the weights for both the Gauss-Legendre and the Chebyshev (Fejér-1) meshes whose size can be set by the user.

New program summary

Program title: SINGQUADCatalogue identifier: AEBR_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBR_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.: 4128No. of bytes in distributed program, including test data, etc.: 25 815Distribution format: tar.gzProgramming language: Fortran 90Computer: Any with a Fortran 90 compilerOperating system: Linux, Windows, MacRAM: Depending on the complexity of the problemClassification: 4.11Nature of problem: Program provides Gauss-Legendre and Chebyshev (Fejér-1) weights for various singular integrands.Solution method: The weights are obtained from the condition that the quadrature of order N must be exact for a polynomial of degree?(N−1). The weights are expressed as moments of the singular kernels associated with Legendre or Chebyshev polynomials. These moments are obtained in analytic form amenable for computation.Additional comments: If the NAGWare f95 compiler is used, the option, “-kind = byte”, must be included in the compile command lines of the Makefile.Running time: The test run supplied with the distribution takes a couple of seconds to execute.     

S@M, a mathematica implementation of the spinor-helicity formalism

In this paper we present the package S@M (Spinors@Mathematica) which implements the spinor-helicity formalism in Mathematica. The package allows the use of complex-spinor algebra along with the multi-purpose features of Mathematica. The package defines the spinor objects with their basic properties along with functions to manipulate them. It also offers the possibility of evaluating the spinorial objects numerically at every computational step. The package is therefore well suited to be used in the context of on-shell technology, in particular for the evaluation of scattering amplitudes at tree- and loop-level.

Program summary

Program title: S@MCatalogue identifier: AEBF_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBF_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.: 14 404No. of bytes in distributed program, including test data, etc.: 77 536Distribution format: tar.gzProgramming language: MathematicaComputer: All computers running MathematicaOperating system: Any system running MathematicaClassification: 4.4, 5, 11.1Nature of problem: Implementation of the spinor-helicity formalismSolution method: Mathematica implementationRunning time: The notebooks provided with the package take only a few seconds to run.     

MontePython: Implementing Quantum Monte Carlo using Python

We present a cross-language C++/Python program for simulations of quantum mechanical systems with the use of Quantum Monte Carlo (QMC) methods. We describe a system for which to apply QMC, the algorithms of variational Monte Carlo and diffusion Monte Carlo and we describe how to implement theses methods in pure C++ and C++/Python. Furthermore we check the efficiency of the implementations in serial and parallel cases to show that the overhead using Python can be negligible.

Program summary

Program title: MontePythonCatalogue identifier: ADZP_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZP_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.: 49 519No. of bytes in distributed program, including test data, etc.: 114 484Distribution format: tar.gzProgramming language: C++, PythonComputer: PC, IBM RS6000/320, HP, ALPHAOperating system: LINUXHas the code been vectorised or parallelized?: Yes, parallelized with MPINumber of processors used: 1-96RAM: Depends on physical system to be simulatedClassification: 7.6; 16.1Nature of problem: Investigating ab initio quantum mechanical systems, specifically Bose-Einstein condensation in dilute gases of ⁸⁷RbSolution method: Quantum Monte CarloRunning time: 225 min with 20 particles (with 4800 walkers moved in 1750 time steps) on 1 AMD Opteron^TM Processor 2218 processor; Production run for, e.g., 200 particles takes around 24 hours on 32 such processors.     

GeodesicViewer - A tool for exploring geodesics in the theory of relativity

The GeodesicViewer realizes exocentric two- and three-dimensional illustrations of lightlike and timelike geodesics in the general theory of relativity. By means of an intuitive graphical user interface, all parameters of a spacetime as well as the initial conditions of the geodesics can be modified interactively. This makes the GeodesicViewer a useful instrument for the exploration of geodesics in four-dimensional Lorentzian spacetimes.

Program summary

Program title: GeodesicViewerCatalogue identifier: AEFP_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFP_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.: 168 868No. of bytes in distributed program, including test data, etc.: 6 076 202Distribution format: tar.gzProgramming language: C++, Qt, Qwt, OpenGLComputer: All platforms with a C++ compiler, Qt, Qwt, OpenGLOperating system: Linux, Mac OS XRAM: 24 MbytesClassification: 1.5External routines:

•

Gnu Scientific Library (GSL) (http://www.gnu.org/software/gsl/)

•

Motion4D (included in the package). The Motion4D library can also be downloaded from CPC. Catalogue identifier: AEEX

•

Qt (http://qt.nokia.com/downloads)

•

Qwt (http://qwt.sourceforge.net/)

•

OpenGL (http://www.opengl.org/)

Nature of problem: Illustrate geodesics in four-dimensional Lorentzian spacetimes.Solution method: Integration of ordinary differential equations. 3D-Rendering via OpenGL.Running time: Interactive. The examples given take milliseconds.     

JJGEN: A flexible program for generating lists of jj-coupled configuration state functions

The success of large scale relativistic multiconfiguration Dirac-Hartree-Fock calculations for atomic systems rely on judiciously chosen configuration expansions. Dependent on the atomic system as well as on the studied properties, various correlation effects need to be considered. Based on the active set approach, this program allows the user to generate general lists of jj-coupled configuration state functions to be used as input to the grasp2K multiconfiguration Dirac-Hartree-Fock package [P. Jönsson, X. He, C. Froese Fischer, I.P. Grant, Comput. Phys. Comm. (2007), in press].

Program summary

Program title: JJGENCatalogue identifier: ADZG_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZG_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 673No. of bytes in distributed program, including test data, etc.: 430 543Distribution format: tar.gzProgramming language: FortranComputer: Intel compatible PCOperating system: Linux, UnixWord size: 32 bitsClassification: 7.3Nature of problem: Generation of lists of jj-coupled configuration state functions to describe different electron correlation effects in many-electron atoms.Solution method: From a set of reference configurations a list of jj-coupled configuration state functions is generated by excitations to an active set of orbitals. Imposing restrictions on the allowed excitations the configuration expansion can be targeted to describe different correlation effects.Restrictions: The complexity of the cases that can be handled is entirely determined by the grasp2K package [P. Jönsson, X. He, C. Froese Fischer, I.P. Grant, Comput. Phys. Comm. (2007), in press] used for the generation of the electronic wave-functions.Running time: CPU time required to execute test cases: few seconds.     

Line-by-line spectroscopic simulations on graphics processing units

We report here on software that performs line-by-line spectroscopic simulations on gases. Elaborate models (such as narrow band and correlated-K) are accurate and efficient for bands where various components are not simultaneously and significantly active. Line-by-line is probably the most accurate model in the infrared for blends of gases that contain high proportions of H₂O and CO₂ as this was the case for our prototype simulation. Our implementation on graphics processing units sustains a speedup close to 330 on computation-intensive tasks and 12 on memory intensive tasks compared to implementations on one core of high-end processors. This speedup is due to data parallelism, efficient memory access for specific patterns and some dedicated hardware operators only available in graphics processing units. It is obtained leaving most of processor resources available and it would scale linearly with the number of graphics processing units in parallel machines. Line-by-line simulation coupled with simulation of fluid dynamics was long believed to be economically intractable but our work shows that it could be done with some affordable additional resources compared to what is necessary to perform simulations on fluid dynamics alone.

Program summary

Program title: GPU4RECatalogue identifier: ADZY_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZY_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.: 62 776No. of bytes in distributed program, including test data, etc.: 1 513 247Distribution format: tar.gzProgramming language: C++Computer: x86 PCOperating system: Linux, Microsoft Windows. Compilation requires either gcc/g++ under Linux or Visual C++ 2003/2005 and Cygwin under Windows. It has been tested using gcc 4.1.2 under Ubuntu Linux 7.04 and using Visual C++ 2005 with Cygwin 1.5.24 under Windows XP.RAM: 1 gigabyteClassification: 21.2External routines: OpenGL (http://www.opengl.org)Nature of problem: Simulating radiative transfer on high-temperature high-pressure gases.Solution method: Line-by-line Monte-Carlo ray-tracing.Unusual features: Parallel computations are moved to the GPU.Additional comments: nVidia GeForce 7000 or ATI Radeon X1000 series graphics processing unit is required.Running time: A few minutes.     

Plato: A localised orbital based density functional theory code

The Plato package allows both orthogonal and non-orthogonal tight-binding as well as density functional theory (DFT) calculations to be performed within a single framework. The package also provides extensive tools for analysing the results of simulations as well as a number of tools for creating input files. The code is based upon the ideas first discussed in Sankey and Niklewski (1989) [1] with extensions to allow high-quality DFT calculations to be performed. DFT calculations can utilise either the local density approximation or the generalised gradient approximation. Basis sets from minimal basis through to ones containing multiple radial functions per angular momenta and polarisation functions can be used. Illustrations of how the package has been employed are given along with instructions for its utilisation.

Program summary

Program title: PlatoCatalogue identifier: AEFC_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFC_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.: 219 974No. of bytes in distributed program, including test data, etc.: 1 821 493Distribution format: tar.gzProgramming language: C/MPI and PERLComputer: Apple Macintosh, PC, Unix machinesOperating system: Unix, Linux and Mac OS XHas the code been vectorised or parallelised?: Yes, up to 256 processors testedRAM: Up to 2 Gbytes per processorClassification: 7.3External routines: LAPACK, BLAS and optionally ScaLAPACK, BLACS, PBLAS, FFTWNature of problem: Density functional theory study of electronic structure and total energies of molecules, crystals and surfaces.Solution method: Localised orbital based density functional theory.Restrictions: Tight-binding and density functional theory only, no exact exchange.Unusual features: Both atom centred and uniform meshes available. Can deal with arbitrary angular momenta for orbitals, whilst still retaining Slater–Koster tables for accuracy.Running time: Test cases will run in a few minutes, large calculations may run for several days.     

HypExp 2, Expanding hypergeometric functions about half-integer parameters

In this article, we describe a new algorithm for the expansion of hypergeometric functions about half-integer parameters. The implementation of this algorithm for certain classes of hypergeometric functions in the already existing Mathematica package HypExp is described. Examples of applications in Feynman diagrams with up to four loops are given.

New version program summary

Program title:HypExp 2Catalogue identifier:ADXF_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXF_v2_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.:106 401No. of bytes in distributed program, including test data, etc.:2 668 729Distribution format:tar.gzProgramming language:MathematicaComputer:Computers running MathematicaOperating system:Linux, Windows, MacRAM:Depending on the complexity of the problemSupplementary material:Library files which contain the expansion of certain hypergeometric functions around their parameters are availableClassification:4.7, 5Does the new version supersede the previous version?:YesNature of problem:Expansion of hypergeometric functions about parameters that are integer and/or half-integer valued.Solution method:New algorithm implemented in Mathematica.Reasons for new version:Expansion about half-integer parameters.Summary of revisions:Ability to expand about half-integer valued parameters added.Restrictions:The classes of hypergeometric functions with half-integer parameters that can be expanded are listed below.Additional comments:The package uses the package HPL included in the distribution.Running time:Depending on the expansion.