TY - GEN A1 - Lang, Jens A1 - Walter, Artur T1 - A Finite Element Method Adaptive in Space and Time for Nonlinear Reaction-Diffusion- Systems. N2 - Large scale combustion simulations show the need for adaptive methods. First, to save computation time and mainly to resolve local and instationary phenomena. In contrast to the widespread method of lines, we look at the reaction- diffusion equations as an abstract Cauchy problem in an appropriate Hilbert space. This means, we first discretize in time, assuming the space problems solved up to a prescribed tolerance. So, we are able to control the space and time error separately in an adaptive approach. The time discretization is done by several adaptive Runge-Kutta methods whereas for the space discretization a finite element method is used. The different behaviour of the proposed approaches are demonstrated on many fundamental examples from ecology, flame propagation, electrodynamics and combustion theory. {\bf Keywords:} initial boundary value problem, Rothe- method, adaptive Runge-Kutta method, finite elements, mesh refinement. {\bf AMS CLASSIFICATION:} 65J15, 65M30, 65M50. T3 - ZIB-Report - SC-92-05 KW - initial boundary value problem KW - Rothe method KW - adaptive Runge-Kutta method KW - finite elements KW - mesh refinement Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-755 ER - TY - GEN A1 - Grötschel, Martin T1 - Discrete Mathematics in Manufacturing. N2 - Manufacturing is a topic that provides rich opportunities for important mathematical contributions to real-world problems. The purpose of this paper is to show, by means of several examples, where and how mathematical problems of a discrete nature arise in manufacturing and to demonstrate the savings and improvements that can be achieved by employing the techniques of combinatorial optimization. The topics covered range from the design phase of a product (e. g.,routing, placement and via minimization in VLSI design), the control of CNC machines (e. g., drilling and plotting), to the management of assembly lines, storage systems and whole factories. We also point out difficulties in the modelling of complex situations and outline the algorithmic methods that are used for the solution of the mathematical problems arising in manufacturing. {\bf Key words:} discrete mathematics , combinatorial optimization, applications to manufacturing. T3 - ZIB-Report - SC-92-03 KW - discrete mathematics KW - combinatorial optimization KW - applications to manufacturing Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-738 ER - TY - GEN A1 - Wakabayashi, Yoshiko T1 - Medians of Binary Relations: Computational Complexity. N2 - Let $\Re$ be the set of all binary relations on a finite set $N$ and $d$ be the symmetric difference distance defined on $\Re$. For a given profile $\Pi = (R_1,...,R_m) \in R^m$, a relation $R* \in \Re $ that minimizes the function $\sum^m_{k=1} d(R_k,R) $ is called a median relation of $\Pi$. A number of problems occuring in the social sciences, in qualitative data analysis and in multicriteria decision making can be modelled as problems of finding medians of a profile of binary relations. In these contexts the profile $\Pi$ represents collected data (preferences, similarities, games) and the objective is that of finding a median relation of $\Pi$ with some special feature (representing e. g., consensus of preferences, clustering of similar objects, ranking of teams, etc.). In this paper we analyse the computational complexity of all such problems in which the median is required to satisfy one or more of the properties: reflexitivity, symmetry, antisymmetry, transitivity and completeness. We prove that whenever transitivity is required (except when symmetry and completeness are also simultaneously required) then the corresponding median problem is $NP$-hard. In some cases we prove that they remain $NP$-hard when the profile $\Pi$ has a fixed number of binary relations. T3 - ZIB-Report - SC-92-04 Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-746 ER - TY - GEN A1 - Schöffel, Klaus T1 - Computational Chemistry Software at ZIB. N2 - The following report intends to provide a survey over the computational chemistry molecular structure software installed on the supercomputers CRAY X-MP/216 and CRAY Y-MP2E/164 at ZIB. It shows what kind of problems can be tackled with the existing chemistry software, which covers a wide range of ab initio, semiempirical, molecular mechanics, and dynamics applications. T3 - ZIB-Report - TR-92-02 Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4879 ER - TY - GEN A1 - Ackermann, Jörg A1 - HELFRICH, K. T1 - Radius of Convergence of the 1/Z-Expansion for Diatomic Molecules: The Ground State of the Isoelectronic H2Sequence. N2 - Using the perturbational-variational Rayleigh-Ritz matrix formalism, the 1/Z-expansion for the ground state of the isoelectronic $H_2$ sequence in the range of the internuclear distance $0.2\le R \le 9.0$ is calculated. Also lower bounds of the radius of convergence, based on Kato's theory of linear operators, are given. The numerical results of the 1/Z-expansion can be compared with the exact results and do not converge in the whole R-range. This behavior is in qualitative agreement with the lower bounds for the radius of convergence and enlights some still open properties of 1/Z- expansions for this sequence in the literature. {\bf PACS:} 31.15 + q; 31.20 Di; 31.20 Tz. T3 - ZIB-Report - SC-92-02 Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-721 ER - TY - GEN A1 - Bornemann, Folkmar A. A1 - Yserentant, Harry T1 - A Basic Norm Equivalence for the Theory of Multilevel Methods. N2 - Subspace decompositions of finite element spaces based on $L2$-like orthogonal projections play an important role for the construction and analysis of multigrid like iterative methods. Recently several authors proved the equivalence of the associated discrete norms with the $H^1$-norm. The present report gives an elementary, self-contained derivation of this result which is based on the use of $ K$-functionals known from the theory of interpolation spaces. {\bf Keywords:} multilevel methods, nonuniform meshes, optimal convergence rates. {\bf AMS(MOS) Subject classifications:} 65N55, 65N30, 65N50. T3 - ZIB-Report - SC-92-01 KW - multilevel methods KW - nonuniform meshes KW - optimal convergence rates Y1 - 1992 U6 - https://doi.org/10.1007/BF01388699 N1 - No preprint available ER - TY - GEN A1 - Nowak, Ulrich A1 - Weimann, Lutz T1 - A Family of Newton Codes for Systems of Highly Nonlinear Equations. N2 - This report presents new codes for the numerical solution of highly nonlinear systems. They realize the most recent variants of affine invariant Newton Techniques due to Deuflhard. The standard method is implemented in the code NLEQ1, whereas the code NLEQ2 contains a rank reduction device additionally. The code NLEQ1S is the sparse version of NLEQ1, i.e. the arising linear systems are solved with sparse matrix techniques. Within the new implementations a common design of the software in view of user interface and internal modularization is realized. Numerical experiments for some rather challenging examples illustrate robustness and efficiency of algorithm and software. T3 - ZIB-Report - TR-91-10 Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4855 ER - TY - GEN A1 - Schöffel, Klaus T1 - Ab initio Quantum Chemical Calculations with GAMESS-UK and GAUSSIAN90 Program Packeges - A Comparison -. N2 - Two commercially available molecular electronic structure software packages GAUSSIAN90 and GAMESS-UK are compared. Basis for this comparison is a benchmark suite which is designed to highlight the typical range of calculations commonly performed by the ab initio computational chemist. T3 - ZIB-Report - TR-92-01 Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4867 ER - TY - GEN A1 - Hohmann, Andreas T1 - An Adaptive Continuation Method for Implicitly Defined Surfaces. N2 - A new method for the numerical aproximation of an implicitly defined surface is presented. It is a generalization of the Euler- Gauss-Newton method for implicitly defined (one- parameter) curves to the case of (two-parameter) surfaces. The basic task in the more general case is an efficient combination of modern CAGD techniques (such as triangular Bernstein-Bzier patches and the nine parameter Hermite interpolant) and the rank deficient Gauss-Newton method. T3 - ZIB-Report - SC-91-20 Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-705 ER - TY - GEN A1 - Godsil, C. D. A1 - Grötschel, Martin A1 - Welsh, D. J. A. T1 - Combinatorics in Statistical Physics T3 - ZIB-Report - SC-91-19 Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-694 ER - TY - GEN A1 - Wulkow, Michael T1 - Adaptive Treatment of Polyreactions in Weighted Sequence Spaces. N2 - Countable systems of ordinary differential equations appear frequently in chemistry, physics, biology and statistics. They can be considered as ordinary differential equations in sequence spaces. In this work, a fully adaptive algorithm for the computational treatment of such systems is developed. The method is based on a time discretization of an abstract Cauchy problem in Hilbert space and a discrete Galerkin approach for the discretization of the arising stationary subproblems. The Galerkin method uses orthogonal functions of a discrete variable, which are generated by certain weight functions. A theory of countable systems in the associated weighted sequence spaces is developed as well as a theory of the Galerkin method. The Galerkin equations are solved adaptively either by use of analytical properties of the orthogonal functions or by an appropriate numerical summation. The resulting algorithm CODEX is applied to examples of technological interest, in particular from polymer chemistry. T3 - ZIB-Report - SC-91-17 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-673 ER - TY - GEN A1 - Ackermann, Jörg A1 - Wulkow, Michael T1 - The Treatment of Macromolecular Processes with Chain-Length-Dependent Reaction Coefficients - An Examplefrom Soot Formation. N2 - The description of chain length distributions in macromolecular reaction kinetics leads to so-called countable systems of differential equations. In particular, when the appearing reaction rate coefficients depend on the chain length of the reacting macromolecules itself, an efficient numerical treatment of these systems is very difficult. Then even the evaluation of the right-hand side of the system can become prohibitively expensive with respect to computing time. In this paper we show how the discrete Galerkin method can be applied to such problems. The existing algorithm CODEX is improved by use of a multiplicative error correction scheme for time discretization and a new type of numerical preprocessing by means of a Gauss summation. Both ideas are exemplary for a wide class of approximation types and are described very briefly here. The new numerical techniques are tested on an example from soot formation, where the coagulation of molecules is modeled in terms of reaction coefficients depending on the surface of the particles and their collision frequency. T3 - ZIB-Report - SC-91-18 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-682 ER - TY - GEN A1 - Schöpf, R. A1 - Deuflhard, Peter T1 - OCCAL: A mixed symbolic- numeric Optimal Control CALculator. N2 - The numerical solution of optimal control problems by indirect methods (such as multiple shooting or collocation) requires a considerable amount of analytic calculation to establish a numerically tractable system. These analytic calculations, though being rather tedious in realistic examples, are nowadays mostly still done by hand-and thus prone to calculation errors. The paper aims at automating this analytic processing to a reasonable extent by means of a modern symbolic manipulation language (here: REDUCE). In its present stage of development the package OCCAL (mnemotechnically for \underline{O}ptimal \underline{C}ontrol \underline{CAL}culator) permits an interactive use, covering tasks like automatic determination of control and, in case of a singular control, of its order. In simpler problems, the present version of OCCAL automatically produces the full subroutine input for a MULtiple shooting code (MULCON) with adaptive numerical CONtinuation. In more complicate problems where singular sub-arcs may occur or where the sequence of sub-arcs of the optimal trajectory is unclear OCCAL is a significant help in reducing analytic pre-processing. Numerical examples illustrate the performance of OCCAL/MULCON. T3 - ZIB-Report - SC-91-13 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-638 ER - TY - GEN A1 - Hoppe, Ronald H. W. A1 - Kornhuber, Ralf T1 - Adaptive Multilevel - Methods for Obstacle Problems. N2 - We consider the discretization of obstacle problems for the Laplacian by piecewise linear finite elements. Assuming that the discrete problems are reduced to a sequence of linear problems by suitable active set strategies, the linear problems are solved iteratively by preconditioned c-g iterations. The proposed preconditioners are treated theoretically as abstract additive Schwarz methods and are implemented as truncated hierarchical basis preconditioners. To allow for local mesh refinement we derive semi-local and local a posteriori error estimates, providing lower and upper estimates for the discretization error. The theoretical results are illustrated by numerical computations. T3 - ZIB-Report - SC-91-16 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-667 ER - TY - GEN A1 - Gatermann, Karin A1 - Hohmann, Andreas T1 - Hexagonal Lattice Dome - Illustration of a Nontrivial Bifurcation Problem. N2 - The deformation of a hexagonal lattice dome under an external load is an example of a parameter dependent system which is equivariant under the symmetry group of a regular hexagon. In this paper the mixed symbolic-numerical algorithm SYMCON is applied to analyze its steady state solutions automatically showing their different symmetry and stability properties. T3 - ZIB-Report - SC-91-08 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-587 ER - TY - GEN A1 - Ziegler, Günter M. T1 - On the Difference Between Real and Complex Arrangements. N2 - If $B$ is an arrangement of linear complex Hyperplanes in $C^d$, then the following can be constructed from knowledge of its intersection lattice: (a) the cohomology groups of the complement [Br], (b) the cohomology algebra of the complement [OS], (c) the fundamental group of the complement, if $d\le2$, (d) the singularity link up to homeomorphism, if $d\le3$, (e) the singularity link up to homotopy type [ZZ]. If $B'$ is, more generally, a 2-arrangement in $ R^{2d}$ (an arrangement of real subspaces of codimension 2 with even-dimensional intersections), then the intersection lattice still determines (a) the cohomology groups of the complement [GM] and (e) the homotopy type of the singularity link [ZZ]. We show, however, that for 2-arrangements the data (b), (c) and (d) are not determined by the intersection lattice. They require the knowledge of extra information on sign patterns, which can be computed as determinants of linear relations, or (equivalently) as linking coefficients in the sense of knot theory. T3 - ZIB-Report - SC-91-14 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-648 ER - TY - GEN A1 - Ziegler, Günter M. A1 - Zivaljevic, R. T. T1 - Homotopy Types of Subspace Arrangements via Diagrams of Spaces. N2 - We prove combinatorial formulas for the homotopy type of the union of the subspaces in an (affine, compactified affine, spherical or projective) subspace arrangement. From these formulas we derive results of Goresky & MacPherson on the homology of the arrangement and the cohomology of its complement. The union of an arrangement can be interpreted as the direct limit of a diagram of spaces over the intersection poset. A closely related space is obtained by taking the homotopy direct limit of this diagram. Our method consists in constructing a combinatorial model diagram over the same poset, whose homotopy limit can be compared to the original one by usual homotopy comparison results for diagrams of spaces. T3 - ZIB-Report - SC-91-15 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-658 ER - TY - GEN A1 - Schmidt, Frank T1 - An Adaptive Approach to the Numerical Solution of Fresnel's Wave Equation. N2 - An adaptive approach to the numerical solution of the wave propagation in integrated optics devices with 1D cross sections is described. First, Fresnel's approximation of the exact wave equation resulting from Maxwell's equations is considered. A criterion to estimate the validity of this approximation is derived. Fresnel's wave equation being formally equivalent to Schroedinger's equation uniquely defines an initial-boundary-value problem, which is solved numerically by a stepwise calculation of the propagating field. Discretization in longitudinal direction first with stepsize control leads to a stationary subproblem for the transversal field distribution, which is then handled by an adaptive finite element method. Thus full adaptivity of the algorithm is realized. The numerical examples are concentrated on taper structures playing an essential role in integrated optics devices for telecommunication systems. T3 - ZIB-Report - SC-91-12 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-628 ER - TY - GEN A1 - Sturmfels, Bernd A1 - Ziegler, Günter M. T1 - Extension Spaces of Oriented Matriods. N2 - We study the space of all extensions of a real hyperplane arrangement by a new pseudo- hyperplane, and, more generally, of an oriented matroid by a new element. The question whether this space has the homotopy type of a sphere is a special case of the "Generalized Baues Problem" of Billera, Kapranov & Sturmfels, via the Bohne-Dress Theorem on zonotopal tilings. We prove that the extension space is spherical for the class of strongly euclidean oriented matroids. This class includes the alternating matroids and all oriented matroids of rank at most 3 or of corank at most 2. In general it is not even known whether the extension space is connected. We show that the subspace of realizable extensions is always connected but not necessarily spherical. T3 - ZIB-Report - SC-91-11 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-618 ER - TY - GEN A1 - Ziegler, Günter M. T1 - Higher Bruhat Orders and Cyclic Hyperplane Arrangements. N2 - We study the higher Bruhat orders $B(n,k)$ of Manin & Schechtman [MaS] and - characterize them in terms of inversion sets, - identify them with the posets $U(C^{n+1,r},n+1)$ of uniform extensions of the alternating oriented matroids $C^{n,r}$ for $r:=n-k$ (that is, with the extensions of a cyclic hyperplane arrangement by a new oriented pseudoplane), - show that $B(n,k)$ is a lattice for $k =1$ and for $r\le 3$, but not in general, - show that $B(n,k)$ is ordered by inclusion of inversion sets for $k=1$ and for $r\le 4$. However, $B(8,3)$ is not ordered by inclusion. This implies that the partial order $B_\subseteq (n,k)$ defined by inclusion of inversion sets differs from $B(n,k)$ in general. We show that the proper part of $B_\subseteq (n,k)$ is homotopy equivalent to $S^{r-2}$. Consequently, - $B(n,k)\simeq S^{r-2}$ for $k=1$ and for $r\le 4$. In contrast to this, we find that the uniform extension poset of an affine hyperplane arrangement is in general not graded and not a lattice even for $r=3$, and that the proper part is not always homotopy equivalent to $S^{r(M)-2}$. T3 - ZIB-Report - SC-91-10 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-607 ER - TY - GEN A1 - Deuflhard, Peter A1 - Nowak, Ulrich A1 - Pöhle, Uwe A1 - Schmidt, B. Ch. A1 - Weyer, J. T1 - Die Ausbreitung von HIV/AIDS in Ballungsgebieten. N2 - The rapidly increasing number of AIDS cases requires a realistic estimation of the future development of the HIV/AIDS disease. For that purpose we develop a large system of coupled nonlinear differential equations describing simultaneously the dynamics of the development of the disease, the population size, the gender and age structure. A set of 1650 coupled equations are linked by balanced parameters. The balancing procedure is described by a set of (formally) 2,178,000 additional algebraic conditions. As the considered system is stiff, it requires new special extrapolation methods combined with techniquees of dynamical sparsing for the solution of sparsely filled systems. According to our simulations we expect 19,0,000 deaths caused by AIDS in the Federal Republic of Germany (former territories) in the year 2000. Such an epidemical spread would tie up about 4-7 percent of the actual health care workers. T3 - ZIB-Report - TR-91-09 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4844 ER - TY - GEN A1 - Bornemann, Folkmar A. T1 - A Sharpened Condition Number Estimate for the BPX Preconditioner of Elliptic Finite Element Problems on Highly Nonuniform Triangulations. N2 - In this paper it is shown that for highly nonuniformly refined triangulations the condition number of the BPX preconditioner for elliptic finite element problems grows at most linearly in the depth of refinement. This is achieved by viewing the computational available version of the BPX preconditioner as an abstract additive Schwarz method with exact solvers. {\bf AMS CLASSIFICATION:} 65F10, 65F35, 65N20, 65N30. T3 - ZIB-Report - SC-91-09 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-596 ER - TY - GEN A1 - Deuflhard, Peter A1 - Potra, Florian T1 - A Refined Gauss-Newton- Mysovskii Theorem. N2 - The present paper contains a generalization of a refinement of the Newton- Mysovskii theorem, recently obtained by the authors, to the case of Gauss-Newton procedures for solving nonlinear least-squares problems with full Jacobians. Invariant sufficient conditions are given that ensure the convergence of the Gauss-Newton iterates towards a solution of the problem, as well as the uniqueness of that solution in an explicitely defined neighborhood. It is shown by a counter- example that the results do not carry over to the rank deficient case. T3 - ZIB-Report - SC-91-04 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-541 ER - TY - GEN A1 - Hoppe, Ronald H. W. A1 - Kornhuber, Ralf T1 - Multilevel Preconditioned CG-Iterations for Variational Inequalities. N2 - We consider such variational inequalities which either describe obstacle problems or result from an implicit time discretization of moving boundary problems of two phase Stefan type. Based on a discretization in space by means of continuous, piecewise linear finite elements with respect to a nested hierarchy of triangulations, in both cases we use iterative processes consisting of inner and outer iterations. The outer iterations are either active set strategies or generalized Newton methods while the inner iterations are preconditioned cg- iterations with multilevel preconditioners. T3 - ZIB-Report - SC-91-06 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-561 ER - TY - GEN A1 - Gajewski(Ed.), Herbert A1 - Deuflhard, Peter A1 - Markowich(ed.), A. T1 - Tagung NUMSIM '91 5.-08. Mai 1991. Collected Abstracts and Papers. N2 - NUMSIM 91 fand vom 6.-08. Mai 1991 im Wissenschaftlichen Kommunikations- & Konferenz- Zentrum der Humboldt-Universität in Gosen bei Berlin statt. Im Mittelpunkt dieses Seminars standen Probleme der numerischen Simulation von Ladungstransport- und Technologieprozessen der Mikro- und Optoelektronik. Es führte Spezialisten der physikalischen Modellierung, der numerischen Mathematik und mathematischen Analysis sowie Nutzer von Simulationsprogrammen aus dem deutschsprachigen Raum zusammen. Ziel war die Vermittlung und der Austausch von Erfahrungen, die Diskussion gemeinsamer Aufgabenstellungen und Projekte und - in Anbetracht der deutschen Vereinigung nicht zuletzt - das gegenseitige Kennen- und Verstehenlernen. T3 - ZIB-Report - TR-91-08 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4838 ER - TY - GEN A1 - Bornemann, Folkmar A. T1 - An Adaptive Multilevel Approach to Parabolic Equations in Two Space Dimensions. N2 - A new adaptive multilevel approach for linear partial differential equations is presented, which is able to handle complicated space geometries, discontinuous coefficients, inconsistent initial data. Discretization in time first (Rothe's method) with order and stepsize control is perturbed by an adaptive finite element discretization of the elliptic subproblems, whose errors are controlled independently. Thus the high standards of solving adaptively ordinary differential equations and elliptic boundary value problems are combined. A theory of time discretization in Hilbert space is developed which yields to an optimal variable order method based on a multiplicative error correction. The problem of an efficient solution of the singularly perturbed elliptic subproblems and the problem of error estimation for them can be uniquely solved within the framework of preconditioning. A Multilevel nodal basis preconditioner is derived, which allows the use of highly nonuniform triangulations. Implementation issues are discussed in detail. Numerous numerical examples in one and two space dimensions clearly show the significant perspectives opened by the new algorithmic approach. Finally an application of the method is given in the area of hyperthermia, a recent clinical method for cancer therapy. T3 - ZIB-Report - TR-91-07 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4821 ER - TY - GEN A1 - Fiedler, Bernold A1 - Scheurle, J. T1 - Discretization of Homoclinic Orbits, Rapid Forcing and "Invisible" Chaos. N2 - One-step discretizations of order $p$ and step size $\varepsilon$ of ordinary differential equations can be viewed as time-$\varepsilon$ maps of \begin{displaymath} \dot{x} (t) = f(\lambda ,x(t)) + \varepsilon^p g(\varepsilon,\lambda,t/\varepsilon,x(t)), x \in R^N,\lambda \in R, \end{displaymath} where $g$ has period $\varepsilon$ in $t$. This is a rapidly forced nonautonomous system. We study the behavior of a homoclinit orbit $\Gamma$ for $\varepsilon = 0, \lambda =0$, under discretization. Under generic assumptions we show that $\Gamma$ becomes transverse for positive $\in$. The transversality effects are estimated from above to be exponentially small in $\in$. For example, the length $l(\varepsilon$) of the parameter interval of $\lambda$ for which $\Gamma$ persists can be estimated by \begin{displaymath} l(\varepsilon)\le Cexp(-2\pi\eta/\varepsilon), \end{displaymath} where $C,\eta$ are positive constants. The coefficient $\eta$ is related to the minimal distance from the real axis of the poles of $\Gamma(t)$ in the complex time domain. Likewise, the region where complicated, "chaotic" dynamics prevail is estimated to be exponentially small, provided $x \in R^2$ and the saddle quantity of the associated equilibrium is nonzero. Our results are visualized by high precision numerical experiments. The experiments show that, due to exponential smallness, homoclinic transversality becomes pratically invisible under normal circumstances, already for only moderately small step size. {\bf Keywords:} Homoclinic orbit, ordinary differential equations, discretization, transversality, averaging, exponential smallness, chaos. {\bf Subject Classifications:} (AMS): 34C15, 34C35, 58F14, 65L60 T3 - ZIB-Report - SC-91-05 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-552 ER - TY - GEN A1 - Lang, Jens A1 - Walter, Artur T1 - An Adaptive Discontinuous Finite Element Method for the Transport Equation. N2 - In this paper we introduce a discontinuous finite element method. In our approach, it is possible to combine the advantages of finite element and finite difference methods. The main ingredients are numerical flux approximation and local orthogonal basis functions. The scheme is defined on arbitrary triangulations and can be easily extended to nonlinear problems. Two different error indicators are derived. Especially the second one is closely connected to our approach and able to handle arbitrary variing flow directions. Numerical results are given for boundary value problems in two dimensions. They demonstrate the performance of the scheme, combined with the two error indicators. {\bf Key words:} neutron transport equation, discontinuous finite element, adaptive grid refinement. {\bf Subject classifications:} AMS(MOS) 65N30, 65M15. T3 - ZIB-Report - SC-91-07 KW - neutron transport equation KW - discontinuous finite element KW - adaptive grid refinement Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-579 ER - TY - GEN A1 - Schöffel, Klaus T1 - Computational Chemistry Software for CRAY X-MP/2. N2 - The following report intends to provide a survey over the computational chemistry molecular structure software installed on CRAY X-MP/24 at ZIB. It shows what kind of problems can be tackled with the existing chemistry software, which covers a wide range of ab initio, semiempirical, molecular mechanics, and dynamics applications. T3 - ZIB-Report - TR-91-06 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4816 ER - TY - GEN A1 - Griewank, Andreas T1 - Sequential Evaluation of Adjoints and Higher Derivative Vectors by Overloading and Reverse Accumulation. N2 - Most nonlinear computations require the evaluation of first and higher derivatives of vector functions defined by computer programs. It is shown here how vectors of such partial derivatives can be obtained automatically and efficiently if the computer language allows overloading (as is or will be the case for C++, PASCAL-XSC, FORTRAN90, and other modern languages). Here, overloading facilitates the extension of arithmetic operations and univariate functions from real or complex arguments to truncated Taylor-series (or other user- defined types), and it generates instructions for the subsequent evaluation of adjoints. Similar effects can be achieved by precompilation of FORTRAN77 programs. The proposed differentiation algorithm yields gradients and higher derivatives at a small multiple of the run-time and RAM requirement of the original function evaluation program. {\bf Keywords:} Automatic Differentiation, Chain Rule, Overloading, Taylor Coefficients, Gradients, Hessians, Reverse Accumulation, Adjoint Equations. {\bf Abbreviated title:} Automatic Differentiation by Overloading. T3 - ZIB-Report - SC-91-03 KW - Automatic Differentiation KW - Chain Rule KW - Overloading KW - Taylor Coefficients KW - Gradients KW - Hessians KW - Reverse Accumulation KW - Adjoint Equations Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-539 ER - TY - GEN T1 - Jahresbericht 1990 T3 - Jahresbericht - 1990 KW - Annual Report Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-9831 SN - 0934-5892 ER - TY - GEN A1 - Gottschewski, Jürgen T1 - Supercomputing During the German Reunification. N2 - The success of the German Reunification requires a unique migration from the traditional Eastern Block system to a modern free society. The successful migration must include a shift from a low state of the art in technology to the high level of science and technology present in the Western Community. The role of supercomputing in this transition will include leading the process of organizational change in research and development. The paper contains a short overview on the current state of this evolution, on the currently used political and administrative tools to keep this evolution on track, and on the influence on the Konrad-Zuse-Zentrum, which is, as a Berlin institution, in the geographical and cultural center of this transformation. T3 - ZIB-Report - TR-91-05 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4803 ER - TY - GEN A1 - Kornhuber, Ralf A1 - Roitzsch, Rainer T1 - Self Adaptive Computation of the Breakdown Voltage of Planar pn-Junctions with Multistep Field Plates. N2 - The breakdown voltage highly depends on the electric field in the depletion area whose computation is the most time consuming part of the simulation. We present a self adaptive Finite Element Method which reduces dramatically the required computation time compared to usual Finite Difference Methods. A numerical example illustrates the efficiency and reliability of the algorithm. T3 - ZIB-Report - SC-91-02 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-525 ER - TY - GEN A1 - Lang, Jens T1 - An Adaptive Finite Element Method for Convection-Diffusion Problems by Interpolation Techniques. N2 - For adaptive solution of convection- difussion problems with the streamline-diffusion finite element method, an error estimator based on interpolation techniques is developed. It can be shown that for correctness of this error estimator a restriction of the maximum angle is to be sufficient. Compared to usual methods, the adaptive process leads to more accurate solutions at much less computational cost. Numerical tests are enclosed. {\bf Keywords: } Adaptive finite elements, convection- diffusion equation, internal and boundary layers, streamline-diffusion. {\bf Subject Classifications:} AMS(MOS): 65N15, 65N30 T3 - ZIB-Report - TR-91-04 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4795 ER - TY - GEN A1 - Noelle, Sebastian W. C. T1 - On the Limits of Operator Splitting: Numerical Experiments for the Complex Burgers Equation. N2 - A high resolution Godunov-type method is developed and applied to a two-dimensional version of the Burgers equation. On one-dimensional testproblems the scheme is second order accurate for smooth flow and resolves discontinuities sharply. On two-dimensional problems, however, large numerical oscillations are introduced when discontinuities traveling obliquely to the grid directions are present. These oscillations are caused by a Strang- type operator splitting, and cannot be eliminated using flux limiters or small timesteps. {\bf Keywords:} High resolution methods, Godunovs method, operator splitting, complex Burgers equation. {\bf AMS:} 35: A40, L45, L60, L65, L, 67. 65: M05, M10. T3 - ZIB-Report - TR-91-03 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4789 ER - TY - GEN A1 - Lügger, Joachim A1 - Dalitz, Wolfgang T1 - Verteilung mathematischer Software mittels elektronischer Netze: Die elektronische Softwarebibliothek eLib. N2 - Mit zunehmender Dezentralisierung im Bereich des Wissenschaftlichen Rechnens steigt auch der Bedarf an Softwarequellen qualitativ hochwertiger mathematischer Algorithmen. Für den Software- und Informationsaustausch wird die Einrichtung öffentlicher, elektronischer Softwarebibliotheken vorgeschlagen. Am Beispiel der eLib des ZIB werden die notwendigen technischen Konzepte, die u.a. auch in den heutigen Gegebenheiten der technischen Kommunikation begründet sind, diskutiert. Die eLib ist fenster- und menüorientiert aufgebaut und folglich im Dialog sehr einfach zu benutzen. Sie kann aber auch mit ihrer Kommandosprache allein auf dem Wege über elektronische Post abgefragt werden. Über die eLib sind, nach dem Verbund mit den Netzbibliotheken NetLib und Reduce-NetLib jetzt etwa 6000 mathematische Programme und allgemeine Dokumente verfügbar. Das Dienstleistungsangebot der eLib setzt Beiträge aus einer aktiven Benutzergemeinschaft voraus. Es besteht die Hoffnung, daß durch den öffentlichen Softwareaustausch sich Entwickler und Anwender wechselseitig in ihrer Arbeit befruchten. T3 - ZIB-Report - TR-91-02 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4774 ER - TY - GEN A1 - Bornemann, Folkmar A. A1 - Erdmann, Bodo A1 - Roitzsch, Rainer T1 - KASKADE - Numerical Experiments. N2 - The C-implementation of KASKADE, an adaptive solver for linear elliptic differential equations in 2D, is object of a set of numerical experiments to analyze the use of resources (time and memory) with respect to numerical accuracy. We study the dependency of the reliability, robustness, and efficiency of the program from the parameters controlling the algorithm. T3 - ZIB-Report - TR-91-01 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4762 ER - TY - GEN A1 - Bornemann, Folkmar A. T1 - An Adaptive Multilevel Approach to Parabolic Equations III. N2 - Part III of the paper is devoted to the construction of an adaptive FEM solver in two spatial dimensions, which is able to handle the singularly perturbed elliptic problems arising from discretization in time. The problems of error estimation and multilevel iterative solution of the linear systems - both uniformly well behaved with respect to the time step - can be solved simultaneously within the framework of preconditioning. A multilevel nodal basis preconditioner able to handle highly nonuniform meshes is derived. As a numerical example an application of the method to the bioheat-transfer equation is included. {\bf AMS CLASSIFICATION:} 65F10, 65F35, 65M50, 65M60, 65N30. T3 - ZIB-Report - SC-91-01 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-519 ER - TY - GEN A1 - Walter, Artur T1 - Sparse Secant Methods for the Iterative Solution of Large Nonsymmetric Linear Systems. N2 - A variety of secant methods has been revisited in view of the construction of iterative solvers for large nonsymmetric linear systems $ Ax = b $ stemming from the discretization of convection diffusion equations. In the first section, we tried to approximate $ A ^{-1} $ directly. Since the sparsity structure of A- is not known, additional storage vectors are needed during the iteration. In the next section, an incomplete factorization $ LU $ of $ A $ is the starting point and we tried to improve this easy invertible approximation of $ A $. The update is constructed in such a way that the sparsity structure of $ L $ and $ U $ is maintained. Two different sparsity preserving updates are investigated from theoretical and practical point of view. Numerical experiments on discretized PDEs of convection diffusion type in 2- D with internal layers and on "arbitrary" matrices with symmetric sparsity structure are given. {\bf Key words:} nonsymmetric linear system, sparse secant method, Broyden's method, incomplete factorization. T3 - ZIB-Report - SC-90-20 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-483 ER - TY - GEN A1 - Wulkow, Michael A1 - Deuflhard, Peter T1 - Towards an Efficient Computational Treatment of Heterogeneous Polymer Reactions. N2 - The discrete Galerkin method developed by the authors has turned out to be an efficient tool for the computational treatment of very large scale ODE systems arising in polyreaction kinetics. Up to now, this approach has been worked out in detail for homogeneous polymer reactions. The present paper deals with one line of possible extensions of the method to the case of so-called heterogeneous processes, which may appear e. g. in smog reactions. The associated mathematical models involve reaction coefficients depending on the chain length of the reacting polymer. The herein suggested extension is worked out in some detail on the basis of the earlier paper. In addition, a numerical example describing polymer degradation is included. T3 - ZIB-Report - SC-90-01 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-292 ER - TY - GEN A1 - Giloi, Wolfgang K. T1 - Konrad Zuses Plankalkül als Vorläufer moderner Programmier-modelle. (Vortrag) T3 - ZIB-Report - TR-90-13 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4759 ER - TY - GEN A1 - Nowak, Ulrich A1 - Weimann, Lutz T1 - GIANT - A Software Package for the Numerical Solution of Very Large Systems of Highly Nonlinear Equations. N2 - This report presents the final realization and implementation of a global inexact Newton method proposed by Deuflhard. In order to create a complete piece of software, a recently developed iterative solver (program GBIT) due to Deuflhard, Freund, Walter is adapted and serves as the standard iterative linear solver. Alternative linear iterative solvers may be adapted as well, e.g. the widely distributed code GMRES. The new software package GIANT (Global Inexact Affine Invariant Newton Techniques) allows an efficient and robust numerical solution of very large scale highly nonlinear systems. Due to the user friendly interface and its modular design, the software package is open for an easy adaptation to specific problems. Numerical experiments for some selected problems illustrate performance and usage of the package. T3 - ZIB-Report - TR-90-11 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4746 ER - TY - GEN A1 - Caprasse, H. A1 - Demaret, J. A1 - Gatermann, Karin A1 - Melenk, Herbert T1 - Power-Law Type Solutions of Fourth-Order Gravity N2 - We study the power-law type solutions of the fourth order field equations derived from a generic quadratic Lagrangian density in the case of multidimensional Bianchi I cosmological models. All the solutions of the system of algebraic equations have been found, using computer algebra, from a search of the Groebner bases associated to it. While, in space dimension $ d = 3 $ , the Einsteinian Kasner metric is still the most general power-law type solution, for $ d > 3 $ , no solution, other than the Minkowski space-time, is common to the three systems of equations associated with the three contributions to the Lagrangian density. In the case of a pure Riemann-squared contribution (suggested by a recent calculation of the effective action for the heterotic string), the possibility exists to realize a splitting of the $ d $-dimensional space into a ( $ d - 3 $)-dimensional internal space and a physical 3- dimensional space, the latter expanding in time as a power bigger than 2 (about 4.5 when $ d = 9 $). T3 - ZIB-Report - SC-90-18 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-466 ER - TY - GEN A1 - Wang, Daoliu T1 - Symplectic Difference Schemes for Perturbed Hamiltonian Systems. N2 - In this paper we consider symplectic difference schemes for perturbed Hamiltonian systems of integrable ones, which can cover many important problems. Symplectic difference schemes for general Hamiltonian systems can also be used to these problems. But the perturbation property has not been paid proper attention to, which is important in the method proposed here. Numerical simulation shows that, for this method the time step size can be taken quite large and the qualitative property , such as preserving invariant tori, is also better than usual symplectic difference schemes. T3 - ZIB-Report - SC-90-17 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-456 ER - TY - GEN A1 - Hege, Hans-Christian A1 - Stüben, Hinnerk T1 - Vectorization and Parallelization of Irregular Problems via Graph coloring. N2 - Efficient implementations of irregular problems on vector and parallel architectures are generally hard to realize. An important class of problems are Gauß-Seidel iteration schemes applied to irregular data sets. The unstructured data dependences arising there prevent restructuring compilers from generating efficient code for vector or parallel machines. It is shown, how to structure the data dependences by decomposing the underlying data set using graph coloring techniques and by specifying a particular execution order already on the algorithm level. Methods to master the irregularities originating from different types of tasks are proposed. An application is given and some open issues and future developments are discussed. T3 - ZIB-Report - SC-90-16 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-443 ER - TY - GEN A1 - Roitzsch, Rainer A1 - Kornhuber, Ralf T1 - BOXES - a Program to Generate Triangulations from a Rectangular Domain Description. N2 - BOXES computes a triangulation from a 2D domain description which consists of an arbitrary set of rectangles. Each rectangle may have attributes to control the triangulating process, define subdomain classes, or specify boundary conditions. The output of the program can be used as a coarse grid for KASKADE or one of its variants. Additional features are extensive checking of the user input, graphical display, and simple editing. T3 - ZIB-Report - TR-90-09 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4738 ER - TY - GEN A1 - Ackermann, Jörg A1 - Wulkow, Michael T1 - MACRON - A Program Package for Macromalecular Reaction Kinetics. N2 - This paper presents the new program package MACRON for the simulation of macromolecular kinetics including standard chemical reactions. Such problems lead to countable (possibly) infinite systems of ordinary differential equations (CODE's), which are numerically treated by the so-called discrete Galerkin method here. By a chemical compiler the required analytical preprocessing is performed, such that the complete reaction system, standard kinetics as well as macromolecular reactions, can be entered in the chemical formalism. Typical macromolecular reaction steps are chain addition, termination, chain transfer and degradation (cracking). In order to ensure efficiency and reliability, high sophisticated numerical routines are built within the package. MACRON can be used without a detailed knowledge of the used numerical methods. As an illustration the application of MACRON to some realistic problems is presented. T3 - ZIB-Report - SC-90-14 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-423 ER - TY - GEN A1 - Wulkow, Michael A1 - Ackermann, Jörg T1 - Numerical Treatment of Polyreactions - Recent Developments. N2 - The mathematical modeling of macromolecular reactions leads to countable (possibly infinite) systems of ordinary differential equations (CODE's). This paper reviews two recent developments of the so-called discrete Galerkin method, which has been developed for the numerical treatment of countable systems, which arise e.g. in polymer chemistry. The first approach can be considered as a method of lines with moving basis functions and has been implemented recently in the program package MACRON. The second type of the Galerkin method is characterized by a so-called outer time discretization of the complete problem and an appropriate and efficient solution of the arising subproblems. This method is realized in the research code CODEX. T3 - ZIB-Report - SC-90-15 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-439 ER - TY - GEN A1 - Wulkow, Michael T1 - Numerical Treatment of Countable Systems of Ordinary Differential Equations. N2 - Countable systems of ordinary differential equations appear frequently in chemistry, physics, biology and medicine. They can be considered as ordinary differential equations in sequence spaces. In this work, a full adaptive algorithm for the computational treatment of such systems is developed. The method combines time discretization with extrapolation in Hilbert spaces with a discrete Galerkin approach as discretization of the stationary subproblems. The Galerkin method is based on orthogonal functions of a discrete variable , which are generated by certain weight functions. A theory of countable systems in the associated weighted sequence spaces is developed as well as a theory of the Galerkin method. The Galerkin equations can be assembled either by use of analytical properties of the orthogonal functions or numerically by a multilevel summation algorithm. The resulting algorithm CODEX is applied to many examples of technological interest, in particular from polymer chemistry. T3 - ZIB-Report - TR-90-08 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4728 ER - TY - GEN A1 - Bornemann, Folkmar A. T1 - An Adaptive Multilevel Approach to Parabolic Equations II. N2 - In continuation of part I this paper develops a variable-order time discretization in Hilbert space based on a multiplicative error correction. Matching of time and space errors as explained in part I allows to construct an adaptive multilevel discretization of the parabolic problem. In contrast to the extrapolation method in time, which has been used in part I, the new time discretization allows to separate space and time errors and further to solve fewer elliptic subproblems with less effort, which is essential in view of the application to space dimension greater than one. Numerical examples for space dimension one are included which clearly indicate the improvement. T3 - ZIB-Report - SC-90-13 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-410 ER - TY - GEN A1 - Walter, Artur T1 - Improvement of Incomplete Factorizations by a Sparse Secant Method. N2 - In the present paper, the improvement of an incomplete factorization of a non-symmetric matrix A is discussed. Starting from the ideas of sparsity preserving quasi-Newton methods, an algorithm is developed which improves the approximation of A by the incomplete factorization maintaining the sparsity structure of the matrices. No renumbering of the unknowns or the admittance of additional fill-in is necessary. The linear convergence of the algorithm is proved under the assumption, that $ L $ and $ U $* have the same sparsity structure and an incomplete factorization with some reasonable approximation property exits. In combination with this algorithm, the method of incomplete factorization and its several modifications are applicable to a wider class of problems with improved convergence qualities. This is shown by a numerical example. {\bf Key Words:} non-symmetric linear system, sparse secant method, incomplete factorization. AMS(MOS) {\bf Subject Classifications:} 65F10, 65N20, 65N30. T3 - ZIB-Report - SC-90-12 KW - non-symmetric linear systems KW - sparse secant method KW - incomplete factorization Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-401 ER - TY - GEN A1 - Kornhuber, Ralf A1 - Roitzsch, Rainer T1 - Self Adaptive FEM Simulation of Reverse Biased pn-Junctions. N2 - The potential distribution of reverse biased pn-junctions can be described by a double obstacle problem for the Laplacian. This problem is solved by a self adaptive Finite Element Method involving automatic termination criteria for the iterative solver, local error estimation and local mesh refinement. Special attention is paid to the efficient resolution of the geometries typically arising in semiconductor device simulation. The algorithm is applied to a reverse biased pn- junction with multi-step field plate and stop- electrode to illustrate its efficiency and reliability. T3 - ZIB-Report - SC-90-10 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-389 ER - TY - GEN A1 - Gatermann, Karin A1 - Hohmann, Andreas T1 - Symbolic Exploitation of Symmetry in Numerical Path-following. N2 - Parameter-dependent systems of nonlinear equations with symmetry are treated by a combination of symbolic and numerical computations. In the symbolic part of the algorithm the complete analysis of the symmetry occurs, and it is here where symmetrical normal forms, symmetry reduced systems, and block diagonal Jacobians are computed. Given a particular problem, the symbolic algorithm can create and compute through the list of possible bifurcations thereby forming a so-called tree of decisions correlated to the different types of symmetry breaking bifurcation points. The remaining part of the algorithm deals with the numerical pathfollowing based on the implicit reparametrisation as suggested and worked out by Deuflhard/Fiedler/Kunkel. The symmetry preserving bifurcation points are computed using recently developed augmented systems incorporating the use of symmetry. {\bf Keywords:} pathfollowing, mixed symbolic-numeric algorithm, parameter-dependent, nonlinear systems, linear representations. T3 - ZIB-Report - SC-90-11 KW - pathfollowing KW - mixed symbolic-numeric algorithm KW - parameter dependent KW - nonlinear systems KW - linear representations Y1 - 1990 N1 - No preprint available; the article appeared in: IMPACT Comput. Sci. Eng. 3 (1991) pp. 330-365 ER - TY - GEN A1 - Deuflhard, Peter T1 - Zuses Werk weiterdenken. (Vortrag) T3 - ZIB-Report - TR-90-07 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4717 ER - TY - GEN A1 - Deuflhard, Peter A1 - Potra, Florian T1 - Asymptotic Mesh Independence of Newton-Galerkin Methods via a Refined Mysovskii Theorem. N2 - The paper presents a theoretical characterization of the often observed asymptotic mesh independence of Newton's method, which means that Newton's method applied to discretized operator equations behaves essentially the same for all sufficiently fine discretizations. The theory does not need any uniform Lipschitz assumptions that were necessary in comparable earlier treatments. The refined Newton-Mysovskii theorem, which will be of interest in a wider context, gives both existence and uniqueness of the solution and quadratic convergence for sufficiently good starting points. Attention is restricted to Galerkin approximations even though similar results should hold for finite difference methods - but corresponding proofs would certainly be more technical. As an illustrative example, adaptive 1-D collocation methods are discussed. T3 - ZIB-Report - SC-90-09 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-379 ER - TY - GEN A1 - Chevalier, C. A1 - Melenk, Herbert A1 - Warnatz, J. T1 - Automatic Generation of Reaction Mechanisms for Description of Oxidation of Higher Hydrocarbons. N2 - Oxidation mechanisms even for rather simple hydrocarbons like heptane consist due to the occurrence of many isomeric structures of thousands of reactions of hundreds of species. The automatic generation of these reaction mechanisms using artificial intelligence means is described. Results are presented for n-heptane-air mixtures, where a hand-written reaction mechanism tested against experimental data is available. T3 - ZIB-Report - SC-90-08 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-361 ER - TY - GEN A1 - Deuflhard, Peter A1 - Hohmann, Andreas T1 - Einführung in die Numerische Mathematik. N2 - Diese Vorlesung ist eine Einführung in die Numerische Mathematik als einem der drei Bereiche (neben den Naturwissenschaften und der Informatik) des am besten mit dem Begriff Scientific Computing charakterisierten Forschungsgebietes. Aufgabe dieser relativ jungen Wissenschaft ist die Entwicklung von Rechenverfahren für Probleme aus den Naturwissenschaften mit Hilfe mathematischer Methoden. T3 - ZIB-Report - TR-90-06 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4704 ER - TY - GEN A1 - Deuflhard, Peter A1 - Nowak, Ulrich A1 - Wulkow, Michael T1 - Recent Developments in Chemical Computing. N2 - The paper surveys three aspects of chemical computing, which seem to play a role in recent developments. First, extrapolation methods for the numerical treatment of differential- algebraic equations are introduced. The associated extrapolation code LIMEX has reached a certain level of sophistication, which makes it a real competitor to the elsewhere widely used multi-step code DASSL of Petzold. Second, adaptive methods of lines for partial differential equations such as those arising in combustion problems are treated. Both static and dynamic regridding techniques are discussed in some detail. Finally, some new ideas about the treatment of the kinetic equations arising from polymer reactions are presented. The new feature of the suggested approach is the application of a Galerkin procedure using sets of orthogonal polynomials over a discrete variable (which, of course, in the case of polymer reactions is the polymer degree). The new approach may open the door to a new reliable low dimensional treatment of complex polymer reactions. T3 - ZIB-Report - SC-90-07 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-359 ER - TY - GEN A1 - Deuflhard, Peter A1 - Freund, R. A1 - Walter, Artur T1 - Fast Secant Methods for the Iterative Solution of Large Nonsymmetric Linear Systems. N2 - A family of secant methods based on general rank-1 updates has been revisited in view of the construction of iterative solvers for large non- Hermitian linear systems. As it turns out, both Broydens "good" and "bad" update techniques play a special role - but should be associated with two different line search principles. For Broydens "bad" update technique, a minimum residual principle is natural - thus making it theorectically comparable with a series of well-known algorithms like GMRES. Broydens "good" update technique, however, is shown to be naturally linked with a minimum "next correction" principle - which asymptotically mimics a minimum error principle. The two minimization principles differ significantly for sufficiently large system dimension. Numerical experiments on discretized PDE's of convection diffusion type in 2-D with internal layers give a first impression of the possible power of the derived "good" Broyden variant. {\bf Key Words:} nonsymmetric linear system, secant method, rank-1 update, Broydens method, line search, GMRES. AMS(MOS) {\bf Subject Classifications:} 65F10, 65N20. T3 - ZIB-Report - SC-90-05 KW - nonsymmetric linear systems KW - secant method KW - rank-1 update KW - Broyden's method KW - line search KW - GMRES Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-330 ER - TY - GEN A1 - Wang, Daoliu T1 - On Symplectic Difference Schemes for Hamiltonian Systems. N2 - Symplectic difference schemes have been shown to be a right formalism for numerical computation of Hamiltonian systems. They are suitable to long time computation and of good qualitative properties. These properties are ensured by the fact that a symplectic difference scheme approximating to a time-independent Hamiltonian system can be regarded as a perturbed time-dependent Hamiltonian system of the original one. That is, a solution of a symplectic difference scheme is a solution of a certain perturbed time dependent Hamiltonian system evaluated at discrete (time) points. This is the main result of the paper. Moreover, linear symplectic difference schemes approximating to a linear time-independent Hamiltonian system can be regarded as a perturbed time-independent Hamiltonian system. So it has all properties that a linear Hamiltonian system has. Based on these results, stochastic webs and chaos in symplectic difference schemes are also discussed. They will appear in numerical simulation for Hamiltonian systems, even with one degree of freedom. T3 - ZIB-Report - SC-90-06 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-345 ER - TY - GEN T1 - Jahresbericht 1989 T3 - Jahresbericht - 1989 KW - Annual Report Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-9847 SN - 0934-5892 ER - TY - GEN A1 - Bornemann, Folkmar A. T1 - An Adaptive Multilevel Approach to Parabolic Equations I. General Theory & 1D-Implementation. N2 - A new adaptive multilevel approach for parabolic PDE's is presented. Full adaptivity of the algorithm is realized by combining multilevel time discretization, better known as extrapolation methods, and multilevel finite element space discretization. In the theoretical part of the paper the existence of asymptotic expansions in terms of time-steps for single-step methods in Hilbert space is established. Finite element approximation then leads to perturbed expansions, whose perturbations, however, can be pushed below a necessary level by means of an adaptive grid control. The theoretical presentation is independent of space dimension. In this part I of the paper details of the algorithm and numerical examples are given for the 1D case only. The numerical results clearly show the significant perspectives opened by the new algorithmic approach. T3 - ZIB-Report - SC-90-04 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-320 ER - TY - GEN A1 - Hege, Hans-Christian T1 - Datenabhängigkeitsanalyse und Programmtransformationen auf CRAY-Rechnern mit dem Fortran-Präprozessor fpp N2 - The FORTRAN preprocessor fpp in the newly introduced Autotasking System of CRAY Research allows automatic vectorization and parallelization on basis of a data dependence analysis. An introduction into data dependence analysis is given, showing how data dependence graphs unveil opportunities for program transformations like vectorization and concurrentization. The report contains a complete description of the preprocessors functionality, its options and directives for increasing the effectiveness of the dependence analyzer and steering the code transformations. Finally, some advice is given for the practical use of fpp on CRAY computers. T3 - ZIB-Report - TR-90-04 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4685 ER - TY - GEN T1 - CRAY-Handbuch. Einführung in die Benutzung der CRAY X-MP unter UNICOS 5.1 N2 - Das CRAY-Handbuch des ZIB beschreibt nach einer grundlegenden bersicht detailliert die Handhabung von UNICOS-Dateien, die Verwendung von UNICOS-Kommandos und schließlich die (für Benutzer allein zugelassene) Ausführung von Batchjobs auf der CRAY. Neben der Erreichbarkeit der CRAY über das TCP/IP-Netz werden die Compiler FORTRAN, PASCAL und C beschrieben und schließlich auf die Optimierung der Rechenzeit ebenso wie auf die Fehlersuche eingegangen. Ein umfangreicher Anhang führt bereitgestellte Programmpakete, weiterführende Literatur, Zahlenbereiche und Zeichendarstellung an der CRAY auf. T3 - ZIB-Report - TR-90-03 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4679 ER - TY - GEN A1 - Gatermann, Karin T1 - Symbolic solution of polynomial equation systems with symmetry. N2 - Systems of polynomial equations often have symmetry. The Buchberger algorithm which may be used for the solution ignores this symmetry. It is restricted to moderate problems unless factorizing polynomials are found leading to several smaller systems. Therefore two methods are presented which use the symmetry to find factorizing polynomials, decompose the ideal and thus decrease the complexitiy of the system a lot. In a first approach projections determine factorizing polynomials as input for the solution process, if the group contains reflections with respect to a hyperplane. Two different ways are described for the symmetric group Sm and the dihedral group Dm. While for Sm subsystems are ignored if they have the same zeros modulo G as another subsystem, for the dihedral group Dm polynomials with more than two factors are generated with the help of the theory of linear representations and restrictions are used as well. These decomposition algorithms are independent of the finally used solution technique. We used the REDUCE package Groebner to solve examples from CAPRASSE, DEMARET and NOONBURG which illustrate the efficiency of our REDUCE program. A short introduction to the theory of linear representations is given. In a second approach problems of another class are transformed such that more factors are found during the computation; these transformations are based on the theory of linear representations. Examples illustrate these approaches. The range of solvable problems is enlarged significantly. T3 - ZIB-Report - SC-90-03 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-311 ER - TY - GEN A1 - Deuflhard, Peter T1 - Global Inexact Newton Methods for Very Large Scale Nonlinear Problems. N2 - Newton methods for nonlinear problems are known to require the solution of a sequence of linear problems of the same type. For very large scale problems, as understood herein, the arising linear systems can only be solved by iterative methods. Then Newtons iteration appears as outer iteration. The question of interest will be to control the accuracy of the inner iteration such that the convergence speed of Newtons method is preserved. The purpose of the paper is to combine the concept of inexact Newton methods with the concept of the affine invariant exact Newton methods - which is important for problems with ill- conditioned Jacobian matrices (such as typical 2-D or 3-D discretized partial differential equations). T3 - ZIB-Report - SC-90-02 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-303 ER - TY - GEN A1 - Maierhöfer, Gerhard A1 - Skorobohatyj, Georg T1 - Implementierung von parallelen Versionen der Gleichungslöser EULEX und EULSIM auf Transputern. N2 - Im vorliegenden Bericht wird die Parallelisierung zweier numerischer Algorithmen zur Lösung gewöhnlicher Differentialgleichungssysteme 1.Ordnung (explizite und semi-implizite Euler- Diskretisierung und $h$-Extrapolation) beschrieben. Implementiert wurden die Algorithmen mit OCCAM2 unter TDS (Transputer Development System) mit bis zu 4 Transputern T800. Meßwerte für die erreichten Beschleunigungen werden anhand mehrer Beispiele von Differentialgleichungs-systemen angegeben. {\bf Schlüsselwörter:} Adaptive, parallele Systeme; OCCAM2; Transputer; numerische Gleichungslöser; Euler-Diskretisierung; $h$-Extrapolation. T3 - ZIB-Report - TR-90-02 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4666 ER - TY - GEN A1 - Gatermann, Karin T1 - Gruppentheoretische Konstruktion von symmetrischen Kubaturformeln. N2 - $G$-invariant cubature formulas for numerical integration over n-dimensional, $G$- invariant integration regions are computed symbolically. The nodes are the common zeros of some $d$-orthogonal polynomials which build an $H$-basis of an ideal. Approaches for these polynomials depending on parameters are made with the help of the theory of linear representations of a group $G$. This theory is also used for the effective computation of necessary conditions which determines the parameters. Another approach uses invariant theory and gröbner bases. T3 - ZIB-Report - TR-90-01 Y1 - 1990 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4652 ER - TY - GEN A1 - Yserentant, Harry T1 - Two Preconditioners Based on the Multi-Level Splitting of Finite Element Spaces. N2 - The hierarchical basis preconditioner and the recent preconditioner of BRAMBLE, PASCIAK and XU are derived and analyzed within a joint framework. This discussion elucidates the close relationship between both methods. Special care is devoted to highly nonuniform meshes; our theory is based exclusively on local properties like the shape regularity of the finite elements. T3 - ZIB-Report - SC-89-09 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-274 ER - TY - GEN A1 - Hoppe, Ronald H. W. T1 - Numerical Solution of Multicomponent Alloy Solidification by Multi-Grid Techniques. N2 - The solidification of an $ N $-component alloy is described by an initial boundary value problem for a system of degenerate parabolic equations modelling heat conduction and mass diffusion. Discretizing implicitly in time and by piecewise linear finite elements in the space variables, at each time step the solution of a system of quasivariational inequalities is required. For the numerical solution of that system, a multi-grid algorithm is developed by making use of game theoretic concepts and duality arguments from convex analysis. Finally, the efficiency of the algorithm is demonstrated by displaying numerical results for a ternary alloy. T3 - ZIB-Report - SC-89-10 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-285 ER - TY - GEN A1 - Maierhöfer, Gerhard T1 - Ein paralleler adaptiver Algorithmus für die numerische Integration. N2 - Der Bericht ist die Ausarbeitung eines Vortrages, der auf dem Transputer Anwender Treffen (TAT'89) im September 1989 gehalten wurde. Es wird die Parallelisierung und Implementierung eines adaptiven Algorithmus zur numerischen Integration (Romberg Quadratur) beschrieben. Ausgewählte Meßergebnisse sind enthalten. {\bf Keywords:} Numerischer Algorithmus, Romberg Quadratur, paralleler adaptiver Algorithmus, dynamische Lastverteilung und Prozessorzahl, lokaler Speicher, Nearest-Neighbour-Architektur, Transputer, TDS, OCCAM2. T3 - ZIB-Report - SC-89-08 KW - numerischer Algorithmus KW - Romberg Quadratur KW - paralleler adaptiver Algorithmus KW - dynamische Lastverteilung und Prozessorzahl KW - lokaler Speicher KW - Near Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-265 ER - TY - GEN A1 - Maierhöfer, Gerhard A1 - Skorobohatyj, Georg T1 - Implementierung des parallelen TRAPEX auf Transputern. N2 - Der folgende Bericht ist eine Ergänzung des ZIB Technical Report TR 88-05. Entsprechend wird hier nicht auf die grundsätzlichen Fragen der Parallelisierbarkeit des sequentiellen TRAPEX eingegangen. Diese sind im TR 88-5 erörtert, die dort beschriebenen Algorithmen (vertikale und horizontale Parallelisierung) werden auch für die Transputerarchitektur verwendet. Meßergebnisse sind im letzten Teil angefügt. T3 - ZIB-Report - TR-89-08 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4649 ER - TY - GEN A1 - Grammel, Martin A1 - Maierhöfer, Gerhard A1 - Skorobohatyj, Georg T1 - Trapex in POOL-Implementierung eines numerischen Algorithmus in einer parallelen objektorientierten Sprache. N2 - Die Parallelisierung eines vorhandenen sequentiellen Programmes erfolgt im allgemeinen in der Weise, daß es auf Grund einer Analyse der Datenabhängigkeiten in mehrere parallel ausführbare Teile zerlegt wird, die ihrerseits sequentiell ablaufen und untereinander Daten austauschen. Eine parallele Programmiersprache sollte diese Vorgehensweise unterstützen, indem sie Sprachmittel zur Verfügung stellt, mittels derer sich die Modularisierung und Kommunikation bequem formulieren läßt. Nun läßt sich prinzipiell jede vorhandene sequentielle Programmiersprache um solche Konstrukte erweitern; günstiger in Hinblick auf die Strukturierung eines parallelen Programmes erscheint es aber möglicherweise, eine solche Programmiersprache zu verwenden, die die genannten Konzepte als elementare zur Verfügung stellt. Dies ist offensichtlich bei objektorientierten Programmiersprachen der Fall. Im vorliegenden Bericht wird an Hand eines Beispiels - TRAPEX -, das in POOL-T implementiert wurde, untersucht, inwieweit eine objektorientierte Sprache zur effektiven Programmierung eines numerischen Algorithmus geeignet ist. {\bf Key Words:} parallele objektorientierte Programmiersprache, Message Passing, Abstract Data Typing, Klassenhierarchi, Modul-Import-Mechanismus, Client/Server-Prinzip, Load Balancing, adaptive numerische Algorithmen, Ordnungs- und Schrittweitensteuerung, Romberg-Quadratur. T3 - ZIB-Report - TR-90-05 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4696 ER - TY - GEN A1 - Bornemann, Folkmar A. T1 - Adaptive multilevel discretization in time and space for parabolic partial differential equations. N2 - The present paper developes an adaptive multilevel approach for parabolic PDE's - as a first step, for one linear scalar equation. Full adaptivity of the algorithm is conceptually realized by simultaneous multilevel discretization in both time and space. Thus the approach combines multilevel time discretization, better known as extrapolation methods, and multilevel finite element space discretization such as the hierarchical basis method. The algorithmic approach is theoretically backed by careful application of fundamental results from semigroup theory. These results help to establish the existence of asymptotic expansions (in terms of time-steps) in Hilbert space. Finite element approximation then leads to perturbed expansions, whose perturbations, however, can be pushed below a necessary level by means of an adaptive grid control. The arising space grids are not required to satisfy any quasi- uniformity assumption. Even though the theoretical presentation is independent of space dimension details of the algorithm and numerical examples are given for the 1-D case only. For the 1-D elliptic solver, which is used, an error estimator is established, which works uniformly well for a family of elliptic problems. The numerical results clearly show the significant perspectives opened by the new algorithmic approach. T3 - ZIB-Report - TR-89-07 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4638 ER - TY - GEN A1 - Huamo, Wu A1 - Shuli, Yang T1 - MmB-A New Class of Accurate High Resolution Schemes for Conservation Laws in Two Dimensions. N2 - In this paper we present the MmB schemes, which preserve the local maximum and minimum bounds of the initial data in the smallest union of mesh elements of previous time step containing the domain of dependence of the solution on the mesh element with center at point $ P $\ under consideration. In 1-D, the MmB schemes are almost identical with TVD schemes. As well-known, there is no second-order TVD scheme in 2-D, nevertheless, we present here two classes of 2-D second-order accurate MmB-schemes. It is proved that 1-D discrete MmB (or TVD) and 1-D semi-discrete TVD schemes may have second-order accuracy at (nonsonic) critical points, but cannot be of uniformly second-order accurate in the whole neighborhood of the critical points. New accurate high resolution flux limiters are suggested. Numerical results for 1-D and 2-D test problems are given. {\bf Keywords:} Difference scheme, TVD, MmB, flux limiter. T3 - ZIB-Report - SC-89-06 KW - difference scheme KW - TVD KW - MmB KW - flux limiter Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-248 ER - TY - GEN A1 - Budde, Uwe A1 - Wulkow, Michael T1 - Computation of Molecular Weight Distributions for Free Radical Polymerization Systems. N2 - Modeling of free radical polymerization leads to very large and usually stiff systems of ordinary differential equations which cannot be solved directly in an efficient way. This paper presents the application of a new approach called discrete Galerkin method to a realistic example - the polymerization of methyl methacrylate(MMA). The method is characterized by a Galerkin approximation on the basis of orthogonal polynomials of a discrete variable which represents the polymer degree. It allows the efficient computation of solutions of complete kinetic schemes with time- or moment-dependent reaction coefficients by reducing the complexity to a few differential equations. The approximation error can be controlled by an error estimation. In the case of MMA polymerization a reduction of computational effort by a factor of about 25 compared to a standard method can be obtained for the quasi-steady-state approximation of the model. In addition solutions of the instationary kinetic scheme can be easily computed. T3 - ZIB-Report - SC-89-07 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-253 ER - TY - GEN A1 - Huamo, Wu T1 - On the Possible Accuracy of TVD Schemes. N2 - The paper presents a detailed analysis of the possible accuracy available for TVD schemes in one dimension with emphasis to the semi-discrete 1-D TVD schemes. The analysis shows that the widely accepted statement [1] of degeneration of accuracy at critical points for TVD schemes should be corrected. We have theorem: TVD schemes using flux limiters $ \varphi $ of the form [1], [2] may be second-order accurate at critical points if $ \varphi $ (3) + $ \varphi $(-1) = 2, but cannot be uniformly second-order accurate in the whole neighborhood of critical point. If $ \varphi $(1) = 1, then the TVD schemes are second-order accurate in the region of smooth solutions sufficiently far from the critical points. Two ways are suggested to improve the accuracy. Numerical example is given. {\bf Keywords:} Semi-discrete schemes, TVD, flux limiter, degeneration of accuracy. T3 - ZIB-Report - SC-89-03 KW - semi-discrete schemes KW - TVD KW - flux limiter KW - degeneration of accuracy Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-212 ER - TY - GEN A1 - Möller, H. Michael T1 - Multivariate Rational Interpolation: Reconstruction of Rational Functions. N2 - In this paper we consider the problem of reconstructing a multivariate rational function, when only its values at sufficiently many points are known. We use for the reconstruction of bivariate rational functions a bivariate rational interpolation operator investigated by Siemaszko [7] and a new one, compare both by examples in a Computer Algebra system, and present their multivariate generalizations. {\bf Keywords:} Multivariate rational interpolation, reconstruction, symbolic computation. T3 - ZIB-Report - SC-89-04 KW - Multivariante rational interpolation KW - reconstruction KW - symbolic computation Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-227 ER - TY - GEN A1 - Kornhuber, Ralf A1 - Roitzsch, Rainer T1 - On Adaptive Grid Refinement in the Presence of Internal or Boundary Layers. N2 - We propose an anisotropic refinement strategy which is specially designed for the efficient numerical resolution of internal and boundary layers. This strategy is based on the directed refinement of single triangles together with adaptive multilevel grid orientation. It is demonstrated by several numerical examples that compared to usual methods, the new anisotropic refinement ends up in more stable and more accurate solutions at much less computational cost. {\bf Keywords:} Adaptive finite elements, directed refinement, adaptive grid orientation, convection diffusion equation, internal and boundary layers. T3 - ZIB-Report - SC-89-05 KW - adaptive finite elements KW - directed refinement KW - adaptive grid orientation KW - convection diffusion equation KW - internal and boundary layers Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-230 ER - TY - GEN A1 - Melenk, Herbert A1 - Neun, Winfried T1 - Implementation of Portable Standard LISP for the SPARC Processor. N2 - The SPARC processor is a RISC (Reduced Instruction Set Computer) microcomputer, built into the SUN4 workstations. Since RISC processors are very well-suited for LISP processing, the implementation of a dialect of LISP (Portable Standard LISP, PSL) boded well for a great speed-up in comparison with other types of microcomputers. A first approach was done at The RAND Corporation in Santa Monica, which was derived from classical processor types like MC68000 or VAX. At the Konrad- Zuse-Zentrum für Informationstechnik Berlin (ZIB) that initial implementation was redesigned in order to adapt PSL to the specific features of the SPARC processor. The present implementation, in some parts, is very close to Cray PSL version also done in ZIB. Some timing informations are given in the appendix. T3 - ZIB-Report - TR-89-06 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4625 ER - TY - GEN A1 - Roitzsch, Rainer T1 - KASKADE Programmer's Manual. N2 - The Programmers Manual contains the documentation of the C-KASKADE implementation. KASKADE is an adaptive finite element code (see SC ). The interfaces of all modules are described. T3 - ZIB-Report - TR-89-05 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4617 ER - TY - GEN A1 - Roitzsch, Rainer T1 - KASKADE User's Manual. N2 - The Users Manual contains an introduction to the usage of the C-KASKADE implementation. KASKADE is an adaptive finite element code (see SC ). The manual includes a description of all commands, a simple recipe to define a new problem by changing the KASKADE source and a complete example. T3 - ZIB-Report - TR-89-04 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4605 ER - TY - GEN T1 - Jahresbericht 1988 T3 - Jahresbericht - 1988 KW - Annual Report Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-9858 SN - 0934-5892 ER - TY - GEN A1 - Zha, Hongyuan T1 - Restricted Singular Value Decomposition of Matrix Triplets. N2 - In this paper we introduce the concept of restricted singular values (RSV's) of matrix triplets. A theorem concerning the RSV's of a general matrix triplet $ (A,B,C) $, where $ A \in C^{m\times n} $, $B\in C^{m\times p} $ and $ C\in C^{q\times n} $, which is called restricted singular value decomposition (RSVD) of matrix triplets, is derived. This result generalizes the wellknown SVD, GSVD and the recently proposed product induced SVD (PSVD). Connection of RSV's with the problem of determination of matrix rank under restricted perturbation is also discussed. {\bf Keywords:} Matrix rank, singular values, generalized singular values, product induced singular values, restricted singular values, matrix decompositions. T3 - ZIB-Report - SC-89-02 KW - Matrix rank KW - singular values KW - generalized singular values KW - product induced singular values KW - restricted singular values KW - matrix decompositions Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-205 ER - TY - GEN A1 - Deuflhard, Peter T1 - Numerik von Anfangswertmethoden für gewöhnliche Differentialgleichungen. N2 - Das voliegende Skriptum entstand aus einer Vorlesung, die ich im WS 87/88 an der Freien Universität Berlin im Fachbereich Mathematik gehalten habe. Mein ursprüngliches Vorlesungsmanuskript wurde von Herrn F. Bornemann in weiten Teilen überarbeitet, reorganisiert und substantiell ergänzt. Der Inhalt stammt größtenteils aus Originalarbeiten jüngeren Datums. Darüberhinaus finden sich zahlreiche Teile, die aus meiner jahrelangen Beschäftigung mit dem Thema entstanden aber unpubliziert geblieben sind. Das Skriptum erhebt nicht den Anspruch, ein Lehrbuch zu sein. Es war zunächst als Ausarbeitung für meinen studentischen Hörerkreis sowie als internes Arbeitspapier für das ZIB bestimmt. Die Kunde von der bloßen Existenz eines solchen Skriptums hat jedoch zu einer derart regen Nachfrage geführt, daß es hiermit als Technischer Report des ZIB einer breiteren ffentlichkeit zugänglich gemacht werden soll. In der vorliegenden Form richtet es sich in erster Linie an Mathematiker; es soll sich jedoch auch für Naturwissenschaftler und Ingenieure eignen, die sich einen Einblick in den theoretischen und algorithmischen Hintergrund der von ihnen verwendeten wissenschaftlichen Software verschaffen wollen. T3 - ZIB-Report - TR-89-02 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4580 ER - TY - GEN A1 - Walter, Artur T1 - Ein Finite-Element-Verfahren zur numerischen Lösung von Erhaltungsgleichungen. N2 - A new approach for the discretisation of hyperbolic conservation laws via a finite element method is developed and analysed. Appropriate forms of the Eulers equation of gas dynamic are considered to employ the algorithm in a reasonable way for this system of nonlinear equations. Both mathematical and physical stability results are obtained. A main part of the paper is devoted to the convergence proof with energy methods under strong regularity of the solution of a scalar nonlinear conservation law. Some hints on the implementation and numerical results for the calculation of transonic gasflow through a Laval nozzle are given. The necessary amount of numerical work is compared to an established finite difference method and the efficiency of the algorithm is shown. A survey on recent literature about finite element methods for hyperbolic problem is included. T3 - ZIB-Report - TR-89-03 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4597 ER - TY - GEN A1 - Zha, Hongyuan T1 - A Numerical Algorithm for Computing the Restricted Singular Value Decomposition of Matrix Triplets. N2 - This paper presents a numerical algorithm for computing the restricted singular value decomposition of matrix triplets (RSVD). It is shown that one can use unitary transformations to separate the regular part from a general matrix triplet. After preprocessing on the regular part, one obtains a matrix triplet consisting of three upper triangular matrices of the same dimensions. The RSVD of this special matrix triplet is computed using the implicit Kogbetliantz technique. The algorithm is well suited for parallel computation. {\bf Keywords:} Restricted singular values, matrix triplets, unitary transformations, implicit Kogbetliantz technique. T3 - ZIB-Report - SC-89-01 KW - restricted singular values KW - matrix triplets KW - unitary transformations KW - implicit Kogbetliantz technique Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-194 ER - TY - GEN A1 - Chow, Shui-Nee A1 - Deng, Bo A1 - Fiedler, Bernold T1 - Homoclinic Bifurcation at Resonant Eigen-values. T3 - ZIB-Report - SC-88-10 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-181 ER - TY - GEN T1 - CRAY-HANDBUCH. Einführung in die Benutzung der CRAY X-MP unter UNICOS. (veraltet) T3 - ZIB-Report - TR-89-01 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4574 ER - TY - GEN A1 - Kornhuber, Ralf A1 - Roitzsch, Rainer T1 - Adaptive Finite-Elemet- Methoden für konvektionsdominierte Randwertprobleme bei partiellen Differentialgleichungen. N2 - Ausgangspunkt bei der Behandlung konvektiv dominierter, elliptischer Probleme sind die bekannten hierarchischen Finite-Element-Methoden für den rein elliptischen Fall. Als stabile Erweiterung des Standard-Galerkin-Verfahrens wird das Stromlinien-Diffusions-Verfahren durch physikalische Überlegungen motiviert und kurz diskutiert. Anschließend zeigen wir, daß diese Methode erst in Verbindung mit einer hier erstmals vorgestellten lokalen Ausrichtung der Kanten wirksam eingesetzt werden kann. Zusammen mit einer ebenfalls neu entwickelten richtungsorientierten Verfeinerungsstrategie erhält man eine erheblich stabilere, genauere und schnellere Auflösung von Grenzschichten als mit herkömmlichen Methoden. T3 - ZIB-Report - SC-88-09 Y1 - 1989 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-170 ER - TY - GEN A1 - Maierhöfer, Gerhard A1 - Skorobohatyj, Georg T1 - Parallel-TRAPEX (Ein paralleler adaptiver Algorithmus zur numerischen Integration; seine Implementierung für SUPRENUM-artige Architekturen mit SUSI). N2 - This paper describes some ways of transforming a sequential adaptive algorithm for numerical evaluation of an integral (Romberg- Quadrature with polynomial Extrapolation method) to a parallel one, such as have been implementad by the authors. We developed an algorithm which preserves the sequential adaptivity and is capable of running on various architectures, dynamically controlling the number of avtive processors depending on the problem. To study the time behaviour, we used the simulator SUSI(SUprenum SImulatur) which is able to simulate SUPRENUM-like architectures. Results are given in part 3. {\bf Keywords:} Romberg quadrature; numerical integration; parallel adaptive algorithm; SUSI; simulation; SUPRENUM; MIMD-Fortran; computer architecture; granularity; load balacing; master-slave-principle. T3 - ZIB-Report - TR-88-05 Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4560 ER - TY - GEN A1 - Hoppe, Ronald H. W. A1 - Kornhuber, Ralf T1 - Multi-Grid Solution of Two Coupled Stefan Equations Arising in Induction Heating of Large Steel Slabs. N2 - Induction heating of large steel slabs can be described by a coupled system of nonlinear evolution equations of Stefan type representing the temporal and spatial distribution of the induced magnetic field and the generated temperature within the slab. Discretizing these equations implicitly in time and by finite differences in space, at each time step the solution of a system of difference inclusions is required. For the solution of that system two multi-grid algorithms are given which combined with a nested iteration type continuation strategy to proceed in time result in computationally highly efficient schemes for the numerical simulation of the induction heating process. {\bf Keywords:} induction heating, system of two coupled Stefan equations, multi-grid algorithms. {\bf Subject Classification:} AMS(MOS): 35K60, 35R35, 65H10, 65N05, 65N20, 78A25, 78A55. T3 - ZIB-Report - SC-88-08 KW - induction heating KW - system of two coupled Stefan equations KW - multi-grid algorithms Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-160 ER - TY - GEN A1 - Melenk, Herbert A1 - Möller, H. Michael A1 - Neun, Winfried T1 - Symbolic Solution of Large Stationary Chemical Kinetics Problems. N2 - The paper presents a new application of computer algebra to the treatment of steady states of reaction systems. The method is based on the Buchberger algorithm. This algorithm was modified such that it can exploit the special structure of the equations derived from reaction systems, so even large systems can be handled. In contrast to numerical approximation techniques, the algebraic solution gives a complete and definite overview of the solution space and it is even applicable when parameter values are unknown or undetermined. The algorithm, its adaptation to the problem class and its application to selected examples are presented. T3 - ZIB-Report - SC-88-07 Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-157 ER - TY - GEN A1 - Deuflhard, Peter A1 - Wulkow, Michael T1 - Computational Treatment of Polyreaction Kinetics by Orthogonal Polynomials of a Discrete Variable. N2 - The paper presents a new approach to the computational treatment of polyreaction kinetics. This approach is characterized by a Galerkin method based on orthogonal polynomials of a discrete variable, the polymer degree (or chain length). In comparison with the known competing approaches (statistical moment treatment, Galerkin methods for continuous polymer models), the suggested method is shown to avoid the disadvantages and preserve the adventages of either of them. The basic idea of the method is the construction of a discrete inner product associated with a reasonably chosen probability density function. For the so-called Schulz-Flory distribution one thus obtains the discrete Laguerre polynomials, whereas the Poisson distribution leads to the Charlier polynomials. Numerical experiments for selected polyreaction mechanisms illustrate the efficiency of the proposed method. T3 - ZIB-Report - SC-88-06 Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-140 ER - TY - GEN A1 - Deuflhard, Peter A1 - Leinen, P. A1 - Yserentant, Harry T1 - Concepts of an Adaptive Hierarchical Finite Element Code. N2 - The paper presents the mathematical concepts underlying the new adaptive finite element code KASKADE, which, in its present form, applies to linear scalar second-order 2-D elliptic problems on general domains. Starting point for the new development is the recent work on hierarchical finite element bases due to Yserentant (1986). It is shown that this approach permits a flexible balance between iterative solver, local error estimator, and local mesh refinement device - which are the main components of an adaptive PDE code. Without use of standard multigrid techniques, the same kind of computational complexity is achieved - independent of any uniformity restrictions on the applied meshes. In addition, the method is extremely simple and all computations are purely local - making the method particularly attractive in view of parallel computing. The algorithmic approach is illustrated by a well-known critical test problem. {\bf Keywords:} finite elements, hierarchical basis, adaptive mesh refinement, preconditioned conjugate gradient methods. T3 - ZIB-Report - SC-88-05 KW - finite elements KW - hierarchical basis KW - adaptive mesh refinement KW - preconditioned conjugate gradient methods Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-131 ER - TY - GEN A1 - Buhtz, Danuta Anna A1 - Langendorf, Jens A1 - Paetsch, Olaf T1 - GRAZIL - 3D. Ein graphisches Anwendungsprogramm zur Darstellung von Kurven- und Funktionsverläufen im räumlichen Koordinatensystem. N2 - GRAZIL - 3D ist ein interaktives Programmpaket zur graphischen Darstellung von dreidimensionalen Kurvenverläufen. Das Aussehen der räumlichen Ausgabe des Bildes wird durch Angabe des Blickwinkels und der Perspektive festgelegt. T3 - ZIB-Report - TR-88-04 Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4556 ER - TY - GEN T1 - Jahresbericht 1987 T3 - Jahresbericht - 1987 KW - Annual Report Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-9860 SN - 0933-7903 ER - TY - GEN A1 - Melenk, Herbert A1 - Neun, Winfried T1 - Parallel Polynomial Operations in the Buchberger Algorithm. N2 - The Buchberger algorithm is a basic tool for the solution of systems of polynomial equations in an environment of computer algebra applications. A model for overlapped processing of different steps of the algorithm is presented, which uses the data structure of the polynomials (distributive representation) for synchronization. The model can be applied for multi processors with fast access to shared data. It is tested with Cray X-MP multi processors based on a parallel version of Portable Standard Lisp (PSL 3.4). T3 - ZIB-Report - SC-88-04 Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-125 ER - TY - GEN A1 - Melenk, Herbert A1 - Neun, Winfried T1 - REDUCE User s Guide for the CRAY 1 / CRAY X-MP Series Running UNICOS. Version 3.3. N2 - This document describes operating procedures for running REDUCE specific to the CRAY 1 and CRAY X-MP computers running the Operating System UNICOS. The document was derived from the corresponding document for Vax/UNIX prepared by A. C. Hearn and L. R. Seward, The Rand Corporation, Santa Monica, (CP85). T3 - ZIB-Report - TR-88-02 Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4539 ER - TY - GEN A1 - Melenk, Herbert A1 - Neun, Winfried T1 - REDUCE Installation Guide for CRAY 1 / X-MP Systems Running UNICOS. Version 3.3. N2 - This guide describes the CRAY/UNICOS REDUCE distribution tape and the procedures for installing, testing and maintaining REDUCE on a CRAY 1 or CRAY X-MP running UNICOS. This document was derived from the corresponding document for Vax/UNIX prepared by A. C. Hearn and L.R. Seward, The Rand Corporation, Santa Monica, publication CP84. T3 - ZIB-Report - TR-88-03 Y1 - 1988 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4541 ER -