A parallel radix‐4 block cyclic reduction algorithm

A conventional block cyclic reduction algorithm operates by halving the size of the linear system at each reduction step, that is, the algorithm is a radix‐2 method. An algorithm analogous to the block cyclic reduction known as the radix‐q partial solution variant of the cyclic reduction (PSCR) method allows the use of higher radix numbers and is thus more suitable for parallel architectures as it requires fever reduction steps. This paper presents an alternative and more intuitive way of deriving a radix‐4 block cyclic reduction method for systems with a coefficient matrix of the form tridiag{ ? I,D, ? I}. This is performed by modifying an existing radix‐2 block cyclic reduction method. The resulting algorithm is then parallelized by using the partial fraction technique. The parallel variant is demonstrated to be less computationally expensive when compared to the radix‐2 block cyclic reduction method in the sense that the total number of emerging subproblems is reduced. The method is also shown to be numerically stable and equivalent to the radix‐4 PSCR method. The numerical results archived correspond to the theoretical expectations. Copyright © 2013 John Wiley & Sons, Ltd.     

基于Beta函数及其偏导数的广义积分的高精度快速计算

考虑Beta函数偏导数的计算以及与此相关的广义积分的高精度快速计算问题.首先将Beta函数B(x,y)的定义扩展到整个复平面上,并建立了在整个复平面上Beta函数B(x,y)的偏导数的递推公式.对许多广义积分我们给出Beta函数偏导数的表示形式,因而利用Beta函数的偏导数计算这些广义积分.数值计算表明,算法无论从计算精度还是计算速度,远好于数值积分.另外,得到了B_(p,q)(x,y)存在闭形式的条件,并给出一些广义积分的闭形式.     

Fast multidimensional localized parallel NMR spectroscopy for the analysis of samples

A parallel localized spectroscopy (PALSY) method is presented to speed up the acquisition of multidimensional NMR (nD) spectra. The sample is virtually divided into a discrete number of nonoverlapping slices that relax independently during consecutive scans of the experiment, affording a substantial reduction in the interscan relaxation delay and the total experiment time. PALSY was tested for the acquisition of three experiments 2D COSY, 2D DQF‐COSY and 2D TQF‐COSY in parallel, affording a time‐saving factor of 3–4. Some unique advantages are that the achievable resolution in any dimension is not compromised in any way: it uses conventional NMR data processing, it is not prone to generate spectral artifacts, and once calibrated, it can be used routinely with these and other combinations of NMR spectra. Copyright © 2010 John Wiley & Sons, Ltd.     

北京大学高平均流强光阴极电子枪驱动激光系统

DC-SC超导光阴极注入器将直流光阴极电子枪和超导加速腔结合在一起,其特点之一是适合于提供高平均流强电子束.文中介绍了直流光阴极电子枪的驱动激光系统的研制.该驱动激光系统的主机为瑞士TBP公司的GE-100-XHP型激光器,是采用可饱和吸收镜技术的被动锁模激光器,其工作晶体为Nd:YVO₄,输出波长为1064nm,脉冲重复频率为81.25MHz,平均功率10W.该驱动激光系统配套的CLX-1100同步控制器可以将激光脉冲与射频参考信号严格锁定,随机抖动不大于1ps.根据使用要求,作者研制了高效的激光倍频器,获得了大于1W的266nm紫外激光.同时,激光到光阴极的传输设计中,采用了傅立叶光学技术,在紫外激光功率650mW左右的条件下,从自行研制的Cs₂Te光阴极上实现了大于500μA的光电流引出     

Optimal Release Times in Single-Stage Manufacturing Systems with Blocking: Optimal Control Perspective

This paper concerns a due-date matching problem in a single-stage manufacturing system. Given a finite sequence of jobs and their service order, and given the delivery due date of each job, the problem is to choose the jobs release (arrival) times so as to match as closely as possible their completion times to their respective due dates. The system is modelled as a deterministic single-server FIFO queue with an output buffer for storing jobs whose service is completed prior to their due dates. The output buffer has a finite capacity; when it is full, the server is being blocked. Associated with each job there is a convex cost function penalizing its earliness as well as tardiness. The due-date matching problem is cast as an optimal control problem, whose objective is to minimize the sum of the above cost functions by the choice of the jobs arrival (release) times. Time-box upper-bound and lower-bound constraints are imposed on the jobs output (delivery) times. The optimal-control setting brings to bear on the development of fast and efficient algorithms having intuitive geometric appeal and potential for online implementation.Communicated by W. B. GongResearch supported in part by the National Science Foundation under Grant ECS-9979693 and by the Georgia Tech Manufacturing Research Center under Grant B01-D06.     

柠檬醛对黄曲霉质膜物理损伤机制的研究

通过反射干涉对比显微技术和快速图像分析技术,捕获细胞质膜外光晕变化和细胞周边边缘激发变化的信息,并采用闪烁本征模式分解测定柠檬醛作用于黄曲霉质膜前后各模式（m=0,1,2）所对应的平均平方振幅<δn2m>、剪切弹性μ和弹性弯曲模量Kc,与理论模式光谱相比,研究柠檬醛损伤质膜的物理机制.结果表明,m=1的平动模式主导了光晕波谱,其平均平方振幅呈双指数衰变,该模式几乎与剪切弹性无关,即μ=0;弹性弯曲模量Kc的最小和最大值分别为9.45×10-18J和5.43×10-17J,与柠檬醛的浓度成正相关性.揭示该醛对黄曲霉质膜的损伤,除改变膜的选择通透性和电导率外,同时使质膜平动模式的平均平方振幅呈双指数衰变,弹性弯曲模量Kc增加,从而降低膜伸缩的弹性.     

三相快速流化床中固体颗粒循环流量的研究

采用快速循环流化床而使粒子、液体处于循环状态操作，粒子不会外溢堆积，可使床层在相对稳定状态下进行工作，从而克服普通三相流化床返混大，转化率低的缺点。     

Correction for Shift-variant Characteristics of Single Photon Emission Computed Tomography Acquisition System Using Sensitivity Functions

The single photon emission computed tomography (SPECT) imaging system has shift-variant characteristics due to non-uniform attenuation of gamma-ray, collimator design, scattered photons, etc. In order to provide quantitatively accurate SPECT images, these characteristics should be compensated in the reconstruction. This paper presents a new method to correct the shift-variant characteristics, which is based on a continuous-discrete mapping model and filtered backprojection (FBP) method, in which the projection data are assumed to be acquired by narrow ray sum beams in the FBP method and the assumed data set is expressed as a linear combination of the actual projection data. Narrow ray sum beams are approximated by a weighted sum of the original sensitivity functions. Thus, at the reconstruction the projection data are first modified using an approximation and the FBP method is then applied to the corrected projection data and a SPECT image is reconstructed. We further propose a technique that requires the inversion of smaller matrices than the conventional algebraic method; the amount of calculation and memory space become smaller and the stability of the calculation is greatly improved as well. The results of the numerical simulations are also demonstrated.     

Fast direct solver for Poisson equation in a 2D elliptical domain

In this article, we extend our previous work 3 for developing some fast Poisson solvers on 2D polar and spherical geometries to an elliptical domain. Instead of solving the equation in an irregular Cartesian geometry, we formulate the equation in elliptical coordinates. The solver relies on representing the solution as a truncated Fourier series, then solving the differential equations of Fourier coefficients by finite difference discretizations. Using a grid by shifting half mesh away from the pole and incorporating the derived numerical boundary value, the difficulty of coordinate singularity can be elevated easily. Unlike the case of 2D disk domain, the present difference equation for each Fourier mode is coupled with its conjugate mode through the numerical boundary value near the pole; thus, those two modes are solved simultaneously. Both second‐ and fourth‐order accurate schemes for Dirichlet and Neumann problems are presented. In particular, the fourth‐order accuracy can be achieved by a three‐point compact stencil which is in contrast to a five‐point long stencil for the disk case. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 72–81, 2004     

Determining the geometry of hydrogen bonds in solids with picometer accuracy by quantum-chemical calculations and NMR spectroscopy.

The structure of multiply hydrogen-bonded systems is determined with picometer accuracy by a combined solid-state NMR and quantum-chemical approach. On the experimental side, advanced 1H-15N dipolar recoupling NMR techniques are capable of providing proton-nitrogen distances of up to about 250 pm with an accuracy level of +/-1 pm for short distances (i.e., around 100 pm) and +/-5 pm for longer ones (i.e., 180 to 250 pm). The experiments were performed under fast magic-angle spinning, which ensures sufficient dipolar decoupling and spectral resolution of the 1H resonance lines. On the quantum-chemical side, the structures of the hydrogen-bonded systems were computationally optimised, yielding complete sets of nitrogen-proton and proton-proton distances, which are essential for correctly interpreting the experimental NMR data. In this way, nitrogen-proton distances were determined with picometer accuracy, so that vibrational averaging effects on dipole-dipole couplings need to be considered. The obtained structures were finally confirmed by the complete agreement of computed and experimental 'H and '5N chemical shifts. This demonstrates that solid-state NMR and quantum-chemical methods ideally complement each other and, in a combined manner, represent a powerful approach for reliable, high-precision structure determination whenever scattering techniques are inapplicable.