Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDecomposition method for block-tridiagonal matrix systems
86
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The decomposition method which makes the parallel solution of the block-tridiagonal matrix systems possible is presented. The performance of the method is analytically estimated based on the number of elementary multiplicative operations for its parallel and serial parts. The computational speedup with respect to the conventional sequential Thomas algorithm is assessed for various types of the application of the method. It is observed that the maximum of the analytical speedup for a given number of blocks on the diagonal is achieved at some finite number of parallel processors. The values of the parameters required to reach the maximum computational speedup are obtained. The benchmark calculations show a good agreement of analytical estimations of the computational speedup and practically achieved results. The application of the method is illustrated by employing the decomposition method to the matrix system originated from a boundary value problem for the two-dimensional integro-differential Faddeev equations. The block-tridiagonal structure of the matrix arises from the proper discretization scheme including the finite-differences over the first coordinate and spline approximation over the second one. The application of the decomposition method for parallelization of solving the matrix system reduces the overall time of calculation up to 10 times.
rate research
Read More
We introduce a numerical method and python package, https://github.com/andillio/CHiMES, that simulates quantum systems initially well approximated by mean field theory using a second order extension of the classical field approach. We call this the field moment expansion method. In this way, we can accurately approximate the evolution of first and second field moments beyond where the mean field theory breaks down. This allows us to estimate the quantum breaktime of a classical approximation without any calculations external to the theory. We investigate the accuracy of the field moment expansion using a number of well studied quantum test problems. Interacting Bosonic systems similar to scalar field dark matter are chosen as test problems. We find that successful application of this method depends on two conditions: the quantum system must initially be well described by the classical theory, and that the growth of the higher order moments be hierarchical.
In this paper, two types of Schur complement based preconditioners are studied for twofold and block tridiagonal saddle point problems. One is based on the nested (or recursive) Schur complement, the other is based on an additive type Schur complement after permuting the original saddle point systems. We discuss different preconditioners incorporating the exact Schur complements. It is shown that some of them will lead to positive stable preconditioned systems. Our theoretical analysis is instructive for devising various exact and inexact preconditioners, as well as iterative solvers for many twofold and block tridiagonal saddle point problems.
Reachability analysis aims at identifying states reachable by a system within a given time horizon. This task is known to be computationally expensive for linear hybrid systems. Reachability analysis works by iteratively applying continuous and discrete post operators to compute states reachable according to continuous and discrete dynamics, respectively. In this paper, we enhance both of these operators and make sure that most of the involved computations are performed in low-dimensional state space. In particular, we improve the continuous-post operator by performing computations in high-dimensional state space only for time intervals relevant for the subsequent application of the discrete-post operator. Furthermore, the new discrete-post operator performs low-dimensional computations by leveraging the structure of the guard and assignment of a considered transition. We illustrate the potential of our approach on a number of challenging benchmarks.
An square matrix is $k$-Toeplitz if its diagonals are periodic sequences of period $k$. We find rational formulas for the determinant, the characteristic polynomial, and the elements of the inverse of a tridiagonal $k$-Toeplitz matrix (in particular, of any tridiagonal matrix) over any commutative unital ring, using only elementary linear algebra.
The displacement field for three dimensional dynamic elasticity problems in the frequency domain can be decomposed into a sum of a longitudinal and a transversal part known as a Helmholtz decomposition. The Cartesian components of both the longitudinal and transverse fields satisfy scalar Helmholtz equations that can be solved using a desingularized boundary element method (BEM) framework. The curl free longitudinal and divergence free transversal conditions can also be cast as additional scalar Helmholtz equations. When compared to other BEM implementations, the current framework leads to smaller matrix dimensions and a simpler conceptual approach. The numerical implementation of this approach is benchmarked against the 3D elastic wave field generated by a rigid vibrating sphere embedded in an infinite linear elastic medium for which the analytical solution has been derived. Examples of focussed 3D elastic waves generated by a vibrating bowl-shaped rigid object with convex and concave surfaces are also considered. In the static zero frequency limit, the Helmholtz decomposition becomes non-unique, and both the longitudinal and transverse components contain divergent terms that are proportional to the inverse square of the frequency. However, these divergences are equal and opposite so that their sum, that is the displacement field that reflects the physics of the problem, remains finite in the zero frequency limit.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
