For two-point boundary value problems in weak formulation with piecewise constant coefficients and piecewise continuous right-hand side functions we derive a representation of its weak solution by local Green's functions. Then we use it to generate exact three-point discretizations by Galerkin's method on essentially arbitrary grids. The coarsest possible grid is the set of points at which the piecewise constant coefficients and the right- hand side functions are discontinuous. This grid can be refined to resolve any solution properties like boundary and interior layers much more correctly. The proper basis functions for the Galerkin method are entirely defined by the local Green's functions. The exact discretizations are of completely exponentially fitted type and stable. The system matrices of the resulting tridiagonal systems of linear equations are in any case irreducible M-matrices with a uniformly bounded norm of its inverse.

For solving systems of grid equations approximating elliptic boundary value problems a method of constructing variable preconditioning procedures is presented. The main purpose is to discuss how an efficient preconditioning iterative procedure can be constructed in the case of elliptic problems with disproportional coefficients, e.g. equations with a large coefficient in the reaction term (or a small diffusion coefficient). The optimality of the suggested technique is based on fictitious space and multilevel decom- position methods. Using an additive form of the preconditioners, we intro- duce factors into the preconditioners to optimize the corresponding conver- gence rate. The optimization with respect to these factors is used at each step of the iterative process. The application of this technique to two-level $p$-hierarchical precondi- tioners and domain decomposition methods is considered too.

A new backward stable, structure preserving method of complexity O(n^3) is presented for computing the stable invariant subspace of a real Hamiltonian matrix and the stabilizing solution of the continuous-time algebraic Riccati equation. The new method is based on the relationship between the invariant subspaces of the Hamiltonian matrix H and the extended matrix /0 H\ and makes use \H 0/ of the symplectic URV-like decomposition that was recently introduced by the authors.

Rank-revealing ULV and URV factorizations are useful tools to determine the rank and to compute bases for null-spaces of a matrix. However, in the practical ULV (resp. URV ) factorization each left (resp. right) null vector is recomputed from its corresponding right (resp. left) null vector via triangular solves. Triangular solves are required at initial factorization, refinement and updating. As a result, algorithms based on these factorizations may be expensive, especially on parallel computers where triangular solves are expensive. In this paper we propose an alternative approach. Our new rank-revealing ULV factorization, which we call ¨top-down¨ ULV factorization ( TDULV -factorization) is based on right null vectors of lower triangular matrices and therefore no triangular solves are required. Right null vectors are easy to estimate accurately using condition estimators such as incremental condition estimator (ICE). The TDULV factorization is shown to be equivalent to the URV factorization with the advantage of circumventing triangular solves.

This paper is concerned with a specific finite element strategy for solving elliptic boundary value problems in domains with corners and edges. First, the anisotropic singular behaviour of the solution is described. Then the finite element method with anisotropic, graded meshes and piecewise linear shape functions is investigated for such problems; the schemes exhibit optimal convergence rates with decreasing mesh size. For the proof, new local interpolation error estimates for functions from anisotropically weighted spaces are derived. Finally, a numerical experiment is described, that shows a good agreement of the calculated approximation orders with the theoretically predicted ones.

This paper deals with parallelization methods for time-dependent problems where the time steps are shared out among the processors. A Full Multigrid technique serves as solution algorithm, hence information of the preceding time step and of the coarser grid is necessary to compute the solution at each new grid level. Applying the usual extrapolation formula to process this information, the parallelization will not be very efficient. We developed another extrapolation technique which causes a much higher parallelization effect. Test examples show that no essential loss of exactness appears, such that the method presented here shall be well-applicable.

The paper gives an overview over local inequalities for anisotropic simplicial Lagrangian finite elements. The main original contributions are the estimates for higher derivatives of the interpolation error, the formulation of the assumptions on admissible anisotropic finite elements in terms of geometrical conditions in the three-dimensional case, and an anisotropic variant of the inverse inequality. An application of anisotropic meshes in the context of a stabilized Galerkin method for a convection-diffusion problem is given.

Nonsteady Navier-Stokes equations represent a differential-algebraic system of strangeness index one after any spatial discretization. Since such systems are hard to treat in their original form, most approaches use some kind of index reduction. Processing this index reduction it is important to take care of the manifolds contained in the differential-algebraic equation (DAE). We investigate for several discretization schemes for the Navier-Stokes equations how the consideration of the manifolds is taken into account and propose a variant of solving these equations along the lines of the theoretically best index reduction. Applying this technique, the error of the time discretisation depends only on the method applied for solving the DAE.

The main goal of the thesis is to give explicit formulas for the computation of the local Picard group of some binoid algebras. In particular the Stanley-Reisner and the general monomial case are covered. In order to do so, we introduce a new topology on the spectrum of the binoid algebra over a field, that we call combinatorial topology, coarser than Zariski topology, that mimics the topology on the spectrum of a pointed monoid. Then we use the tools of cohomology of sheaves on schemes of pointed monoids in order to prove some formulas for computing the cohomology of the sheaf of units on the punctured spectrum of a simplicial pointed monoid. We prove that the Picard group of a Stanley-Reisner algebra is trivial with the Zariski topology and, finally, we use these results and the combinatorial topology to obtain the explicit formulas. Lastly, we extend this result to any quotient of the polynomial ring over a monomial ideal.

Diese Dissertation behandelt ein Problem aus der Deformationsquantisierung: Nachdem man die Quantisierung eines klassischen Systems konstruiert hat, würde man gerne ihre mathematischen Eigenschaften verstehen (sowohl die des klassischen Systems als auch die des Quantensystems). Falls beide Systeme durch *-Algebren über dem Körper der komplexen Zahlen beschrieben werden, bedeutet dies dass man die Eigenschaften bestimmter *-Algebren verstehen muss: Welche Darstellungen gibt es? Was sind deren Eigenschaften? Wie können die Zustände in diesen Darstellungen beschrieben werden? Wie kann das Spektrum der Observablen beschrieben werden? Um eine hinreichend allgemeine Behandlung dieser Fragen zu ermöglichen, wird das Konzept von abstrakten O*-Algebren entwickelt. Dies sind im Wesentlichen *-Algebren zusammen mit einem Kegel positiver linearer Funktionale darauf (z.B. die stetigen positiven linearen Funktionale wenn man mit einer *-Algebra startet, die mit einer gutartigen Topologie versehen ist). Im Anschluss daran wird dieser Ansatz dann auf zwei Beispiele aus der Deformationsquantisierung angewandt, die im Detail untersucht werden. / This thesis discusses and proposes a solution for one problem arising from deformation quantization: Having constructed the quantization of a classical system, one would like to understand its mathematical properties (of both the classical and quantum system). Especially if both systems are described by ∗-algebras over the field of complex numbers, this means to understand the properties of certain ∗-algebras: What are their representations? What are the properties of these representations? How can the states be described in these representations? How can the spectrum of the observables be described? In order to allow for a sufficiently general treatment of these questions, the concept of abstract O ∗-algebras is introduced. Roughly speaking, these are ∗ -algebras together with a cone of positive linear functionals on them (e.g. the continuous ones if one starts with a ∗-algebra that is endowed with a well-behaved topology). This language is then applied to two examples from deformation quantization, which will be studied in great detail.