PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Parallel Algorithms”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18Linked to original sources

Phase field model for three-dimensional dendritic growth with fluid flow.

We study the effect of fluid flow on three-dimensional (3D) dendrite growth using a phase-field model on an adaptive finite-element grid. In order to simulate 3D fluid flow, we use an averaging method for the flow problem coupled to the phase-field method and the semi-implicit approximated projection method (SIAPM). We describe a parallel implementation for the algorithm, using the CHARM++ FEM framework, and demonstrate its efficiency. We introduce an improved method for extracting dendrite tip position and tip radius, facilitating accurate comparison to theory. We benchmark our results for 2D dendrite growth with solvability theory and previous results, finding them to be in good agreement. The physics of dendritic growth with fluid flow in three dimensions is very different from that in two dimensions, and we discuss the origin of this behavior.

Algorithms↗

Three-dimensional finite-difference bidomain modeling of homogeneous cardiac tissue on a data-parallel computer.

In this paper, a data-parallel computer is used to provide the memory and reduction in computer time for solving large finite-difference bidomain problems. The finite-difference grid is mapped effectively to the processors of the parallel computer, simply by mapping one node to one (virtual) processor. Implemented on the connection machines (CM's) CM-200 and CM-5, the data-parallel finite-difference algorithm has allowed the solution of finite-difference bidomain problems with over 2 million nodes. Details on the algorithm are presented together with computational performance results.

Algorithms↗

Adaptive imaging using the generalized coherence factor.

Sound-velocity inhomogeneities degrade both spatial and contrast resolutions. This paper proposes a new adaptive imaging technique that uses the generalized coherence factor (GCF) to reduce the focusing errors resulting from the sound-velocity inhomogeneities. The GCF is derived from the spatial spectrum of the received aperture data after proper receive delays have been applied. It is defined as the ratio of the spectral energy within a prespecified low-frequency range to the total energy. It is demonstrated that the low-frequency component of the spectrum corresponds to the coherent portion of the received data, and that the high-frequency component corresponds to the incoherent portion. Hence, the GCF reduces to the coherence factor defined in the literature if the prespecified low-frequency range is restricted to dc only. In addition, the GCF is also an index of the focusing quality and can be used as a weighting factor for the reconstructed image. The efficacy of the GCF technique is demonstrated for focusing errors resulting from the sound-velocity inhomogeneities. Simulations and real ultrasound data are used to evaluate the efficacy of the proposed GCF technique. The characteristics of the GCF, including the effects of the signal-to-noise ratio and the number of channels, are also discussed. The GCF technique also is compared with the correlation-based technique and the parallel adaptive receive compensation algorithm; the improvement in image quality obtained with the proposed technique rivals that of the latter technique. In the presence of a displaced phase screen, this proposed technique also outperforms the correlation-based technique. Computational complexity and implementation issues also are addressed.

Algorithms↗

Spatial transformations for eye-hand coordination.

Eye-hand coordination is complex because it involves the visual guidance of both the eyes and hands, while simultaneously using eye movements to optimize vision. Since only hand motion directly affects the external world, eye movements are the slave in this system. This eye-hand visuomotor system incorporates closed-loop visual feedback but here we focus on early feedforward mechanisms that allow primates to make spatially accurate reaches. First, we consider how the parietal cortex might store and update gaze-centered representations of reach targets during a sequence of gaze shifts and fixations. Recent evidence suggests that such representations might be compared with hand position signals within this early gaze-centered frame. However, the resulting motor error commands cannot be treated independently of their frame of origin or the frame of their destined motor command. Behavioral experiments show that the brain deals with the nonlinear aspects of such reference frame transformations, and incorporates internal models of the complex linkage geometry of the eye-head-shoulder system. These transformations are modeled as a series of vector displacement commands, rotated by eye and head orientation, and implemented between parietal and frontal cortex through efficient parallel neuronal architectures. Finally, we consider how this reach system might interact with the visually guided grasp system through both parallel and coordinated neural algorithms.

Animals↗

Adaptive imaging using an optimal receive aperture size.

Sidelobe contribution from off-axis targets degrades image quality in a coherent array imaging system. In ultrasound imaging, focusing errors resulting from sound-velocity inhomogeneities in human tissue--also known as phase aberrations reduce the coherence of the received signals and elevate the sidelobe level. This paper proposes an adaptive receive-aperture technique based on thresholding of the coherence factor (CF). The CF describes the coherence of the received array signals and can be used as an index of focusing quality. This paper demonstrates that thresholding of the CF allows the mainlobe-dominated signals to be distinguished from the sidelobe-dominated signals, after which the receive-aperture size at each imaging position can be optimally determined so as to enhance the mainlobe-dominated signals and suppress the sidelobe-dominated signals. Thus, image quality degradation resulting from sound-velocity inhomogeneities can be reduced. Simulations and measured ultrasound data are used to evaluate the efficacy of the proposed technique. The characteristics of the proposed technique including the effects of the signal-to-noise ratio (SNR) and the transmit focal depth, and speckle reduction are discussed. The proposed technique is also compared with the parallel adaptive receive compensation algorithm and shown to produce a better improvement in image quality.

Algorithms↗

Adaptive aberration correction of abdominal images using PARCA.

Clinical images of the abdomen may be degraded by aberrations due to propagation through the abdominal wall, resulting in loss of contrast and reduction in overall image quality. The Parallel Adaptive Receive Compensation Algorithm (PARCA) uses multiple receive beams to compensate for these aberrations. Previous studies have demonstrated that PARCA can restore image quality in tissue-equivalent phantoms scanned with controlled phase and amplitude fluctuations. In this study, PARCA is used to improve clinical images of the abdomen. PARCA processed images show both improved contrast of structures and improved overall image quality.

Abdomen↗

Systematic design process for slanted grating couplers.

We have developed a systematic design process for recently proposed slanted grating couplers (SGCs) that operates in the strong coupling regime. Based on rigorous analysis of SGC properties, this design process utilizes the k-vector diagram and a rigorous grating leaky-mode solver to intentionally enforce the phase-match and Bragg conditions. We demonstrate that the resultant SGC designs have performance similar to those obtained by parallel micro-genetic algorithm (muGA) optimization with the two-dimensional finite-difference time-domain (2D FDTD) method. Only two 2D FDTD simulations are necessary in the later stages of our systematic design process. Therefore the time saving is tremendous compared to a muGA 2D FDTD design tool, which can require thousands of individual 2D FDTD simulation runs. We illustrate the utility of our new systematic design process with an embedded slanted grating coupler example.

Journal Article↗

Spectrofluorometric resolution of closely overlapping drug mixtures using chemometrics methods.

A modified parallel factors analysis (PARAFAC) algorithm with a penalty diagonalization error (PDE), newly proposed by the present authors, was utilized to simultaneously resolve drug mixtures of propranolol (PRO), dipyridamole (DIP) and amiloride (AMI) without any loss of sensitivity. The analyses were performed in aqueous solution. The experimental results demonstrated that the profiles of the spectra and the concentrations could be accurately resolved using the PDE algorithm with a high sensitivity and stable repeatability. That is to say, the closely overlapping problem of the spectra could be easily solved. Furthermore, simultaneous determinations of three kinds of tablets, which contain PRO, AMI and DIP, respectively, were successfully performed with satisfactory results.

Algorithms↗

Visualization of intraaneurysmal flow patterns with transluminal flow images of 3D MR angiograms in conjunction with aneurysmal configurations.

BACKGROUND AND PURPOSE: How the complex flow phenomena generated within unruptured cerebral aneurysms relate to the corresponding aneurysmal geometry is unknown. To estimate the interaction between flow patterns and morphologic features of unruptured cerebral aneurysms, we developed a method to visualize intraanuerysmal flow patterns with transluminal flow imaging of 3D MR angiograms in conjunction with aneurysmal configurations. METHODS: Transluminal images of the vessel lumen were reconstructed with use of a parallel volume-rendering algorithm by selecting information on the margin of lumina from the volume data sets of 3D time-of-flight MR angiograms. Transluminal flow images were then created by superimposing flow-related intraluminal information onto transluminal images. Intraaneurysmal flow patterns were evaluated in three cases of unruptured cerebral aneurysms, based on the animated display of transluminal flow images with stepwise extracted intraluminal volume data of signal intensity, in conjunction with the corresponding aneurysmal configurations depicted on 3D MR angiograms. RESULTS: Transluminal flow images showed 3D visualization of flow-related signal intensity distribution obtained from volume data of MR angiograms, so that qualitative information regarding intraaneurysmal flow patterns could be estimated with respect to morphologic features of cerebral aneurysms. CONCLUSION: Transluminal flow images of 3D MR angiograms allowed feasible visualization of intraaneurysmal flow patterns that were studied. More work is required to validate the technique and clarify the significance of being able to visualize intraaneurysmal flow patterns.

Aged↗

[Evaluation of hyperacute cerebral ischemic lesions with transparent 3D visualization of echo-planar diffusion-weighted magnetic resonance imaging].

Transparent three-dimensional (3D) visualization of hyperacute cerebral ischemic lesions was reconstructed using parallel volume-rendering algorithm by selecting the high signal intensity area for ischemic lesion and the low signal intensity area for parenchymal margin, respectively, from the volume dataset of echo-planar diffusion-weighted magnetic resonance imaging (DWI). The 3D DWI can depict the ischemic lesions at the hyperacute stage of cerebral ischemia, transparently through the contours of brain and ventricular structures. Application of 3D DWI and evaluation of ischemic lesions were discussed in three cases of hyperacute cerebral infarction.

Acute Disease↗

A parallel computing approach to genetic sequence comparison: the master-worker paradigm with interworker communication.

We have implemented a parallel version of a dynamic programming biological sequence comparison algorithm to study the potential applicability of using parallel computers for genetic sequence comparisons. Our parallel program is built using C-Linda, a machine-independent parallel programming language, and was tested on both a 10 CPU Sequent Symmetry and a 64 CPU Intel Hypercube. C-Linda implements a shared associative memory model, "tuple space," through which multiple processes can communicate and coordinate control. In our master-worker (MW) parallel implementation, a master process creates several worker processes, extracts a test sequence and multiple library sequences from a database and stores them in tuple space. Each worker reads the test sequence and then repeatedly extracts library strings from tuple space, performs pairwise sequence comparison using a local comparison algorithm to generate a similarity score, and returns the similarity scores to tuple space. The master collects the scores from tuple space and identifies the best match over all library sequences. We also implemented a method of global interworker communication to reduce the total search time by stopping those string comparisons that had no chance of improving on the current best match. Comparisons of the total run time, speedup, and efficiency were made for parallel and sequential versions of a basic MW implementation as well as versions with the global abort threshold.

Algorithms↗

Optimal design of optical reference signals by use of a genetic algorithm.

A new technique for the generation of optical reference signals with optimal properties is presented. In grating measurement systems a reference signal is needed to achieve an absolute measurement of the position. The optical signal is the autocorrelation of two codes with binary transmittance. For a long time, the design of this type of code has required great computational effort, which limits the size of the code to approximately 30 elements. Recently, the application of the dividing rectangles (DIRECT) algorithm has allowed the automatic design of codes up to 100 elements. Because of the binary nature of the problem and the parallel processing of the genetic algorithms, these algorithms are efficient tools for obtaining codes with particular autocorrelation properties. We design optimum zero reference codes with arbitrary length by means of a genetic algorithm enhanced with a restricted search operator.

Algorithms↗

Utilization of cross-plane rays for three-dimensional reconstruction by filtered back-projection.

Present popular computed tomography (CT) algorithms reconstruct an object from the ray measurements lying on a set of parallel planes. This paper presents an algorithm that can also utilize "cross-plane" rays (i.e., rays that cross through many planes) to reconstruct the object. In this reconstruction algorithm, the ray measurements are grouped into two-dimensional projections, filtered, and stored. The filtered projections can then be back-projected onto a three-dimensional matrix or any plane through the three-dimensional volume. General theoretical aspects are presented and then applied to the special case in which ray measurements have been made in all directions. The algorithm is tested using computer-generated data. Expressions for the noise power spectrum and the variance in the reconstruction are derived. It is shown that the noise-to-signal ratio per detected photon for this reconstruction method is close to a theoretical limit, as it also is for normal CT. The ability to use ray measurements that cross many planes is especially useful in emission CT, where order-of-magnitude improvements in image quality per unit dose can be achieved.

Computers↗

Scalable partitioning and exploration of chemical spaces using geometric hashing.

Virtual screening (VS) has become a preferred tool to augment high-throughput screening(1) and determine new leads in the drug discovery process. The core of a VS informatics pipeline includes several data mining algorithms that work on huge databases of chemical compounds containing millions of molecular structures and their associated data. Thus, scaling traditional applications such as classification, partitioning, and outlier detection for huge chemical data sets without a significant loss in accuracy is very important. In this paper, we introduce a data mining framework built on top of a recently developed fast approximate nearest-neighbor-finding algorithm(2) called locality-sensitive hashing (LSH) that can be used to mine huge chemical spaces in a scalable fashion using very modest computational resources. The core LSH algorithm hashes chemical descriptors so that points close to each other in the descriptor space are also close to each other in the hashed space. Using this data structure, one can perform approximate nearest-neighbor searches very quickly, in sublinear time. We validate the accuracy and performance of our framework on three real data sets of sizes ranging from 4337 to 249 071 molecules. Results indicate that the identification of nearest neighbors using the LSH algorithm is at least 2 orders of magnitude faster than the traditional k-nearest-neighbor method and is over 94% accurate for most query parameters. Furthermore, when viewed as a data-partitioning procedure, the LSH algorithm lends itself to easy parallelization of nearest-neighbor classification or regression. We also apply our framework to detect outlying (diverse) compounds in a given chemical space; this algorithm is extremely rapid in determining whether a compound is located in a sparse region of chemical space or not, and it is quite accurate when compared to results obtained using principal-component-analysis-based heuristics.

Journal Article↗

Prediction of RNA base pairing probabilities on massively parallel computers.

We present an implementation of McCaskill's algorithm for computing the base pair probabilities of an RNA molecule for massively parallel message passing architectures. The program can be used to routinely fold RNA sequences of more than 10,000 nucleotides. Applications to complete viral genomes are discussed.

Algorithms↗

Visualization of RNA secondary structures using highly parallel computers.

Results of RNA secondary structure prediction algorithm are usually given as a set of hydrogen bonds between bases. However, we cannot know the precise structure of an RNA molecule by only knowing which bases form hydrogen bonds. One way to understand the structure of an RNA molecule is to visualize it using a planar graph so that we can easily know the geometric relations among the substructures such as stacking regions and loops. To do this, we consider bases to be particles on a plane and introduce a repulsive force and an attractive force among these particles and determine their positions according to these forces. A naive algorithm requires O(N2) time but we can reduce it to O(NlogN) with an approximation algorithm which is often used in the area of N-body simulation. Our program is written in parallel object-oriented language 'Schematic' which is recently developed. Efficiency of our implementation on a parallel computer and results of visualization of secondary structures are presented using cadang-cadang coconut viroid as an example.

Algorithms↗

Unsupervised reconstruction of a three-dimensional left ventricular strain from parallel tagged cardiac images.

A new algorithm, called the Unsupervised Tag ExTraction and Heart strain(E) Reconstruction (UNTETHER) algorithm, is presented for quantifying three-dimensional (3D) myocardial strain in tagged cardiac MR images. Five human volunteers and five postinfarct patients were imaged. 3D strains measured by UNTETHER and a user-supervised technique were compared. Each study was analyzed in 49 +/- 8 min with UNTETHER, compared to approximately 4 hr with the user-supervised technique. For pooled human data, the correlation coefficient between the two methods for circumferential shortening (E(cc)) was r = 0.91 at the mid-wall (P < 0.0005). UNTETHER is capable of measuring wall motion abnormalities resulting from coronary artery disease, and has the potential to overcome the main limitations (time and user-supervision requirements) to routine clinical use of tagged cardiac MRI.

Algorithms↗

Fast space-filling molecular graphics using dynamic partitioning among parallel processors.

We present a novel algorithm for the efficient generation of high-quality space-filling molecular graphics that is particularly appropriate for the creation of the large number of images needed in the animation of molecular dynamics. Each atom of the molecule is represented by a sphere of an appropriate radius, and the image of the sphere is constructed pixel-by-pixel using a generalization of the lighting model proposed by Porter (Comp. Graphics 1978, 12, 282). The edges of the spheres are antialiased, and intersections between spheres are handled through a simple blending algorithm that provides very smooth edges. We have implemented this algorithm on a multiprocessor computer using a procedure that dynamically repartitions the effort among the processors based on the CPU time used by each processor to create the previous image. This dynamic reallocation among processors automatically maximizes efficiency in the face of both the changing nature of the image from frame to frame and the shifting demands of the other programs running simultaneously on the same processors. We present data showing the efficiency of this multiprocessing algorithm as the number of processors is increased. The combination of the graphics and multiprocessor algorithms allows the fast generation of many high-quality images.

Algorithms↗