LIBERI: Library for numerical evaluation of electron-repulsion integrals

We provide a C library, called LIBERI, for numerical evaluation of four-center electron repulsion integrals, based on successive reduction of integral dimension by using Fourier transforms. LIBERI enables us to compute the integrals for numerically defined basis functions within ¹⁰−5 Hartree accuracy as well as their derivatives with respect to the atomic nuclear positions. Damping of the Coulomb interaction can also be imposed to take account of screening effect.

Program summary

Program title: LIBERICatalogue identifier: AEGG_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGG_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.: 44 091No. of bytes in distributed program, including test data, etc.: 1 692 085Distribution format: tar.gzProgramming language: CComputer: allOperating system: any Unix-like systemRAM: 5-10 MbClassification: 7.4External routines: Lapack (http://www.netlib.org/lapack/), Blas (http://www.netlib.org/blas/), FFTW3 (http://www.fftw.org/)Nature of problem: Numerical evaluation of four-center electron-repulsion integrals.Solution method: Four-center electron-repulsion integrals are computed for given basis function set, based on successive reduction of integral dimension using Fourier transform.Running time: 0.5 sec for the demo program supplied with the package.     

PArthENoPE: Public algorithm evaluating the nucleosynthesis of primordial elements

We describe a program for computing the abundances of light elements produced during Big Bang Nucleosynthesis which is publicly available at http://parthenope.na.infn.it/. Starting from nuclear statistical equilibrium conditions the program solves the set of coupled ordinary differential equations, follows the departure from chemical equilibrium of nuclear species, and determines their asymptotic abundances as function of several input cosmological parameters as the baryon density, the number of effective neutrino, the value of cosmological constant and the neutrino chemical potential. The program requires commercial NAG library routines.

Program summary

Program title: PArthENoPECatalogue identifier: AEAV_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAV_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 033No. of bytes in distributed program, including test data, etc.: 46 002Distribution format: tar.gzProgramming language: Fortran 77Computer: PC-compatible running Fortran on Unix, MS Windows or LinuxOperating system: Windows 2000, Windows XP, LinuxClassification: 1.2, 1.9, 17.8External routines: NAG LibrariesNature of problem: Computation of yields of light elements synthesized in the primordial universe.Solution method: BDF method for the integration of the ODEs, implemented in a NAG routine.Running time: 90 sec with default parameters on a Dual Xeon Processor 2.4 GHz with 2 GB RAM.     

Calculation of four-particle harmonic-oscillator transformation brackets

A procedure for precise calculation of the three- and four-particle harmonic-oscillator (HO) transformation brackets is presented. The analytical expressions of the four-particle HO transformation brackets are given. The computer code for the calculations of HO transformation brackets proves to be quick, efficient and produces results with small numerical uncertainties.

Program summary

Program title: HOTBCatalogue identifier: AEFQ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFQ_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.: 1247No. of bytes in distributed program, including test data, etc.: 6659Distribution format: tar.gzProgramming language: FORTRAN 90Computer: Any computer with FORTRAN 90 compilerOperating system: Windows, Linux, FreeBSD, True64 UnixRAM: 8 MBClassification: 17.17Nature of problem: Calculation of the three-particle and four-particle harmonic-oscillator transformation brackets.Solution method: The method is based on compact expressions of the three-particle harmonics oscillator brackets, presented in [1] and expressions of the four-particle harmonics oscillator brackets, presented in this paper.Restrictions: The three- and four-particle harmonic-oscillator transformation brackets up to the e=28.Unusual features: Possibility of calculating the four-particle harmonic-oscillator transformation brackets.Running time: Less than one second for the single harmonic-oscillator transformation bracket.References:

[1]

G.P. Kamuntavi?ius, R.K. Kalinauskas, B.R. Barret, S. Mickevi?ius, D. Germanas, Nuclear Physics A 695 (2001) 191.

     

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).     

Symbolic computation of the Hartree-Fock energy from a chiral EFT three-nucleon interaction at NLO

We present the first of a two-part Mathematica notebook collection that implements a symbolic approach for the application of the density matrix expansion (DME) to the Hartree-Fock (HF) energy from a chiral effective field theory (EFT) three-nucleon interaction at N²LO. The final output from the notebooks is a Skyrme-like energy density functional that provides a quasi-local approximation to the non-local HF energy. In this paper, we discuss the derivation of the HF energy and its simplification in terms of the scalar/vector-isoscalar/isovector parts of the one-body density matrix. Furthermore, a set of steps is described and illustrated on how to extend the approach to other three-nucleon interactions.

Program summary

Program title: SymbHFNNNCatalogue identifier: AEGC_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGC_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.: 96 666No. of bytes in distributed program, including test data, etc.: 378 083Distribution format: tar.gzProgramming language: Mathematica 7.1Computer: Any computer running Mathematica 6.0 and later versionsOperating system: Windows Xp, Linux/UnixRAM: 256 MbClassification: 5, 17.16, 17.22Nature of problem: The calculation of the HF energy from the chiral EFT three-nucleon interaction at N²LO involves tremendous spin-isospin algebra. The problem is compounded by the need to eventually obtain a quasi-local approximation to the HF energy, which requires the HF energy to be expressed in terms of scalar/vector-isoscalar/isovector parts of the one-body density matrix. The Mathematica notebooks discussed in this paper solve the latter issue.Solution method: The HF energy from the chiral EFT three-nucleon interaction at N²LO is cast into a form suitable for an automatic simplification of the spin-isospin traces. Several Mathematica functions and symbolic manipulation techniques are used to obtain the result in terms of the scalar/vector-isoscalar/isovector parts of the one-body density matrix.Running time: Several hours     

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     

JuNoLo - Jülich nonlocal code for parallel post-processing evaluation of vdW-DF correlation energy

Nowadays the state of the art Density Functional Theory (DFT) codes are based on local (LDA) or semilocal (GGA) energy functionals. Recently the theory of a truly nonlocal energy functional has been developed. It has been used mostly as a post-DFT calculation approach, i.e. by applying the functional to the charge density calculated using any standard DFT code, thus obtaining a new improved value for the total energy of the system. Nonlocal calculation is computationally quite expensive and scales as N² where N is the number of points in which the density is defined, and a massively parallel calculation is welcome for a wider applicability of the new approach. In this article we present a code which accomplishes this goal.

Program summary

Program title: JuNoLoCatalogue identifier: AEFM_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFM_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 980No. of bytes in distributed program, including test data, etc.: 2 126 072Distribution format: tar.gzProgramming language: Fortran 90Computer: any architecture with a Fortran 90 compilerOperating system: Linux, AIXHas the code been vectorised or parallelized?: Yes, from 1 to 65536 processors may be used.RAM: depends strongly on the problem's size.Classification: 7.3External routines:• FFTW (http://www.tw.org/)• MPI (http://www.mcs.anl.gov/research/projects/mpich2/ or http://www.lam-mpi.org/)Nature of problem: Obtaining the value of the nonlocal vdW-DF energy based on the charge density distribution obtained from some Density Functional Theory code.Solution method: Numerical calculation of the double sum is implemented in a parallel F90 code. Calculation of this sum yields the required nonlocal vdW-DF energy.Unusual features: Binds to virtually any DFT program.Additional comments: Excellent parallelization features.Running time: Depends strongly on the size of the problem and the number of CPUs used.     

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.     

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.     

Algorithmic derivation of Dyson-Schwinger equations

We present an algorithm for the derivation of Dyson-Schwinger equations of general theories that is suitable for an implementation within a symbolic programming language. Moreover, we introduce the Mathematica package DoDSE¹ which provides such an implementation. It derives the Dyson-Schwinger equations graphically once the interactions of the theory are specified. A few examples for the application of both the algorithm and the DoDSE package are provided.

Program summary

Program title: DoDSECatalogue identifier: AECT_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECT_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.: 105 874No. of bytes in distributed program, including test data, etc.: 262 446Distribution format: tar.gzProgramming language: Mathematica 6 and higherComputer: all on which Mathematica is availableOperating system: all on which Mathematica is availableClassification: 11.1, 11.4, 11.5, 11.6Nature of problem: Derivation of Dyson-Schwinger equations for a theory with given interactions.Solution method: Implementation of an algorithm for the derivation of Dyson-Schwinger equations.Unusual features: The results can be plotted as Feynman diagrams in Mathematica.Running time: Less than a second to minutes for Dyson-Schwinger equations of higher vertex functions.     

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.     

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.     

A cascadic monotonic time-discretized algorithm for finite-level quantum control computation

A computer package (CNMS) is presented aimed at the solution of finite-level quantum optimal control problems. This package is based on a recently developed computational strategy known as monotonic schemes.Quantum optimal control problems arise in particular in quantum optics where the optimization of a control representing laser pulses is required. The purpose of the external control field is to channel the system's wavefunction between given states in its most efficient way. Physically motivated constraints, such as limited laser resources, are accommodated through appropriately chosen cost functionals.

Program summary

Program title: CNMSCatalogue identifier: ADEB_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADEB_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.: 770No. of bytes in distributed program, including test data, etc.: 7098Distribution format: tar.gzProgramming language: MATLAB 6Computer: AMD Athlon 64 × 2 Dual, 2:21 GHz, 1:5 GB RAMOperating system: Microsoft Windows XPWord size: 32Classification: 4.9Nature of problem: Quantum controlSolution method: IterativeRunning time: 60-600 sec     

REACH: A program for coarse-grained biomolecular simulation

REACH (Realistic Extension Algorithm viaCovariance Hessian) is a program package for residue-scale coarse-grained biomolecular simulation. The program calculates the force constants of a residue-scale elastic network model in single-domain proteins using the variance-covariance matrix obtained from atomistic molecular dynamics simulation. Secondary-structure dependence of the force constants is integrated. The method involves self-consistent, direct mapping of atomistic simulation results onto a coarse-grained force field in an efficient automated procedure without requiring iterative fits and avoiding system dependence.

Program summary

Program title: REACHCatalogue identifier: AEDA_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDA_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.: 42 244No. of bytes in distributed program, including test data, etc.: 3 682 118Distribution format: tar.gzProgramming language: FORTRAN 77Computer: x86 PCOperating system: GNU/Linux, SUSE and Red HatRAM: Depends on the system size to be calculatedWord size: 32 or 64 bitsClassification: 3External routines: LAPACK, BLASNature of problem: A direct calculation of force field for residue-scale coarse-grained biomolecular simulation derived from atomistic molecular dynamics trajectory.Solution method: A variance-covariance matrix and the associated Hessian (second-derivative) matrix are calculated from an atomistic molecular dynamics trajectory of single-domain protein internal motion and the off-diagonal Hessian matrix is fitted to that of a residue-scale elastic network model. The resulting force constants for the residue pair interactions are expressed as model functions as a function of pairwise distance.Running time: Depends on the system size and the number of MD trajectory frames used. The test run provided with the distribution takes only a few seconds to execute.     

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.     

BOKASUN: A fast and precise numerical program to calculate the Master Integrals of the two-loop sunrise diagrams

We present the program BOKASUN for fast and precise evaluation of the Master Integrals of the two-loop self-mass sunrise diagram for arbitrary values of the internal masses and the external four-momentum. We use a combination of two methods: a Bernoulli accelerated series expansion and a Runge-Kutta numerical solution of a system of linear differential equations.

Program summary

Program title: BOKASUNCatalogue identifier: AECG_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECG_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.: 9404No. of bytes in distributed program, including test data, etc.: 104 123Distribution format: tar.gzProgramming language: FORTRAN77Computer: Any computer with a Fortran compiler accepting FORTRAN77 standard. Tested on various PC's with LINUXOperating system: LINUXRAM: 120 kbytesClassification: 4.4Nature of problem: Any integral arising in the evaluation of the two-loop sunrise Feynman diagram can be expressed in terms of a given set of Master Integrals, which should be calculated numerically. The program provides a fast and precise evaluation method of the Master Integrals for arbitrary (but not vanishing) masses and arbitrary value of the external momentum.Solution method: The integrals depend on three internal masses and the external momentum squared p². The method is a combination of an accelerated expansion in 1/p² in its (pretty large!) region of fast convergence and of a Runge-Kutta numerical solution of a system of linear differential equations.Running time: To obtain 4 Master Integrals on PC with 2 GHz processor it takes 3 μs for series expansion with pre-calculated coefficients, 80 μs for series expansion without pre-calculated coefficients, from a few seconds up to a few minutes for Runge-Kutta method (depending on the required accuracy and the values of the physical parameters).     

Milne, a routine for the numerical solution of Milne's problem

The routine Milne provides accurate numerical values for the classical Milne's problem of neutron transport for the planar one speed and isotropic scattering case. The solution is based on the Case eigen-function formalism. The relevant X functions are evaluated accurately by the Double Exponential quadrature. The calculated quantities are the extrapolation distance and the scalar and the angular fluxes. Also, the H function needed in astrophysical calculations is evaluated as a byproduct.

Program summary

Program title: MilneCatalogue identifier: AEGS_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGS_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.: 701No. of bytes in distributed program, including test data, etc.: 6845Distribution format: tar.gzProgramming language: Fortran 77Computer: PC under Linux or WindowsOperating system: Ubuntu 8.04 (Kernel version 2.6.24-16-generic), Windows-XPClassification: 4.11, 21.1, 21.2Nature of problem: The X functions are integral expressions. The convergence of these regular and Cauchy Principal Value integrals are impaired by the singularities of the integrand in the complex plane. The DE quadrature scheme tackles these singularities in a robust manner compared to the standard Gauss quadrature.Running time: The test included in the distribution takes a few seconds to run.