ALBERT

All Library Books, journals and Electronic Records Telegrafenberg

Your email was sent successfully. Check your inbox.

An error occurred while sending the email. Please try again.

Proceed reservation?

Export
Filter
  • Articles  (42)
  • Articles: DFG German National Licenses  (42)
  • CR: G1.3  (42)
  • 2015-2019
  • 1985-1989  (42)
  • 1945-1949
  • Mathematics  (42)
  • Architecture, Civil Engineering, Surveying
Collection
  • Articles  (42)
Source
  • Articles: DFG German National Licenses  (42)
Publisher
Years
Year
Topic
  • Mathematics  (42)
  • Architecture, Civil Engineering, Surveying
  • 1
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 50 (1986), S. 263-271 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We study semi-iterative (two step interative) methods of the form $$x_{n + 1} = T^* T(\alpha x_n + \gamma x_{n - 1} ) + \beta x_n + (1 - \beta )x_{n - 1} - (\alpha + \gamma )T^* y$$ for the approximate solution of ill-posed or ill-conditioned linear equationsTx=y in (infinite or finite dimensional) Hilbert spaces. We present results on convergence, convergence rates, the influence of perturbed data, and on the comparison of different methods.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 2
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 52 (1987), S. 241-250 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 15A12 ; 49D15 ; 65F35 ; CR: G1.3 ; G1.6
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We derive lower bounds for the ∞-condition number of then×n-Vandermonde matrixV n(x) in the cases where the node vectorx T=[x1, x2,...,xn] has positive elements or real elements located symmetrically with respect to the origin. The bounds obtained grow exponentially inn. withO(2n) andO(2n/2), respectively. We also compute the optimal spectral condition numbers ofV n(x) for the two node configurations (including the optimal nodes) and compare them with the bounds obtained.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 3
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 54 (1989), S. 591-599 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05, 65G05, 15A06 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary In this paper a Gauss-Jordan algorithm with column interchanges is presented and analysed. We show that, in contrast with Gaussian elimination, the Gauss-Jordan algorithm has essentially differing properties when using column interchanges instead of row interchanges for improving the numerical stability. For solutions obtained by Gauss-Jordan with column interchanges, a more satisfactory bound for the residual norm can be given. The analysis gives theoretical evidence that the algorithm yields numerical solutions as good as those obtained by Gaussian elimination and that, in most practical situations, the residuals are equally small. This is confirmed by numerical experiments. Moreover, timing experiments on a Cyber 205 vector computer show that the algorithm presented has good vectorisation properties.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 4
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 55 (1989), S. 265-280 
    ISSN: 0945-3245
    Keywords: AMS(MOS) ; 65H10 ; 58C99 ; 55M25 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The homotopy method can be used to solve eigenvalue-eigenvector problems. The purpose of this paper is to report the numerical experience of the homotopy method of computing eigenpairs for real symmetric tridiagonal matrices together with a couple of new theoretical results. In practice, it is rerely of any interest to compute all the eigenvalues. The homotopy method, having the order preserving property, can provide any specific eigenvalue without calculating any other eigenvalues. Besides this advantage, we note that the homotopy algorithm is to a large degree a parallel algorithm. Numerical experimentation shows that the homotopy method can be very efficient especially for graded matrices.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 5
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 52 (1987), S. 279-300 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary In this paper we compare several implementations of Kogbetliantz's algorithm for computing the SVD on sequential as well as on parallel machines. Comparisons are based on timings and on operation counts. The numerical accuracy of the different methods is also analyzed.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 6
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 53 (1988), S. 571-593 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The Chebyshev and second-order Richardson methods are classical iterative schemes for solving linear systems. We consider the convergence analysis of these methods when each step of the iteration is carried out inexactly. This has many applications, since a preconditioned iteration requires, at each step, the solution of a linear system which may be solved inexactly using an “inner” iteration. We derive an error bound which applies to the general nonsymmetric inexact Chebyshev iteration. We show how this simplifies slightly in the case of a symmetric or skew-symmetric iteration, and we consider both the cases of underestimating and overestimating the spectrum. We show that in the symmetric case, it is actually advantageous to underestimate the spectrum when the spectral radius and the degree of inexactness are both large. This is not true in the case of the skew-symmetric iteration. We show how similar results apply to the Richardson iteration. Finally, we describe numerical experiments which illustrate the results and suggest that the Chebyshev and Richardson methods, with reasonable parameter choices, may be more effective than the conjugate gradient method in the presence of inexactness.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 7
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 54 (1988), S. 19-32 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F20 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary A direct method is developed for solving linear least squares problems $$\mathop {\min \left\| {Ax - b} \right\|_2 }\limits_x $$ , whereA is large and sparse and the solution is subject to lower and upper boundsl≦x≦u. The problem is initially transformed to upper triangular form by a sparseQR-factorization. An active set algorithm is then used. The key step is the stable updating of theR-factor associated with the columns ofA corresponding to the free variables, when theQ-factor is not available. For this a new method is developed, which uses the semi-normal equations and iterative refinement.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 8
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 627-633 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We propose a “one-sided” or “implicit” variant of the Jacobi diagonalization algorithm for positive definite matrices. The variant is based on a previous Cholesky decomposition and currently uses essentially one square array which, on output, contains the matrix of eigenvectors thus reaching the storage economy of the symmetric QL algorithm. The current array is accessed only columnwise which makes the algorithm attractive for various parallelized and/or vectorized implementations. Even on a serial computer our algorithm shows improved efficiency, in particular if the Cholesky step is made with diagonal pivoting. On matrices of ordern=25–200 our algorithm is about 2–3.5 times slower than QL thus being almost on the halfway between the standard Jacobi and QL algorithms. The previous Cholesky decomposition can be performed with higher precision without extra time or storage costs thus offering considerable gains in accuracy with highly conditioned input matrices.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 9
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 721-734 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The acceleration by Tchebychev iteration for solving nonsymmetric eigenvalue problems is dicussed. A simple algorithm is derived to obtain the optimal ellipse which passes through two eigenvalues in a complex plane relative to a reference complex eigenvalue. New criteria are established to identify the optimal ellipse of the eigenspectrum. The algorithm is fast, reliable and does not require a search for all possible ellipses which enclose the spectrum. The procedure is applicable to nonsymmetric linear systems as well.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 10
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 54 (1989), S. 639-654 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 35P15, 47A55, 65M15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary LetA, B be essentially self-adjoint and positive definite differential operators defined inL 2(G). Using Svirskij's construction of the base operator and some results from the analytic perturbation theory of linear operators a formula providing eigenvalue lower bounds of the problemAu=λBu is derived. In this formula a rough lower bound of some higher eigenvalue and the residual convergence of the Rayleigh-Ritz eigenfunction approximations are needed. Some numerical results are presented.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 11
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 55 (1989), S. 463-476 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05, 65N30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Description / Table of Contents: Résumé On considère la méthode de dissections emboîtées basée sur des théorèmes de séparation introduite par Gilbert-Tarjan et Roman utilisée pour la résolution par élimination de Gauss de grands systèmes linéaires creux. Plus précisemment, on étudie une structure de données par blocs similaire à celle proposée par George dans le cadre des graphes en grille, et on démontre les propriétés suivantes: d'une part, pour des familles de graphes à degré borné admettant unn σ-théorème de séparation, 1/2≦σ〈1, le stockage secondaire de la structure par blocs contenant la matrice factorisée est linéaire par rapport à la taille du système; d'autre part, en rajoutant une hypothèse non restrictive sur la manière d'effectuer la séparation, la structure peut être construite en temps linéaire par une factorisation logique par blocs. Des exemples numériques illustrent ces résultats théoriques.
    Notes: Summary We consider the nested dissection method based on separator theorems introduced by Gilbert-Tarjan and Roman used for solving large sparse systems of linear equations. More precisely, we study a block storage scheme such as proposed by George for regular square grids and we prove the following results: first, for families of graphs of bounded degree withn σ-separator theorem, 1/2≦σ〈1, the overhead storage of the block data structure for the factored matrix is linear in system dimension; on the other hand, by adding a non restrictive assumption on the separation, this structure can be constructed by a block symbolic factorization which runs in time linear in matrix dimension. Numerical experiments illustrating these theoretical results are provided.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 12
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 55 (1989), S. 667-684 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F20 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: 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.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 13
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 215-227 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05 ; 65N20 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Discretization of the Poisson equation on a rectangle by finite differences using the standard five-point stencil yields a linear system of algebraic equations, which can be solved rapidly by the cyclic reduction method. In this method a sequence of tridiagonal linear systems is solved. The matrices of these systems commute, and we investigate numerical aspects of their ordering. In particular, we present two new ordering schemes that avoid overflow and loss of accuracy due to underflow. These ordering schemes improve the numerical performance of the subroutine HWSCRT of FISHPAK. Our orderings are also applicable to the solution of Helmholtz's equation by cyclic reduction, and to related numerical schemes, such as FACR methods.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 14
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 179-213 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; 65H10 ; CR: G1.3 ; G1.5
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We first discuss the solution of a fixed point equationxΦx of a Fréchet differentiable self-mapping Φ by iterative methods of the general form $$x_m : = \left[ {\Phi \left( {\sum\limits_{i = 0}^{m - 1} {x_i \gamma _{i,m - 1} } } \right) - \sum\limits_{j = 0}^{m - 1} {x_j \beta _{j,m} } } \right]\frac{1}{{\beta _{m,m} }},m = 1,2, \ldots ,$$ defined by two infinite nonsingular upper triangular matricesB=(β j, m ) andC=(γ j, m ) with column sums 1. We show that two such methods, defined byB, C and $$\tilde B,\tilde C$$ , respectively, are in a certain sense equivalent if and only if $$C^{ - 1} B = \tilde C^{ - 1} \tilde B$$ . In particular, $$\hat B: = C^{ - 1} $$ and the unit matrix (replacingC) define an equivalent so-called semiiterative method. We introduce (k, l)-step methods as those whereB andC have upper bandwidthk andl−1, respectively. They require storing max {k, l} previous iterates only. For stationary methodsB andC have Toeplitz structure except for their first row, which is chosen such that the column sum condition holds. An Euler method, which may require to store all iterates, is equivalent to a (stationary) (k, l)-step method if and only if the underlying conformal mapg is a rational function of the formg(w)=w(μ 0+...+μ k −k )/ (v 0+...+v l−1 w −+1). By choosingg withg(1)=1 such that for some ρ〈1 it maps |w|〉ρ onto the exterior of some continuumS known to contain the eigen values of the Fréchet derivative of Φ, one obtains a feasible procedure for designing a locally converging stationary (k, l)-step method custom-made for a set of problems. In the case Φx≔Tx+d, with a linear operatorT, where one wants to solve a linear system of equations, we show that the residual polynomials of a stationary semiiterative method are generalized Faber polynomials with respect to a particular weight function. Using another weight function leads to what we call almost stationary methods. (The classical Chebyshev iteration is an example of such a method.) We define equivalent almost stationary (k, l)-step methods and give a corresponding convergence result.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 15
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 625-626 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary A-posteriori bound is given for the computed eigenpair ( $$\hat \lambda ,\hat x$$ ), of the eigenvalue problemAx=λx, which is shown to be more realistic than the available one. A simple expression is further presented for calculating the backward error.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 16
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 52 (1988), S. 365-375 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary For each λ in some domainD in the complex plane, letF(λ) be a linear, compact operator on a Banach spaceX and letF be holomorphic in λ. Assuming that there is a ξ so thatI−F(ξ) is not one-to-one, we examine two local methods for approximating the nonlinear eigenvalue ξ. In the Newton method the smallest eigenvalue of the operator pencil [I−F(λ),F′(λ)] is used as increment. We show that under suitable hypotheses the sequence of Newton iterates is locally, quadratically convergent. Second, suppose 0 is an eigenvalue of the operator pencil [I−F(ξ),I] with algebraic multiplicitym. For fixed λ leth(λ) denote the arithmetic mean of them eigenvalues of the pencil [I−F(λ),I] which are closest to 0. Thenh is holomorphic in a neighborhood of ξ andh(ξ)=0. Under suitable hypotheses the classical Muller's method applied toh converges locally with order approximately 1.84.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 17
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 52 (1988), S. 665-682 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The paper deals with the problems of fast inversion of matricesA=T+H, whereT is Toeplitz andH is Hankel. Several algorithms are presented and compared, among them algorithms working for arbitrary strongly nonsingular matricesA=T+H.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 18
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 55 (1989), S. 501-519 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F50 ; 65K05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We study the computation of sparse null bases of equilibrium matrices in the context of the force method in structural optimization. Two classes of structural problems are considered. For the class of rigid jointed skeletal structures, we use a partitioning method suggested by Henderson and Maunder to partition the problem into a set of independent null basis computations. For the class of structures represented by a continuum, we compute a sizable fraction of the null vectors in a basis from a consideration of the finite element formulation and the bipartite graph of the equilibrium matrix. The remaining null vectors are computed by the triangular algorithm in [6]. The new algorithms find sparser bases than the triangular algorithm and are also faster.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 19
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 51 (1987), S. 143-180 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65N20 ; 65F10 ; 60K25 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Markovian queueing networks having overflow capacity are discussed. The Kolmogorov balance equations result in a linear homogeneous system, where the right null-vector is the steady-state probability distribution for the network. Preconditioned conjugate gradient methods are employed to find the null-vector. The preconditioner is a singular matrix which can be handled by separation of variables. The resulting preconditioned system is nonsingular. Numerical results show that the number of iterations required for convergence is roughly constant independent of the queue sizes. Analytic results are given to explain this fast convergence.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 20
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 51 (1987), S. 429-439 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary In this paper we consider an extension to nonlinear algebraic systems of the class of algorithms recently proposed by Abaffy, Broyden and Spedicato for general linear systems. We analyze the convergence properties, showing that under the usual assumptions on the function and some mild assumptions on the free parameters available in the class, the algorithm is locally convergent and has a superlinear rate of convergence (per major iteration, which is computationally comparable to a single Newton's step). Some particular algorithms satisfying the conditions on the free parameters are considered.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 21
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 53 (1988), S. 351-366 
    ISSN: 0945-3245
    Keywords: AMS(MOS)65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary New techniques for fast Toeplitz QR decomposition are presented. The methods are based on the shift invariance property of a Toeplitz matrix. The numerical properties of the algorithms are discussed and some comparisons are made with two other fast Toeplitz orthogonalization methods.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 22
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 53 (1988), S. 485-498 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: 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.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 23
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 47 (1985), S. 483-504 
    ISSN: 0945-3245
    Keywords: AMS(MOS) ; 65F30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary LetA be a realm×n matrix with full row rankm. In many algorithms in engineering and science, such as the force method in structural analysis, the dual variable method for the Navier-Stokes equations or more generally null space methods in quadratic programming, it is necessary to compute a basis matrixB for the null space ofA. HereB isn×r, r=n−m, of rankr, withAB=0. In many instancesA is large and sparse and often banded. The purpose of this paper is to describe and test a variation of a method originally suggested by Topcu and called the turnback algorithm for computing a banded basis matrixB. Two implementations of the algorithm are given, one using Gaussian elimination and the other using orthogonal factorization by Givens rotations. The FORTRAN software was executed on an IBM 3081 computer with an FPS-164 attached array processor at the Triangle Universities Computing Center and on a CYBER 205 vector computer. Test results on a variety of structural analysis problems including two- and three-dimensional frames, plane stress, plate bending and mixed finite element problems are discussed. These results indicate that both implementations of the algorithm yielded a well-conditioned, banded, basis matrixB whenA is well-conditioned. However, the orthogonal implementation yielded a better conditionedB for large, illconditioned problems.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 24
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 50 (1987), S. 613-632 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65G05, 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary A forward error analysis is presented for the Björck-Pereyra algorithms used for solving Vandermonde systems of equations. This analysis applies to the case where the points defining the Vandermonde matrix are nonnegative and are arranged in increasing order. It is shown that for a particular class of Vandermonde problems the error bound obtained depends on the dimensionn and on the machine precision only, being independent of the condition number of the coefficient matrix. By comparing appropriate condition numbers for the Vandermonde problem with the forward error bounds it is shown that the Björck-Pereyra algorithms introduce no more uncertainty into the numerical solution than is caused simply by storing the right-hand side vector on the computer. A technique for computing “running” a posteriori error bounds is derived. Several numerical experiments are presented, and it is observed that the ordering of the points can greatly affect the solution accuracy.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 25
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 73-91 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Linear convergence of the row cyclic Jacobi and Kogbetliantz methods can be guaranteed if certain constraints concerning the angles of rotations are implemented. Unlike the results of Forsythe and Henrici, convergence can be achieved without under-rotations.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 26
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 157-177 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; 65N20 ; 65N30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary A recursive way of constructing preconditioning matrices for the stiffness matrix in the discretization of selfadjoint second order elliptic boundary value problems is proposed. It is based on a sequence of nested finite element spaces with the usual nodal basis functions. Using a nodeordering corresponding to the nested meshes, the finite element stiffness matrix is recursively split up into two-level block structures and is factored approximately in such a way that any successive Schur complement is replaced (approximated) by a matrix defined recursively and thereform only implicitely given. To solve a system with this matrix we need to perform a fixed number (v) of iterations on the preceding level using as an iteration matrix the preconditioning matrix already defined on that level. It is shown that by a proper choice of iteration parameters it suffices to use $$v 〉 \left( {1 - \gamma ^2 } \right)^{ - \tfrac{1}{2}} $$ iterations for the so constructedv-foldV-cycle (wherev=2 corresponds to aW-cycle) preconditioning matrices to be spectrally equivalent to the stiffness matrix. The conditions involve only the constant λ in the strengthened C.-B.-S. inequality for the corresponding two-level hierarchical basis function spaces and are therefore independent of the regularity of the solution for instance. If we use successive uniform refinements of the meshes the method is of optimal order of computational complexity, if $$\gamma ^2〈 \tfrac{8}{9}$$ . Under reasonable assumptions of the finite element mesh, the condition numbers turn out to be so small that there are in practice few reasons to use an accelerated iterative method like the conjugate gradient method, for instance.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 27
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 283-289 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Comparison results for weak regular splittings of monotone matrices are derived. As an application we get upper and lower bounds for the convergence rate of iterative procedures based on multisplittings. This yields a very simple proof of results of Neumann-Plemmons on upper bounds, and establishes lower bounds, which has in special cases been conjectured by these authors.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 28
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 269-282 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65H10 ; 65W05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Linear multisplitting methods are known as parallel iterative methods for solving a linear systemAx=b. We extend the idea of multisplittings to the problem of solving a nonlinear system of equationsF(x)=0. Our nonlinear multisplittings are based on several nonlinear splittings of the functionF. In a parallel computing environment, each processor would have to calculate the exact solution of an individual nonlinear system belonging to ‘his’ nonlinear multisplitting and these solutions are combined to yield the next iterate. Although the individual systems are usually much less involved than the original system, the exact solutions will in general not be available. Therefore, we consider important variants where the exact solutions of the individual systems are approximated by some standard method such as Newton's method. Several methods proposed in literature may be regarded as special nonlinear multisplitting methods. As an application of our systematic approach we present a local convergence analysis of the nonlinear multisplitting methods and their variants. One result is that the local convergence of these methods is determined by an induced linear multisplitting of the Jacobian ofF.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 29
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 365-395 
    ISSN: 0945-3245
    Keywords: AMS(MOS) 65F, 65G ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Part I of this work deals with the forward error analysis of Gaussian elimination for general linear algebraic systems. The error analysis is based on a linearization method which determines first order approximations of the absolute errors exactly. Superposition and cancellation of error effects, structure and sparsity of the coefficient matrices are completely taken into account by this method. The most important results of the paper are new condition numbers and associated optimal component-wise error and residual estimates for the solutions of linear algebraic systems under data perturbations and perturbations by rounding erros in the arithmetic floating-point operations. The estimates do not use vector or matrix norms. The relative data and rounding condition numbers as well as the associated backward and residual stability constants are scaling-invariant. The condition numbers can be computed approximately from the input data, the intermediate results, and the solution of the linear system. Numerical examples show that by these means realistic bounds of the errors and the residuals of approximate solutions can be obtained. Using the forward error analysis, also typical results of backward error analysis are deduced. Stability theorems and a priori error estimates for special classes of linear systems are proved in Part II of this work.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 30
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 69-83 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Let $$\mathfrak{A}$$ be a real irreduciblen×n interval matrix. Then a necessary and sufficient condition is given for the sequence $$\{ \mathfrak{A}^k \} $$ of the powers of an interval matrix $$\mathfrak{A}$$ to converge to a matrix $$\mathfrak{A}^\infty $$ which is not the null matrix. In addition a criterion for $$\mathfrak{A}$$ is proved to decide whether the limit matrix $$\mathfrak{A}^\infty $$ satisfies the condition of symmetry $$\mathfrak{A}^\infty = - \mathfrak{A}^\infty $$ .
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 31
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 48 (1986), S. 349-368 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The Symmetric Tridiagonal Eigenproblem has been the topic of some recent work. Many methods have been advanced for the computation of the eigenvalues of such a matrix. In this paper, we present a divide-and-conquer approach to the computation of the eigenvalues of a symmetric tridiagonal matrix via the evaluation of the characteristic polynomial. The problem of evaluation of the characteristic polynomial is partitioned into smaller parts which are solved and these solutions are then combined to form the solution to the original problem. We give the update equations for the characteristic polynomial and certain auxiliary polynomials used in the computation. Furthermore, this set of recursions can be implemented on a regulartree structure. If the concurrency exhibited by this algorithm is exploited, it can be shown that thetime for computation of all the eigenvalues becomesO(nlogn) instead ofO(n 2) as is the case for the approach where the order is increased by only one at every step. We address the numerical problems associated with the use of the characteristic polynomial and present a numerically stable technique for the eigenvalue computation.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 32
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Ann×n complex matrixB is calledparacontracting if ‖B‖2≦1 and 0≠x∈[N(I-B)]⊥⇒‖Bx‖2〈‖x‖2. We show that a productB=B k B k−1 ...B 1 ofk paracontracting matrices is semiconvergent and give upper bounds on the subdominant eigenvalue ofB in terms of the subdominant singular values of theB i 's and in terms of the angles between certain subspaces. Our results here extend earlier results due to Halperin and due to Smith, Solomon and Wagner. We also determine necessary and sufficient conditions forn numbers in the interval [0, 1] to form the spectrum of a product of two orthogonal projections and hence characterize the subdominant eigenvalue of such a product. In the final part of the paper we apply the upper bounds mentioned earlier to provide an estimate on the subdominant eigenvalue of the SOR iteration matrix ℒω associated with ann×n hermitian positive semidefinite matrixA none of whose diagonal entries vanish.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 33
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 339-349 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05, 65B05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Description / Table of Contents: Zusammenfassung In dieser Arbeit wird die folgende Fragestellung behandelt: Wenn man weiß, daß gewisse TeilsystemeA 1·x 1=b 1, ...,A m ·x m =b m eines linearen GleichungssystemsA·x=b eindeutig lösbar sind, was kann man über die Lösbarkeit des SystemsA·x=b sagen und wie kann man gegebenenfallsx ausx 1,...,x m berechnen? Die Lösung dieses Problems erlaubt die rekursive Berechnung gewisser linearer oder quasilinearer Folgentransformationen [7, 13, 15].
    Notes: Summary This note is concerned with the following problem: Given a systemA·x=b of linear equations and knowing that certains of its subsystemsA 1·x 1=b 1, ...,A m ·x m =b m can be solved uniquely what can be said about the regularity ofA and how to find the solutionx fromx 1, ...,x m ? This question is of particular interest for establishing methods computing certain linear or quasilinear sequence transformations recursively [7, 13, 15].
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 34
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 255-268 
    ISSN: 0945-3245
    Keywords: AMS(MOS) 65F05 ; 65F25 ; 65G99 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The accumulation of rounding errors in a method used to compute the solution of an underdetermined system of linear equations at the least distance from a given point is being studied. The method used is based on orthogonal matrix decompositions.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 35
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 443-454 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The DAP architecture brings into consideration the Jacobi method where several non-interacting rotations can be performed in parallel. However the design of the algorithm is crucial in a parallel environment. In this paper we shall consider two techniques specifically designed to reduce the organisation and the arithmetic components of the parallel Jacobi method.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 36
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 479-491 
    ISSN: 0945-3245
    Keywords: AMS(MOS):65F30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary If the columns of a matrix are orthonormal and it is partitioned into a 2-by-1 block matrix, then the singular value decompositions of the blocks are related. This is the essence of the “CS decomposition”. The computation of these related SVD's requires some care. Stewart has given an algorithm that uses the LINPACK SVD algorithm together with a Jacobitype “clean-up” operation on a cross-product matrix. Our technique is equally stable and fast but avoids the cross product matrix. The simplicity of our technique makes it more amenable to parallel computation on systolic-type computer architectures. These developments are of interest because a good way to compute the generalized singular value decomposition of a matrix pair (A, B) is to compute the CS decomposition of a certain orthogonal column matrix related toA andB.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 37
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 48 (1986), S. 239-249 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We show that the greedy algorithm introduced in [1] and [5] to perform the parallel QR decomposition of a dense rectangular matrix of sizem×n is optimal. Then we assume thatm/n 2 tends to zero asm andn go to infinity, and prove that the complexity of such a decomposition is asymptotically2n, when an unlimited number of processors is available.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 38
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 48 (1986), S. 525-542 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10, 41A50, 41A17 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We investigate approximation problems of the form (P) $$\mathop {\min }\limits_{p:p(0) = 1} \mathop {\max |p(\lambda )|}\limits_{\lambda \in S} $$ withp ranging over all polynomials of degree ≦k andS any line segment of the complex plane. Problems of this type appear in connection with error bounds for the minimum residual algorithm for solving linear systems with coefficient matrixA=e iθ(dI−N), whereN=−N H andd∈ℝ. Besides a few general results on (P), we present the explicit solution of (P) for segmentsS which are parallel to the imaginary and symmetric to the real axis. This result leads to an optimal semi-iterative method for the class of matricesA=I−N, whereN=−N T.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 39
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 50 (1986), S. 17-26 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; 65G10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary In classical numerical analysis the asymptotic convergence factor (R 1-factor) of an iterative processx m+1=Axm+b coincides with the spectral radius of then×n iteration matrixA. Thus the famous Theorem of Stein and Rosenberg can at least be partly reformulated in terms of asymptotic convergence factor. Forn×n interval matricesA with irreducible upper bound and nonnegative lower bound we compare the asymptotic convergence factor (α T ) of the total step method in interval analysis with the factorα S of the corresponding single step method. We derive a result similar to that of the Theorem of Stein and Rosenberg. Furthermore we show thatα S can be less than the spectral radius of the real single step matrix corresponding to the total step matrix |A| where |A| is the absolute value ofA. This answers an old question in interval analysis.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 40
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 47 (1985), S. 175-190 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 65F05 ; 65N30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Description / Table of Contents: Summary A ≪nested dissection≫ ordering is given for solving any system of linear equationsAX=B, for the family ofsparse symmetric positive definite matrices corresponding to the class of undirected graphs withbounded degree, and satisfyinga $$\sqrt n $$ -separator theorem. If we compare the methods and results presented by Rose [9], our algorithm is more simple, but the complexity results are less general because of the restriction of bounded degree. Besides, our proofs use continually the arborescent structure on the family of separators, which is, by another way, a partition of the set of vertices for the initial graph. Then, the general implementation scheme in the finite element package MODULEF, fortwo-dimensional finite element problems, is presented, and numerical comparisons between our ordering and the standard envelope method are given.
    Notes: Résumé Nous présentons une numérotation de type ≪Nested Dissection≫ des inconnues d'un système linéaireAX=B pour des ensembles de matricescreuses symétriques définies positives correspondant à des famille de graphes non orientés,à degré borné, et admettant un $$\sqrt n $$ -thérème de séparation. Comparativement aux méthodes et résultats de Rose [9], l'algorithme présenté est plus simple, mais les théorèmes de complexité moins généraux, en raison de l'hypothèse restrictive de degré borné. En outre, les démonstrations font appel en permanence à la structure d'arbre sur la famille des séparateurs qui constitue, par ailleurs, une partition de l'ensemble des sommets du graphe initial. Nous présentons ensuite le schéma général d'implémentation dans le cadre du code d'éléments finis MODULEF pourdes problèmes plans d'éléments finis, et nous donnons quelques mesures comparatives avec la numérotation plus classique qui tend à minimiser le profil de la matrice.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 41
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 49 (1986), S. 81-94 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary This paper presents a new algorithm for computing theQR factorization of anm×n Toeplitz matrix inO(mn) operations. The algorithm exploits the procedure for the rank-1 modification and the fact that both principal (m−1)×(n−1) submatrices of the Toeplitz matrix are identical. An efficient parallel implementation of the algorithm is possible.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 42
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 48 (1986), S. 499-523 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We derive new estimates for the rate of convergence of the conjugate gradient method by utilizing isolated eigenvalues of parts of the spectrum. We present a new generalized version of an incomplete factorization method and compare the derived estimates of the number of iterations with the number actually found for some elliptic difference equations and for a similar problem with a model empirical distribution function.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
Close ⊗
This website uses cookies and the analysis tool Matomo. More information can be found here...