Send mail to:             for the digests or bakeoff
            for comments or help
 Current editor:  Craig Douglas       
Anonymous ftp repository: (

World Wide Web: (access through the CS Department)

Today's editor:  Craig Douglas (

Volume 4, Number 6 (June 30, 1994)

Today's topics:

     Kaskade FEM Programs
     mgghat version 1.1
     Paper by Petr Mejzl\'{\i}k for DDM7
     Preprint by Jung and Ruede
     Copper Mountain 1995
     Bib entries from Konrad-Zuse Zentrum


Date: Thu, 16 Jun 1994 17:11:41 +0200
From: roitzsch@sc.ZIB-Berlin.DE (Rainer Roitzsch)
Subject: Kaskade FEM Programs

    Editor's Note: All of this is now in mgnet/kaskade.

KASKADE version 2.1

This release contains bug fixes, internal reorganizations,
and some extensions which we needed for the time-dependent
code kastio.

New features:

- 3D,

- Iterative solvers gmres, cgs,
- boundary conditions for systems,
  Cauchy function gets additional parameter,
- assemble the convection term,

- temperature plot, color tables,
- colorizing the materials,
- more robust write command (ignoring non-used points etc.),
- simplified input formats.

Internal reorganizations:

- memory management (returning unused memory),
- new events DeleteTriangle and NewGlobSol,
- X11, MacOS, Postscript drivers changed.

The manuals are stored in the directory Manuals.

KASTIO version 2.1

First release (2D/3D) of the adaptive multilevel code for
parabolic equations of Folkmar Bornemann.

Rainer Roitzsch (
Bodo Erdmann (
Jens Lang (

Konrad-Zuse-Zentrum fuer Informationstechnik (ZIB)
Heilbronner Strasse 10
W-1000 Berlin 31


Date: Mon, 20 Jun 94 15:50:32 EDT
From: (William F. Mitchell 301-975-3808)
Subject: mgghat version 1.1


I have an update to MGGHAT.  I just placed it in mgnet/incoming on casper.
The file is mgghat1.1.tgz which is, of course, a gzipped tar file.  Following
is an announcement for the MGNET digest.


    Editor's Note: Per Bill's later request, version 1.0 has been deleted.
    -------------  The new files are all in mgnet/mgghat.

I am pleased to announce the availability of Version 1.1 of MGGHAT.

MGGHAT (MultiGrid Galerkin Hierarchical Adaptive Triangles) is a
FORTRAN program for the solution of second order linear elliptic
partial differential equations of the form

    d    du      d    du
 - -- (p --)  - -- (q --) + ru = f   in D
   dx    dx     dy    dy

with boundary conditions of the form

                             u = g   on d1

        du  dy      du dx
      p --  --  - q -- -- + cu = g   on d2
        dx  ds      dy ds

where p>0, q>0, r, f, c and g are functions of x and y, D
is a polygonal domain in R^2 (possibly with holes), d1 U d2 is
the boundary of D, and d/ds is differentiation with respect to a
counterclockwise parameterization of the boundary (x(s),y(s)).
The second form of the boundary condition, called the natural
boundary condition, reduces to the Neuman boundary condition
when p=q=1.

MGGHAT uses a finite element method with linear, quadratic or
cubic elements (user selectable) over triangles.  The adaptive
refinement via newest vertex bisection and the multigrid iteration
are both based on a hierarchical basis formulation.  Run time and
a posteriori graphical displays are made with gnuplot.

MGGHAT is available from netlib.  It can be obtained, for example,
by sending the message "send index from pdes/mgghat" (without the
quotes) to

MGGHAT can also be obtained from MGNET (multigrid network) by
annonymous ftp to in directory mgnet/mgghat.

A hypertext version of the User's Guide can be accessed by your
favorite WWW brower at the URL

The major changes from version 1.0 to version 1.1 are:

o Added the capability to save the program state and evaluate
  the solution from the saved state.  This is useful for systems
  of equations, time dependent problems, nonlinear equations, etc.

o Added example user files for a system of equations and for a
  time dependent problem.

o Custom run-time graphics replaced with gnuplot.

o Added an (optional) Tcl/Tk based menu for addition/deletion of
  graphical displays during run time, and rotation of 3D plots.

o Repackaged to be more user oriented.

o Now have makefiles for AIX, Convex (not vector), Cray (not vector),
  HPUX, IRIX, Linux, and SunOS.  Makefiles are more user oriented,
  and, when available, have an auto-double option.

o Default parameter values are set in block data rather than the
  user main program.

o New User's Guide written.

o Several minor bug fixes and robustness improvements.

William F. Mitchell                            |
Applied and Computational Mathematics Division |
National Institute of Standards and Technology |
Gaithersburg, MD 20899                         |


Date: 20 June 1994
From: Petr Mejzlik 
Subject: Paper by Petr Mejzl\'{\i}k for DDM7

A Bisection Method to Find All Solutions of a System of Nonlinear Equations

                             Petr Mejzl\'{\i}k
              Institute of Computer Science, Masaryk University
                  Buresova 20, 60200 Brno, Czech Republic

This paper describes an algorithm for the solution of a system of nonlinear 
equations F(x) = 0, where F = (f_1, ..., f_n): D \subset R^m --> R^n 
and D is a compact domain, given that any of the functions f_i is
monotonic when restricted to any single variable at an arbitrary point. The
algorithm finds an approximation of the solutions as a union of
m-dimensional intervals. The computation is based on reduction of the box
containing all the solutions, its bisection, and elimination of subintervals
which do not contain a solution. The algorithm does not require computation
of partial derivatives or their approximations. Its use is illustrated on a
model case.

Petr Mejzlik                           Internet:
Institute of Computer Science          FAX : +42-5-41212747
Masaryk University                     Tel.: +42-5-41321237 ext. 471
Buresova 20                                  +42-5-41213125
60200 Brno, Czech Republic

    Editor's Note: Look in mgnet/DDM7/mejzlik.abs and mgnet/DDM7/


Date: Fri, 1 Jul 1994 12:50:32 +0200
From: Ulrich Ruede 
Subject: Preprint by Jung and Ruede

I have downloaded a new preprint to mgnet:

Implicit Extrapolation Methods for Multilevel Finite Element Computations
        Theory and Applications

        Michael Jung  and  Ulrich Ruede

        Fakultaet fuer Mathematik,
        Technische Universitaet Chemnitz-Zwickau,
        D-09009 Chemnitz

        SPC-Preprint 94_11, Juni 1994

The submission consists of the abstract
in an asci-file and the full paper in a compressed postscript file

For Europeans it may be more convenient to access the paper
through the ftp server here at Chemnitz
in directory
This and in the other subdirectories of pub/Local/mathematik
contain Tech Reports of the various groups here, in particular the
Research Group on Scientific Parallel Computing. For those who use
the World Wide Web (mosaic), the most convenient access to my preprints may be

Uli Ruede

    Editor's Note: These files are in mgnet/papers/Jung-Ruede.


Date: Sat, 11 Jun 94 20:08:29 -0600
From: Steve McCormick 
Subject:  Copper Mountain 1995

                          CALL FOR PAPERS
                         MULTIGRID METHODS
                     Copper Mountain, Colorado
                      March 28 - April 2, 1995
                      ORGANIZING INSTITUTIONS  
                     The University of Colorado 
              Front Range Scientific Computations, Inc.
         The Society for Industrial and Applied Mathematics  
                        Department of Energy
            National Aeronautics and Space Administration
    Tom Manteuffel and Steve McCormick,
    University of Colorado 
     Joel Dendy, Los Alamos National Laboratory
     Van Henson, Naval Postgraduate School
     Jan Mandel, University of Colorado at Denver
     Duane Melson, NASA Langley
     Seymour Parter, University of Wisconsin
     Joseph Pasciak, Brookhaven National Laboratory
     John Ruge, University of Colorado at Denver
     Klaus Stueben, Gesellschaft f. Math. u. Datenverarbeitung
     Pieter Wesseling, Delft University
     Olof Widlund, Courant Institute 
     Paul Frederickson, RIACS
     Craig Douglas, IBM/CERFACS/Yale
     Circus: forum for late-breaking results 
     Workshops: informal discussions on topics of interest
     Student Paper Competition: travel support and special session
     Proceedings: NASA Conference Publication
     Special Journal Publication of Proceedings: to be arranged 
     Student Papers                     Dec.  1, 1994
     Abstracts                          Jan.  1, 1995
     Lodging Reservations               Feb. 28, 1995
     Early Registration                 Feb. 28, 1995
     Papers for Proceedings             Mar. 28, 1995 (at conference)
Anyone who wants to give a talk and will not be in the regular
session may do so at one of the Circus sessions.  Each session will
begin by asking who wishes to speak and setting the Circus schedule
accordingly.  To provide a more informal atmosphere for discussion,
there will also be time set aside for Workshops on specific topics
of interest.  The schedule for the Circus and Workshops will be
determined each day according to demand.
Travel and lodging assistance will be awarded to students judged to
have the best research papers.  Papers with original research due
mainly to the student must be received before December 1, 1994. They
must be singly authored and no more than 10 pages in length. Sending
only an abstract is unacceptable. A panel of judges will determine
the best papers. Authors of the three best papers will be given a
travel stipend, free lodging, and registration, and will be expected
to present their papers in a special session at the conference.
Several other students who participate in the competition may be
awarded free lodging and registration, depending on availability of
There will be morning and late afternoon lectures of approximately
25 minutes each.  Evenings will be reserved for the Circus and
Workshops.  Afternoons will be open for more informal workshops and
Potential contributors should submit an abstract of no more than
200 words by hard copy regular mail or, preferrably, by postscript
or ASCII email to:
          University of Colorado
          Program in Applied Math
          CB 526
          Boulder CO  80309-0526

If you submit your abstract by email, please do NOT send a TeX or 
LaTeX file. It will be returned if you do. Instead, please send
either a ps file or an ASCII file.

Authors will be notified of the status of their abstract by the end
of January. Those notified of acceptance may then submit their paper
to the proceedings by bringing hard copy to the conference. Authors
who wish to make their paper generally available will be encouraged
to submit a ps file to mgnet.
                   LOGISTICS, REGISTRATION, ETC.
                      MAY BE DIRECTED TO:                         


                     YOU MAY ALSO CONTACT:

                    Steve McCormick
                    University of Colorado
                    Program in Applied Math
                    CB 526
                    Boulder CO  80309-0526
                     FEBRUARY BY ACCESSING:

Tuesday, March 28
     7:00-9:00 p.m.      Registration and Reception
Wednesday, March 29
     8:00-12:00          Session 1
     4:30-6:30 p.m.      Session 2
     7:30-9:30 p.m.      Circus/Workshops
Thursday, March 30
     8:00-12:00          Session 3
     4:30-6:30 p.m.      Session 4
     7:30-9:30 p.m.      Circus/Workshops
Friday, March 31
     8:00-12:00          Session 5
     4:30-6:30 p.m.      Session 6
     7:30-9:30 p.m.      Banquet
Saturday, April 1
     8:00-12:00          Session 7
     4:30-6:30 p.m.      Session 8
     7:30-9:30 p.m.      Circus/Workshops
Sunday, April 2
     8:00-12:00          Session 9
                    Conference Registration Form
          Copper Mountain Conference on Multigrid Methods
                       March 28 - April 2, 1995
Please fill out a COMPLETE MAILING ADDRESS so that we can
be sure you receive further information!!
Street Address____________________________________________________
Office Phone_____________________Home Phone_______________________

Email Address_____________________________________________________
Advanced Registration Fee (Postmarked by February 28, 1995)
    ( )   SIAM Member $215
    ( )   Non-member  $245
    ( )   Student     $ 50
Registration Fee (Postmarked after February 28, 1995)
    ( )   SIAM Member $255
    ( )   Non-member  $285
    ( )   Student     $ 50
Please Indicate method of payment:
    ( )   Check or Money Order enclosed (Checks made payable to:  
    ( )   Master Card or Visa: (NO American Express)
Authorized Signature______________________________
Send registration form by regular mail or email to: 
     University of Colorado
     Program in Applied Math
     CB 526
     Boulder CO  80309-0526

Date: Thu, 16 Jun 1994 11:00:37 +0200
From: roitzsch@sc.ZIB-Berlin.DE (Rainer Roitzsch)
Subject: Bib entries from Konrad-Zuse Zentrum

Dear Craig Douglas,

I attach a selection of references to the work on adaptive fem at
the Konrad-Zuse Center (ZIB). We hope that the references can be included
in the MG-bibliography database.

Rainer Roitzsch

    Editor's Note: These will be added to the database in the next few days
    -------------  in the standard MGNet style.

AUTHOR = {Ackermann, Joerg and Roitzsch, Rainer},
TITLE = {A two-dimensional multilevel adaptive finite element method for
the time-independent {S}chroedinger equation},
JOURNAL = {Chem. Phys. Letters},
YEAR = {1993},
VOLUME = {214},
PAGES = {109-117} }

AUTHOR = {Bornemann, Folkmar},
TITLE = {Adaptive multilevel solutions of parabolic equations in {3D}},
YEAR = {1992},
PAGES = {75-81},
EDITOR = {Vichnevetsky, Robert and Knight, Doyle and Richter, Gerard},
BOOKTITLE = {Advances in Computer Methods for Partial Differential Equations --

AUTHOR = {Bornemann, Folkmar},
TITLE = {An adaptive multilevel approach to parabolic equations {III}. {2D}
error estimation and multilevel preconditioning},
JOURNAL = {IMPACT Comput. Sci. Engrg.},
YEAR = {1992},
VOLUME = {4},
PAGES = {1-45} }

AUTHOR = {Bornemann, Folkmar},
TITLE = {An adaptive multilevel approach to parabolic equations {I}.
General theory and {1D}-implementation},
JOURNAL = {IMPACT Comput. Sci. Engrg.},
YEAR = {1990},
VOLUME = {2},
PAGES = {279-317} }

AUTHOR = {Bornemann, Folkmar},
TITLE = {An adaptive multilevel approach to parabolic equations II.
variable-order time discretization based on a multiplicative error
JOURNAL = {IMPACT Comput. Sci. Engrg.},
YEAR = {1991},
VOLUME = {3},
PAGES = {93-122} }

AUTHOR = {Bornemann, Folkmar and Erdmann, Bodo and Kornhuber, Ralf},
TITLE = {Adaptive multilevel methods in three space dimensions},
JOURNAL = {J. Numer. Meth. Engrg.},
YEAR = {1993},
VOLUME = {36},
PAGES = {3187-3203} }

AUTHOR = {Bornemann, Folkmar and Yserentant, Harry},
TITLE = {A basic norm equivalence in the theory of multilevel methods},
JOURNAL = {Numer. Math.},
YEAR = {1993},
VOLUME = {64},
PAGES = {455-476} }

AUTHOR = {Deuflhard, Peter},
TITLE = {Cascadic conjugate gradient methods for elliptic partial
differential equations. {A}lgorithm and results},
BOOKTITLE = {Proc. 7th Domain Decompostion Meeting},
YEAR = {1993},
PUBLISHER = {Penn State Univ.} }

AUTHOR = {Deuflhard, Peter and Freund, R. and Walter, Artur},
TITLE = {Fast secant methods for the iterative solution of large
nonsymmetric linear systems},
JOURNAL = {IMPACT Comput. Sci. Engrg.},
YEAR = {1990},
VOLUME = {2},
PAGES = {244-276} }

AUTHOR = {Deuflhard, Peter and Leinen, Peter and Yserentant, Harry},
TITLE = {Concepts of an adaptive hierarchical finite element code},
JOURNAL = {IMPACT Comput. Sci. Engrg.},
YEAR = {1989},
VOLUME = {1},
PAGES = {3-35} }

AUTHOR = {Deuflhard, Peter and Potra, Florian A.},
TITLE = {Asymptotic mesh independence of {N}ewton-{G}alerkin methods via a
refined {M}ysovskii theorem},
JOURNAL = {SIAM J. Numer. Anal.},
YEAR = {1992},
VOLUME = {29},
NUMBER = {5},
PAGES = {1395-1412} }

AUTHOR = {Erdmann, Bodo and Frei, M. and Hoppe, R.H.W. and Kornhuber, Ralf
and West, U.},
TITLE = {Adaptive finite element methods for variational inequalities},
JOURNAL = {East-West J. Numer. Math.},
YEAR = {1993},
VOLUME = {1},
NUMBER = {3},
PAGES = {165-197} }

AUTHOR = {Hoppe, R.H.W. and Kornhuber, Ralf},
TITLE = {Adaptive multilevel--methods for obstacle problems},
JOURNAL = {SIAM J. Numer. Anal.},
YEAR = {1994},
VOLUME = {31},
NUMBER = {2},
PAGES = {301-323} }

AUTHOR = {Hoppe, R.H.W. and Kornhuber, Ralf},
TITLE = {Multigrid solution of two coupled {S}tefan equations arising in
induction heating of large steel slabs}, JOURNAL = {J. Numer. Meth. Engrg.},
YEAR = {1990},
VOLUME = {30},
PAGES = {779-801} }

AUTHOR = {Kornhuber, Ralf and Roitzsch, Rainer},
TITLE = {On adaptive grid refinement in the presence of internal or
boundary layers},
JOURNAL = {IMPACT Comput. Sci. Engrg.},
YEAR = {1990},
VOLUME = {2},
PAGES = {40-72} }

AUTHOR = {Kornhuber, Ralf and Roitzsch, Rainer},
TITLE = {Self adaptive finite element simulation of bipolar, strongly
reverse biased pn-junctions},
JOURNAL = {Comm. Num. Meth. in Engrg.},
YEAR = {1993},
VOLUME = {9},
PAGES = {243-250} }

AUTHOR = {Kornhuber, Ralf and Roitzsch, Rainer},
TITLE = {Self adaptive computation of the breakdown voltage of planar
pn--junctions with multistep field plates},
YEAR = {1991},
PAGES = {535-543},
EDITOR = {Fichtner et al, W.},
BOOKTITLE = {Proceedings of the 4th Conference of Somulation of Semiconductor
Devices and Processes},
PUBLISHER = {Hartung--Gorre} }

AUTHOR = {Kornhuber, Ralf and Wittum, Gabriel},
TITLE = {Discretization and preconditioning of convection diffusion equations},
YEAR = {1993},
PAGES = {67-77},
EDITOR = {Hackbusch et al, W.},
BOOKTITLE = {Incomplete Decomposition (ILU) -- Algorithms, Theory, and
PUBLISHER = {Vieweg} }

AUTHOR = {Lang, Jens and Walter, Artur},
TITLE = {A finite element method adaptive in space and time for nonlinear
reaction-diffusion systems},
JOURNAL = {IMPACT Comput. Sci. Engrg.},
YEAR = {1992},
VOLUME = {4},
PAGES = {269-314} }

AUTHOR = {Lang, Jens and Walter, Artur},
TITLE = {An adaptive {R}othe method for nonlinear reaction-diffusion systems},
JOURNAL = {Applied Numer. Math.},
YEAR = {1993},
VOLUME = {13},
PAGES = {135-146} }


End of MGNet Digest