共查询到20条相似文献,搜索用时 31 毫秒
1.
M.T. Darvishi F. Khani S. Hamedi-Nezhad B. Zheng 《Journal of Computational and Applied Mathematics》2008
In this article, we develop symmetric block successive overrelaxation (S-block-SOR) methods for finding the solution of the rank-deficient least squares problems. We propose an S2-block-SOR and an S3-block-SOR method for solving such problems and the convergence of these two methods is studied. The comparisons between the S2-block and the S3-block methods are presented with some numerical examples. 相似文献
2.
Tommy Elfving 《Numerische Mathematik》1980,35(1):1-12
Summary We shall in this paper consider the problem of computing a generalized solution of a given linear system of equations. The matrix will be partitioned by blocks of rows or blocks of columns. The generalized inverses of the blocks are then used as data to Jacobi- and SOR-types of iterative schemes. It is shown that the methods based on partitioning by rows converge towards the minimum norm solution of a consistent linear system. The column methods converge towards a least squares solution of a given system. For the case with two blocks explicit expressions for the optimal values of the iteration parameters are obtained. Finally an application is given to the linear system that arises from reconstruction of a two-dimensional object by its one-dimensional projections. 相似文献
3.
Efficient subroutines for dense matrix computations have recently been developed and are available on many high-speed computers.
On some computers the speed of many dense matrix operations is near to the peak-performance. For sparse matrices storage and
operations can be saved by operating only and storing only nonzero elements. However, the price is a great degradation of
the speed of computations on supercomputers (due to the use of indirect addresses, to the need to insert new nonzeros in the
sparse storage scheme, to the lack of data locality, etc.).
On many high-speed computers a dense matrix technique is preferable to sparse matrix technique when the matrices are not large,
because the high computational speed compensates fully the disadvantages of using more arithmetic operations and more storage.
For very large matrices the computations must be organized as a sequence of tasks in each of which a dense block is treated.
The blocks must be large enough to achieve a high computational speed, but not too large, because this will lead to a large
increase in both the computing time and the storage. A special “locally optimized reordering algorithm” (LORA) is described,
which reorders the matrix so that dense blocks can be constructed and treated with some standard software, say LAPACK or NAG.
These ideas are implemented for linear least-squares problems. The rectangular matrices (that appear in such problems) are
decomposed by an orthogonal method. Results obtained on a CRAY C92A computer demonstrate the efficiency of using large dense
blocks. 相似文献
4.
Summary Recently, special attention has been given in the literature, to the problems of accurately computing the least-squares solution of very largescale over-determined systems of linear equations which occur in geodetic applications. In particular, it has been suggested that one can solve such problems iteratively by applying the block SOR (Successive Overrelaxation) and EGS1 (Extrapolated Gauss Seidel 1) plus semi-iterative methods to a linear system with coefficient matrix 2-cyclic or 3-cyclic. The comparison of 2-block SOR and 3-block SOR was made in [1] and showed that the 2-block SOR is better. In [6], the authors also proved that 3-block EGS1-SI is better than 3-block SOR. Here, we first show that the 2-block DJ (Double Jacobi)-SI, GS-SI and EGS1-SI methods are equivalent and all of them are equivalent to the 3-block EGS1-SI method; then, we prove that the composite methods and 2-block SOR have the same asymptotic rate of convergence, but the former has a better average rate of convergence; finally, numerical experiments are reported, and confirm that the 3-block EGS1-SI is better than the 2-block SOR. 相似文献
5.
In this paper a new iterative method is given for solving large sparse least squares problems and computing the minimum norm
solution to underdetermined consistent linear systems. The new scheme is called the generalized successive overrelaxation
(GSOR) method and is shown to be convergent ifA is full column rank. The GSOR method involves a parameter ρ and an auxiliary matrixP. One can choose matrix P so that the GSOR method only involves matrix and vector operations; therefore the GSOR method is
suitable for parallel computations. Besides, the GSOR method can be combined with preconditioning techniques, and therefore
can be expected to be more effective.
This author's work was supported by Natural Science Foundation of Liaoning Province, China. 相似文献
6.
In this paper, we consider a class of Uzawa-SOR methods for saddle point problems, and prove the convergence of the proposed methods. We solve a lower triangular system per iteration in the proposed methods, instead of solving a linear equation Az=b. Actually, the new methods can be considered as an inexact iteration method with the Uzawa as the outer iteration and the SOR as the inner iteration. Although the proposed methods cannot achieve the same convergence rate as the GSOR methods proposed by Bai et al. [Z.-Z. Bai, B.N. Parlett, Z.-Q. Wang, On generalized successive overrelaxation methods for augmented linear systems, Numer. Math. 102 (2005) 1-38], but our proposed methods have less workloads per iteration step. Experimental results show that our proposed methods are feasible and effective. 相似文献
7.
Summary We study the augmented system approach for the solution of sparse linear least-squares problems. It is well known that this method has better numerical properties than the method based on the normal equations. We use recent work by Arioli et al. (1988) to introduce error bounds and estimates for the components of the solution of the augmented system. In particular, we find that, using iterative refinement, we obtain a very robust algorithm and our estimates of the error are accurate and cheap to compute. The final error and all our error estimates are much better than the classical or Skeel's error analysis (1979) indicates. Moreover, we prove that our error estimates are independent of the row scaling of the augmented system and we analyze the influence of the Björck scaling (1967) on these estimates. We illustrate this with runs both on large-scale practical problems and contrived examples, comparing the numerical behaviour of the augmented systems approach with a code using the normal equations. These experiments show that while the augmented system approach with iterative refinement can sometimes be less efficient than the normal equations approach, it is comparable or better when the least-squares matrix has a full row, and is, in any case, much more stable and robust.This author was visiting Harwell and was funded by a grant from the Italian National Council of Research (CNR), Istituto di Elaborazione dell'Informazione-CNR, via S. Maria 46, I-56100 Pisa, ItalyThis author was visiting Harwell from Faculty of Mathematics and Computer Science of the University of Amsterdam 相似文献
8.
Richard E. Ewing 《BIT Numerical Mathematics》1989,29(4):850-866
The simulation of large-scale fluid flow applications often requires the efficient solution of extremely large nonsymmetric linear and nonlinear sparse systems of equations arising from the discretization of systems of partial differential equations. While preconditioned conjugate gradient methods work well for symmetric, positive-definite matrices, other methods are necessary to treat large, nonsymmetric matrices. The applications may also involve highly localized phenomena which can be addressed via local and adaptive grid refinement techniques. These local refinement methods usually cause non-standard grid connections which destroy the bandedness of the matrices and the associated ease of solution and vectorization of the algorithms. The use of preconditioned conjugate gradient or conjugate-gradient-like iterative methods in large-scale reservoir simulation applications is briefly surveyed. Then, some block preconditioning methods for adaptive grid refinement via domain decomposition techniques are presented and compared. These techniques are being used efficiently in existing large-scale simulation codes. 相似文献
9.
In the finite element method, a standard approach to mesh tying is to apply Lagrange multipliers. If the interface is curved, however, discretization generally leads to adjoining surfaces that do not coincide spatially. Straightforward Lagrange multiplier methods lead to discrete formulations failing a first-order patch test [T.A. Laursen, M.W. Heinstein, Consistent mesh-tying methods for topologically distinct discretized surfaces in non-linear solid mechanics, Internat. J. Numer. Methods Eng. 57 (2003) 1197–1242]. 相似文献
10.
Xiaoxia Zhou Yongzhong SongLi Wang Qingsheng Liu 《Journal of Computational and Applied Mathematics》2009
In this paper, we present the preconditioned generalized accelerated overrelaxation (GAOR) method for solving linear systems based on a class of weighted linear least square problems. Two kinds of preconditioning are proposed, and each one contains three preconditioners. We compare the spectral radii of the iteration matrices of the preconditioned and the original methods. The comparison results show that the convergence rate of the preconditioned GAOR methods is indeed better than the rate of the original method, whenever the original method is convergent. Finally, a numerical example is presented in order to confirm these theoretical results. 相似文献
11.
Summary.
In this paper we introduce a class of robust multilevel
interface solvers for two-dimensional
finite element discrete elliptic problems with highly
varying coefficients corresponding to geometric decompositions by a
tensor product of strongly non-uniform meshes.
The global iterations convergence rate is shown to be of
the order
with respect to the number of degrees
of freedom on the single subdomain boundaries, uniformly upon the
coarse and fine mesh sizes, jumps in the coefficients
and aspect ratios of substructures.
As the first approach, we adapt the frequency filtering techniques
[28] to construct robust smoothers
on the highly non-uniform coarse grid. As an alternative, a multilevel
averaging procedure for successive coarse grid correction is
proposed and analyzed.
The resultant multilevel coarse grid
preconditioner is shown to have (in a two level case) the condition
number independent
of the coarse mesh grading and
jumps in the coefficients related to the coarsest refinement level.
The proposed technique exhibited high serial and parallel
performance in the skin diffusion processes modelling [20]
where the high dimensional coarse mesh problem inherits a strong geometrical
and coefficients anisotropy.
The approach may be also applied to magnetostatics problems
as well as in some composite materials simulation.
Received December 27, 1994 相似文献
12.
In this paper we study and compare some preconditioned conjugate gradient methods for solving large-scale higher-order finite element schemes approximating two- and three-dimensional linear elasticity boundary value problems. The preconditioners discussed in this paper are derived from hierarchical splitting of the finite element space first proposed by O. Axelsson and I. Gustafsson. We especially focus our attention to the implicit construction of preconditioning operators by means of some fixpoint iteration process including multigrid techniques. Many numerical experiments confirm the efficiency of these preconditioners in comparison with classical direct methods most frequently used in practice up to now. 相似文献
13.
Standard Galerkin finite element methods or finite difference methods for singular perturbation problems lead to strongly unsymmetric matrices, which furthermore are in general notM-matrices. Accordingly, preconditioned iterative methods such as preconditioned (generalized) conjugate gradient methods, which have turned out to be very successful for symmetric and positive definite problems, can fail to converge or require an excessive number of iterations for singular perturbation problems.This is not so much due to the asymmetry, as it is to the fact that the spectrum can have both eigenvalues with positive and negative real parts, or eigenvalues with arbitrary small positive real parts and nonnegligible imaginary parts. This will be the case for a standard Galerkin method, unless the meshparameterh is chosen excessively small. There exist other discretization methods, however, for which the corresponding bilinear form is coercive, whence its finite element matrix has only eigenvalues with positive real parts; in fact, the real parts are positive uniformly in the singular perturbation parameter.In the present paper we examine the streamline diffusion finite element method in this respect. It is found that incomplete block-matrix factorization methods, both on classical form and on an inverse-free (vectorizable) form, coupled with a general least squares conjugate gradient method, can work exceptionally well on this type of problem. The number of iterations is sometimes significantly smaller than for the corresponding almost symmetric problem where the velocity field is close to zero or the singular perturbation parameter =1.The 2
nd
author's research was sponsored by Control Data Corporation through its PACER fellowship program.The 3
rd
author's research was supported by the Netherlands organization for scientific research (NWO).On leave from the Institute of Mathematics, Academy of Science, 1090 Sofia, P.O. Box 373, Bulgaria. 相似文献
14.
E. J. van Asselt 《BIT Numerical Mathematics》1985,25(2):380-385
Globally convergent nonlinear relaxation methods are considered for a class of nonlinear boundary value problems (BVPs), where the discretizations are continuousM-functions.It is shown that the equations with one variable occurring in the nonlinear relaxation methods can always be solved by Newton's method combined with the Bisection method. The nonlinear relaxation methods are used to get an initial approximation in the domain of attraction of Newton's method. Numerical examples are given. 相似文献
15.
We survey multilevel iterative methods applied for solving large sparse systems with matrices, which depend on a level parameter, such as arise by the discretization of boundary value problems for partial differential equations when successive refinements of an initial discretization mesh is used to construct a sequence of nested difference or finite element meshes.We discuss various two-level (two-grid) preconditioning techniques, including some for nonsymmetric problems. The generalization of these techniques to the multilevel case is a nontrivial task. We emphasize several ways this can be done including classical multigrid methods and a recently proposed algebraic multilevel preconditioning method. Conditions for which the methods have an optimal order of computational complexity are presented.On leave from the Institute of Mathematics, and Center for Informatics and Computer Technology, Bulgarian Academy of Sciences, Sofia, Bulgaria. The research of the second author reported here was partly supported by the Stichting Mathematisch Centrum, Amsterdam. 相似文献
16.
In this paper, we propose two variants of the additive Schwarz method for the approximation of second order elliptic boundary
value problems with discontinuous coefficients, on nonmatching grids using the lowest order Crouzeix-Raviart element for the
discretization in each subdomain. The overall discretization is based on the mortar technique for coupling nonmatching grids.
The convergence behavior of the proposed methods is similar to that of their closely related methods for conforming elements.
The condition number bound for the preconditioned systems is independent of the jumps of the coefficient, and depend linearly
on the ratio between the subdomain size and the mesh size. The performance of the methods is illustrated by some numerical
results.
This work has been supported by the Alexander von Humboldt Foundation and the special funds for major state basic research
projects (973) under 2005CB321701 and the National Science Foundation (NSF) of China (No.10471144)
This work has been supported in part by the Bergen Center for Computational Science, University of Bergen 相似文献
17.
Summary For solving second order elliptic problems discretized on a sequence of nested mixed finite element spaces nearly optimal iterative methods are proposed. The methods are within the general framework of the product (multiplicative) scheme for operators in a Hilbert space, proposed recently by Bramble, Pasciak, Wang, and Xu [5,6,26,27] and make use of certain multilevel decomposition of the corresponding spaces for the flux variable. 相似文献
18.
Galerkin-wavelet methods for two-point boundary value problems 总被引:7,自引:0,他引:7
Summary Anti-derivatives of wavelets are used for the numerical solution of differential equations. Optimal error estimates are obtained in the applications to two-point boundary value problems of second order. The orthogonal property of the wavelets is used to construct efficient iterative methods for the solution of the resultant linear algebraic systems. Numerical examples are given.This work was supported by National Science Foundation 相似文献
19.
Summary Algebraic multilevel analogues of the BEPS preconditioner designed for solving discrete elliptic problems on grids with local refinement are formulated, and bounds on their relative condition numbers, with respect to the composite-grid matrix, are derived. TheV-cycle and, more generally,v-foldV-cycle multilevel BEPS preconditioners are presented and studied. It is proved that for 2-D problems theV-cycle multilevel BEPS is almost optimal, whereas thev-foldV-cycle algebraic multilevel BEPS is optimal under a mild restriction on the composite cell-centered grid. For thev-fold multilevel BEPS, the variational relation between the finite difference matrix and the corresponding matrix on the next-coarser level is not necessarily required. Since they are purely algebraically derived, thev-fold (v>1) multilevel BEPS preconditioners perform without any restrictionson the shape of subregions, unless the refinement is too fast. For theV-cycle BEPS preconditioner (2-D problem), a variational relation between the matrices on two consecutive grids is required, but there is no restriction on the method of refinement on the shape, or on the size of the subdomains. 相似文献
20.
Summary We consider the numerical solution of indefinite systems of linear equations arising in the calculation of saddle points. We are mainly concerned with sparse systems of this type resulting from certain discretizations of partial differential equations. We present an iterative method involving two levels of iteration, similar in some respects to the Uzawa algorithm. We relate the rates of convergence of the outer and inner iterations, proving that, under natural hypotheses, the outer iteration achieves the rate of convergence of the inner iteration. The technique is applied to finite element approximations of the Stokes equations.The work of this author was supported by the Office of Naval Research under contract N00014-82K-0197, by Avions Marcel Dassault, 78 Quai Marcel Dassault, 92214 St Cloud, France, and by Direction des Recherches Etudes et Techniques, 26 boulevard Victor, F-75996 Paris Armées, FranceThe work of this author was supported by Avions Marcel Daussault-Breguet Aviation, 78 quai Marcel Daussault, F-92214 St Cloud, France and by Direction des Recherches Etudes et Techniques, 26 boulevard Victor, F-75996 Paris Armées, FranceThe work of this author was supported by Konrad-Zuse-Zentrum für Informationstechnik Berlin, Federal Republic of Germany 相似文献