A lanczos algorithm with a non-direct product basis was used to compute energy levels of H+ (3), H2D+, D2H+, D+ (3), and T+ (3) with J values as large as 46, 53, 66, 66, and 81. The energy levels are based on a modifi...
详细信息
A lanczos algorithm with a non-direct product basis was used to compute energy levels of H+ (3), H2D+, D2H+, D+ (3), and T+ (3) with J values as large as 46, 53, 66, 66, and 81. The energy levels are based on a modified potential surface of M. Pavanello etal. that is better adapted to the ab initio energies near the dissociation limit.
The simulation of electronic circuits involves the numerical solution of very large-scale, sparse, in general nonlinear, systems of differential-algebraic equations. Often, the size of these systems can be reduced con...
详细信息
The simulation of electronic circuits involves the numerical solution of very large-scale, sparse, in general nonlinear, systems of differential-algebraic equations. Often, the size of these systems can be reduced considerably by replacing the equations corresponding to linear subcircuits by approximate models of much smaller state-space dimension. In this paper, we describe the use of Krylov-subspace methods for generating such reduced-order models of linear subcircuits. Particular emphasis is on reduced-order modeling techniques that preserve the passivity of linear RLC subcircuits. (C) 2000 Elsevier Science B.V. All rights reserved.
We show that the time-evolution operator of kicked quantum systems, although a full matrix of size N x N, can be diagonalized with the help of a new method based on a suitable combination of fast Fourier transform and...
详细信息
We show that the time-evolution operator of kicked quantum systems, although a full matrix of size N x N, can be diagonalized with the help of a new method based on a suitable combination of fast Fourier transform and lanczos algorithm in just O(N-2 ln N) operations. It allows the diagonalization of matrices of sizes up to N approximate to 10(6) going far beyond the possibilities of standard diagonalization techniques which need O(N-3) operations. We have applied this method to the kicked Harper model revealing its intricate spectral properties. (C) 1999 Elsevier Science B.V All rights reserved.
A new algorithm for information retrieval is described. It is a vector space method with automatic query expansion. The original user query is projected onto a Krylov subspace generated by the query and the term-docum...
详细信息
A new algorithm for information retrieval is described. It is a vector space method with automatic query expansion. The original user query is projected onto a Krylov subspace generated by the query and the term-document matrix. Each dimension of the Krylov space is generated by a simple vector space search, using first the user query and then new queries generated by the algorithm and orthogonal to the previous query vectors. The new algorithm is closely related to latent semantic indexing (LSI), but it is a local algorithm that works on a new subspace of very low dimension for each query. This makes it faster and more flexible than LSI. No preliminary computation of the singular value decomposition (SVD) is needed, and changes in the data base cause no complication. Numerical tests on both small (Cranfield) and larger (Financial Times data from the TREC collection) data sets are reported. The new algorithm gives better precision at given recall levels than simple vector space and LSI in those cases that have been compared.
This work proposes a new method for RLCG interconnect model-order reductions in consideration with the adjoint network. Relationships between an original MNA network and its corresponding adjoint MNA network will be e...
详细信息
This work proposes a new method for RLCG interconnect model-order reductions in consideration with the adjoint network. Relationships between an original MNA network and its corresponding adjoint MNA network will be explored first. It will be shown that the congruence transformation matrix used in the one-sided projection can be constructed by using the bi-orthogonal bases developed from the Lanezos-type algorithms. In particular, if the multi-port driving-point impedance of RLCG interconnect circuits is the main concern, the transfer functions and system moments of the adjoint network can be directly calculated from those of the original RLCG interconnect network by exploring symmetric properties of the MNA formulation. Therefore, the cost of constructing the congruence transformation matrix can be simplified by up to 50% of the previous methods. Comparative studies among various standard methods and the proposed methods are also investigated. Experimental results on large-scale RLCG interconnect circuits will demonstrate the accuracy and the efficiency of the proposed method.
The Krylov subspace projection approach is a well-established tool for the reduced-order modeling of dynamical systems in the time domain. In this paper, we address the main issues obstructing the application of this ...
详细信息
The Krylov subspace projection approach is a well-established tool for the reduced-order modeling of dynamical systems in the time domain. In this paper, we address the main issues obstructing the application of this powerful approach to the time-domain solution of exterior wave problems. We use frequency-independent perfectly matched layers to simulate the extension to infinity. Pure imaginary stretching functions based on Zolotarev's optimal rational approximation of the square root are implemented leading to perfectly matched layers with a controlled accuracy over a complete spectral interval of interest. A new Krylov-based solution method via stability-corrected operator exponents is presented which allows us to construct reduced-order models (ROMs) that respect the delicate spectral properties of the original scattering problem. The ROMs are unconditionally stable and are based on a renormalized bi-lanczos algorithm. We give a theoretical foundation of our method and illustrate its performance through a number of numerical examples in which we simulate two-dimensional electromagnetic wave propagation in unbounded domains, including a photonic waveguide example. The new algorithm outperforms the conventional finite-difference time-domain method for problems on large time intervals.
For large-scale matrices, there is no practical algorithm to compute the smallest singular value with a satisfied relative accuracy. The widely used bidiagonalization lanczos method can compute the largest singular va...
详细信息
For large-scale matrices, there is no practical algorithm to compute the smallest singular value with a satisfied relative accuracy. The widely used bidiagonalization lanczos method can compute the largest singular value with good relative accuracy, but not the smallest one. In this paper we transform the smallest singular value of matrix A to the largest eigenvalue of (A(T)A)(-1), and use Rayleigh-Ritz method, which is referred as Inverse-Rayleigh-Ritz (IRR) method. The technique computing quadratic Form plays an: important role in IRR. IRR takes no more hop cost and storage than lanczos-like Krylov methods on A and gives more accurate results.
An algorithm to reduce a symmetric matrix to a similar semiseparable one of semiseparability rank 1, using orthogonal similarity transformations, is proposed in this paper. It is shown that, while running to completio...
详细信息
An algorithm to reduce a symmetric matrix to a similar semiseparable one of semiseparability rank 1, using orthogonal similarity transformations, is proposed in this paper. It is shown that, while running to completion, the proposed algorithm gives information on the spectrum of the similar initial matrix. In fact, the proposed algorithm shares the same properties of the lanczos method and the Householder reduction to tridiagonal form. Furthermore, at each iteration, the proposed algorithm performs a step of the QR method without shift to a principal submatrix to retrieve the semiseparable structure. The latter step can be considered a kind of subspace-like iteration method, where the size of the subspace increases by one dimension at each step of the algorithm. Hence, when during the execution of the algorithm the Ritz values approximate the dominant eigenvalues closely enough, diagonal blocks will appear in the semiseparable part where the norm of the corresponding subdiagonal blocks goes to zero in the subsequent iteration steps, depending on the corresponding gap between the eigenvalues. A numerical experiment is included, illustrating the properties of the new algorithm.
We report calculated microwave and infrared rovibrational transitions of the van der Waals complex He-2 - OCS. The calculations were done using a product basis, a lanczos eigensolver, and potentials built from He-OCS,...
详细信息
We report calculated microwave and infrared rovibrational transitions of the van der Waals complex He-2 - OCS. The calculations were done using a product basis, a lanczos eigensolver, and potentials built from He-OCS, and He-He potential functions taken from the literature. All five of the large amplitude coordinates are treated exactly and calculations are done for J values up to five. All rovibrational levels are converged to 0.001 cm(-1) by using basis sets with as many as 87 million funcions. Good agreement is found with previously reported experimental results. Although we assume that the dipole moment is along the OCS axis, we find transitions with appreciable intensity between different torsion states.
We present exact methods for the numerical integration of the Wannier-Stark system in a many-body scenario including two Bloch bands. Our ab initio approaches allow for the treatment of a few-body problem with bosonic...
详细信息
We present exact methods for the numerical integration of the Wannier-Stark system in a many-body scenario including two Bloch bands. Our ab initio approaches allow for the treatment of a few-body problem with bosonic statistics and strong interparticle interaction. The numerical implementation is based on the lanczos algorithm for the diagonalization of large, but sparse symmetric Floquet matrices. We analyze the scheme efficiency in terms of the computational time, which is shown to scale polynomially with the size of the system. The numerically computed eigensystem is applied to the analysis of the Floquet Hamiltonian describing our problem. We show that this allows, for instance, for the efficient detection and characterization of avoided crossings and their statistical analysis. We finally compare the efficiency of our lanczos diagonalization for computing the temporal evolution of our many-body system with an explicit fourth order Runge-Kutta integration. Both implementations heavily exploit efficient matrix vector multiplication schemes. Our results should permit an extrapolation of the applicability of exact methods to increasing sizes of generic many-body quantum problems with bosonic statistics. (C) 2014 Elsevier B.V. All rights reserved.
暂无评论