Two efficient methods for solving generalized Lyapunov equations and their implementations in FORTRAN 77 are presented. The first one is a generalization of the Bartels--Stewart method and the second is an extension of Hammarling's method to generalized Lyapunov equations. Our LAPACK based subroutines are implemented in a quite flexible way. They can handle the transposed equations and provide scaling to avoid overflow in the solution. Moreover, the Bartels--Stewart subroutine offers the optional estimation of the separation and the reciprocal condition number. A brief description of both algorithms is given. The performance of the software is demonstrated by numerical experiments.

The investigation of new parallel algorithms for MIMD computers requires some postprocessing facilities for quickly evaluating the behavior of those algorithms We present two kinds of visualization tool implementations for 2D and 3D finite element applications to be used on a parallel computer and a host workstation.

Some boundary value problems yield anisotropic solutions, e.g. solutions with boundary layers. If such problems are to be solved with the finite element method (FEM), anisotropically refined meshes can be advantageous. In order to construct these meshes or to control the error one aims at reliable error estimators. For \emph{isotropic} meshes many estimators are known, but they either fail when used on \emph{anisotropic} meshes, or they were not applied yet. For rectangular (or cuboidal) anisotropic meshes a modified error estimator had already been found. We are investigating error estimators on anisotropic tetrahedral or triangular meshes because such grids offer greater geometrical flexibility. For the Poisson equation a residual error estimator, a local Dirichlet problem error estimator, and an $L_2$ error estimator are derived, respectively. Additionally a residual error estimator is presented for a singularly perturbed reaction diffusion equation. It is important that the anisotropic mesh corresponds to the anisotropic solution. Provided that a certain condition is satisfied, we have proven that all estimators bound the error reliably.

In the paper, the parallelization of multi-grid methods for solving second-order elliptic boundary value problems in two-dimensional domains is discussed. The parallelization strategy is based on a non-overlapping domain decomposition data structure such that the algorithm is well-suited for an implementation on a parallel machine with MIMD architecture. For getting an algorithm with a good paral- lel performance it is necessary to have as few communication as possible between the processors. In our implementation, communication is only needed within the smoothing procedures and the coarse-grid solver. The interpolation and restriction procedures can be performed without any communication. New variants of smoothers of Gauss-Seidel type having the same communication cost as Jacobi smoothers are presented. For solving the coarse-grid systems iterative methods are proposed that are applied to the corresponding Schur complement system. Three numerical examples, namely a Poisson equation, a magnetic field problem, and a plane linear elasticity problem, demonstrate the efficiency of the parallel multi- grid algorithm.

CHAPTER I. Basic Equations. Fundamental Matrices. Thermo-Radiation Conditions 1. Basic differential equations of thermoelasticity theory 2. Fundamental matrices 3. Thermo-radiating conditions. Somigliana type integral representations CHAPTER II. Formulation of Boundary Value and Interface Problems 4. Functional spaces 5. Formulation of basic and mixed BVPs 6. Formulation of crack type problems 7. Formulation of basic and mixed interface problems CHAPTER III. Uniqueness Theorems 8. Uniqueness theorems in pseudo-oscillation problems 9. Uniqueness theorems in steady state oscillation problems CHAPTER IV. Potentials and Boundary Integral Operators 10. Thermoelastic steady state oscillation potentials 11. Pseudo-oscillation potentials CHAPTER V. Regular Boundary Value and Interface Problems 12. Basic BVPs of pseudo-oscillations 13. Basic exterior BVPs of steady state oscillations 14. Basic interface problems of pseudo-oscillations 15. Basic interface problems of steady state oscillations CHAPTER VI. Mixed and Crack Type Problems 16. Basic mixed BVPs 17. Crack type problems 18. Mixed interface problems of steady state oscillations 19. Mixed interface problems of pseudo-oscillations

The parallel version of precondition techniques is developed for matrices arising from the Galerkin boundary element method for two-dimensional domains with Dirichlet boundary conditions. Results were obtained for implementations on a transputer network as well as on an nCUBE-2 parallel computer showing that iterative solution methods are very well suited for a MIMD computer. A comparison of numerical results for iterative and direct solution methods is presented and underlines the superiority of iterative methods for large systems.

The parallel version of precondition iterative techniques is developed for matrices arising from the panel boundary element method for three-dimensional simple connected domains with Dirichlet boundary conditions. Results were obtained on an nCUBE-2 parallel computer showing that iterative solution methods are very well suited also in three-dimensional case for implementation on a MIMD computer and that they are much more efficient than usual direct solution techniques.

Extrapolation methods for the solution of partial differential equations are commonly based on the existence of error expansions for the approximate solution. Implicit extrapolation, in the contrast, is based on applying extrapolation indirectly, by using it on quantities like the residual. In the context of multigrid methods, a special technique of this type is known as \034 -extrapolation. For finite element systems this algorithm can be shown to be equivalent to higher order finite elements. The analysis is local and does not use global expansions, so that the implicit extrapolation technique may be used on unstructured meshes and in cases where the solution fails to be globally smooth. Furthermore, the natural multilevel structure can be used to construct efficient multigrid and multilevel preconditioning techniques. The effectivity of the method is demonstrated for heat conduction problems and problems from elasticity theory.

The parallel version of precondition iterative techniques is developed for matrices arising from the panel boundary element method for three-dimensional simple connected domains with Dirichlet boundary conditions. Results were obtained on an nCUBE-2 parallel computer showing that iterative solution methods are very well suited also in three-dimensional case for implementation on a MIMD computer and that they are much more efficient than usual direct solution techniques.

The two-dimensional interface crack problem is investigated for anisotropic bodies in the Comninou formulation. It is established that, as in the isotropic case, properly incorporating contact zones at the crack tips avoids contradictions connected with the oscillating asymptotic behaviour of physical and mechanical characteristics leading to the overlapping of material. Applying the special integral representation formulae for the displacement field the problem in question is reduced to the scalar singular integral equation with the index equal to -1. The analysis of this equation is given. The comparison with the results of previous authors shows that the integral equations corresponding to the interface crack problems in the anisotropic and isotropic cases are actually the same from the point of view of the theoretical and numerical analysis.