Converting source code from one programming language to another is a problem that occurs regularly in real life, but has attracted limited attention and has not been investigated systematically. This paper presents th...
详细信息
Building k-nearest neighbor(kNN) graphs is a necessary step in such areas as data mining and machine learning. So in this paper, we attempt to study the kNN furthermore, we first propose a parallel algorithm for appro...
详细信息
Building k-nearest neighbor(kNN) graphs is a necessary step in such areas as data mining and machine learning. So in this paper, we attempt to study the kNN furthermore, we first propose a parallel algorithm for approximate kNN graph construction and then apply the kNN graph to the application of clustering. Experiments show that our MPI/OpenMP mixed mode codes can make the construction of approximate kNN graph faster and make the parallelization and implementation easier. Finally, we compare the results of agglomerative clustering methods by using our parallel algorithm to illustrate the applicability of this method.
Traditional photography focuses on the optimization of lenses for a perfect imaging ***, with the great developments of computational resources and optical modulation devices, we can achieve more powerful imaging abil...
详细信息
Traditional photography focuses on the optimization of lenses for a perfect imaging ***, with the great developments of computational resources and optical modulation devices, we can achieve more powerful imaging abilities with concise *** photography is such an emerging interdisciplinary field by incorporating computational strategy in traditional imaging system to break the limitations in various dimensions such as spectrum, time and space. Recent advances in different aspects have aroused great interests and introduced tremendous applications in biology, material science and computer vision.
Modeling of complex physical systems with Modelica usually leads to the high-index differential algebraic equation system(DAE), index reduction is an important part of solving the high-index DAE. The structure index r...
详细信息
Modeling of complex physical systems with Modelica usually leads to the high-index differential algebraic equation system(DAE), index reduction is an important part of solving the high-index DAE. The structure index reduction algorithm is one of the popular methods, but in special cases, it fails. Combinatorial relaxation algorithm can detect and correct the breakdown situation. And the maximum weight matching of bipartite graph is an important part of the combinatorial relaxation algorithm. In order to choose the proper method for the large-scale, dense bipartite graph, this paper provides three implementations of the Hungarian algorithm. The experiment results and the theory show that the BFS single-augmented method is better than others.
This paper proposes an improved KM algorithm to computing the structural index of linear time-invariant Differential Algebraic Equation(DAE) systems. The problem is of practical significance in index reduction based o...
详细信息
This paper proposes an improved KM algorithm to computing the structural index of linear time-invariant Differential Algebraic Equation(DAE) systems. The problem is of practical significance in index reduction based on structural index of DAE system and combinatorial relaxation theory. This improved KM algorithm combines greedy idea and classical KM algorithm. It first computes matches as much as possible using greedy technology, and then call KM algorithm to search the matches for the unmatched vertices during the step of greedy technology. The improved KM algorithm reduces the running time bound by a factor of, the number of matches searched using greedy algorithm. Generally, the time complexity is, the optimal time is.
In this paper, we consider second order elliptic ODE eigenproblems on general grids. We construct an efficient algorithm for computing the eigenvalue by using weighted mean combination of the linear finite element met...
详细信息
In this paper, we consider second order elliptic ODE eigenproblems on general grids. We construct an efficient algorithm for computing the eigenvalue by using weighted mean combination of the linear finite element method and corresponding 2nd-order finite difference *** first take the arithmetic mean of the two methods. Then we compute the quasi-optimal combined parameters for different eigenvalues to improve our efficient algorithm. The algorithm we construct convergence faster and have higher accuracy than the linear finite element method and corresponding 2nd-order finite difference method. Some numerical examples tested on both uniform meshes and nonuniform meshes are given to illustrate the computational cost of different numerical methods for solving eigenvalue *** efficiency, all the matrices use sparse storage in our algorithm.
Under investigation in this paper is the Whitham-Broer-Kaup (WBK) system, which describes the dispersive long wave in shallow water. Through a variable transformation, the WBK system is casted into a general Broer-Kau...
详细信息
Under investigation in this paper is the Whitham-Broer-Kaup (WBK) system, which describes the dispersive long wave in shallow water. Through a variable transformation, the WBK system is casted into a general Broer-Kaup system whose Lax pair can be derived by the Ablowitz-Kaup-Newell-Segur technology. With symbolic computation, based on the aforementioned Lax pair, the N-fold Darboux transformation is constructed with a gauge transformation and the multi-soliton solutions are obtained. Finally, the elastic interactions of the two-soliton solutions (including the head-on and overtaking collisions) for the WBK system are graphically studied. Those multi-soliton collisions can beused to illustrate the bidirectional propagation of the waves in shallow water.
Many fields of science and engineering, such as astronomy, medical imaging, seismology and spectroscopy, have been revolutionized by Fourier methods. The fast Fourier transform (FFT) is an efficient algorithm to compu...
详细信息
Many fields of science and engineering, such as astronomy, medical imaging, seismology and spectroscopy, have been revolutionized by Fourier methods. The fast Fourier transform (FFT) is an efficient algorithm to compute the discrete Fourier transform (DFT) and its inverse. The emerging class of high performance computing architectures, such as GPU, seeks to achieve much higher performance and efficiency by exposing a hierarchy of distinct memories to programmers. However, the complexity of GPU programming poses a significant challenge for programmers. In this paper, based on the Kronecker product form multi-dimensional FFTs, we propose an automatic performance tuning framework for various OpenCL GPUs. Several key techniques of GPU programming on AMD and NVIDIA GPUs are also identified. Our OpenCL FFT library achieves up to 1.5 to 4 times, 1.5 to 40 times and 1.4 times the performance of clAmdFft 1.0 for 1D, 2D and 3D FFT respectively on an AMD GPU, and the overall performance is within 90% of CUFFT 4.0 on two NVIDIA GPUs.
This paper is to investigate a variable-coefficient modified Kortweg-de Vries (vc-mKdV) model, which describes some situations from fluid mechanics, ocean dynamics, and plasma mechanics. By the AblowRz-Kaup-NewellSe...
详细信息
This paper is to investigate a variable-coefficient modified Kortweg-de Vries (vc-mKdV) model, which describes some situations from fluid mechanics, ocean dynamics, and plasma mechanics. By the AblowRz-Kaup-NewellSegur procedure and symbolic computation, the Lax pair of the vc-MKdV model is derived. Then, based on the aforementioned Lax pair, the Darboux transformation is constructed and a new one-soliton-like solution is obtained as weft Features of the one-soliton-like solution are analyzed and graphically discussed to illustrate the influence of the variable coefficients in the solitonlike propagation.
Register pressure problem has been a known problem for compiler because of the mismatch between the infinite number of pseudo registers and the finite number of hard registers. Too heavy register pressure may results ...
详细信息
暂无评论