Forward Error Analysis of Gaussian Elimination. Part I: Error and Residual Estimates.
We introduce a new way of the analysis of iterative aggregation-disaggregation methods for computing stationary probability distribution vectors of stochastic matrices. This new approach is based on the Fourier transform of the error propagation matrix. Exact formula for its spectrum can be obtained if the stochastic matrix is circulant. Some examples are presented.
The error analysis of preconditioned waveform relaxation iterations for differential systems is presented. This analysis extends and refines previous results by Burrage, Jackiewicz, Nørsett and Renaut by incorporating all terms in the expansion of the error of waveform relaxation iterations in the Laplace transform domain. Lower bounds for the size of the window of rapid convergence are also obtained. The theory is illustrated for waveform relaxation methods applied to differential systems resulting...
A unified approach to the theory and construction of direct methods is presented. The approach is based on the idea of the transfer of conditions. In examples it is shown how to obtain a particular method from the general algorithm.
In this paper we present a nonsingularity result which is a generalization of Nekrasov property by using two different permutations of the index set. The main motivation comes from the following observation: matrices that are Nekrasov matrices up to the same permutations of rows and columns, are nonsingular. But, testing all the permutations of the index set for the given matrix is too expensive. So, in some cases, our new nonsingularity criterion allows us to use the results already calculated...
This paper addresses the derivation of new second-kind Fredholm combined field integral equations for the Krylov iterative solution of tridimensional acoustic scattering problems by a smooth closed surface. These integral equations need the introduction of suitable tangential square-root operators to regularize the formulations. Existence and uniqueness occur for these formulations. They can be interpreted as generalizations of the well-known Brakhage-Werner [A. Brakhage and P. Werner, Arch....
We present algorithms for the determination of polynomials orthogonal with respect to a positive weight function multiplied by a polynomial with simple roots inside the interval of integration. We apply these algorithms to search for and calculate all possible sequences of imbedded quadratures of maximal polynomials order of precision for the generalized Laguerre and Hermite weight functions.
The integration to steady state of many initial value ODEs and PDEs using the forward Euler method can alternatively be considered as gradient descent for an associated minimization problem. Greedy algorithms such as steepest descent for determining the step size are as slow to reach steady state as is forward Euler integration with the best uniform step size. But other, much faster methods using bolder step size selection exist. Various alternatives are investigated from both theoretical and practical...
The purpose of our work is to develop an automatic shape optimization tool for runner wheel blades in reaction water turbines, especially in Kaplan turbines. The fluid flow is simulated using an in-house incompressible turbulent flow solver based on recently introduced isogeometric analysis (see e.g. J. A. Cotrell et al.: Isogeometric Analysis: Toward Integration of CAD and FEA, Wiley, 2009). The proposed automatic shape optimization approach is based on a so-called hybrid optimization which combines...
A method of characterizing all eigenvalues of a preconditioned discretized scalar diffusion operator with Dirichlet boundary conditions has been recently introduced in Gergelits, Mardal, Nielsen, and Strakoš (2019). Motivated by this paper, we offer a slightly different approach that extends the previous results in some directions. Namely, we provide bounds on all increasingly ordered eigenvalues of a general diffusion or elasticity operator with tensor data, discretized with the conforming finite...