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  (15)
  • CR: G1.3  (15)
  • Springer  (15)
  • American Meteorological Society
  • Blackwell Publishing Ltd
  • Copernicus
  • Hindawi
  • Institute of Electrical and Electronics Engineers
  • Molecular Diversity Preservation International
  • Springer Science + Business Media
  • 2020-2022
  • 2010-2014
  • 1985-1989  (15)
  • 1960-1964
  • 1989  (15)
  • Mathematics  (15)
  • Geography
  • Architecture, Civil Engineering, Surveying
Collection
  • Articles  (15)
Publisher
  • Springer  (15)
  • American Meteorological Society
  • Blackwell Publishing Ltd
  • Copernicus
  • Hindawi
  • +
Years
  • 2020-2022
  • 2010-2014
  • 1985-1989  (15)
  • 1960-1964
Year
Topic
  • Mathematics  (15)
  • Geography
  • Architecture, Civil Engineering, Surveying
  • 1
    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 ...
  • 2
    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 ...
  • 3
    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 ...
  • 4
    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 ...
  • 5
    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 ...
  • 6
    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 ...
  • 7
    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 ...
  • 8
    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 ...
  • 9
    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 ...
  • 10
    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 ...
Close ⊗
This website uses cookies and the analysis tool Matomo. More information can be found here...