A program for the predictor-corrector numerov method

Title of program: PCNUM Catalogue number: AARJ Program ontainable from: CPC Program Library, Queen's University of Belfast, N. Ireland (see application form in this issue) Computer: IBM 370/3031; Installation: University of Windsor Computer Centre Operating system: OSMVT under VM Programming language used: FORTRAN IV High speed storage required: 10 Kwords (for 8 channels) No. of bits in a word: 32 Overlay structure: none No. of magnetic tapes required: none Other peripherals used: card reader, line printer No. of cards in combined program and test deck: 829 Card punching code: EBCDIC     

Efficient evolution of unpolarized and polarized parton distributions with QCD-Pegasus

The Fortran package QCD-Pegasus is presented. This program provides fast, flexible and accurate solutions of the evolution equations for unpolarized and polarized parton distributions of hadrons in perturbative QCD. The evolution is performed using the symbolic moment-space solutions on a one-fits-all Mellin inversion contour. User options include the order of the evolution including the next-to-next-to-leading order in the unpolarized case, the type of the evolution including an emulation of brute-force solutions, the evolution with a fixed number nf of flavors or in the variable-nf scheme, and the evolution with a renormalization scale unequal to the factorization scale. The initial distributions are needed in a form facilitating the computation of the complex Mellin moments.

Program summary

Title of program: QCD-PegasusVersion: 1.0Catalogue identifier: ADVNProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVNProgram obtainable from: CPC Program Library Queen's University of Belfast, N. IrelandLicense: GNU Public LicenseComputers: allOperating systems: allProgram language:Fortran 77 (using the common compiler extension of procedure names with more than six characters)Memory required to execute: negligible (<1 MB)Other programs called: noneExternal files needed: noneNumber of lines in distributed program, including test data, etc.: 8157Number of bytes in distributed program, including test data, etc.: 240 578Distribution format: tar.gzNature of the physical problem: Solution of the evolution equations for the unpolarized and polarized parton distributions of hadrons at leading order (LO), next-to-leading order and next-to-next-to-leading order of perturbative QCD. Evolution performed either with a fixed number nf of effectively massless quark flavors or in the variable-nf scheme. The calculation of observables from the parton distributions is not part of the present package.Method of solution: Analytic solution in Mellin space (beyond LO in general by power-expansion around the lowest-order expansion) followed by a fast Mellin inversion to x-space using a fixed one-fits-all contour.Restrictions on complexity of the problem: The initial distributions for the evolution are required in a form facilitating an efficient calculation of their complex Mellin moments. The ratio of the renormalization and factorization scales μr/μ has to be a fixed number.Typical running time: One to ten seconds, on a PC with a 2.0 GHz Pentium-IV processor, for performing the evolution of 200 initial distributions to 500 (x,μ) points each. For more details see Section 6.     

ARVO: A Fortran package for computing the solvent accessible surface area and the excluded volume of overlapping spheres via analytic equations

In calculating the solvation energy of proteins, the hydration effects, drug binding, molecular docking, etc., it is important to have an efficient and exact algorithms for computing the solvent accessible surface area and the excluded volume of macromolecules. Here we present a Fortran package based on the new exact analytical methods for computing volume and surface area of overlapping spheres. In the considered procedure the surface area and volume are expressed as surface integrals of the second kind over the closed region. Using the stereographic projection the surface integrals are transformed to a sum of double integrals which are reduced to the curve integrals. MPI Fortran version is described as well. The package is also useful for computing the percolation probability of continuum percolation models.

Program summary

Title of program: ARVOCatalogue identifier: ADULProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADULProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandOperating system under which the program has been tested: LINUX system and Windows systemProgramming language used: FORTRANComputer: PC Pentium; SPP'2000;Number of bytes in distributed program, including test data, etc.: 322 633Number of lines in distributed program, including test data, etc.: 5051Distribution format: tar.gzCard punching code: ASCIINature of physical problem: Molecular mechanics computations, continuum percolations.Method of solution: Numerical algorithm based on the analytical formulas, after using the stereographic transformation.Restriction on the complexity of the problem: The program does not account explicitly for cavities inside the molecule.Typical running time: Depends on the size of the molecule under consideration.Unusual features of the program: No     

Parallel implementation of an adaptive and parameter-free N-body integrator

Previously, Pruett et al. (2003) [3] described an N-body integrator of arbitrarily high order M with an asymptotic operation count of O(M²N²). The algorithm's structure lends itself readily to data parallelization, which we document and demonstrate here in the integration of point-mass systems subject to Newtonian gravitation. High order is shown to benefit parallel efficiency. The resulting N-body integrator is robust, parameter-free, highly accurate, and adaptive in both time-step and order. Moreover, it exhibits linear speedup on distributed parallel processors, provided that each processor is assigned at least a handful of bodies.

Program summary

Program title: PNB.f90Catalogue identifier: AEIK_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEIK_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.: 3052No. of bytes in distributed program, including test data, etc.: 68 600Distribution format: tar.gzProgramming language: Fortran 90 and OpenMPIComputer: All shared or distributed memory parallel processorsOperating system: Unix/LinuxHas the code been vectorized or parallelized?: The code has been parallelized but has not been explicitly vectorized.RAM: Dependent upon NClassification: 4.3, 4.12, 6.5Nature of problem: High accuracy numerical evaluation of trajectories of N point masses each subject to Newtonian gravitation.Solution method: Parallel and adaptive extrapolation in time via power series of arbitrary degree.Running time: 5.1 s for the demo program supplied with the package.     

A correction to the POTLIB Library described in “POTLIB 2001: A potential energy surface library for chemical systems” : [Comput. Phys. Comm. 144 (2002) 169-187]

Program summary

Program Title: POTLIB 2001, version 1.0Catalogue identifier: ADPJProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADPJProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland     

Applying linear interaction energy method for binding affinity calculations of podophyllotoxin analogues with tubulin using continuum solvent model and prediction of cytotoxic activity

Podophyllotoxin and its analogues have important therapeutic value in the treatment of cancer, due to their ability to induce apoptosis in cancer cells in a proliferation-independent manner. These ligands bind to colchicine binding site of tubulin near the α- and β-tubulin interface and interfere with tubulin polymerization. The binding free energies of podophyllotoxin-based inhibitors of tubulin were computed using a linear interaction energy (LIE) method with a surface generalized Born (SGB) continuum solvation model. A training set of 76 podophyllotoxin analogues was used to build a binding affinity model for estimating the free energy of binding for 36 inhibitors (test set) with diverse structural modifications. The average root mean square error (RMSE) between the experimental and predicted binding free energy values was 0.56 kcal/mol which is comparable to the level of accuracy achieved by the most accurate methods, such as free energy perturbation (FEP) or thermodynamic integration (TI). The squared correlation coefficient between experimental and SGB–LIE estimates for the free energy for the test set compounds is also significant (R² = 0.733). On the basis of the analysis of the binding energy, we propose that the three-dimensional conformation of the A, B, C and D rings is important for interaction with tubulin. On the basis of this insight, 12 analogues of varying ring modification were taken, tested with LIE methodology and then validated with their experimental potencies of tubulin polymerization inhibition. Low levels of RMSE for the majority of inhibitors establish the structure-based LIE method as an efficient tool for generating more potent and specific inhibitors of tubulin by testing rationally designed lead compounds based on podophyllotoxin derivatization.     

Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

Here we develop simple numerical algorithms for both stationary and non-stationary solutions of the time-dependent Gross-Pitaevskii (GP) equation describing the properties of Bose-Einstein condensates at ultra low temperatures. In particular, we consider algorithms involving real- and imaginary-time propagation based on a split-step Crank-Nicolson method. In a one-space-variable form of the GP equation we consider the one-dimensional, two-dimensional circularly-symmetric, and the three-dimensional spherically-symmetric harmonic-oscillator traps. In the two-space-variable form we consider the GP equation in two-dimensional anisotropic and three-dimensional axially-symmetric traps. The fully-anisotropic three-dimensional GP equation is also considered. Numerical results for the chemical potential and root-mean-square size of stationary states are reported using imaginary-time propagation programs for all the cases and compared with previously obtained results. Also presented are numerical results of non-stationary oscillation for different trap symmetries using real-time propagation programs. A set of convenient working codes developed in Fortran 77 are also provided for all these cases (twelve programs in all). In the case of two or three space variables, Fortran 90/95 versions provide some simplification over the Fortran 77 programs, and these programs are also included (six programs in all).

Program summary

Program title: (i) imagetime1d, (ii) imagetime2d, (iii) imagetime3d, (iv) imagetimecir, (v) imagetimesph, (vi) imagetimeaxial, (vii) realtime1d, (viii) realtime2d, (ix) realtime3d, (x) realtimecir, (xi) realtimesph, (xii) realtimeaxialCatalogue identifier: AEDU_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDU_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.: 122 907No. of bytes in distributed program, including test data, etc.: 609 662Distribution format: tar.gzProgramming language: FORTRAN 77 and Fortran 90/95Computer: PCOperating system: Linux, UnixRAM: 1 GByte (i, iv, v), 2 GByte (ii, vi, vii, x, xi), 4 GByte (iii, viii, xii), 8 GByte (ix)Classification: 2.9, 4.3, 4.12Nature of problem: These programs are designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in one-, two- or three-space dimensions with a harmonic, circularly-symmetric, spherically-symmetric, axially-symmetric or anisotropic trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Solution method: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation, in either imaginary or real time, over small time steps. The method yields the solution of stationary and/or non-stationary problems.Additional comments: This package consists of 12 programs, see “Program title”, above. FORTRAN77 versions are provided for each of the 12 and, in addition, Fortran 90/95 versions are included for ii, iii, vi, viii, ix, xii. For the particular purpose of each program please see the below.Running time: Minutes on a medium PC (i, iv, v, vii, x, xi), a few hours on a medium PC (ii, vi, viii, xii), days on a medium PC (iii, ix).

Program summary (1)

Title of program: imagtime1d.FTitle of electronic file: imagtime1d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 1 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in one-space dimension with a harmonic trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems.

Program summary (2)

Title of program: imagtimecir.FTitle of electronic file: imagtimecir.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 1 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in two-space dimensions with a circularly-symmetric trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems.

Program summary (3)

Title of program: imagtimesph.FTitle of electronic file: imagtimesph.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 1 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in three-space dimensions with a spherically-symmetric trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems.

Program summary (4)

Title of program: realtime1d.FTitle of electronic file: realtime1d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in one-space dimension with a harmonic trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems.

Program summary (5)

Title of program: realtimecir.FTitle of electronic file: realtimecir.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in two-space dimensions with a circularly-symmetric trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems.

Program summary (6)

Title of program: realtimesph.FTitle of electronic file: realtimesph.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77Typical running time: Minutes on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in three-space dimensions with a spherically-symmetric trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems.

Program summary (7)

Title of programs: imagtimeaxial.F and imagtimeaxial.f90Title of electronic file: imagtimeaxial.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Few hours on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in three-space dimensions with an axially-symmetric trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems.

Program summary (8)

Title of program: imagtime2d.F and imagtime2d.f90Title of electronic file: imagtime2d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 2 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Few hours on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in two-space dimensions with an anisotropic trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems.

Program summary (9)

Title of program: realtimeaxial.F and realtimeaxial.f90Title of electronic file: realtimeaxial.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 4 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time Hours on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in three-space dimensions with an axially-symmetric trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems.

Program summary (10)

Title of program: realtime2d.F and realtime2d.f90Title of electronic file: realtime2d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 4 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Hours on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in two-space dimensions with an anisotropic trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems.

Program summary (11)

Title of program: imagtime3d.F and imagtime3d.f90Title of electronic file: imagtime3d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum RAM memory: 4 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Few days on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in three-space dimensions with an anisotropic trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in imaginary time over small time steps. The method yields the solution of stationary problems.

Program summary (12)

Title of program: realtime3d.F and realtime3d.f90Title of electronic file: realtime3d.tar.gzCatalogue identifier:Program summary URL:Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandDistribution format: tar.gzComputers: PC/Linux, workstation/UNIXMaximum Ram Memory: 8 GByteProgramming language used: Fortran 77 and Fortran 90Typical running time: Days on a medium PCUnusual features: NoneNature of physical problem: This program is designed to solve the time-dependent Gross-Pitaevskii nonlinear partial differential equation in three-space dimensions with an anisotropic trap. The Gross-Pitaevskii equation describes the properties of a dilute trapped Bose-Einstein condensate.Method of solution: The time-dependent Gross-Pitaevskii equation is solved by the split-step Crank-Nicolson method by discretizing in space and time. The discretized equation is then solved by propagation in real time over small time steps. The method yields the solution of stationary and non-stationary problems.     

Oscillation theorems for second-order nonlinear neutral differential equations

The purpose of this paper is to study the oscillation of the second-order neutral differential equations of the form

(a(t)[z^′(t)]^γ)^′+q(t)x^β(σ(t))=0,     

MERADGEN 1.0: Monte Carlo generator for the simulation of radiative events in parity conserving doubly-polarized Møller scattering

The Monte Carlo generator MERADGEN 1.0 for the simulation of radiative events in parity conserving doubly-polarized Møller scattering has been developed. Analytical integration wherever it is possible provides rather fast and accurate generation. Some numerical tests and histograms are presented.

Program summary

Program title: MERADGEN 1.0Catalogue identifier:ADYM_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYM_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneProgramming language: FORTRAN 77Computer(s) for which the program has been designed: allOperating system(s) for which the program has been designed: LinuxRAM required to execute with typical data: 1 MBNo. of lines in distributed program, including test data, etc.:2196No. of bytes in distributed program, including test data, etc.:23 501Distribution format:tar.gzHas the code been vectorized or parallelized?: noNumber of processors used: 1Supplementary material: noneExternal routines/libraries used: noneCPC Program Library subprograms used: noneNature of problem: Simulation of radiative events in parity conserving doubly-polarized Møller scattering.Solution method: Monte Carlo method for simulation within QED, analytical integration wherever it is possible that provides rather fast and accurate generation.Restrictions: noneUnusual features: noneAdditional comments: noneRunning time: The simulation of 10⁸ radiative events for itest:=1 takes up to 45 seconds on AMD Athlon 2.80 GHz processor.     

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.     

QCDMAPT: Program package for Analytic approach to QCD

A program package, which facilitates computations in the framework of Analytic approach to QCD, is developed and described in detail. The package includes both the calculated explicit expressions for relevant spectral functions up to the four-loop level and the subroutines for necessary integrals.

Program summary

Program title: QCDMAPTCatalogue identifier: AEGP_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGP_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.: 2579No. of bytes in distributed program, including test data, etc.: 180 052Distribution format: tar.gzProgramming language: Maple 9 and higherComputer: Any which supports Maple 9Operating system: Any which supports Maple 9Classification: 11.1, 11.5, 11.6Nature of problem: Subroutines helping computations within Analytic approach to QCD.Solution method: A program package for Maple is provided. It includes both the explicit expressions for relevant spectral functions and the subroutines for basic integrals used in the framework of Analytic approach to QCD.Running time: Template program running time is about a minute (depends on CPU).     

HPL, a Mathematica implementation of the harmonic polylogarithms

In this paper, we present an implementation of the harmonic polylogarithm of Remiddi and Vermaseren [E. Remiddi, J.A.M. Vermaseren, Int. J. Modern Phys. A 15 (2000) 725, hep-ph/9905237] for Mathematica. It contains an implementation of the product algebra, the derivative properties, series expansion and numerical evaluation. The analytic continuation has been treated carefully, allowing the user to keep the control over the definition of the sign of the imaginary parts. Many options enables the user to adapt the behavior of the package to his specific problem.

Program summary

Program title: HPLCatalogue identifier:ADWXProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADWXProgram obtained from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions:noneProgramming language: MathematicaNo. of lines in distributed program, including test data, etc.:13 310No. of bytes in distributed program, including test data, etc.: 1 990 584Distribution format: tar.gzComputer:all computers running MathematicaOperating systems:operating systems running MathematicaNature of problem: Computer algebraic treatment of the harmonic polylogarithms which appear in the evaluation of Feynman diagramsSolution method: Mathematica implementation     

The orbifolder: A tool to study the low-energy effective theory of heterotic orbifolds

The orbifolder is a program developed in C++ that computes and analyzes the low-energy effective theory of heterotic orbifold compactifications. The program includes routines to compute the massless spectrum, to identify the allowed couplings in the superpotential, to automatically generate large sets of orbifold models, to identify phenomenologically interesting models (e.g. MSSM-like models) and to analyze their vacuum configurations.Program summaryProgram title: orbifolderCatalogue identifier: AELR_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AELR_v1_0.htmlProgram obtainable from: CPC Program Library, Queen?s University, Belfast, N. IrelandLicensing provisions: GNU General Public License version 3No. of lines in distributed program, including test data, etc.: 145 572No. of bytes in distributed program, including test data, etc.: 930 517Distribution format: tar.gzProgramming language: C++Computer: Personal computerOperating system: Tested on Linux (Fedora 15, Ubuntu 11, SuSE 11)Word size: 32 bits or 64 bitsClassification: 11.1External routines: Boost (http://www.boost.org/), GSL (http://www.gnu.org/software/gsl/)Nature of problem: Calculating the low-energy spectrum of heterotic orbifold compactifications.Solution method: Quadratic equations on a lattice; representation theory; polynomial algebra.Running time: Less than a second per model.     

A recursive method to calculate UV-divergent parts at one-loop level in dimensional regularization

A method is introduced to calculate the UV-divergent parts at one-loop level in dimensional regularization. The method is based on the recursion, and the basic integrals are just the scaleless integrals after the recursive reduction, which involve no other momentum scales except the loop momentum itself. The method can be easily implemented in any symbolic computer language, and a implementation in Mathematica is ready to use.Program summaryProgram title: UVPartCatalogue identifier: AELY_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AELY_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 361No. of bytes in distributed program, including test data, etc.: 412 084Distribution format: tar.gzProgramming language: MathematicaComputer: Any computer where the Mathematica is running.Operating system: Any capable of running Mathematica.Classification: 11.1External routines: FeynCalc (http://www.feyncalc.org/), FeynArts (http://www.feynarts.de/)Nature of problem: To get the UV-divergent part of any one-loop expression.Solution method: UVPart is a Mathematica package where the recursive method has been implemented.Running time: In general it is below one second.     

A fast algorithm for voxel-based deterministic simulation of X-ray imaging

Deterministic method based on ray tracing technique is known as a powerful alternative to the Monte Carlo approach for virtual X-ray imaging. The algorithm speed is a critical issue in the perspective of simulating hundreds of images, notably to simulate tomographic acquisition or even more, to simulate X-ray radiographic video recordings. We present an algorithm for voxel-based deterministic simulation of X-ray imaging using voxel-driven forward and backward perspective projection operations and minimum bounding rectangles (MBRs). The algorithm is fast, easy to implement, and creates high-quality simulated radiographs. As a result, simulated radiographs can typically be obtained in split seconds with a simple personal computer.

Program summary

Program title: X-rayCatalogue identifier: AEAD_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAD_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.: 416 257No. of bytes in distributed program, including test data, etc.: 6 018 263Distribution format: tar.gzProgramming language: C (Visual C++)Computer: Any PC. Tested on DELL Precision 380 based on a Pentium D 3.20 GHz processor with 3.50 GB of RAMOperating system: Windows XPClassification: 14, 21.1Nature of problem: Radiographic simulation of voxelized objects based on ray tracing technique.Solution method: The core of the simulation is a fast routine for the calculation of ray-box intersections and minimum bounding rectangles, together with voxel-driven forward and backward perspective projection operations.Restrictions: Memory constraints. There are three programs in all.

•

A. Program for test 3.1(1): Object and detector have axis-aligned orientation;

•

B. Program for test 3.1(2): Object in arbitrary orientation;

•

C. Program for test 3.2: Simulation of X-ray video recordings.

1.

Program A Memory required to execute with typical data: 207 Megabytes, depending on the size of the input file. Typical running time: 2.30 s. (Tested in release mode, the same below.)

2.

Program B (the main program) Memory required to execute with typical data: 114 Megabytes, depending on the size of the input file. Typical running time: 1.60 s.

3.

Program C Memory required to execute with typical data: 215 Megabytes, depending on the size of the input file. Typical computation time: 27.26 s for cast-5, 101.87 s for cast-6.

     

Slavnov-Taylor1.0: A Mathematica package for computation in BRST formalism

Slavnov-Taylor is a Mathematica package which allows us to perform automatic symbolic computation in BRST formalism. This article serves as a self-contained guide to prospective users, and indicates the conventions and approximations used.

Program summary

Title of program:Slavnov-TaylorCatalogue identifier:ADSSProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSSProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandProgramming Language: Mathematica 4.0Platform: Any platform supporting Mathematica 4.0Computers tested on: Pentium PCOperating systems under which the program has been tested: LinuxMemory required to execute: Minimal: 1.254.784 bytes, Standard: 1.281.248 bytesNo. of bytes in distributed program, including test data, etc.: 8368 bytesDistribution format: tar gzip fileKeywords: Slavnov-Taylor, BRST, MathematicaNature of physical problem: Symbolic computation in the Slavnov-Taylor formalism for gauge theories in 4-dimensional space-time based on a semi-simple compact Lie group for a general BRS transformations.Restrictions on the complexity of the problem: Only matter in the adjoint is allowed.Typical running time: less than one second     

-XSummer- Transcendental functions and symbolic summation in Form

Harmonic sums and their generalizations are extremely useful in the evaluation of higher-order perturbative corrections in quantum field theory. Of particular interest have been the so-called nested sums, where the harmonic sums and their generalizations appear as building blocks, originating for example, from the expansion of generalized hypergeometric functions around integer values of the parameters. In this paper we discuss the implementation of several algorithms to solve these sums by algebraic means, using the computer algebra system Form.

Program summary

Title of program:XSummerCatalogue identifier:ADXQ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXQ_v1_0Program obtainable from:CPC Program Library, Queen's University of Belfast, N. IrelandLicense:GNU Public License and Form LicenseComputers:allOperating system:allProgram language:FormMemory required to execute:Depending on the complexity of the problem, recommended at least 64 MB RAMNo. of lines in distributed program, including test data, etc.:9854No. of bytes in distributed program, including test data, etc.:126 551Distribution format:tar.gzOther programs called:noneExternal files needed:noneNature of the physical problem:Systematic expansion of higher transcendental functions in a small parameter. The expansions arise in the calculation of loop integrals in perturbative quantum field theory.Method of solution:Algebraic manipulations of nested sums.Restrictions on complexity of the problem:Usually limited only by the available disk space.Typical running time:Dependent on the complexity of the problem.