On the distribution of errors in the iterative solution of a system of linear equations. II

1975 ◽  
Vol 25 (4) ◽  
pp. 461-463 ◽  
Author(s):  
Tetsuro Yamamoto
2020 ◽  
Author(s):  
Mikhail Kruglyakov ◽  
Alexey Kuvshinov

<p> In this contribution, we present novel global 3-D electromagnetic forward solver based on a numerical solution of integral equation (IE) with contracting kernel. Compared to widely used x3dg code which is also based on IE approach, new solver exploits alternative (more efficient and accurate) numerical algorithms to calculate Green’s tensors, as well as an alternative (Galerkin) method to construct the system of linear equations (SLE). The latter provides guaranteed convergence of the iterative solution of SLE. The solver outperforms x3dg in terms of accuracy, and, in contrast to (sequential) x3dg, it allows for efficient parallel computations, meaning that the code has practically linear scalability up to the hundreds of processors.</p>


Author(s):  
A. C. Aitken

SynopsisThe convergence of customary processes of iteration for solving linear equations, in particular simple and Seidelian iteration, is studied from the standpoint of matrices. A new variant of Seidelian iteration is introduced. In the positive definite case it always converges, the characteristic roots of its operator being real and positive and less than unity.


Author(s):  
David Ek ◽  
Anders Forsgren

AbstractThe focus in this paper is interior-point methods for bound-constrained nonlinear optimization, where the system of nonlinear equations that arise are solved with Newton’s method. There is a trade-off between solving Newton systems directly, which give high quality solutions, and solving many approximate Newton systems which are computationally less expensive but give lower quality solutions. We propose partial and full approximate solutions to the Newton systems. The specific approximate solution depends on estimates of the active and inactive constraints at the solution. These sets are at each iteration estimated by basic heuristics. The partial approximate solutions are computationally inexpensive, whereas a system of linear equations needs to be solved for the full approximate solution. The size of the system is determined by the estimate of the inactive constraints at the solution. In addition, we motivate and suggest two Newton-like approaches which are based on an intermediate step that consists of the partial approximate solutions. The theoretical setting is introduced and asymptotic error bounds are given. We also give numerical results to investigate the performance of the approximate solutions within and beyond the theoretical framework.


2019 ◽  
Vol 2019 (1) ◽  
Author(s):  
A. Khalid ◽  
M. N. Naeem ◽  
P. Agarwal ◽  
A. Ghaffar ◽  
Z. Ullah ◽  
...  

AbstractIn the current paper, authors proposed a computational model based on the cubic B-spline method to solve linear 6th order BVPs arising in astrophysics. The prescribed method transforms the boundary problem to a system of linear equations. The algorithm we are going to develop in this paper is not only simply the approximation solution of the 6th order BVPs using cubic B-spline, but it also describes the estimated derivatives of 1st order to 6th order of the analytic solution at the same time. This novel technique has lesser computational cost than numerous other techniques and is second order convergent. To show the efficiency of the proposed method, four numerical examples have been tested. The results are described using error tables and graphs and are compared with the results existing in the literature.


Sign in / Sign up

Export Citation Format

Share Document