Intel Fortran Research Articles

Patterns of Nanoporous Spherical Packing Emerging under Influence of Curvature and Confinement

Nanoporous membranes are popular in nanotechnology due to biomedical and industrial applications. During the past decade, experimental, theoretical and computational research into porous membranes and soft materials has opened up new mathematical dimensions. In bulk, diblock copolymers exhibit ordered morphologies such as parallel matrices of lamellae, bicontinuous matrices of gyroids, hexagonal matrices of cylinders and body-centred cubic matrices of spheres. In melt, confinement plays an essential role in tuning the frustration of the diblock copolymer system to predict more nanostructures. These nanostructures depend on the composition of the copolymers, their confining geometries and the degree of structural frustration. An isotropic 9-point stencil for Laplacian is constructed. The discrete finite-difference technique is used in polar grids to discretize the macromolecule of the diblock copolymer system to study spherical patterns to study the effect of curvature and confinement with a well-known and efficient cell dynamic simulation model. Intel FORTRAN (IFORT) codes are generated to run the CDS model and visualisation of simulation results is observed with the help of OPENDX. A comparison of the proposed study with existing experimental and computational studies is also presented.

FINITE ELEMENT MODELLING AND ANALYSIS OF FREE VIBRATIONS OF AXIALLY FUNCTIONALY GRADED BEAMS WITH NON-UNIFORM CROSS-SECTIONS

The finite element package Abaqus is used for modeling and analysis of free vibrations of axially functional graded beams with non-uniform cross-sections. Considerable attention has been paid to this aspect of the geometry, as it has a significant effect on the behavior of the beam as a structural component. User-defined material model subroutines (UMAT) were developed in a unified way for one- and three-dimensional models of the specified beams to provide numerical implementation of functional material gradients within each finite element. UMAT procedures were programmed in FORTRAN in the MS Visual Studio environment and compiled with the main program package using the Intel Fortran Compiler. In the procedures, material heterogeneity was assigned to each material integration point of appropriate standard one- and three-dimensional finite elements, eliminating the need to design graded finite elements via user-defined element subroutine. Frequencies and mode shapes of free vibrations of one- and three-dimensional axially functional graded beams with non-uniform cross-sections were found using the Lanczos method in Abaqus. An analysis of the accuracy and effectiveness of the proposed modeling approach was carried out, comparing the obtained results with data known in the literature. The presented modeling technique based on user-defined material subroutines provides valuable opportunities for scientists and engineers engaged in the dynamic analysis of structural components made of functionally gradient materials and having variable geometry along the axial direction. Keywords: axially functionally graded beams; non-uniform cross-section beams; free vibrations; graded finite element; Abaqus user defined subroutines

An Interactive System Based on the IASP91 Earth Model for Earthquake Data Processing

System software for interactive human–computer data processing based on the IASP91 Earth model was designed. An interactive data processing system for visualizing earthquake data was designed and implemented via the Intel Fortran platform. The system reads and processes broadband seismic data acquired by field stations, mainly including the reading and import of raw data, pre-processing, identification of seismic phases and inter-correlation traveltimes picking. In the data processing step, shortcomings have been improved and functions have been gradually refined and enhanced, making it easier and faster to process data. It has already processed more than 1000 large seismic events received by the station from 2013 to 2018. The practical application shows that the human–computer interaction system is easy to operate, accurate, fast and flexible, and is an effective tool for processing seismic data.

Advanced exergy analysis of a PEM fuel cell with hydrogen energy storage integrated with organic Rankine cycle for electricity generation

In this study, a system consisting of PEMFC and ORC is modeled and analyzed from the thermodynamic aspects. Accordingly, to have a better understanding of the performance of the system, advanced exergy analysis is applied. The analysis is developed in Intel Fortran Compiler within the Microsoft Visual Studio. For modeling all dependent components, separate subroutines in the Microsoft Visual Studio platform were made by implementing the first and second law of thermodynamics. Based on the conventional exergy destruction rate, for the hybrid power system, Ėx,D,total is 484.4 kW. The results show that the exergy destruction rate in the PEMFC stack is the highest among other elements. The PEM fuel cell exergy efficiency is 46.89% according to the conventional exergy study. Moreover, the PEM fuel cell and overall electrical efficiencies are 39.06% and 44.81%, respectively. According to the advanced exergy analysis, except for the condenser and compressor, the avoidable part of exergy destruction rates is higher than those unavoidable. The largest endogenous and exogenous exergy destruction rates belong to PEMFC and compressor correspondingly. In addition, the compressor has a poor performance since its endogenous exergy destruction is 34.24 kW, and its unavoidable one is 19.55 kW. Besides, a large portion of its unavoidable (i.e., 57.09%) cannot be wiped out even we use the best technology. Moreover, the enhancement of components performance breeds eliminating the percentage of exogenous available exergy destruction rate of the turbine (57.66%) with the highest value.In the present work, the parametric study is conducted by increasing the current density as an essential design parameter to analyze its effects on some critical parameters such as the overall efficiency, the net generated power, and so on. For implementing the parametric study, the Engineering Equation Solver (EES) was used. The combination of Fortran and EES helped to obtain the required outcomes from the overall system. Besides, the overall efficiency increases until it reaches a maximum value of 47.31%and drops thereafter.

UnDiFi-2D: An unstructured discontinuity fitting code for 2D grids

UnDiFi-2D, an open source (free software) Unstructured-grid, Discontinuity Fitting code, is presented. The aim of UnDiFi-2D is to model gas-dynamic discontinuities in two-dimensional (2D) flows as if they were true discontinuities of null thickness that bound regions of the flow-field where a smooth solution to the governing PDEs exists. UnDiFi-2D therefore needs to be coupled with an unstructured CFD solver that is used to discretize the governing PDEs within the smooth regions of the flow-field. Two different, in-house developed, CFD solvers are also included in the current distribution.The main features of the UnDiFi-2D software can be summarized as follows:Programming languageUnDiFi-2D is written in standard Fortran 77/95; its design is highly modular in order to enhance simplicity of use, maintenance and allow coupling with virtually any existing CFD solver;Usability, maintenance and enhancement In order to improve the usability, maintenance and enhancement of the code also the documentation has been carefully taken into account. The git distributed versioning system has been adopted to facilitate collaborative maintenance and code development;CopyrightsUnDiFi-2D is a free software that anyone can use, copy, distribute, change and improve under the GNU Public License version 3.The present paper is a manifesto of the first public release of the UnDiFi-2D code. It describes the currently implemented features, which are the result of more than a decade of still ongoing CFD developments. This work is focused on the computational techniques adopted and a detailed description of the main characteristics is reported. UnDiFi-2D capabilities are demonstrated by means of examples test cases. The design of the code allows to easily include existing CFD codes and is aimed at ease code reuse and readability. Program summaryProgram title: UnDiFi-2DCPC Library link to program files:https://doi.org/10.17632/5hwssmc2mx.1Licensing provisions: GNU General Public License, version 3Programming language: Fortran; developed and tested with Intel Fortran Compiler v. 18.0.3 and GNU gfortran.External routines: The code depends on several libraries and third-party packages which are detailed in the corpus of the text.Nature of problem: Numerical computation of flows with discontinuities.Solution method: Shock-fitting technique.Additional comments including restrictions and unusual features:•At present, UnDiFi-2D is validated for inviscid steady and unsteady two-dimensional flows without changes in the number of discontinuity lines and interaction points.•UnDiFi-2D implements a shock-fitting algorithm and can be coupled with unstructured cell-vertex solvers, with an Arbitrary Lagrangian-Eulerian (ALE) formulation.•UnDiFi-2D project adopts git [1], a free and open source distributed version control system. A public repository dedicated to UnDiFi-2D project [2] has been created on github [3], a web-based hosting service for software development projects using git versioning system. Finally, a comprehensive documentation is provided in the form of user manual developed in Pandoc [4].

The diagnostic capability of the Cray, gfortran, Intel, Nag and Oracle Fortran compilers

This version has been updated with details of the Intel 19.1 release and the Nag 7.0 release.

The diagnostic capability of the Cray, gfortran, Intel, Nag and Oracle Fortran compilers

We use several compilers for the development of the examples in our books. The following table summarises the compilers we have used:

GenASiSBasics: Object-oriented utilitarian functionality for large-scale physics simulations (Version 3)

GenASiS Basics provides Fortran 2003 classes furnishing extensible object-oriented utilitarian functionality for large-scale physics simulations on distributed memory supercomputers. This functionality includes physical units and constants; display to the screen or standard output device; message passing; I/O to disk; and runtime parameter management and usage statistics. This revision---Version 3 of Basics---includes a significant name change, some minor additions to functionality, and a major addition to functionality: infrastructure facilitating the offloading of computational kernels to devices such as GPUs.

A coupled model for simulating water flow and solute transport in furrow irrigation

For optimal water and fertilizer management under furrow irrigation, it is important to understand the water and solute dynamics on the land surface and in the subsurface. An efficient mathematical tool is required to describe these dynamic processes. We propose a coupled model in which surface water flow and solute transport are described using the zero-inertia equation and the average cross-sectional convection-dispersion equation, respectively, while the two-dimensional Richards equation and the convection-dispersion equation are used to simulate water flow and solute transport in soils, respectively. Solutions are computed numerically using finite differences for surface water flow and finite volumes for solute transports in furrow. Subsurface water flow and solute transport equations are solved using the CHAIN_2D code. An iterative method is used to couple computations of surface and subsurface processes. Both surface and subsurface water flow and solute transport modules are coded in program subroutines and functions in the Intel FORTRAN environment. The coupled model was validated by comparing its simulation results with measured data. Results showed that simulated water front advances in the furrow and water contents in the soil agreed with the observations reasonably well. Good simulations can be achieved with a relatively fine temporal resolution. Numerical oscillations can be eliminated by adopting appropriate time steps. As compared with the traditional furrow irrigation model, the proposed model can better quantify soil water and solute dynamics by considering interactions between surface and subsurface water flow and solute transport processes. The proposed model can be used as a decision tool to design and manage furrow irrigation.

Modelling fracture in heterogeneous materials on HPC systems using a hybrid MPI/Fortran coarray multi-scale CAFE framework

A 3D multi-scale cellular automata finite element (CAFE) framework for modelling fracture in heterogeneous materials is described. The framework is implemented in a hybrid MPI/Fortran coarray code for efficient parallel execution on HPC platforms. Two open source BSD licensed libraries developed by the authors in modern Fortran were used: CGPACK, implementing cellular automata (CA) using Fortran coarrays, and ParaFEM, implementing finite elements (FE) using MPI. The framework implements a two-way concurrent hierarchical information exchange between the structural level (FE) and the microstructure (CA). MPI to coarrays interface and data structures are described. The CAFE framework is used to predict transgranular cleavage propagation in a polycrystalline iron round bar under tension. Novel results enabled by this CAFE framework include simulation of progressive cleavage propagation through individual grains and across grain boundaries, and emergence of a macro-crack from merging of cracks on preferentially oriented cleavage planes in individual crystals. Nearly ideal strong scaling up to at least tens of thousands of cores was demonstrated by CGPACK and by ParaFEM in isolation in prior work on Cray XE6. Cray XC30 and XC40 platforms and CrayPAT profiling were used in this work. Initially the strong scaling limit of hybrid CGPACK/ParaFEM CAFE model was 2000 cores. After replacing all-to-all communication patterns with the nearest neighbour algorithms the strong scaling limit on Cray XC30 was increased to 7000 cores. TAU profiling on non-Cray systems identified deficiencies in Intel Fortran 16 optimisation of remote coarray operations. Finally, coarray synchronisation challenges and opportunities for thread parallelisation in CA are discussed.

Numerical and Analytical Methods of Study of Stochastic Systems with Delay

This paper is a brief review of methods of qualitative and quantitative study of various classes of stochastic systems with aftereffect. We describe schemes of analysis of stochastic ordinary and partial differential and integrodifferential equations, which are also applicable for their deterministic analogs. These numerical and analytical schemes are implemented in the software package Mathematica and in the language Intel Fortran.

Symbolic and numeric scheme for solution of linear integro-differential equations with random parameter uncertainties and Gaussian stochastic process input

• A system of linear stochastic integro-differential equations is under consideration. • Gaussian nonstationary noises and random model uncertainties excite the system. • The problem is to find the mean and covariance functions of the state vector process. • Green matrix-valued functions are used to solve the problem. • As example, a response of a thermoviscoelastic system is estimated by our scheme. The paper describes a theoretical apparatus and an algorithmic part of application of the Green matrix-valued functions for time-domain analysis of systems of linear stochastic integro-differential equations. It is suggested that these systems are subjected to Gaussian nonstationary stochastic noises in the presence of model parameter uncertainties that are described in the framework of the probability theory. If the uncertain model parameter is fixed to a given value, then a time-history of the system will be fully represented by a second-order Gaussian vector stochastic process whose properties are completely defined by its conditional vector-valued mean function and matrix-valued covariance function. The scheme that is proposed is constituted of a combination of two subschemes. The first one explicitly defines closed relations for symbolic and numeric computations of the conditional mean and covariance functions, and the second one calculates unconditional characteristics by the Monte Carlo method. A full scheme realized on the base of Wolfram Mathematica and Intel Fortran software programs, is demonstrated by an example devoted to an estimation of a nonstationary stochastic response of a mechanical system with a thermoviscoelastic component. Results obtained by using the proposed scheme are compared with a reference solution constructed by using a direct Monte Carlo simulation.

GenASiS Mathematics : Object-oriented manifolds, operations, and solvers for large-scale physics simulations

The large-scale computer simulation of a system of physical fields governed by partial differential equations requires some means of approximating the mathematical limit of continuity. For example, conservation laws are often treated with a ‘finite-volume’ approach in which space is partitioned into a large number of small ‘cells,’ with fluxes through cell faces providing an intuitive discretization modeled on the mathematical definition of the divergence operator. Here we describe and make available Fortran 2003 classes furnishing extensible object-oriented implementations of simple meshes and the evolution of generic conserved currents thereon, along with individual ‘unit test’ programs and larger example problems demonstrating their use. These classes inaugurate the Mathematics division of our developing astrophysics simulation code GenASiS (Gen eral A strophysical Si mulation S ystem), which will be expanded over time to include additional meshing options, mathematical operations, solver types, and solver variations appropriate for many multiphysics applications. Program SummaryProgram Title: SineWave, SawtoothWave, RiemannProblem, RayleighTaylor, SedovTaylor, and FishboneMoncrief (fluid dynamics example problems illustrating GenASiSMathematics)Program Files doi:http://dx.doi.org/10.17632/mzvxngwtw6.1Licensing provisions: GPLv3Programming language: Fortran 2003 (tested with gfortran 6.2.0, Intel Fortran 16.0.3, Cray Compiler 8.5.3)External routines/libraries: MPI [1] and Silo [2]Nature of problem: By way of illustrating GenASiSMathematics functionality, solve example fluid dynamics problems.Solution method: Finite-volume discretization; second-order slope-limited reconstruction; HLL Riemann Solver; Runge–Kutta time integration.Additional comments including Restrictions and Unusual features: The example problems named above are not ends in themselves, but serve to illustrate the functionality available though GenASiSMathematics. In addition to these more substantial examples, we provide individual unit test programs for the classes comprised by GenASiSMathematics.GenASiSMathematics is available in the CPC Program Library and also at https://github.com/GenASiS. [1]http://www.mcs.anl.gov/mpi/[2]https://wci.llnl.gov/simulation/computer-codes/silo

OpenMP GNU and Intel Fortran programs for solving the time-dependent Gross–Pitaevskii equation

We present Open Multi-Processing (OpenMP) version of Fortran 90 programs for solving the Gross-Pitaevskii (GP) equation for a Bose-Einstein condensate in one, two, and three spatial dimensions, optimized for use with GNU and Intel compilers. We use the split-step Crank-Nicolson algorithm for imaginary- and real-time propagation, which enables efficient calculation of stationary and non-stationary solutions, respectively. The present OpenMP programs are designed for computers with multi-core processors and optimized for compiling with both commercially-licensed Intel Fortran and popular free open-source GNU Fortran compiler. The programs are easy to use and are elaborated with helpful comments for the users. All input parameters are listed at the beginning of each program. Different output files provide physical quantities such as energy, chemical potential, root-mean-square sizes, densities, etc. We also present speedup test results for new versions of the programs.

Heat-coupled Gaussian continuous-wave double-pass optical parametric oscillator: thermally induced phase mismatching for periodically poled MgO:LiNbO3 crystal.

In this paper, a model describing the thermal effects on the optical parametric oscillator (OPO) of Gaussian continuous waves in double-pass cavities is presented. Eight equations, including forward and backward nonlinear waves, heat, and thermally induced phase mismatching equations, were coupled and solved simultaneously to investigate the effect of heat generation on the OPO's efficiency. The model was applied for a periodically poled MgO:LiNbO3 crystal with an excellent agreement with experimental data. The numerical calculations have been carried out by a homemade code written in Intel FORTRAN, which used a finite difference method.

GenASiS Basics: Object-oriented utilitarian functionality for large-scale physics simulations (Version 2)

GenASiSBasics provides Fortran 2003 classes furnishing extensible object-oriented utilitarian functionality for large-scale physics simulations on distributed memory supercomputers. This functionality includes physical units and constants; display to the screen or standard output device; message passing; I/O to disk; and runtime parameter management and usage statistics. This revision–Version 2 of Basics–makes mostly minor additions to functionality and includes some simplifying name changes. New version program summaryProgram Title: SineWaveAdvection, SawtoothWaveAdvection, and RiemannProblem (fluid dynamics example problems illustrating GenASiSBasics); ArgonEquilibrium and ClusterFormation (molecular dynamics example problems illustrating GenASiSBasics)Program Files doi:http://dx.doi.org/10.17632/6w9ygpygmc.1Licensing provisions: GPLv3Programming language: Fortran 2003 (tested with gfortran 6.1.0, Intel Fortran 16.0.3, Cray Compiler 8.5.3)Journal reference of previous version: Computer Physics Communications, 196 (2015) 506Does the new version supersede the previous version?: YesReasons for the new version: This version makes mostly minor additions to functionality and includes some simplifying name changes.Summary of revisions: Several additions to functionality are minor. Two new singleton objects are KIND_SMALL and KIND_TINY, for smaller sized numbers than those specified by the previously available KIND_DEFAULT. The class MeasuredValueForm can now handle some more complicated cases of unit string processing. The numerical values in the CONSTANT singleton have been updated to 2016 values [3], and CONSTANT and UNIT contain a few additional members.A new class TimerForm can be used to track the wall time occupied by various segments of code. The PROGRAM_HEADER singleton now contains an array member of this new class. With calls like [Display omitted] the user can initialize their own timers; on return iMyTimer contains the index of the newly initialized timer. The calls [Display omitted] and [Display omitted] should surround the block of code to be timed. The information displayed by calling the ShowStatistics method of PROGRAM_HEADER includes data from all initialized timers, including one for overall execution time which is present by default.The code now expects to be compiled with OpenMP, typically by applying compiler flags. Strictly speaking this is only required for the PROGRAM_HEADER singleton, which queries the number of threads via a library call. In GenASiSBasics, OpenMP directives (which appear as comments as far as the Fortran 2003 standard is concerned) are only used in the Clear and Copy commands.There have been a number of name changes, mostly for simplification and consistency. These include the classes in ArrayArrays, where for example ArrayInteger_1D_Form is now simply Integer_1D_Form. Similar streamlining changes have been made to MessagePassing classes: IncomingMessageArrayRealForm is now MessageIncoming_1D_R_Form, for instance. The class VariableGroupArrayMetadata is now VariableGroup_1D_Form. The name ParametersStreamForm has been changed by one character (deletion of an s) to ParameterStreamForm. The member Selected of VariableGroupForm has been changed to iaSelected, where the prefix ia is a conventional prefix we use for an array of array indices.The interface and functionality of the SetGrid member of StructuredGridImageForm have been modified so as not to include boundary cells exterior to the computational domain, which prevented display of the computational domain in 3D plots with VisIt [4] unless a Box operator was applied. See the fluid dynamics examples for the modified usage.Finally, version 4.10 of the Silo library [2] introduced an include file named silo_f9x.inc, which the FileSystem classes of GenASiSBasics now expect to be available instead of silo.inc.Nature of problem: By way of illustrating GenASiSBasics functionality, solve example fluid dynamics and molecular dynamics problems.Solution method: For fluid dynamics examples, finite-volume. For molecular dynamics examples, leapfrog and velocity-Verlet integration.External routines/libraries: MPI [1] and Silo [2]Additional comments including Restrictions and Unusual features:The example problems named above are not ends in themselves, but serve to illustrate our object-oriented approach and the functionality available though GenASiSBasics. In addition to these more substantial examples, we provide individual unit test programs for each of the classes comprised by GenASiSBasics.GenASiSBasics is available in the CPC Program Library and also at https://github.com/GenASiS. [1]http://www.mcs.anl.gov/mpi/[2]https://wci.llnl.gov/simulation/computer-codes/silo[3]C. Patrignani et al. (Particle Data Group), Chin. Phys. C 40 (2016) 100001[4]https://wci.llnl.gov/simulation/computer-codes/visit

FEMSIM + HRMC: Simulation of and structural refinement using fluctuation electron microscopy for amorphous materials

FEMSIM, a Fortran code, is used to simulate the fluctuation electron microscopy signal, the variance, V(k), from a model atomic structure. FEMSIM has been incorporated into a hybrid-reverse Monte Carlo code that combines an embedded atom or Finnis–Sinclair potential with the deviation between simulated and experimental V(k) data to refine an atomic model with structure constrained by both the potential and experimental data. The resulting models have experimentally-derived medium-range order. Program summaryProgram Title : FEMSIM + HRMCProgram Files doi : http://dx.doi.org/10.17632/8dt6jvtjpf.1Licensing provisions : MIT Open SourceProgramming language: Intel FORTRAN 2000Nature of problem: Simulation of fluctuation microscopy experimental data; atomic modeling using empirical potentials and experimental dataSolution method: Monte Carlo and Simulated AnnealingAdditional comments: An Intel Fortran compiler mpif90 (OpenMPI) is recommended

Finite element model for femtosecond laser pulse heating using dual phase lag effect

The ultrashort pulse laser in material processing of micro/nanoscale has emerged as a useful application due to its capability of high precision control on heating and heat affected zone. The non-Fourier effect in heat conduction becomes more predominant as the duration of heating becomes extremely small which is the case of ultrashort pulse laser. To analyze the heat transfer process by an ultrashort pulse laser, a three dimensional finite element based numerical model is developed assuming finite speed of thermal wave propagation with two phase lags. The present work also describes the finite element formulation using phase lag boundary conditions. The numerical model is developed indigenously using an Intel fortran compiler. The transient heating and cooling phenomena are analyzed for the effect of two relaxation times, variation of pulse width, and multiple pulses. The model results are validated with ultrashort pulse laser heating of a thin gold film. The accuracy of the computed results shows promising application of the developed numerical model for differential influence of process parameters in micro/nano scale heat transport by ultrashort pulse laser processing.

OpenMP Fortran and C programs for solving the time-dependent Gross–Pitaevskii equation in an anisotropic trap

We present new version of previously published Fortran and C programs for solving the Gross-Pitaevskii equation for a Bose-Einstein condensate with contact interaction in one, two and three spatial dimensions in imaginary and real time, yielding both stationary and non-stationary solutions. To reduce the execution time on multicore processors, new versions of parallelized programs are developed using Open Multi-Processing (OpenMP) interface. The input in the previous versions of programs was the mathematical quantity nonlinearity for dimensionless form of Gross-Pitaevskii equation, whereas in the present programs the inputs are quantities of experimental interest, such as, number of atoms, scattering length, oscillator length for the trap, etc. New output files for some integrated one- and two-dimensional densities of experimental interest are given. We also present speedup test results for the new programs.

RNGAVXLIB: Program library for random number generation, AVX realization

Abstract The library RNGSSELIB for random number generators (RNGs) based upon the SSE2 command set is presented. The library contains realization of a number of modern and most reliable generators. Usage of SSE2 command set allows to substantially improve performance of the generators. Three new RNG realizations are also constructed. We present detailed analysis of the speed depending on compiler usage and associated optimization level, as well as results of extensive statistical testing for all generators using available test packages. Fast SSE implementations produce exactly the same output sequence as the original algorithms. Program summary Program title: RNGSSELIB Catalogue identifier: AEIT_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEIT_v1_0.html Program obtainable from: CPC Program Library, Queenʼs University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 4177 No. of bytes in distributed program, including test data, etc.: 21 228 Distribution format: tar.gz Programming language: C. Computer: PC. Operating system: UNIX, Windows. RAM: 1 Mbytes Classification: 4.13. Nature of problem: Any calculation requiring uniform pseudorandom number generator, in particular, Monte Carlo calculations. Solution method: The library contains realization of a number of modern and reliable generators: mt19937 , mrg32k3a and lfsr113 . Also new realizations for the method based on parallel evolution of an ensemble of dynamical systems are constructed: GM19 , GM31 and GM61 . The library contains both usual realizations and realizations based on SSE command set. Usage of SSE commands allows the performance of all generators to be substantially improved. Restrictions: For SSE realizations of the generators, Intel or AMD CPU supporting SSE2 command set is required. In order to use the realization lfsr113sse , CPU must support SSE4 command set. Running time: Running time is of the order of 20 sec for generating 109 pseudorandom numbers with a PC based on Intel Core i7-940 CPU. Running time is analysed in detail in Section 5 of the paper.