Implementation in Fpgas of Jacobi Method to Solve the Eigenvalue and Eigenvector Problem

Author(s):  
Ignacio Bravo ◽  
Pedro Jimenez ◽  
Manuel Mazo ◽  
Jose Luis Lazaro ◽  
Alfredo Gardel
Author(s):  
Gaoxin Zhou ◽  
Zhi Gang

In recent years, high order harmonic (or eigenvector) of neutron diffusion equation has been widely used in on-line monitoring system of reactor power. There are two kinds of calculation method to solve the equation: corrected power iteration method and Krylov subspace methods. Fu Li used the corrected power iteration method. When solving for the ith harmonic, it tries to eliminate the influence of the front harmonics using the orthogonality of the harmonic function. But its convergence speed depends on the occupation ratio. When the dominant ratios equal to 1 or close to 1, convergence speed of fixed source iteration method is slow or convergence can’t be achieved. Another method is the Krylov subspace method, the main idea of this method is to project the eigenvalue and eigenvector of large-scale matrix to a small one. Then we can solve the small matrix eigenvalue and eigenvector to get the large ones. In recent years, the restart Arnoldi method emerged as a development of Krylov subspace method. The method uses continuous reboot Arnoldi decomposition, limiting expanding subspace, and the orthogonality of the subspace is guaranteed using orthogonalization method. This paper studied the refined algorithms, a method based on the Krylov subspace method of solving eigenvalue problem for large sparse matrix of neutron diffusion equation. Two improvements have been made for a restarted Arnoldi method. One is that using an ingenious linear combination of the refined Ritz vector forms an initial vector and then generates a new Krylov subspace. Another is that retaining the refined Ritz vector in the new subspace, called, augmented Krylov subspace. This way retains useful information and makes the resulting algorithm converge faster. Several numerical examples are the new algorithm with the implicitly restart Arnoldi algorithm (IRA) and the implicitly restarted refined Arnoldi algorithm (IRRA). Numerical results confirm efficiency of the new algorithm.


2004 ◽  
Vol 19 (11) ◽  
pp. 863-870 ◽  
Author(s):  
S. I. MUSLIH

Multi-Hamiltonian systems are investigated by using the Hamilton–Jacobi method. Integration of a set of total differential equations which includes the equations of motion and the action integral function is discussed. It is shown that this set is integrable if and only if the total variations of the Hamiltonians vanish. Two examples are studied.


2015 ◽  
Vol 2015 ◽  
pp. 1-15 ◽  
Author(s):  
Jun Tao ◽  
Peng Wang ◽  
Haitang Yang

In the framework of the deformed quantum mechanics with a minimal length, we consider the motion of a nonrelativistic particle in a homogeneous external field. We find the integral representation for the physically acceptable wave function in the position representation. Using the method of steepest descent, we obtain the asymptotic expansions of the wave function at large positive and negative arguments. We then employ the leading asymptotic expressions to derive the WKB connection formula, which proceeds from classically forbidden region to classically allowed one through a turning point. By the WKB connection formula, we prove the Bohr-Sommerfeld quantization rule up toOβ2. We also show that if the slope of the potential at a turning point is too steep, the WKB connection formula is no longer valid around the turning point. The effects of the minimal length on the classical motions are investigated using the Hamilton-Jacobi method. We also use the Bohr-Sommerfeld quantization to study statistical physics in deformed spaces with the minimal length.


Sign in / Sign up

Export Citation Format

Share Document