In situ ray tracing and computational steering for interactive blood flow simulation

Recent algorithm and hardware developments have significantly improved our capability to interactively visualise time-varying flow fields. However, when visualising very large dynamically varying datasets interactively there are still limitations in the scalability and efficiency of these methods. Here we present a rendering pipeline which employs an efficient in situ ray tracing technique to visualise flow fields as they are simulated. The ray casting approach is particularly well suited for the visualisation of large and sparse time-varying datasets, where it is capable of rendering fluid flow fields at high image resolutions and at interactive frame rates on a single multi-core processor using OpenMP. The parallel implementation of our in situ visualisation method relies on MPI, requires no specialised hardware support, and employs the same underlying spatial decomposition as the fluid simulator. The visualisation pipeline allows the user to operate on a commodity computer and explore the simulation output interactively. Our simulation environment incorporates numerous features that can be utilised in a wide variety of research contexts.     

Library of sophisticated functions for analysis of nuclear spectra

In the paper we present compact library for analysis of nuclear spectra. The library consists of sophisticated functions for background elimination, smoothing, peak searching, deconvolution, and peak fitting. The functions can process one- and two-dimensional spectra. The software described in the paper comprises a number of conventional as well as newly developed methods needed to analyze experimental data.

Program summary

Program title: SpecAnalysLib 1.1Catalogue identifier: AEDZ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDZ_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 154No. of bytes in distributed program, including test data, etc.: 2 379 437Distribution format: tar.gzProgramming language: C++Computer: Pentium 3 PC 2.4 GHz or higher, Borland C++ Builder v. 6. A precompiled Windows version is included in the distribution packageOperating system: Windows 32 bit versionsRAM: 10 MBWord size: 32 bitsClassification: 17.6Nature of problem: The demand for advanced highly effective experimental data analysis functions is enormous. The library package represents one approach to give the physicists the possibility to use the advanced routines simply by calling them from their own programs. SpecAnalysLib is a collection of functions for analysis of one- and two-parameter γ-ray spectra, but they can be used for other types of data as well. The library consists of sophisticated functions for background elimination, smoothing, peak searching, deconvolution, and peak fitting.Solution method: The algorithms of background estimation are based on Sensitive Non-linear Iterative Peak (SNIP) clipping algorithm. The smoothing algorithms are based on the convolution of the original data with several types of filters and algorithms based on discrete Markov chains. The peak searching algorithms use the smoothed second differences and they can search for peaks of general form. The deconvolution (decomposition - unfolding) functions use the Gold iterative algorithm, its improved high resolution version and Richardson-Lucy algorithm. In the algorithms of peak fitting we have implemented two approaches. The first one is based on the algorithm without matrix inversion - AWMI algorithm. It allows it to fit large blocks of data and large number of parameters. The other one is based on the calculation of the system of linear equations using Stiefel-Hestens method. It converges faster than the AWMI, however it is not suitable for fitting large number of parameters.Restrictions: Dimensionality of the analyzed data is limited to two.Unusual features: Dynamically loadable library (DLL) of processing functions users can call from their own programs.Running time: Most processing routines execute interactively or in a few seconds. Computationally intensive routines (deconvolution, fitting) execute longer, depending on the number of iterations specified and volume of the processed data.     

ALTDSE: An Arnoldi–Lanczos program to solve the time-dependent Schrödinger equation

We describe a general ab initio and non-perturbative method to solve the time-dependent Schrödinger equation (TDSE) for the interaction of a strong attosecond laser pulse with a general atom. While the field-free Hamiltonian and the dipole matrices may be generated using an arbitrary primitive basis, they are assumed to have been transformed to the eigenbasis of the problem before the solution of the TDSE is propagated in time using the Arnoldi–Lanczos method. Probabilities for survival of the ground state, excitation, and single ionization can be extracted from the propagated wavefunction.

Program summary

Program title: ALTDSECatalogue identifier: AEDM_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDM_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.: 2154No. of bytes in distributed program, including test data, etc.: 30 827Distribution format: tar.gzProgramming language: Fortran 95. [A Fortran 2003 call to “flush” is used to simplify monitoring the output file during execution. If this function is not available, these statements should be commented out.].Computer: Shared-memory machinesOperating system: Linux, OpenMPHas the code been vectorized or parallelized?: YesRAM: Several Gb, depending on matrix size and number of processorsSupplementary material: To facilitate the execution of the program, Hamiltonian field-free and dipole matrix files are provided.Classification: 2.5External routines: LAPACK, BLASNature of problem: We describe a computer program for a general ab initio and non-perturbative method to solve the time-dependent Schrödinger equation (TDSE) for the interaction of a strong attosecond laser pulse with a general atom [1,2]. The probabilities for survival of the initial state, excitation of discrete states, and single ionization due to multi-photon processes can be obtained.Solution method: The solution of the TDSE is propagated in time using the Arnoldi–Lanczos method. The field-free Hamiltonian and the dipole matrices, originally generated in an arbitrary basis (e.g., the flexible B-spline R-matrix (BSR) method with non-orthogonal orbitals [3]), must be provided in the eigenbasis of the problem as input.Restrictions: The present program is restricted to a ¹S^e initial state and linearly polarized light. This is the most common situation experimentally, but a generalization is straightforward.Running time: Several hours, depending on the number of threads used.References: [1] X. Guan, O. Zatsarinny, K. Bartschat, B.I. Schneider, J. Feist, C.J. Noble, Phys. Rev. A 76 (2007) 053411. [2] X. Guan, C.J. Noble, O. Zatsarinny, K. Bartschat, B.I. Schneider, Phys. Rev. A 78 (2008) 053402. [3] O. Zatsarinny, Comput. Phys. Comm. 174 (2006) 273.     

LED光谱特性及其在中间视觉下的应用

为了研究中间视觉下LED的光度特性,运用LED光色电综合分析仪测量了不同电流下五种颜色多个LED的光谱,然后计算了不同中间视觉亮度下的光通量,并分析和研究了中间视觉修正系数Rm随电流和亮度的变化。研究表明:电流在12～24mA范围内,同一颜色LED的相对光谱整体变化小于8%,因此在普通照明要求下,可以用平均光谱来描述LED的相对光谱。在中间视觉下,绿、白这两种颜色LED的中间视觉修正系数Rm的对数与亮度的对数存在线性关系,红、黄、蓝LED的Rm的对数与亮度的对数存在二次曲线关系。其中,蓝、绿、白LED的Rm随亮度增大而减小,红、黄LED的Rm随亮度增大而增大。该研究为中间视觉下LED灯具的设计提供了理论基础。     

Study of low work function materials for hot cavity resonance ionization laser ion sources

The selectivity of a hot cavity resonance ionization laser ion source (RILIS) is most often limited by contributions from competing surface ionization of the hot walls of the ionization cavity. In this article we present investigations on the properties of designated high temperature, low work function materials regarding their performance and suitability as cavity material for RILIS. Tungsten test cavities, impregnated with a mixture of barium oxide and strontium oxide (BaOSrO on W), or alternatively gadolinium hexaboride (GdB₆) were studied in comparison to a standard tungsten RILIS cavity as being routinely used for hot cavity laser ionization at ISOLDE. Measurement campaigns took place at the off-line mass separators at ISOLDE/CERN, Geneva and RISIKO/University of Mainz.     

New measurement and evaluation of the excitation function of Ni(p,n) reaction for the production of Cu

One of the radioisotopes for which a growing interest exists in nuclear medicine is ⁶⁴Cu. Its branched decay makes it suitable for both diagnostic and therapeutic purposes. Activation cross sections of the proton induced reaction on enriched ⁶⁴Ni have been studied using the stacked foil technique up to 24 MeV. The experimental cross sections are compared with values available from literature. Thick target yields, based on the discrete measured values of the cross sections are calculated and allow a better estimation of the optimum production parameters.     

Interaction of intense short laser pulses with positronium

The ionization of the positronium is studied in the over-the-barrier regime. The Volkov and momentum-space strong-field approximation are applied to describe the ionization of the hydrogen-like atomic systems. Classical trajectory Monte Carlo calculations were also performed and a good agreement with quantum-mechanical calculations were found. The scaling of the ionization probability density as function of the reduced mass of the target is analyzed. We find that the maxima of the ionization probability density is determined by laser field parameters, and it’s width can be characterised by the charge of the core and reduced mass of the target.     

Development of a high intensity Al beam from SiC targets for accelerated beam experiments at ISAC

An intense beam of ^26gAl has been developed for accelerated beam experiments at TRIUMF’s ISAC facility. Studies of the on-line production of Al radionuclides from thick silicon carbide targets have been performed as part of a program of beam development for astrophysical reaction studies at ISAC. While the release of short-lived Al nuclides from SiC was found to be slow, development of new target material forms and high-power target containers has allowed operation of SiC targets with proton currents of up to 70 μA on target. In addition, operation with the TRIUMF resonant ionization laser ion source (TRILIS) has produced ^26gAl beam intensities of 5.1 × 10¹⁰ s⁻¹.     

Parallel dynamic and kinetic regular triangulation in three dimensions

A parallel algorithm for regular triangulations is presented. For the purpose of fully dynamic and kinetic particle simulations it allows vertex insertion, deletion, movement, and weight changes. We describe new algorithms for incremental construction of regular triangulations, parallel vertex deletion and insertion. Finally, a parallel Lawson flip algorithm for vertex displacements is presented. The performance analysis demonstrates a significant parallel efficiency for various system sizes and performed changes.