MODIFIED ENGAGEMENT METHOD FOR MATRIX OPERATION

We describe a modified engagement method for matrix operation based on a two-dimensional crossed-ring interconnection network, Our method incorporates fewer steps than that reported by Bocker et al. [Appl. Opt. 22, 804 (1983)], and its performance is found to be the most efficient (minimum steps) in comparison with other systolic and/or engagement methods for matrix operation. Thus, it may be helpful for other optical and electronic implementations of matrix operations. One compact optoelectronic integrity approach for implementing the modified engagement method is briefly described. (C) 1995 Optical Society of America

Detection of hepatitis B virus markers using a biosensor based on imaging ellipsometry

A biosensor based on imaging ellipsometry (BIE) has been developed and validated in 169 patients for detecting five markers of hepatitis B virus (HBV) infection. The methodology has been established to pave the way for clinical diagnosis, including ligand screening, determination of the sensitivity, set-up of cut-off values (CoVs) and comparison with other clinical methods. A matrix assay method was established for ligand screening. The CoVs of HBV markers were derived with the help of receiver operating characteristic curves. Enzyme-linked immunosorbent assay (ELISA) was the reference method. Ligands with high bioactivity were selected and sensitivities of 1 ng/mL and 1 IU/mL for hepatitis B surface antigen (HBsAg) and surface antibody (anti-HBs) were obtained respectively. The CoVs of HBsAg, anti-HBs, hepatitis B e antigen, hepatitis B e antibody and core antibody were as follows: 15%, 18%, 15%, 20% and 15%, respectively, which were the percentages over the values of corresponding ligand controls. BIE can simultaneously detect up to five markers within 1 h with results in acceptable agreement with ELISA, and thus shows a potential for diagnosing hepatitis B with high throughput.

氧碘激光器相位延迟片的制备与性能研究

多层介质反射镜在非正入射的时候，两个不同的偏振态之间会产生不同的相移。利用矩阵法，根据菲涅耳公式和电磁场边界条件，推导出p，s波的相移。通过优化设计．入射角为54°，在1285～1345nm之间p，s波获得了270°±1°的相移，同时也使反射率在99．5％以上。用离子束溅射技术制备相位延迟膜，用分光光度计测试了光谱特性和用椭偏仪测试了相位特性，在相应波段获得了262．4°±1．8°的相移，同时也使反射率在99．6％以上。误差的主要来源是离子源工作特性会产生不均匀的过渡层和最外层会吸收一些水气、灰尘等也产生

SpMV的自动性能优化实现技术及其应用研究

在科学计算中,稀疏矩阵向量乘(SpMV)是一个十分重要且经常被大量调用的计算内核.由于SpMV一般实现算法的浮点计算和存储访问次数比率非常低,且其存储访问模式极为不规则,其实际运行性能往往很低.通过采用寄存器分块算法和启发式分块大小选择算法,将稀疏矩阵分成小的稠密分块,重用保存在寄存器中向量x元素,可以提高该计算内核的性能.剖析和总结了OSKI软件包所采用的若干关键优化技术,并进行了实际应用性能测试.测试表明,在实际应用这些优化技术的过程中,应用程序对SpMV的调用次数要达到上百次的量级,才能抵消由于应用这些性能优化技术所带来的额外时间开销,取得性能加速效果.在Pentium4和AMD Athlon平台上,测试了10个矩阵,其平均加速比分别达到了1.69和1.48.

RAM(h)模型下SpMV存储访问复杂度的分析

稀疏矩阵向量乘(SpMV)采取压缩行存储格式的算法性能非常差,而寄存器分块算法可以使得数据尽量在靠近处理器的存储层次中访问而提高性能.利用RAM(h)模型进行分析和比较不同算法形式的存储访问复杂度,可以比较两种算法的优劣.通过RAM(h)分析SpMV两种实现形式的存储访问复杂度,同时在奔腾四平台上,测试了7个稀疏矩阵的SpMV性能,并统计了这两种算法中L1,L2,和TLB的缺失率,实验结果与模型分析的数据一致.

OpenMP循环调度算法及SpMV多核并行化研究

OpenMP是一种支持Fortran，C/C++的共享存储并行编程标准。它基于fork-join的并行执行模型，将程序划分为并行区和串行区。近几年来，OpenMP在SMP(Symmetric Multi-Processing)和多核体系结构的并行编程中得到了广泛的应用。随着多核处理器的发展，实际的应用程序如何充分利用多个处理器核来提高运算效率也成为研究的热点。 在科学计算中，循环结构是最核心的并行对象之一。考虑到负载平衡、调度开销、同步开销等多方面因素，OpenMP标准制定了Static调度、Dynamic调度、Guided调度和Runtime调度等不同策略。针对Guided调度策略不适合递减型循环结构的缺点，本文提出了一种改进的new_guided调度策略，并在OMPi编译器上加以实现。New_guided调度策略的主要思想是对前半部分的循环采用Static调度，后半部分的循环采用Guided调度。此外，本文针对不同的循环结构，在多核处理器上对不同的调度策略进行了评测。测试结果表明，在一般情况下，OpenMP默认的Static策略的调度性能最差；对于规则的循环结构和递增的循环结构，Dynamic调度策略、Guided调度策略和new_guided策略的性能差别不大；对于递减型的循环结构，Dynamic调度策略和new_guided策略的性能相当，要优于Guided调度策略；对于求解Mandelbrot集合这类计算量集中在中间的随机循环结构，Dynamic调度策略优于其它策略，new_guided策略的性能介于Dynamic调度和Guided调度之间。 随着多核处理器的问世和发展，多线程程序设计也已经成为一个不可回避的问题。稀疏矩阵向量乘(SpMV, Sparse Matrix-Vector Multiplication)是一个十分重要且经常被大量调用的科学计算内核。SpMV的存储访问一般都极不规则，导致现有的SpMV算法效率都比较低。目前，多核处理器芯片上的内核数量正在逐步增加。这使得在多核处理器上对SpMV进行并行化加速变得非常重要。本文介绍了稀疏矩阵的两种常用的存储格式CSR和BCSR，并采用OpenMP实现了SpMV的多核并行化。此外，本文还讨论了寄存器分块算法、压缩列索引等优化技术，以及不同调度策略对多线程并行后的SpMV的影响。在曙光天阔服务器S4800A1上的测试表明，大部分矩阵都取得了可扩展、甚至是超线性的加速比，但是对于部分规模较大的矩阵，加速效果并不明显。在我们的测试中，与基于CSR实现的多线程SpMV相比，采用寄存器分块算法优化后的SpMV运算速度平均提高了28.09%。在基于CSR实现的多线程SpMV中，采用列索引优化技术后的程序比优化前的速度平均提高了13.05%。此外，本文实现了一种基于非零元个数的调度策略。在该策略中，每个线程处理几乎相同数量的非零元。我们将它和OpenMP标准提供的三种调度策略进行了测试和分析。测试结果表明：与OpenMP提供的调度策略相比，基于非零元个数的调度策略能取得更好的负载平衡；Dynamic调度和Guided调度在多线程SpMV中的性能基本相当，均优于Static调度策略。

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

Performance of an embedded optical vector matrix multiplication processor architecture

An embedded architecture of optical vector matrix multiplier (OVMM) is presented. The embedded architecture is aimed at optimising the data flow of vector matrix multiplier (VMM) to promote its performance. Data dependence is discussed when the OVMM is connected to a cluster system. A simulator is built to analyse the performance according to the architecture. According to the simulation, Amdahl's law is used to analyse the hybrid opto-electronic system. It is found that the electronic part and its interaction with optical part form the bottleneck of system.

Transfer of paralytic shellfish toxins via marine food chains: A simulated experiment

Objective To study the transfer of paralytic shellfish toxins (PST) using four simulated marine food chains: dinoflagellate Alexandrium tamarense -> Arterriia Artemia salina -> Mysid shrimp Neomysis awatschensis; A. tamarense-N. awatschensis: A. taniarense A. salina -> Perch Lateolabrax japonicus; and A. tamarense -> L. japonicus. Methods The ingestion of A. tamarense, a producer of PST, by L. japonicus, N. awatschensis, and A. salina was first confirmed by microscopic observation of A. tamarense cells in the intestine samples of the three different organisms, and by the analysis of Chl.a levels iii the samples. Toxin accumulation in L. japonicus and N. awatschensis directly from the feeding on A. tamarense or indirectly ibrough the vector of A. salina was then studied. The toxicity of samples was measured using the AOAC mouse bioassay method, and the toxin content and profile of A. tamarense were analyzed by the HPLC method. Results Both A. salina and N. awatschensis could ingest A. tamarense cells. However, the ingestion capability of A. salina exceeded that of N. awatschensis. After the exposure to the culture of A. tamarense (2 000 cells(.)mL(-1)) for 70 minutes, the content of ChLa in A. salina and N. awatschensis reached 0.87 and 0.024 mu g-mg(-1), respectively. Besides, A. tamarense cells existed in the intestines of L. japonicus, N. awatschensis and A. salina by microscopic observation. Therefore, the three organisms could ingest A. tamarense cells directly. A. salina could accumulate high content of PST, and the toxicity of A. salina in samples collected on days 1, 4, and 5 of the experiment was 2.18, 2.6, and 2.1 MU(.)g(-1), respectively. All extracts from the samples could lead to death of tested mice within 7 minutes, and the toxin content in arternia sample collected on the 1st day was estimated to be 1.65x10(-5) pg STX equa Vindividual. Toxin accumulation in L. japonicus and N. awatschensis directly from the feeding on A. tamarense or indirectly froin the vector of A. salina was also studied. The mice injected with extracts from L. japonicus and N. awatschensis samples that accumulated PST either directly or indirectly showed PST intoxication symptoms, indicating that low levels of PST existed in these samples. Conclusion Paralytic shellfish toxins can be transferred to L. japonicus, N. awatschensis, and A. salina from A. taniarense directly or indirectly via the food chains.

High-Strength Composite Fibers: Realizing True Potential of Carbon Nanotubes in Polymer Matrix through Continuous Reticulate Architecture and Molecular Level Couplings

Carbon nanotubes have unprecedented mechanical properties as defect-free nanoscale building blocks, but their potential has not been fully realized in composite materials due to weakness at the interfaces. Here we demonstrate that through load-transfer-favored three-dimensional architecture and molecular level couplings with polymer chains, true potential of CNTs can be realized in composites as Initially envisioned. Composite fibers with reticulate nanotube architectures show order of magnitude improvement in strength compared to randomly dispersed short CNT reinforced composites reported before. The molecular level couplings between nanotubes and polymer chains results in drastic differences in the properties of thermoset and thermoplastic composite fibers, which indicate that conventional macroscopic composite theory falls to explain the overall hybrid behavior at nanoscale.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

Determination of 266 pesticide residues in apple juice by matrix solid-phase dispersion and gas chromatography-mass selective detection

A macro matrix solid-phase dispersion (MSPD) method was developed to extract 266 pesticides from apple juice samples prior to gas chromatography-mass selective detection (GC-MSD) determination. A 10 g samples was mixed with 20 g diatomaceous earth. The mixture was transferred into a glass column. Pesticide residues were leached with a 160 mL hexane-dichloromethane (1:1) at 5 mL/min. Two hundred and sixty-six pesticides were divided into three groups and detected by GC-MSD under selective ion monitoring. The proposed method takes advantage of both liquid-liquid extraction and conventional MSPD methods. Application was illustrated by the analysis of 236 apple juice samples produced in Shaanxi province China mainland this year. (C) 2004 Elsevier B.V. All rights reserved.

DTADH and quantum critical phenomena caused by anisotropy and external magnetic field for spin-1/2 Heisenberg diamond chains

Using the transfer matrix renormalization group (TMRG) method, we study the connection between the first derivative of the thermal average of driving-term Hamiltonian (DTADH) and the trace of quantum critical behaviors at finite temperatures. Connecting with the exact diagonalization method, we give the phase diagrams and analyze the properties of each phase for both the ferromagnetic and anti-ferromagnetic frustrated J(3) anisotropy diamond chain models. The finite-temperature scaling behaviors near the critical regions are also investigated. Further, we show the critical behaviors driven by external magnetic field, analyze the formation of the 1/3 magnetic plateau and the influence of different interactions on those critical points for both the ferrimagnetic and anti-ferromagnetic distorted diamond chains.