NA Digest Sunday, May 7, 1995 Volume 95 : Issue 19

Today's Editor:
Cleve Moler
The MathWorks, Inc.

Submissions for NA Digest:

Mail to

Information about NA-NET:

Mail to

URL for the World Wide Web: -------------------------------------------------------

From: Gene Golub <golub@sccm.Stanford.EDU>
Date: Fri, 5 May 95 22:58:05 PDT
Subject: Happy News

Congratulations to

Nancy Nichols and Mike Green


Daniela Calvetti and Lothar Reichel

on their recent marriages.

Gene Golub


From: Des Higham <>
Date: Tue, 2 May 95 16:58:34 BST
Subject: Seeking 1956 De Vogelaere Report on Hamiltonian Systems

Prof. Ron Mitchell recently told me about a reference concerning
numerical methods for Hamiltonian problems:

Methods of integration which preserve the contact transformation property of
the Hamiltonian equations,
by Rene De Vogelaere,
Report 4, Department of Mathematics, University of Notre Dame, 1956.

The Physical Sciences Librarian at the University of Notre Dame was able
to confirm that the report was written between January and September 1956
and is of 14 pages in length. However, no copies are available.
The report would seem to be a very early reference in an area that is currently
extremely active--does anyone have a copy of the report, or information about
whether it was ever published?

Des Higham
University of Dundee


From: Simon Chamlian <>
Date: Tue, 2 May 95 13:24:02 EDT
Subject: C Program for Gauss Kronrod Quadrature


Does anyone know where I can find a Gauss Kronrod adaptive quadrature
subroutine (similar to the one given in: "Numerical methods and
software", by Kahaner, Moler and Nash) written in C or C++
(not converted to by f2c).

Thanks in advance.



From: Thomas Turnbull <>
Date: Thu, 4 May 1995 13:48:02 -0400
Subject: Parallel BiMMeR Matrix Multiplication Routines

We would like to announce a new release of the BiMMeR parallel matrix
multiplication routines. These routines were developed by
Huss-Lederman, Jacobson, Tsao, Turnbull, and Zhang at SRC. They use a
virtual 2D torus wrap data layout and the Broadcast-Multiple-Roll
algorithm. The code has been recently converted to MPI, the message
passing interface standard, and can be executed on a wide range of
parallel machines. The following operations are now supported:

C = alpha*A*B+beta*C,
C = alpha*A'*B+beta*C, (A' is A transpose)
C = alpha*A*B'+beta*C,

for any number of nodes with an implied 2D topology.

For a detailed description of the algorithms used, see Matrix
Multiplication on the Intel Touchstone Delta by Huss-Lederman,
Jacobson, Tsao and Zhang, Concurrency: Practice and Experience, Vol. 6
(7), Oct. 1994, pp. 571-594. The code and accompanying reports are
available via anonymous ftp from in pub/prism
or via a WWW browser from (The BiMMeR
code is in file bimmer.tar.Z, and the referenced article is in file


From: Rossana Vermiglio <>
Date: Fri, 5 May 1995 11:16:28 MET-DST
Subject: Change of Address of Rossana Vermiglio

My Department moved in a new building and my new address
is the following

Rossana Vermiglio
Dipartimento di matematica e Informatica
Universita' di Udine

via delle Scienze, 208 (loc. RIZZI)


phone number +39 (432) 558423
fax number +39 (432) 558499


From: William J. Stewart <>
Date: Thu, 4 May 1995 06:55:53 -0500
Subject: New Book on Markov Chains

William J. Stewart
Princeton University Press


It is often possible to represent the behavior of a physical system
by describing all the different states which it can occupy and
by indicating how it moves from one state to another in time.
If the future evolution of the system depends only on its current state,
the system may be represented by a Markov process. When the state space
is discrete, the term ``Markov Chain'' is employed.

In this book, the first to offer a systematic and detailed treatment
of the numerical solution of Markov chains, William Stewart provides
scientists on many levels with the power to put this theory to
use in the actual world, where it has application in areas as
diverse as engineering, economics and education. His efforts make
essential reading in a rapidly growing field.

Here Stewart explores all aspects of numerically computing
solutions of Markov chains, especially when the state space is
huge. He provides extensive background to both discrete-time
and continuous-time Markov chains and examines many different
numerical computing methods --- direct, single and multi-vector
iterative, and projection methods. Additionally, he considers
recursive methods often used when the structure of the Markov
chain is upper Hessenberg; iterative aggregation/disaggregation
methods that are particularly appropriate when it is NCD (nearly
completely decomposable) and reduced schemes for cases in which
the chain is periodic. There are chapters on methods for computing
transient solutions, on stochastic automata networks, and finally
on currently available software. Throughout, Stewart draws on
numerous examples and comparisons among the methods he so
thoroughly explains.


ISBN: 0-691-03699-3

Orders by telephone 1-800-777-4726, MUST supply the key code P5081.

Orders by fax must be sent to Princeton University to the
attention of PAUL MOONEY. The fax number is 609-258-1335

Orders by mail must be sent to PAUL MOONEY c/o Princeton
University Press, 41 William Street, Princeton, NJ 08540.


From: John R. Rice <>
Date: Fri, 05 May 1995 11:52:57 -0500
Subject: IFIP Conference on Quality of Numerical Software

Working Group 2.5 on Numerical Software

8-12 July 1996, Oxford, UK

First Announcement and Call for Interest

A working conference on The Quality of Numerical Software:
Assessment and Enhancement, will be arranged in the summer of
1996 in Oxford, England. This will be the seventh working
conference organized by the IFIP Working Group on Numerical
Software (WG 2.5) on behalf of the IFIP Technical Committee on
Programming (TC 2).

This message is an invitation to you from the conference
program committee to contribute your ideas for this conference.
At present we are very interested in suggestions regarding the
scope of the conference, who to invite to attend and who to
invite as speakers, and on which specific topics, and whether you
are interested in participating. The initial scope for the
conference is given below.

Previous working conferences arranged by WG 2.5 have focused
on performance evaluation of numerical software, numerical compu-
tation and programming languages, software for partial differen-
tial equations, problem solving environments, aspects of computa-
tion on asynchronous parallel processors, and scientific problem
solving environments. The number of participants is deliberately
kept small in order to preserve informality and allow substantial
time for interaction.

The conference program committee consists of Bo Einarsson
(Linkoping, Sweden), Brian Ford (Oxford, United Kingdom), Eric
Grosse (Murray Hill, New Jersey), Elias Houstis (Patras, Greece),
John Rice (West Lafayette, Indiana), Mladen Vouk (Raleigh, North
Carolina), and Margaret Wright (Murray Hill, New Jersey).

For further information you are invited to contact the Pro-
gram Committee co-chairs:

Brian Ford, NAG Ltd, Wilkinson House, Jordan Hill Road,
Oxford 0X2 8DR ENGLAND, UK. Telephone +44 865 511245
(Office) or +44 99 868 343 (Home). Electronic mail:

John Rice, Department of Computer Science, Purdue Univer-
sity, 1398 CS Bldg, West Lafayette, IN 47907-1398. Tele-
phone 317-494-6003 (Office) or 317-463-6507. Electronic

We hope that you will join us in organizing a stimulating
and enjoyable Working Conference on this challenging and timely


Brian Ford John R. Rice


From: Andrew Ilin <>
Date: Mon, 01 May 1995 13:47:21 -0500 (CDT)
Subject: Conference on Spectral and High Order Methods

Monday -- Friday, June 5-9, 1995
The Wyndham Warwick Hotel, Houston, Texas

Featured Speakers:

Guy Battle, Texas A&M U Christine Bernardi, U of Paris
Monique Dauge, U of Rennes Benqi Guo, U of Manitoba
Wilhelm Heinrichs, Heinrich-Heine U George Karniadakis, Brown U
Yvon Maday, U of Paris Ralph Metcalfe, U of Houston
J. Tinsley Oden, U of Texas

LATE-BREAKING COMMUNICATIONS: Contributed abstracts for short oral and poster
communications can be submitted until May 15. Abstracts should be 2 to 3 pages
in length. Authors of accepted abstracts will be notified no later than May 15.
Interested parties should contact the Conference Coordinator, Susan Owens, at
FAX: 713-743-3505, Phone: 713-743-8688, E-mail:

SUBJECTS: Spectral methods, high order finite differences and finite elements,
h-p version of finite elements and wavelet methods for partial differential
equations; applications in electronic chip manufacturing processes, fluid and
structural dynamics problems and biotechnology; efficient solvers and efficient
algorithms for spectral and high order methods; parallel computation aspects of
spectral and high order methods. Participants will include physicists,
mathematicians, and scientists and engineers from many other areas involving
spectral and high order methods.

ICOSAHOM 1995 is sponsored by the National Science Foundation, the University
of Houston, and the Texas Center for Advanced Molecular Computation.

LIMITED SUPPORT of up to $550 per person is available for six individuals.
A statement of research interests and need (1 page max) is enclosed.
Deadline: May 15. Only limited support is available.
Women and minorities are encouraged to apply.

ADDITIONAL INFORMATION AND REGISTRATION: Look at URL for the World Wide Web: or send email to


From: Trini Flores <>
Date: Mon, 01 May 95 14:36:21 EST
Subject: Symposium on Geophysical Inverse Problems

Symposium on Inverse Problems: Geophysical Applications
Conducted by SIAM with the cooperation of Gesellschaft fur
Angewandte Mathematik und Mechanik (GAMM)

December 16-19, 1995
Marriott Tenaya Lodge at Yosemite
Fish Camp, California


This symposium is the second of a series of symposia on
inverse problems and their applications. The first symposium held
June 27-July 1, 1994 in St. Wolfgang, Austria focused on applications
in diffusion processes. The second symposium will focus on the study of
connections of gravitational, seismic, electromagnetic, and thermal fields, with
the structure of our planet. These fields may be of natural origin, such as
oscillations caused by earthquakes and telluric currents, or they
may be induced as in seismic prospecting and geophysical prospecting by
electrical means. The interpretation of geophysical data leads to inverse
problems, i.e. from physical laws and the data of field measurements determine
the characteristics of the medium.

Symposium Themes

Groundwater Flow
Geophysical Prospecting
Electromagnetic Waves

Organizing Committee

Heinz W. Engl (Co-chair)
William Rundell (Co-chair)
David L. Colton
Alfred Louis

(One-hour plenary talks)

Underground Imaging of Electrically Conducting Plumes
James G. Berryman

A Geometrical Analysis of Duality Methods for the Inversion of Seismic Data
Guy Chavent

Inverse Problems for Groundwater Contamination and Petroleum Applications
Richard E. Ewing

Inverse Problems in Geodesy
Willi Freeden

Inverse Problems for Model-Data Synthesis in the Atmospheric and Oceanic Science
Michael Ghil

(Title to be determined)
Kurt J. Marfurt

The Inversion of Body Wave Attributes Derived from Seismic Refraction Data
Robert L. Nowack

The Mathematics of Velocity Analysis
William W. Symes

The symposium program and registration material will be available
in September, 1995.

For more information, contact:

3600 University City Science Center
Philadelphia, PA 19104-2688 U.S.A.
Telephone: 215-382-9800
Fax: 215-386-7999


From: Timothy Tautges <>
Date: Mon, 1 May 95 17:03:51 MDT
Subject: Meshing Roundtable


October 16-17, 1995
Albuquerque, New Mexico
Sponsored by Sandia National Laboratories

Mesh generation is an enabling technology for the finite element,
finite difference, and finite volume analysis methods. Automated
meshing techniques have the potential to significantly simplify the
pre-processing task, and as such, their development is of significant
interest to the finite element community. The 4th International
Meshing Roundtable will consist of technical presentations by
researchers working in the field of automated meshing algorithms, and
a panel discussion among leading researchers in this field. This
Roundtable will also include a meeting room with exhibits by mesh
generation software vendors.

This announcement and further information on this Roundtable is being
made available on the World Wide Web, at URL This page will contain
up-to-date information on Roundtable speakers, agenda, and other
useful information.


We are soliciting technical presentations from researchers working in
the field of mesh generation. Formal papers from these presentations
will be published in a conference proceedings. The topics to be covered
at this Roundtable include but are not limited to:

Unstructured or structured hexahedral meshing algorithms
Tetrahedral meshing algorithms
Adaptive meshing algorithms and experience
Geometry decomposition techniques
Grid relaxation techniques
Dirty geometry management


As a new feature of this year's Meshing Roundtable we are soliciting
vendors to exhibit and/or demonstrate their software. Vendor
representatives should contact Tim Tautges (
for more information.


For general information, contact:
Tammy Wilson
Conference Coordinator
MS 0833
Sandia National Labs
Albuquerque NM 87185-0833
( =
phone: (505) 844-0180
fax: (505) 844-8251


From: Jan Verwer <>
Date: Tue, 2 May 1995 10:14:46 GMT
Subject: Workshop on Innovative Time Integrators


This workshop takes place at CWI (Center for Mathematics and Computer Science)
in AMSTERDAM, The Netherlands. The focus is on nonstandard ODE and PDE
time-stepping techniques. Keywords for lectures are classified into two
main categories:

1. Aspects of special purpose time integration methods for systems
of multi-space dimensional time-dependent PDE problems.
2. Long time integration, dynamical systems and multiple time-scales.

The workshop is organized by CHUS SANZ-SERNA, Universidad de Valladolid, Spain
and JAN VERWER from CWI, with financial support provided by the EC/HCM program
'The Equations of Fluid Mechanics and Related Subjects'. This support covers
local expenses of participants. Part of the participants will belong to research
groups from the HCM network and part from outside the network. The number of
attendees is limited to approximately 35 of which about 3/4 will be from

The workshop is open to anyone interested, but, since there is limited space,
participation is by invitation only. Rules for participation can be found
in the FULL ANNOUNCEMENT which can be obtained from Mrs. Simone Panka-van der
Wolff (, subject HCM Workshop). Deadline for application, through
submitting a 3-page Latex abstract, is Nov. 1, 1995. This abstract has to be
mailed both to and (subject HCM Workshop).


From: Michael Heroux <>
Date: Thu, 4 May 95 23:55:37 CDT
Subject: Position at Cray Research

The Math/Algorithms Group at Cray Research has an open position
for an individual with extensive experience in scientific computing
to work with Cray applications analysts and customers on efficient
design and implementation of scientific and engineering applications.

Training and experience should include:

- A Ph.D. in Mathematics, Computer Science or related field.

- A strong background in numerical PDEs and numerical linear algebra,
especially direct and iterative methods for sparse linear systems.

- A knowledge of finite element, finite volume and spectral methods.

- Experience in efficient implementation of large scale scientific
applications on distributed and shared memory parallel/vector machines.

- The ability to work as a member of a team and consultant.

- Experience in a scientific or engineering area.

To apply please send CV and references to:

Mike Heroux
Cray Research, Inc.
655 Lone Oak Drive
Eagan, MN 55121

or by email:


From: SIAM <>
Date: Fri, 05 May 95 13:43:35 EST
Subject: Contents, SIAM Mathematical Analysis

SIAM J. Math. Anal.
Volume 24, Number 3

Elliptic-Parabolic Equations with Hysteresis Boundary Conditions
Ulrich Hornung and R. E. Showalter

Regularity for the Interfaces of Evolutionary p-Laplacian Functions
Hi Jun Choe and Jongsik Kim

Energy Estimates Relating Different Linear Elastic Models of a Thin
Cylindrical Shell II: The Case of Free Boundary
Jyrki Piila and Juhani Pitkaranta

A Global Existence and Uniqueness Theorem for a Model Problem in
Dynamic Elasto-Plasticity with Isotropic Strain-Hardening
A. Nouri and M. Rascle

On Coupled Integral H-Like Equations of Chandrasekhar
Jonq Juang

Analyse spectrale d'une bande acoustique multistratifiee I:
Principe d'absorption limite pour une stratification simple
Elisabeth Croc et Yves Dermenjian

Pyramidal Algorithms for Littlewood--Paley Decompositions
M. A. Muschietti and B. Torresani

Semiclassical Asymptotics beyond All Orders for Simple Scattering Systems
Alain Joye and Charles-Edouard Pfister

Bifurcation of Spatial Central Configurations from Planar Ones
Richard Moeckel and Carles Simo

A Mathematical Model of Traffic Flow on a Network of Unidirectional Roads
Helge Holden and Nils Henrik Risebro

Characterization of Lp-Solutions for the Two-Scale Dilation Equations
Ka-Sing Lau and Jianrong Wang

Interval Oscillation Conditions for Difference Equations
Q. Kong and A. Zettl

Nontensor Product Wavelet Packets in L2(Rs)
Zuowei Shen

Asymptotic Regularity of Compactly Supported Wavelets
Hans Volkmer


End of NA Digest