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

Today's editor:  Craig Douglas (

Volume 3, Number 3 (March 31, 1993)

Today's topics:

     Copper Mountain Papers
     Approximate Schemes for Dense Matrix Algebra
     Announcement for the Amsterdam Multigrid Conference
     List of Contributions for the Amsterdam Multigrid Conference
     MGGHAT now available
     New paper on embedding Multilayer methods onto MIMD machines


Date: Wed, 31 Mar 93 11:30:00 EST
From: "Craig C. Douglas" 
Subject: Copper Mountain Papers

This is just a reminder that the papers are due during the conference and will
be lightly refereed beginning Wednesday, April 7, 1993.

Please remember to download your files to mgnet before, during, or after the
conference.  The NASA macros are in the directory mgnet/CopperMtn/Macros.

There are now a number of papers (in mgnet/CopperMtn) and abstracts
(in mgnet/CopperMtn/Abstracts) available via anonymous ftp from


Date: Thu, 4 Mar 93 10:32:29 -0500
From: "John D. McCalpin" 
Subject: Approximate Schemes for Dense Matrix Algebra

We are looking for one or more mathematicians interested in
collaborating on a project involving the application of new approximate
schemes for dense matrix algebra (multi-grid or wavelet transform
based) to a problem in optimal data assimilation using various versions
of the Kalman filter.   

The matrices are very large, O(10**4) to O(10**5), and dense, but
hierarchically band dominant.   For the mathematician, the project
would involve the derivation of error bounds for the approximate
scheme(s) for a class of matrices that differ slightly from those
studied in the recent literature, and assistance in developing
efficient parallel algorithms.   From the dynamical side, we wish to
investigate the properties of the Kalman filter (i.e. its loss of
optimality) for various types and degrees of degradation of the
information content of the error covariance array.

We are interested in pursuing funding through both the NSF
math/geosciences collaboration initiation program (NSF 92-127) and 
through the new round of the NSF HPCC program (as part of a larger
group of investigators).

John D. McCalpin, University of Delaware 
Ichiro Fukumori, JPL/NASA
John D. McCalpin              
Assistant Professor           
College of Marine Studies, U. Del.


Date: Fri, 12 Mar 93 9:57:44 MET
From: Piet Wesseling (NW)  and
Subject: Announcement for the Amsterdam Multigrid Conference

                             July 6-9, 1993
                            The Netherlands

                         Call for Registration

Following earlier conferences in 1981, 1985 and 1990 EMG'93, the fourth 
European Multigrid Conference will be held in Amsterdam, July 6-9, 1993. 
The conference aims at providing a forum for the presentation and 
discussion of recent research in the general area of multilevel methods. 
Participation will be worldwide.
Themes of the conference will range from basic research to industrial 
applications, including but not limited to the following subjects:
*       computational fluid dynamics, 
*       reservoir engineering, 
*       statistical physics, 
*       parallel computing, 
*       adaptive computing, 
*       numerical analysis of multigrid methods.

The program of the fourth European Multigrid Conference will consist of 
invited lectures in plenary sessions and some 50 selected presentations 
in 2 parallel sessions. Program chairmen are P.W. Hemker (CWI/University 
of Amsterdam) and P. Wesseling (Delft University of Technology).

Invited lectures will be given by:
*       D. Braess (Bochum, Germany), Multigrid for nonconforming finite 
*       A. Brandt (Rehovot, Israel), Multiscale computations in physics 
        and chemistry
*       C. Douglas (New Haven, USA), A generalized multigrid theory
*       L. Fuchs (Stockholm, Sweden), Multigrid methods for calculation 
        of flows in multi-components geometries
*       W. Hackbush (Kiel, Germany), The frequency-decomposition 
        multigrid method
*       P. Lauwers (Bonn, Germany), Multigrid algorithms in lattice 
        gauge theories
*       K. Stueben (Sankt Augustin, Germany), Adaptive multigrid on 
        distributed memory computers
*       S. Vandewalle (Leuven, Belgium), Multicomputer solution of 
        parabolic problems
*       C.H. Venner (Enschede, The Netherlands), Multilevel solution of 
        integral and integro-differential problems
*       G. Wittum (Heidelberg, Germany), On robust multigrid methods
The program of contributed papers will be made available through 
anonymous ftp from "" (Conferences/EMG93) and mgnet 
( If you don't have ftp-facilities you may contact the 
conference administrator (Beau) at fax +31-20-5924199 to get the program.

Invited lectures and a number of selected presentations will be published 
as full papers in conference proceedings, to be published by Birkhauser. 
Other presentations can be published as full papers by CWI. A booklet 
with abstracts of all presentations will be available at the conference.

Famous for its canals, gabled houses, bridges and museums, Amsterdam 
attracts millions of tourists each year. The old inner city is the most 
compact in Europe with almost 7000 National Trust buildings dating from 
the 16th century, when it was Europe's most important commercial center, 
to the 20th century. Through the centuries, Amsterdam has been 
characterized by its international orientation and its tolerance of diverse 
political and religious opinions. Amsterdam is the largest city in the 
Netherlands and the country's nominal capital. It is a major international 
center of commerce, banking, press and publishing, and ranks high on 
the list of cities favored with international conventions. At the beginning of 
July weather in Amsterdam is mild, with daily highs usually in the 18-24 C 
range (65-75 F) and occasional showers.

Amsterdam is a major international airport served by most international 
carriers and linked directly to over a hundred cities worldwide. The 
Schiphol International Airport is about 20 kilometers from Amsterdam's city 
center. There are frequent trains to Amsterdam's central railway station 
(Centraal Station, running 24 hours per day) as well as to Amsterdam-
Zuid/WTC station, in the vicinity of EMG's hotel accommodation.

EMG'93 will be held in the "Oudemanhuispoort", which is part of the 
University of Amsterdam and beautifully located in the old inner city of 
Amsterdam (street address: Oudezijds Achterburgwal 229). Lunches will 
be served in the university's restaurant within easy walking distance of the 
meeting rooms.
The "Oudemanhuispoort" can be easily reached with public transport from 
all parts of the city.

Register and make your hotel reservation in advance using the Registr-
ation and Accommodation Form. Mail or email the form to the Conference 
Administrator. No registration or reservation will be effective until full 
payment is received.

The conference fee for EMG'93 is NLG 400 for advance registration (to be 
paid by 1 May 1993) and NLG 500 after 1 May 1993. The fee includes ad-
mission to all sessions, lunches and coffee breaks, the conference dinner 
and a copy of the book af abstracts.

A limited number of inexpensive hotelrooms have been reserved at the 
Guesthouse of the Free University. Accommodation is available at NLG 
40.-- (rooms) and NLG 82.50 (apartments) per night, excl. breakfast. All 
accommodation is for single use.
To book accommodation, fill out the hotel reservations part of the Registr-
ation and Accommodation Form. The full hotel costs will have to be paid in 
advance. After having received your reservation you will be invoiced for 
these costs.

All payments, net of all charges, are to be made in Dutch Guilders (NLG) 
and can be made by either
*       Banker's draft (bank check), made payable to "Stichting Wiskunde 
        en Informatica Conferenties", and mailed to the conference 
*       Money transfer to account # of "Stichting 
        Wiskunde en Informatica Conferenties" at the RABO-bank, 
        Middenweg 88, Amsterdam (Postal giro of bank: 187744)
Make sure that payments mention "EMG93" and your name.

P.W. Hemker
F. Snijders
P. Wesseling
S. van der Wolff

EMG'93 - European Multigrid Conference 1993
c/o CWI
Ms. Simone van der Wolff
P.O. Box 4079
1009 AB Amsterdam
The Netherlands
Tel. +31-20-5924189
Fax. +31-20-5924199

IBM Nederland NV
KNAW - Royal Academy of Sciences
University of Amsterdam

Please type or print

Name (last):
Postal Code:
Special Requests:

registers for EMG'93.


Please reserve a hotelroom in price category:
[]      NLG 40 
[]      NLG 82.50
Date of arrival:
Date of departure:
Number of nights:

[] I enclose a banker's draft
[] I transferred the registration fee to your bankaccount

Email to:
                EMG'93 Secretariat
                CWI/Ms. Simone van der Wolff
                P.O. Box  4079
                NL 1009 AB Amsterdam


Date: Fri, 12 Mar 93 10:08:52 MET
From: Piet Wesseling (NW)   and
Subject: List of Contributions for the Amsterdam Multigrid Conference

% This LaTeX file containes the list of
% authors and titles of contributions to the
%     4th European Multigrid Conference
% to be held 6th - 9th July 1993
% at the University of Amsterdam

% Sponsored by:
%     CWI
%     AKZO N.V.
%     IBM N.V.
%     KNAW
%     University of Amsterdam



{\bf Preliminary version, March 8th, 1993}


{\bf Invited papers}
\item D. Braess (Bochum, Germany)
      {\em Multigrid methods for nonconforming finite elements}   
\item A. Brandt (Rehovot, Israel)
      {\em Multiscale Computations in Physics and Chemistry} 
\item Craig Douglas (New Haven, USA)
      {\em A generalized multigrid theory}
\item L. Fuchs (Stockholm, Sweden) {\em Multi-grid methods for
      calculation of flows in multi-components geometries}
\item W. Hackbusch (Kiel, Germany)
      {\em The frequency decomposition multigrid method}
\item P. Lauwers (Bonn, Germany)
      {\em Multigrid Algorithms in lattice gauge theories}
\item K. Stueben (Sankt Augustin, Germany)
      {\em Adaptive multigrid on distributed memory computers} 
\item S. Vandewalle (Leuven, Belgium)
      {\em Multicomputer Solution of Parabolic Problems}
\item C. Venner (Enschede, The Netherlands)
      {\em Multilevel solution of Integral and Integro-differential problems} 
\item G. Wittum (Heidelberg, Germany)
      {\em On Robust Multi-grid Methods}          
{\bf Contributed papers}
\item Y. Luh {\em  GMD Bonn, BRD } \\ 
A study on narrow stencil discretisations and multigrid solvers for the 
generalized Stokes equations 
\item E. Katzer {\em  Univ.Kiel, BRD } \\ 
A subspace decomposition twogrid method for hyperbolic equations 
\item P.Ekstrand and L.Fuchs {\em  RIT Stockholm, Sweden } \\ 
A pressure correction method for compressible viscous flow for both low and high Mach numbers 
\item X.S.Bai and L.Fuchs {\em  RIT Stockholm, Sweden } \\ 
A multigrid method applied to calculation of turbulent reacting flows in 
\item E.Dick and K.Riemslagh {\em  Univ.Gent, Belgie } \\ 
Multigrid methods for steady Euler equations based on multistage Jacobi-
\item K.Riemslagh and E.Dick {\em  Univ.Gent, Belgie } \\ 
Implementation aspects of the multigrid formulation of flux-difference
splitting algorithms for steady Euler equations on unstructured adaptive grids 
\item H.Kuerten and B.Geurts {\em  UT, Enschede } \\ 
A multigrid multiblock solver for compressible turbulent flow 
\item N.G.Wright and P.A.Sleigh {\em  Univ. London + Leeds } \\ 
Different approaches to discretisation within a multilevel algorithm 
\item K. de Cock {\em  NLR, Amsterdam } \\ 
Multigrid acceleration of the 2D Euler equations applied to high lift systems 
\item R. Stevenson {\em  RUU, Utrecht } \\ 
Robustness of multigrid applied to anisotropic boundary value problems on 
general domains 
\item B.Gustafsson and P.Lotstedt {\em  Uppsala + SAAB, Sweden } \\ 
Solution of the stationary Euler equations by an improved GMRES algorithm and 
multigrid acceleration 
\item Kloppmann, Schwamborn and Singh {\em  DLR Gottingen, BRD } \\ 
Multigrid solution of the 2D Navier-Stokes equations for transonic internal and 
external flows 
\item vdPloeg, Botta and Wubs {\em  RUU, Groningen } \\  
Grid-independent convergence based on preconditioning techniques 
\item Liu, Liu and McCormick {\em  Univ.Colorado } \\ 
Multigrid methods for flow transition in a 3D channel 
\item M. Griebel {\em  TU Muenchen, BRD } \\ 
Parallel point oriented multilevel methods 
\item Costiner and Ta'asan {\em  Weizmann Inst., Israel } \\  
Separation techniques for multilevel eigenvalue solvers applied to
a non-linear Schroedinger eigenvalue problem 
\item P.Johansson and L.Davidson {\em  TU Gothenburgh, Sweden } \\ 
A full multigrid method applied to turbulent flow using
the SIMPLEC algorithm together with a collocated grid arrangement 
\item He Jun and Kang Li-shan {\em  Univ.Wuhan, China } \\ 
Parallel multigrid method based on wavelet multiresolution 
\item W.Heinrichs {\em  Univ. Dusseldorf, BRD } \\ 
Spectral multigrid methods for the reformulated Stokes equations 
\item Brandt, Mikulinsky and Chaimson {\em  Weizmann Instititute  } \\ 
Multigrid parallelization of parabolic equations 
\item R.H.W. Hoppe {\em  TU Muenchen } \\ 
Adaptive multilevel finite element methods for variational inequalities 
\item T. Neunhoeffer {\em  TU Muenchen } \\ 
Multigrid methods for mixed finite element discretizations of variational 
\item A. Reusken {\em  TU Eindhoven } \\  
Multigrid with matrix dependent transfer operators for convection diffusion 
\item S. Ta'asan {\em  Weizmann + NASA Langley } \\ 
Optimal multigrid solvers for steady state inviscid flow problems 
\item Z. Haras and  S. Ta'asan {\em  Weizmann + NASA Langley } \\ 
LDS Method for linear evolution equations 
\item A.Borzi and  A. Koubek {\em  SISSA, Trieste, IT } \\ 
A multigrid method for the resolution of thermodynamic Bethe Ansatz equation 
\item A.Borzi  {\em  SISSA, Trieste, IT } \\ 
Burgers equation and multigrid techniques 
\item Keyzer and  D.Roose {\em  Univ. Leuven } \\ 
Parallel steady Euler equations calculations using multigrid methods
and adaptive irregular meshes 
\item B.Bialecki and  M.Dryja {\em  Univ.Kentucky + Warsaw } \\ 
Multigrid and multilevel methods for orthogonal spline collocation problems 
\item A.I. Zhmakin {\em  Acad.Sci, Petersburg }  \\
A multigrid method for low Mach number compressible Navier-Stokes
equations on 2D unstructured grids 
\item G. Golubovici and  C. Popa {\em  Univ. Constanta, Roum } \\ 
Interpolation and related coarsening techniques
for the algebraic multigrid method 
\item J.E. Dendy {\em  Los Alamos, USA } \\ 
Multigrid methods for SIMD machines 
\item C. Srinivasa {\em  NAL, Bangalore, India } \\ 
Solution of full potential equation on an airfoil by multigrid technique 
\item Catalano, de Palma, Napolitano and Pascazio
 {\em  Univ. Bari, IT } \\ 
A multidimensional solution adaptive multigrid solver for inviscid cascades 
\item J. Fuhrmann {\em  Acad.Sci. Berlin } \\  
On matrix data structures and the stability of multigrid algorithms 
\item Koobus, Carre Lallemand and Dervieux
 {\em  INRIA, Sophia, FR } \\ 
Agglomeration unstructured multigrid for diffusion problems 
\item S. Turek {\em  Univ. Heidelberg, DE } \\ 
Multigrid techniques for a class of discretely divergence-free finite element 
\item M.G. Edwards and C. Rogers {\em  BP Research, GB } \\ 
Renormalisation and Multigrid for reservoir simulation 
\item Tom Fogwell {\em  IT Corp Martinez } \\ 
A multigrid method for multiphase fluid flow in porous media 
\item Zhu, Lacor and Hirsch {\em  Vrije Univ. Brussel } \\ 
A new residual smoothing method for multigrid acceleration applied to the  
Navier-Stokes equations 
\item U. Ruede {\em  TU Muenchen } \\ 
Multilevel, extrapolation and sparse grid methods 
\item Oosterlee and Wesseling {\em  TU Delft } \\ 
A multigrid waveform relaxation method for time dependent 
incompressible Navier-Stokes equations 
\item Yunqing Huang and Yaoxing Li {\em  Xiangtan Univ., China } \\
The convergence of multigrid methods for solving finite element equations
in the presence of singularities 
\item P.M. de Zeeuw and Hemker (?) {\em  CWI, Amsterdam } \\ 
Asymmetric prolongations and restrictions by equation reconstruction 
\item H.T.M. van der Maarel {\em  CWI, Amsterdam } \\  
A-posteriori estimation of the local discretisation error
for conservation laws and adaptive local grid refinement
\item O.Dorok, F.Schieweck and L.Tobiska {\em  UT Magdeburg } \\  
A Multigrid method for Solving the Bousinesq approximation of the
Navier-Stokes Equations 
\item Zeng and Wesseling {\em  TU Delft } \\ 
Multigrid for the incompressible Navier-Stokes equations 
on staggered grids in general coordinates 
Jan Janssen {\em KU Leuven} \\
Multigrid waveform relaxation with spatial finite element discretisation
L.R.~Matheson and R.E.~Tarjan {\em Princeton Univ., USA} \\
Critical analysis of concurrent and standard
multigrid methods on massively parallel computers
Z.P.~Novak and M.~Salcudean {\em Univ.~B.C., Canada} \\
Fluid flow calculations by a nonlinear multigrid method
L.~Giraud and R.S.~Tuminaro {\em CERFACS, Toulouse} \\
A domain decomposition algorithm for drift-diffusion equations
M.Lemke, K.~Witsch and D.~Quinlan
{\em Univ.~D\"usselforf and Univ.~Colorado} \\
Parallel multilevel methods with adaptive mesh refinement -
Results for singular perturbation problems
\item Xu Changfa {\em Huazhong Univ., Wuhan, China} \\
B-wavelet MG-methods



Date: Wed, 24 Mar 93 10:47:46 EST
From: (Bill Mitchell)
Subject: MGGHAT now available

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

   -( p * u  )  -( q * u  )  +  r * u  = f    in omega
           x  x         y  y

                                    u  = g  on boundary 1 of omega
                           u  + c * u  = g  on boundary 2 of omega

where p,q,r,f,c and g are functions of x and y, omega is a polygonal
domain in R2, n is the outward unit normal, and boundaries 1 and 2
form a partition of the boundary of omega.

This program is now available for public use.  It uses a finite element method
with linear, quadratic or cubic elements over triangles.  The adaptive
refinement via newest vertex bisection and the multigrid iteration are both
based on a hierarchical basis formulation.  Visualization is available at run
time through an X Windows display, and a posteriori through output files that
can be used as GNUPLOT input.

The user supplies subroutines containing the coefficient functions p, q, r,
f, c, and g, a subroutine with the inital (coarse grid) triangulation to
define the domain, and a main program that sets program parameters and calls
subroutine mgghat.  The program is mostly FORTRAN 77, but uses C for the
X11 graphics (can be disabled).

MGGHAT is available from netlib.  Send the message 
"send index from pdes/mgghat" to

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

-- Bill

William F. Mitchell                 |
GE Advanced Technology Labs         |
Bldg. 145-2 Moorestown Corp. Center |  Voice: (609) 866-6536
Moorestown, NJ  08057               |  Fax:   (609) 866-6397


Date: Wed, 31 Mar 93 16:53:11 BST
Subject: New paper on embedding Multilayer methods onto MIMD machines

  I have downloaded a new paper 

         Embedding a ``Tree-Code" on a MIMD parallel computer 
              using a domain decomposition paradigm. 


  "Tree-codes" and ``Multigrid" methods are both hierarchical, recursive 
algorithms for solving vector field equations, faster than traditional methods.
If the methods are non-adaptive, then both Tree-codes and Multigrid methods
require identical parallel communication strategies since they share the same 
internal structure.

  This paper contains a description of Tree-codes followed by an overview of one
Tree-code in particular, namely the Greengard-Rokhlin Fast Multipole Method 
in 2 dimensions.  The parallel version of this method is then discussed at 
length with the presentation of two communication strategies; An "All-to-all" 
broadcast and a routine which enables neighbouring transputers to swap 
information in a predetermined sequence. 

  The computer used is the Meiko Computing Surface, a MIMD parallel computer 
employing T800 transputers running the CSTools communication harness using a
SPMD (Single Program, Multiple Data), domain decomposition paradigm. 

  The 3D versions of the communication algorithms have also been implemented
and are  simple extensions of the 2D routines.

Keywords: Fast Multipole Method, Multigrid, message passing.`

  Gavin J. Pringle
  Napier University
  EH14 1DJ
  Tel : +44 31 455 4631
  Email :

    Editor's Note: mgnet/papers/Pringle/Embed.dvi.Z,
    -------------  mgnet/papers/Pringle/,  and


End of MGNet Digest