Solving Systems of Linear Algebraic Equations with the Use of an Selective Regularization

Author(s):  
Vladimir N. Lutay

The solution of systems of linear algebraic equations, which matrices can be poorly conditioned or singular is considered. As a solution method, the original matrix is decomposed into triangular components by Gauss or Chole-sky with an additional operation, which consists in increasing the small or zero diagonal terms of triangular matrices during the decomposition process. In the first case, the scalar products calculated during decomposition are divided into two positive numbers such that the first is greater than the second, and their sum is equal to the original one. In further operations, the first number replaces the scalar product, as a result of which the value of the diagonal term increases, and the second number is stored and used after the decomposition process is completed to correct the result of calculations. This operation increases the diagonal elements of triangular matrices and prevents the appearance of very small numbers in the Gauss method and a negative root expression in the Cholesky method. If the matrix is singular, then the calculated diagonal element is zero, and an arbitrary positive number is added to it. This allows you to complete the decomposition process and calculate the pseudo-inverse matrix using the Greville method. The results of computational experiments are presented.

2021 ◽  
Vol 28 (3) ◽  
pp. 234-237
Author(s):  
Gleb D. Stepanov

This article describes an algorithm for obtaining a non-negative basic solution of a system of linear algebraic equations. This problem, which undoubtedly has an independent interest, in particular, is the most time-consuming part of the famous simplex method for solving linear programming problems.Unlike the artificial basis Orden’s method used in the classical simplex method, the proposed algorithm does not attract artificial variables and economically consumes computational resources.The algorithm consists of two stages, each of which is based on Gaussian exceptions. The first stage coincides with the main part of the Gaussian complete exclusion method, in which the matrix of the system is reduced to the form with an identity submatrix. The second stage is an iterative cycle, at each of the iterations of which, according to some rules, a resolving element is selected, and then a Gaussian elimination step is performed, preserving the matrix structure obtained at the first stage. The cycle ends either when the absence of non-negative solutions is established, or when one of them is found.Two rules for choosing a resolving element are given. The more primitive of them allows for ambiguity of choice and does not exclude looping (but in very rare cases). Use of the second rule ensures that there is no looping.


Author(s):  
Marina L. Mozgaleva ◽  
Pavel A. Akimov ◽  
Taymuraz B. Kaytukov

he distinctive paper is devoted to so-called multigrid (particularly two-grid) method of structural analysis based on discrete Haar basis (one-dimensional, two-dimensional and three-dimensional problems are under consideration). Approximations of the mesh functions in discrete Haar bases of zero and first levels are described (the mesh function is represented as the sum in which one term is its approximation of the first level, and the second term is so-called complement (up to the initial state) on the grid of the first level). Special projectors are constructed for the spaces of vector functions of the original grid to the space of their approximation on the first-level grid and its complement (the refinement component) to the initial state. Basic scheme of the two-grid method is presented. This method allows solution of boundary problems of structural mechanics with the use of matrix operators of significantly smaller dimension. It should be noted that discrete analogue of the initial operator equation is a system of linear algebraic equations which is constructed with the use of finite element method or finite difference method. Block Gauss method can be used for direct solution.


1966 ◽  
Vol 10 (01) ◽  
pp. 25-48
Author(s):  
Richard P. Bernicker

A linearized two-dimensional theory is presented for high-speed hydrofoils near the free surface. The "direct" problem (hydrofoil shape specified) is attacked by replacing the actual foil with vortex and source sheets. The resulting integral equation for the strength of the singularity distribution is recast into an infinite set of linear algebraic equations relating the unknown constants in a Glauert-type vorticity expansion to the boundary condition on the foil. The solution is achieved using a matrix inversion technique and it is found that the matrix relating the known and unknown constants is a function of depth of submergence alone. Inversion of this matrix at each depth allows the vorticity constants to be calculated for any arbitrary foil section by matrix multiplication. The inverted matrices have been calculated for several depth-to-chord ratios and are presented herein. Several examples for specific camber and thickness distributions are given, and results indicate significant effects in the force characteristics at depths less than one chord. In particular, thickness effects cause a loss of lift at shallow submergences which may be an appreciable percentage of the total design lift. The second part treats the "indirect" problem of designing a hydrofoil sectional shape at a given depth to achieve a specified pressure loading. Similar to the "direct" problem treated in the first part, integral equations are derived for the camber and thickness functions by replacing the actual foil by vortex and source sheets. The solution is obtained by recasting these equations into an infinite set of linear algebraic equations relating the constants in a series expansion of the foil geometry to the known pressure boundary conditions. The matrix relating the known and unknown constants is, again, a function of the depth of submergence alone, and inversion techniques allow the sectional shape to be determined for arbitrary design pressure distributions. Several examples indicate the procedure and results are presented for the change in sectional shape for a given pressure loading as the depth of submergence of the foil is decreased.


2020 ◽  
Vol 28 (2) ◽  
pp. 149-159
Author(s):  
Jiří Kopal ◽  
Miroslav Rozložník ◽  
Miroslav Tůma

AbstractThe problem of solving large-scale systems of linear algebraic equations arises in a wide range of applications. In many cases the preconditioned iterative method is a method of choice. This paper deals with the approximate inverse preconditioning AINV/SAINV based on the incomplete generalized Gram–Schmidt process. This type of the approximate inverse preconditioning has been repeatedly used for matrix diagonalization in computation of electronic structures but approximating inverses is of an interest in parallel computations in general. Our approach uses adaptive dropping of the matrix entries with the control based on the computed intermediate quantities. Strategy has been introduced as a way to solve di cult application problems and it is motivated by recent theoretical results on the loss of orthogonality in the generalized Gram– Schmidt process. Nevertheless, there are more aspects of the approach that need to be better understood. The diagonal pivoting based on a rough estimation of condition numbers of leading principal submatrices can sometimes provide inefficient preconditioners. This short study proposes another type of pivoting, namely the pivoting that exploits incremental condition estimation based on monitoring both direct and inverse factors of the approximate factorization. Such pivoting remains rather cheap and it can provide in many cases more reliable preconditioner. Numerical examples from real-world problems, small enough to enable a full analysis, are used to illustrate the potential gains of the new approach.


The problem of diffraction of a vertical electric dipole field on a spiral conductive sphere and a cone has been solved. By the method of regularization of the matrix operator of the problem, an infinite system of linear algebraic equations of the second kind with a compact matrix operator in Hilbert space $\ell_2$ is obtained. Some limiting variants of the problem statement are considered.


Author(s):  
A. Biloshchytskyi ◽  
А. Kuchansky ◽  
Yu. Andrashko ◽  
S. Biloshchytska

A modification of the PageRank method based on link ranking is proposed to evaluate the research results of subjects of the scientific space, taking into account self-citation. The method of reducing the influence of self-citation on the final evaluation of the results of research activity of subjects of the scientific space is described. The evaluation of the results of research is calculated using the modified PR-q method, taking into account self-citation as a solution of a system of linear algebraic equations, matrix of which consists of coefficients determined by the number of citations of publications of one scientist in the publications of another scientist. The described method can be used for the task of evaluating the activity of the components of the scientific space: scientists, higher education institutions and their structural units. For the task of evaluating the research activity of subjects of the scientific space, a method based on link ranking (PageRank method for web pages) and taking into account the self-citation of scientists is proposed. The latter allows for an adequate assessment, taking into account the abuses associated with the authors excessive self-citation. The essence of the constructed method lies in the construction of a system of linear algebraic equations, whose coefficients of the matrix reflect the citations of some scientists by others in the citation network of scientific publications. The value of the coefficients of the matrix of such a system of linear algebraic equations is subject to certain restrictions, which allow to reduce the influence of the factor of excessive self-citation of the author on his overall assessment of research activity. The described method can be used to calculate the complex evaluation of the components of the scientific space: the scientist, the institution of higher education and its separate structural units. Evaluating research results provides an opportunity to verify the relevance of the research process to the goals identified at the planning stage and, if necessary, to adjust the progress of those studies. Also, the calculation of research evaluations of the components (objects and entities) of the scientific space is a powerful tool for managing research projects.


2019 ◽  
pp. 112-115
Author(s):  
M. Z. Benenson

The  article  discusses  the  use  of  graphics  processing  units  for  solving  large  system  of  linear  algebraic  equations  (SLAE).  A heterogeneous multiprocessor computing platform produced by the NIIVK, whose architecture allows the integration of general‑ purpose microprocessor modules with graphics processor modules was used as an equipment for solving SLAEs. The description  of the SLAE solution program, developed on the basis of the CUBLAS CUDA software interface library, is given. A method is proposed for increasing the accuracy of calculations of linear systems based on the use of a modified Gauss method. It has been  established that the use of the modified Gauss method practically does not increase the program operation time with a significant  increase in the accuracy of calculations. It is concluded that the use of graphics processors for solving SLAEs allows processing  matrices of a larger size compared to the use of general‑purpose microprocessors.


Author(s):  
Petr Denisov ◽  
◽  
Anna Balaban ◽  

The article proposes the modification of a technique for assessing the magnetization of permanent magnets from the known field pattern. The identification method is based on solving an ill-conditioned system of linear algebraic equations by the Tikhonov regularization method. The method of boundary integral equations based on scalar potentials is used to compile the matrix of coefficients. The article presents the algorithm that uses parallel computations when performing the most time-consuming operations to reduce the time for solving the inverse problem. In order to check the proposed method, a program was developed that allows to simulate the measurement process: to calculate the direct problem and find the magnetic induction at the points of the air gap, then introduce the error into the "measurement results" and solve the inverse problem. The results of nu-merical experiments that allow us to evaluate the advantages of parallel implementation using the capabilities of modern multi-core processors are presented.


2021 ◽  
Vol 28 (4) ◽  
pp. 434-451
Author(s):  
Gleb D. Stepanov

The article considers a method for solving a linear programming problem (LPP), which requires finding the minimum or maximum of a linear functional on a set of non-negative solutions of a system of linear algebraic equations with the same unknowns. The method is obtained by improving the classical simplex method, which when involving geometric considerations, in fact, generalizes the Gauss complete exclusion method for solving systems of equations. The proposed method, as well as the method of complete exceptions, proceeds from purely algebraic considerations. It consists of converting the entire LPP, including the objective function, into an equivalent problem with an obvious answer. For the convenience of converting the target functional, the equations are written as linear functionals on the left side and zeros on the right one. From the coefficients of the mentioned functionals, a matrix is formed, which is called the LPP matrix. The zero row of the matrix is the coefficients of the target functional, $a_{00}$ is its free member. The algorithms are described and justified in terms of the transformation of this matrix. In calculations the matrix is a calculation table. The method under consideration by analogy with the simplex method consists of three stages. At the first stage the LPP matrix is reduced to a special 1-canonical form. With such matrices one of the basic solutions of the system is obvious, and the target functional on it is $ a_{00}$, which is very convenient. At the second stage the resulting matrix is transformed into a similar matrix with non-positive elements of the zero column (except $a_{00}$), which entails the non-negativity of the basic solution. At the third stage the matrix is transformed into a matrix that provides non-negativity and optimality of the basic solution. For the second stage the analog of which in the simplex method uses an artificial basis and is the most time-consuming, two variants without artificial variables are given. When describing the first of them, along the way, a very easy-to-understand and remember analogue of the famous Farkas lemma is obtained. The other option is quite simple to use, but its full justification is difficult and will be separately published.


Sign in / Sign up

Export Citation Format

Share Document