Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aeea_v1_0.tar.gz(3730 Kbytes)
Manuscript Title: 2DRMP: a suite of two-dimensional R-matrix propagation codes
Authors: N.S. Scott, M.P. Scott, P.G. Burke, T. Stitt, V. Faro-Maza a, C. Denis, A. Maniopoulou
Program title: 2DRMP:
Catalogue identifier: AEEA_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 180(2009)2424
Programming language: Fortran 95, MPI.
Computer: Tested on CRAY XT4 [1]; IBM eServer 575 [2]; Itanium II cluster [3].
Operating system: Tested on UNICOS/lc [1]; IBM AIX [2]; Red Hat Linux Enterprise AS [3].
Has the code been vectorised or parallelized?: Yes. 16 cores were used for small test run.
Keywords: atomic collision processes, electron impact excitation, electron impact ionization, FARM, intermediate energies, IERM, R-matrix, RMPS, two-dimensional R-matrix propagation.
PACS: 34.80.Dp.
Classification: 2.4.

External routines: BLAS, LAPACK, PBLAS, ScaLAPACK.

Subprograms used:
Cat Id Title Reference
ADAZ_v1_1 FARM_SDRMP CPC 180(2009)2450

Nature of problem:
2DRMP is a suite of programs aimed at creating virtual experiments on high performance architectures to enable the study of electron scattering from H-like atoms and ions at intermediate energies.

Solution method:
Two-dimensional R-matrix propagation theory. The (r1, r2) space of the internal region is subdivided into a number of subregions. Local R-matrices are constructed within each subregion and used to propagate a global R-matrix, ℜ, across the internal region. On the boundary of the internal region ℜ is transformed onto the IERM target state basis. Thus, the two-dimensional R-matrix propagation technique transforms an intractable problem into a series of tractable problems enabling the internal region to be extended far beyond that which is possible with the standard one-sector codes. A distinctive feature of the method is that both electrons are treated identically and the R-matrix basis states are constructed to allow for both electrons to be in the continuum. The subregion size is flexible and can be adjusted to accommodate the number of cores available.

The implementation is currently restricted to electron scattering from H-like atoms and ions.

Additional comments:
The programs have been designed to operate on serial computers and to exploit the distributed memory parallelism found on tightly coupled high performance clusters and supercomputers. 2DRMP has been systematically and comprehensively documented using ROBODoc [4] which is an API documentation tool that works by extracting specially formatted headers from the program source code and writing them to documentation files.

Running time:
The wall clock running time for the small test run using 16 cores and performed on [3] is as follows: bp (7 secs); rint2 (34 secs); newrd (32 secs); diag (21 secs); amps (11 secs); prop (24 secs).

[1] HECToR, CRAY XT4 running UNICOS/lc, http://www.hector.ac.uk/ (accessed 22 July, 2009).
[2] HPCx, IBM eServer 575 running IBM AIX , http://www.hpcx.ac.uk/ (accessed 22 July, 2009).
[3] HP Cluster,Itanium II cluster running Red Hat Linux Enterprise AS, Queen s University Belfast, http://www.qub.ac.uk/directorates/InformationServices/ Research/HighPerformanceComputing/Services/Hardware/HPResearch/ (accessed 22 July, 2009).
[4] Automating Software Documentation with ROBODoc, http://www.xs4all.nl/~rfsber/Robo/ (accessed 22 July, 2009).