NA Digest Monday, August 12, 1996 Volume 96 : Issue 30

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: Cleve Moler <>
Date: Mon, 5 Aug 1996 08:40:55 -0400 (EDT)
Subject: Small Correction for Last Week's Digest

If you are archiving these Digests, please change the first line
after the mail header in last week's Digest to have the correct date
and issue number: Sunday, August 4, 1996 Volume 96 : Issue 29

-- Cleve


From: Allison Bogardo <>
Date: Thu, 08 Aug 96 15:53:33 EST
Subject: Call for Nominations for Wilkinson Prize

SIAM is seeking nominations for the 1997 James H. Wilkinson Prize
in Numerical Analysis and Scientific Computing. The prize is
awarded every four years for research in, or other contributions
to, numerical analysis and scientific computing during the six
years preceding the award. The purpose of the prize is to
stimulate younger contributors and to help them in their careers.

The prize consists of a $1,000 cash award and an engraved plaque
and will be awarded at the 1997 SIAM Annual Meeting to be held at
Stanford University, July 14-18.


Please send nominations to:

Allison Bogardo
3600 University City Science Center
Philadelphia, PA 19104-2688
Telephone: (215) 382-9800
Fax: (215) 386-7999


From: Joel Storch <>
Date: Mon, 5 Aug 1996 09:12:00 -0700
Subject: Change of Address for Joel Storch

I have recently left Draper Laboratory and am now at The Aerospace
Corporation in Los Angeles, CA.

Joel Storch
Phone: (310)336-5932


From: Zbigniew Woznicki <>
Date: Fri, 9 Aug 1996 11:19:57 +0200
Subject: Numerical Linear Algebra Preprint Available

Numerical Linear Algebra Preprint Available

The following paper in a LaTeX input file is available each week
from Monday to Friday via anonymous ftp on (,
file: pub/woznicki/axel.tex

Z. I. Woznicki: "Axelsson's results versus the splitting principles."



From: Wei Liu <>
Date: Mon, 5 Aug 96 14:01:46 -0400
Subject: Spheroidal Wave Functions

Is there any code dealing with spheroidal wave functions?
Many thanks for you information.



From: Peter Knipp <>
Date: Wed, 7 Aug 1996 17:17:03 -0400
Subject: Evaluation of an Exponential Integral

I am interested in evaluating the following function F for many
different values of A and B:

F(A,B)=definite integral (from x=0 to x=1) of [exp(Bx-A/x)]/x.

(The function F is used in Ewald's method for evaluating a two
dimensional lattice sum of Green's functions for Helmholtz'
equation.) The positive quantities A and B are often on the
order of unity (which obviously causes the integration to be
straightforward), but sometimes A can be quite small (which
causes a logarithmic divergence in F), and often A and/or B
are/is on the order of 1000.

My first question: Is F(A,B) "special" in the sense that it can
be related simply to a function named after a dead mathematician?

My second question: Does anyone have tips available for evaluating
F(A,B) in a reasonably efficient and accurate fashion?

Pete Knipp
Christopher Newport University


From: Werner C. Rheinboldt <>
Date: Mon, 05 Aug 1996 14:31:43 -0500 (EST)
Subject: MANPAK and DAEPAK, FORTRAN Packages

Prompted by many inquiries, I have decided to make available my
FORTRAN 77 packages for computations on implicitly defined manifolds
and for the solution of certain nonlinear differential algebraic
equations. They can now be downloaded from the netlib directory

MANPAK is a package of utility programs for computations with
submanifolds of R^n that are implicitly defined by a system of
nonlinear equations F(x) = 0. Here F is a mapping from R^n to
R^m, m < n, such that rank DF(x) = m for x in M = {u in E; F(u)=0}
whence M is a submanifold of R^n of dimension d = n - m. There are
routines for establishing several types of local parametrizations
(coordinate systems) and -- once available -- for computing points
on M with given local coordinates. In addition, there are routines
for computing the first and second derivatives of these local
parametrizations and some other quantities, such as sensitivity
measures and the second fundamental tensor.

DAEPAK is a package of subroutines based on MANPAK for solving
six types of algebraically explicit differential algebraic equations
(DAEs); that is, of DAEs in which either the algebraic equations
and/or the algebraic variables are explicitly specified. There are
initial value problem solvers for the
index 1 DAEs F(u,u',w,t) = 0,
index-1 DAEs A(u)u' = G(u), F(u) = 0,
with optional facility for handling impasse points
index-2 DAEs G(u,u',w,t) = 0, F(u,t) = 0
index-2 DAEs A(u,t)u' + B(u,t)w = G(u,t), F(u,t) = 0
index-3 DAEs A(u,u',t)u" + B(u,u',t)w = G(u,u',t), F(u,t) = 0
index-3 DAEs M(u,t)u" + D_uF(u,t)^T w = G(u,u',t), F(u,t) = 0
(Euler Lagrange problems)

Both packages utilize a set of support routines collected in a
separate package MANAUX. In addtion, sample drivers for the
DAE-solvers are provided. A README file gives further details.

I would welcome any feedback and especially any suggestions for
additions and improvements.


From: Alan Harvey <harvey@csli.Stanford.EDU?
Date: Mon, 5 Aug 1996 17:22:06 -0700
Subject: New Book on Dynamical Systems and Numerical Analysis

I'm pleased to announce the publication of a new book of interest to NA-NET:

Dynamical Systems and Numerical Analysis
A.M. Stuart and A.R. Humphries

[Cambridge Monographs on Applied and Computational Mathematics 2]
published by Cambridge University Press

Dynamical systems are pervasive in the modelling of naturally occurring
phenomena. Most of the models arising in practice cannot be completely
solved by analytic techniques; thus, numerical simulations are of
fundamental importance in gaining an understanding of dynamical systems.
It is therefore crucial to understand the behaviour of numerical
simulations of dynamical systems in order to interpret the data obtained
from such simulations and to facilitate the design of algorithms which
provide correct qualitative information without being unduly expensive.
These two concerns lead to the study of the convergence and stability
properties of numerical methods for dynamical systems.

The first three chapters of this book contain the elements of the theory of
dynamical systems and the numerical solution of initial-value problems. In
the remaining chapters, numerical methods are formulated as dynamical
systems, and the convergence and stability properties of the methods are
examined. Topics studied include the stability of numerical methods for
contractive, dissipative, gradient, and Hamiltonian systems together with
the convergence properties of equilibria, phase portraits, periodic
solutions, and strange attractors under numerical approximation.

This book will be an invaluable tool for graduate students and researchers
in the fields of numerical analysis and dynamical systems.

1. Finite Dimensional Maps
2. Ordinary Differential Equations
3. Numerical Methods for Initial Value Problems
4. Numerical Methods as Dynamical Systems
5. Global Stability
6. Convergence of Invariant Sets
7. Global Properties and Attractors Under Discretization
8. Hamiltonian and Conservative Systems

Published July 1996 700 pages 28 line diagrams 258 exercises
ISBN 0-521-49672-1 (Hardback)

For more information, consult the online catalogs at and, or send enquiries to


From: Heinz W. Engl <>
Date: Wed, 07 Aug 1996 22:10:53 EDT
Subject: New Book on Inverse Problems

The following book has just appeared:

Heinz W. Engl, Martin Hanke, Andreas Neubauer
Regularization of Inverse Problems
Kluwer, Dordrecht, 1996
ISBN: 0-7923-4157-0

The table of contents can be found on my WWW-Page (see below).

Heinz W. Engl, Linz (Austria)

Prof.Dr.Heinz W. Engl E-Mail:
Industriemathematik or
Institut fuer Mathematik secretary:
Johannes-Kepler-Universitaet Phone: +43-(0)732-2468; ext.9219 or 693,
Altenbergerstrasse 69 secretary: 9220; home: +43-(0)732-245518
A-4040 Linz Fax: +43-(0)732-2468855
Oesterreich / Austria Telex: 2-2323 uni li a
World Wide Web:


From: Bernd Fischer <>
Date: Thu, 8 Aug 1996 15:12:19 +0200
Subject: New Book on Iterative Methods

I am pleased to announce that my book:


Wiley-Teubner Series: Advances in Numerical Mathematics
ISBN: 0471 96796 3

is now available.

This book provides a systematic study for the solution of symmetric linear
systems by iterative methods with an emphasis on problems encountered when
solving indefinite systems.
One particular feature which distinguishes this book from the usual
literature on polynomial based iteration methods is its emphasis on the
properties of the underlying polynomials rather than more conventional
``matrix manipulations''. In particular, a development and discussion of
the properties of orthogonal polynomials leads to a unified and complete
analysis of state-of-the-art methods
for solving symmetric linear systems. These include the well-known CG, CR,
STOD, MCR, SYMMLQ, MINRES, LSQR, CRAIG, and Chebyshev iteration techniques.
For each of the numerical methods described a MATLAB implementation is provided,
the codes for which are available on the Internet.
In addition, the material covered includes schemes for computing the
asymptotic convergence rate of indefinite systems, for the estimation and
approximation of the eigenvalue distribution and for the approximation of
the A-norm of the error.
Finally, as a template for the occurrence of symmetric but indefinite
linear systems, a mixed finite element discretization of the Stokes problem
is discussed. The elementary derivations make this chapter accessible for
the non-specialist.
As the book is nearly self-contained, it can be used as an introduction to
the theory and implementation of polynomial based iteration methods.

Bernd Fischer
Institute of Mathematics
Medical University of Luebeck


From: Ruediger Weiss <>
Date: Fri, 9 Aug 1996 13:07:35 +0200 (CES)
Subject: New Book on Parameter-Free Iterative Linear Solvers

New book:
Parameter-Free Iterative Linear Solvers
Ruediger Weiss
published by Akademie Verlag

This book presents a unifying theoretical framework for iterative solvers
for linear systems. The focus is on parameter-free methods, i. e. methods
that can be applied without a priori knowledge of inner properties of
the systems. Despite the generality of the approach convergence estimates
and geometrical properties can be derived.

State-of-the-art solvers are surveyed and classified. Many new relations
between these methods are furnished leading to a deeper insight. Existing
gaps in the theory are pointed out and some visions for the design of more
efficient and robust solvers are given.

A valuation of established methods based on the theoretical results and
on numerical tests is given. Recipes for the application of particular
techniques are derived. The criteria are efficiency, robustness and the
possibility of parallelization.

This book has a triple aim: to simplify the analysis by a unified theory,
to survey and valuate state-of-the-art solvers and to give an introduction
to iterative solvers for non-experts in the field.

For more information turn to the web-site:

Dr. Ruediger Weiss
Universitaet Karlsruhe
Postfach 6980
76128 Karlsruhe
Tel: +49 721 608 4034
Fax: +49 721 32550


From: Biswa Datta <>
Date: Thu, 8 Aug 96 09:23:48 CDT
Subject: Reprint of Numerical Linear Algebra Book


I am pleased to announce that my book Numerical Linear Algebra and
Applications has now been reprinted ( third print) by my publisher
Brooks/Cole Publishing Company.
The MATLAB based software MATCOM that implements most of the algorithms
in the book has also been revised and updated.
I would like to take this opportunity to express my sincere gratitude
to those colleagues and students who helped me with corrections,
suggestions, etc., for this print.

The book can be used in both undergraduate and beginning level
graduate courses in Numerical Linear Algebra, Numerical Analysis,
and Applied Linear Algebra for Mathematics, Computer Science , and
Engineering courses.
If you are planning to use the book for your course this fall or
next spring, I urge you to obtain your complimentary copy from
the publisher, even if you have copies from the previous prints.

Here are the required information on the book :

Author : Biswa Nath Datta ( E-mail :
Publisher : Brooks/Cole Publishing Company (A Division of
International Thomson Publishing Inc.), 511 Forest Lodge Road,
Pacific Grove, California 93950, USA.
Tel : 1-800-423-0563 (for Review)
Fax : (415) 592-3342 (for Review)
Tel: 1-800-354-9706 (for Purchasing)
Fax : (606) 525-0978 (for Purchasing)


From: Leonid P. Lebedev <>
Date: Mon, 12 Aug 96 11:39:09 +0300
Subject: Seeking Collaborator or Position

I am seeking a collaborator or/and a research and teaching position
in Latin America or other countries. The area of interests is
mathematical problems in mechanics; in particular, justification
of numerical methods, the theorems of existence in nonlinear shell
theory, industry applications. I published over 20 papers. I am
coauthor of the book
Functional Analysis; Applications in Mechanics and Inverse Problems
(published by Kluwer Acad. Publishers, 1996).
Additional information will be given by request.

Yours sincerely, Leonid P.Lebedev
Rostov State University, Russia


From: Jerzy Wasniewski <>
Date: Sat, 10 Aug 1996 15:00:13 +0200 (METDST)
Subject: PARA96, The Final Call for the Registration and Hotel Reservation

Organized by:
UNI-C (The Danish Computing Centre for Research and Education) and
IMM (The Department of Mathematical Modelling of the Technical
University of Denmark) in collaboration with
Professor Jack Dongarra from the University of Tennessee and
Oak Ridge National Laboratory

Dates: August 18, 1996 - 2 tutorials:
1) Wavelets/Signal and Image Processing and
2) Programming Parallel Computers
August 19 - 21, 1996 - Workshop sessions
Place: UNI-C, The Danish Computing Centre for Research and Education,
The Technical University of Denmark, Bldg. 304, DK-2800 Lyngby,

The PARA96 meeting will be followed by a linear programming workshop,
August 21-24, 1996, organized by the Department of Mathematical Modelling.

31 speakers from several countries for the Workshop and Tutorial are
invited. 45 contributed talks from 20 countries in 15 parallel sessions
will be presented. Several experts from the computer manufactures are

The PARA96 workshop information are available by website:

More information are available from:
Jerzy Wasniewski, UNI-C, DTU, Bldg. 304
DK-2800 Lyngby, Denmark
tel: +4535878967, fax: +4535878990


From: Yuusuke Iso <>
Date: Mon, 5 Aug 96 18:18:53 JST
Subject: Two Conferences in Kyoto, Japan

"The Carleman Estimate and Inverse Problems"
December 11 - December 13, 1996
Department of Mathematics, Kyoto University (Kyoto, Japan)
Chairmen: Masahiro Yamamoto(Univ.of Tokyo), Yuusuke Iso(Kyoto Univ.)
Sponsor: Sanwa Systems Development Co.,Ltd.

Symposium topics include: uniqueness of solutions to
inverse problems for PDE's; conditional stability
estimates; mathematical and numerical treatment of
ill-posed problems.

Keynote speakers: A.Boukhgueim(Russia), V.Isakov(USA)
D.Tataru(USA), M.Kubo(Japan)

For further information: e-mail;

"Numerical and Mathematical Treatments for Elasticity"
January 8 - January 10, 1997
Research Institute for Mathematical Science
Kyoto University (Kyoto, Japan)
Chairman: Yuusuke Iso (Kyoto University)

The aim of the conference is to give an opportunity to
discuss some problems in elasticity by both mathematicians
and engineers. Especially, mathematical analysis for
anisotropic elasticity, numerical analysis for crack
problems and inverse problems in elasticity are
important topics in the conference.

Keynote speakers: V. Romanov(Russia), etc.

Call for Papers: Please submit an extended abstract
(not more than 2 pages) by November 15, 1996 to
to Prof. Y.Iso (Fax: +81-75-753-3711). The abstract
should be refereed. Proceedings of the conference will
be published by the domestic publication.

For further information: e-mail;


From: James Nagy <>
Date: Tue, 6 Aug 1996 12:05:57 -0500 (CDT)
Subject: Report on NA Workshop in Rousse, Bulgaria

JUNE 24-27, 1996

We attended the First Workshop on Numerical Analysis and Applications
in Rousse, Bulgaria, and would like to summarize our enjoyable
experience. Although there were participants from 22 countries,
the conference had a distinctly eastern European flavour which was
interesting, stimulating and fun. There were many fruitful
conversations with our colleagues, and we look forward to reading the
book of the proceedings, to be published by Springer-Verlag.

The main tracks of the workshop were:

1. Numerical Linear Algebra
2. Numerical Methods for Differential Equations
3. Numerical Modelling
4. High Performance Scientific Computing

Summary of Talks:

The programme consisted of distinguished key speakers who treated us
to expositions of their work of a seminal nature in many topics.

* Tracks 1 and 4: Numerical Linear Algebra and
High Performance Scientific Computing

Bisseling - Numerical linear algebra on bulk synchronous parallel computers
Dimov - Iterative monte carlo algorithms for linear algebra problems
Godunov - Spectral portraits of matrices and criteria for
spectrum splitting
Szulc - Convex combinations of matrices - nonsingularity and
Schur stability
Tyrtyshnikov - Mosaic ranks and skeletons
Vassilevski - Stabilized hierarchical basis preconditioning methods

In addition, there were 27 other talks, including a minisymposium on
Computations in Image Reconstruction and Restoration. Attendees of
these tracks were rewarded with many interesting talks on a wide spectrum
of topics.

* Tracks 2 and 3: Differential Equations and
Numerical Modelling

Brugnano - BVP methods for ODE's
D'yakonov - Operator problems in strengthened Sobolev spaces and numerical
methods for them
Griewank - High order stiff ODE solvers using automatic differentiation
Jovanovic - Interpolation techniques & convergence estimates for fdm
Karageorghis - Conforming spectral decomposition schemes
Kaschiev - Newton's method for solution of one complex eigenvalue
Lebedev - Explicit fdms with variable time steps for stiff systems
Puzynin - Iterative schemes and adaptive algorithms on the base of
continuous analogue of Newton's method for solving nonlinear equations
Shishkin - Diffusion flux for a class of singularly perturbed equations
Vabishchevich - Convection/diffusion problems
Zlatev - Large air pollution models

There were also 38 other contributed papers covering a wide spectrum of
related topics. The topic of stiffness seemed to be a recurring one. This
was seen in talks on Runge Kutta, Multi-step, variable step, automatic
differentiation and several other methods. The problem of stiffness appeared
in many of the numerical modelling applications that were presented.

We all enjoyed the "get together" dinner, the excursion to Veliko Tarnovo
and the official dinner in Gabrovo held in a warm and friendly Bulgaria.

P.S. Williams. London Guildhall University. U.K.
J.G. Nagy. Southern Methodist University. Dallas, Texas.


From: Trini Flores <>
Date: Wed, 07 Aug 96 13:05:19 EST
Subject: SIAM Conference on Sparse Matrices

Second SIAM Conference on Sparse Matrices
October 9-11, 1996
Coeur d'Alene Resort
Coeur d'Alene, Idaho

The conference program will be mailed in mid-August. A condensed
version of it, including information on transportation, hotel,and
registration is now available on the World Wide Web. Point your web
browser to

Important dates:

Monday, September 9, 1996 - deadline for hotel reservation
Wednesday, September 25, 1996 - deadline for advance registration

For additional information, contact SIAM. Telephone: 215-382-9800;
E-mail: For hotel reservation, contact Coeur
d'Alene. Telephone: 208-765-4000; 800-688-5253 (toll free US).


From: Tony Skjellum <tony@Aurora.CS.MsState.Edu>
Date: Thu, 8 Aug 1996 08:55:43 -0500
Subject: Scalable Parallel Libraries Conference

Scalable Parallel Libraries Conference III (SPLC96)
October 24-25, 1996
National Science Foundation Engineering Research Center
for Computational Field Simulation,
Mississippi State, Mississippi

The definitive workshop/conference on scalable libraries and
high performance middleware.

Revised Short Announcement and Call for Papers & Posters


The paper deadlines (see below) are extended.

Paper Deadline: September 15, 1996 REVISED
Notification of acceptance: October 1, 1996 REVISED
Final papers to publisher: October 30, 1996 REVISED

Poster Deadline: September 15, 1996
Notification of acceptance: October 1, 1996

It is STRONGLY preferred that proposal to present be sent by e-mail to

October 23, 1996
Pre-Conference Tutorial: (MPI: The Message Passing Interface)

A practical prerequisite for building scalable libraries is an understanding
of the MPI system. This advanced tutorial covers MPI uses for building
parallel software assuming basic MPI knowledge.


From: Lothar Reichel <>
Date: Thu, 8 Aug 96 16:38:28 EDT
Subject: Conference Honoring Bill Gragg

"ODE to Linear Algebra and Rational Approximation"
A conference on the occasion of William B. Gragg's 60th birthday

A conference is scheduled on Friday and Saturday, November 1-2, 1996, at
the Naval Postgraduate School in Monterey, California, on the occasion of
Bill Gragg's 60th birthday. The scientific program will focus on those areas
of mathematics and scientific computing where Gragg has made contributions
over the years. These include the numerical solution of ordinary differential
equations, linear algebra with applications, and rational approximation. Many
experts in these areas have already agreed to present talks. However, a few
slots may still be available. Contact Lothar Reichel ( if
you would like to join Bill Gragg's colleagues and friends on this occasion.

A special issue of the Journal of Computational and Applied Mathematics
dedicated to Bill Gragg is planned. Deadline for submission is November 2,
1996. Submitted papers will undergo the usual refereeing process. Contact
Greg Ammar ( or Daniela Calvetti
( if you would like to contribute a
paper to the special issue.

Information about the meeting will be posted at the web site

The organizing committee

Greg Ammar
Carlos Borges
Daniela Calvetti
Lothar Reichel


From: Gerardo Toraldo <>
Date: Fri, 9 Aug 1996 18:20:21 -0600
Subject: Conference on Optimization Software

Short Conference on
STATUS AND PERSPECTIVES''. Ischia, Italy 4 - 6 June 1997
H P S N O 97

The Centro di ricerche per il calcolo Parallelo e i Supercal-
colatori (CPS), a joint research center of the CNR (Consiglio
Nazionale delle Ricerche) and the University of Naples
``Federico II'', will host a short conference entitled ``High
Performance Software for Nonlinear Optimization on 4-6 June
1997 in Ischia (which is one of the islands in the bay of Naples),
The Conference will follow the HPSNO95 Conference that was held
in Capri in 95, whose main contributions were published in a
special issue of Computational Optimization and Applications.
The focus of the conference is the current state-of-the-art in
optimization software and, in particular, optimization
software for high performance computers. The conference will
provide an authoritative overview of the nonlinear optimization
field, including algorithms, software evaluation, implementation
issues, applications and future areas of research.
TOPICS OF INTEREST include: (Parallel) Computational experiments,
Large-Scale Constrained and Unconstrained Problems, Global Optimi-
zation, Quadratic Programming, Automatic Differentiation, Linear
and Nonlinear Least Squares, Semidefinite Programming, Numerical
Linear Algebra problems arising in Optimization, Solution of Large
Nonlinear Systems, Linear Complementary Problems, Nonsmooth
Optimization, Combinatorial Optimization, Applications in Sciences,
Management and Engineering.

INVITED SPEAKERS (Preliminary list)
Chris H. Bischof (Argonne National Laboratory, USA)
Nick Gould (Rutherford Appleton Laboratory)
Jorge J. More' (Argonne National Laboratory, USA)
Panos Pardalos (University of Florida, Gainesville, USA)
Mauricio G.C. Resende (AT&T Research, Murray Hill, USA)
Margaret Wright (AT&T Bell Laboratories, Murray Hill, USA)

A selected number of contributed presentation and posters will be
selected for the conference.
A subsequent announcement with deadlines and more details about
the conference will be made available by September 15.
for the Organizing Committee:
Prof. Almerico Murli or to Prof. Gerardo Toraldo
Universita` di Napoli "Federico II", Dip. di Matematica e Appl.,
Complesso Monte S.Angelo, Via Cintia 80126, Napoli - ITALY.


From: Ivan Slapnicar <>
Date: Sun, 11 Aug 1996 11:39:18 +0200 (MET DST)
Subject: Report on Workshop on Accurate Eigensolving and Applications

``International Workshop on Accurate Eigensolving and Applications''

Split, Croatia, July 12-16, 1996

A small group of specialists gathered together in the fascinating
city of Split on the Croatian Coast. The goal was to share their
knowledge of methods to solve the problem $Ax = \lambda Bx$
(or its singular value pendant) to as much accuracy as the data
warrants. In particular, there is a need to describe the structure
of matrices $A,B$, which determine their small eigenvalues (the
ones of interest in most, but not all applications) to high relative
accuracy. For example, consider the stiffness and mass matrices,
usually called $K$, $M$, obtained by finite differences
approximation. For one-dimensional problems the unassembled form
$NDN^t$ does determine the eigensolution $(NDN^t - \lambda M)v = 0$
to high relative accuracy, but more general finite element models,
especially those in high dimensions, are still undecided.

While it was already known that some factorizations like QR
may essentially improve the computed accuracy of the singular
values, we have now learned that the same is true of completely
pivoted LU - a deeply non-unitary process. Also the highly accurate
singular values of a bidiagonal were shown to tune the subsequent
inverse iterations to keep both their high speed and the necessary

The organizers, Professors Barlow, Slapnicar and Veselic, thought
hard about the format of the meeting. The participants fell into
2 groups, about 10 senior researchers active on these problems
(Demmel, Drmac, Hari, Ipsen, Parlett, Ruhe, and the organizers)
and about 15 graduate students who are writing dissertations in this
area. All participants stayed in one of the student dormitories of
the University of Split, the formal host of the Workshop. After
breakfast a small caravan of cars made its way to the Department
of ME, EE and Naval Arch. There were just 3 one hour talks in the
morning and at most one in the afternoon. This left time for
group discussions on Open Problems, Software and a possible

The local organizer, Ivan Slapnicar, had arranged an impressive
program of entertainment, including among others, a guided tour of
Diocletian's palace (from which the city grew in the 5th century),
a full day on a wonderful Zlatni Rat Beach on the island of Brac
and attending an open air performance of Verdi's Nabucco. The weather
was as usual, superb. The participants were looked after very well
by local organizers.

How can you tell a workshop from a regular meeting? Regular meetings
are dominated by overhead projectors, workshops make use of the
blackboard and all talks are peppered with questions from senior
researchers while graduate students can learn the habit of

Although there were only 15 talks in 4 working days there was a
consensus that the remaining time must allow more space for
informal technical work.
The list of talks and the abstracts are available at the web site

At the end of the Workshop, Jesse Barlow announced the follow-up
meeting at the Pennsylvania State University, State College, PA,
in two years.

Beresford Parlett


From: Margret Rothfuss <SciCom@IWR.Uni-Heidelberg.De>
Date: Mon, 5 Aug 1996 13:35:33 +0100
Subject: Position at University of Heidelberg

Interdisciplinary Center for Scientific Computing (IWR)
University of Heidelberg/Germany

The IWR offers two positions within the project jointly funded
by the Federal Ministry of Education, Science, Research and
Technology (BMBF), BASF AG and Hoechst AG. The work
involves the development of optimization methods and algorithms
for optimal experiment design in chemical kinetics and chemical
engineering. The algorithms are designed and implemented in
cooperation with scientists from BASF, Hoechst and
Fachhochschule Frankfurt.

The positions are available for a period of

two years.

The salary is regulated by the German BAT contract:


The ideal candidates have experience in the specification and
implementation of user friendly application software and a solid
background in numerical techniques (SQP) for nonlinear
constrained optimization problems (position 1) respectively
discretization methods for differential-algebraic equations
(position 2). Fluency in German or English is requested, a
doctorate degree in mathematics is desirable.

Application including a detailed CV, list of publications,
photograph and copies of diploma and an abstract of the PhD
thesis should be sent to

Dr. Johannes Schloeder
Keyword: Optimal Experiment Design
Universitaet Heidelberg
Im Neuenheimer Feld 368
D-69120 Heidelberg


From: Willard Miller <>
Date: Mon, 5 Aug 1996 08:09:17 CST
Subject: Position at University of Minnesota

University of Minnesota

The Board of Governors of the Institute for Mathematics and its
Applications (IMA) and the University of Minnesota seek a new Director
of the IMA for an appointment beginning September 1997. The new
Director will be offered a Professorship in the School of Mathematics
of the University of Minnesota. Candidates should have the
qualifications to provide scientific and administrative leadership to
the IMA. Among the duties are program development, fund raising and
preparation of proposals, collaboration with industry, and formulation
of plans for evolution of the Institute. Distinguished academic
credentials and a record of scientific leadership are required. Salary
and term as Director of the IMA are negotiable.

The IMA was established in 1982 with the support of the National
Science Foundation. The Institute has a two-fold mission: To identify
problems and areas of mathematical research needed in other sciences;
and to encourage the participation of mathematicians in these areas of
application by providing settings conducive to the solution of such
problems, and by demonstrating that first-rate mathematics can make a
real impact in the sciences. The IMA scientific programs allow
mathematicians and other scientists to share a stimulating research
environment. The IMA runs annual and summer theme programs with many
workshops and nearly a thousand participants a year, a Postdoctoral
Training Program in the area of the theme year, an ongoing Industrial
Postdoctoral Program, and an Industrial Problems Seminar where
industrial scientists are invited to present problems. There are 13
members in the IMA Participating Corporations Program, a formalized
relationship between the IMA and industrial scientists; and 29 IMA
Participating Institutes, a consortium of universities who provide
valuable support and guidance.

Nominations and applications should be sent to Professor Willard
Miller, Chair, IMA Director Search Committee, IMA, 514 Vincent Hall,
University of Minnesota, Minneapolis, MN 55455. Consideration of
applications will begin December 15, 1996.

The University of Minnesota is an Equal Opportunity Educator and

Willard Miller, Jr.
Professor and Associate Dean for Finance & Planning
Institute of Technology, University of Minnesota
105 Walter Library, 117 Pleasant Street S.E.
Minneapolis, Minnesota 55455 612-624-2006 FAX: 612-624-2841


From: Jamie Washburn <>
Date: Tue, 6 Aug 1996 15:40:19 -0400
Subject: Position at Cornell Theory Center

The Cornell Theory Center, a National Science Foundation center for high
performance computing in science and engineering, seeks candidates for the
position of Scientific Applications Consultant.

This person will provide highly specialized and complex consulting and
scientific applications support to the Theory Center's national and
international scientific and research community; provide in-depth
consulting and advice to researchers in enabling scientific applications on
parallel architectures, analyzing program execution, and optimizing
performance; work closely with users in analyzing problems and developing
solutions; and investigate, resolve or refer reported problems.

Minimum Position Requirements:
BS in a scientific or engineering discipline, or equivalent combination of
education, research, and experience. Advanced degree strongly preferred.
3-5 years FORTRAN or C programming experience in a scientific environment.
Demonstrated ability to work with researchers using high performance
computing. Scientific applications, parallel programming, and UNIX
experience highly desirable. Excellent interpersonal, communication, and
organizational skills. For consideration, send letter and resume to Julia
Addy, Theory Center, Frank H.T. Rhodes Hall, Cornell University, Ithaca, NY
14853-3801 AA/EOE

For additional information about the Theory Center, please refer to our
World Wide Web home page located at:

Jamie Washburn
Human Resources
Theory Center, Cornell Univ.
538 Frank H.T. Rhodes Hall
Ithaca, NY 14853-2601
Phone: (607) 254-8701 Fax: (607) 254-8888


From: Susanne Steinborn <>
Date: Fri, 9 Aug 1996 14:34:26 +0200 (MET DST)
Subject: Positions at Technical University in Clausthal, Germany


Research positions are available immediately at the Institute of Welding
Technology (ISAF) of the Technical University in Clausthal.

Successful candidates have experience in one or more of the areas
welding technology, materials science, automation of processes, mathematics
and information sciences.

Work will be done in research projects of fusion welding and cladding with
automated welding techniques including the numerical simulation and control
of processes.

The salary is regulated by the German BAT contract.
The appointment will be nonpermanent.

Women and handicapped persons are especially encouraged to apply.

Application, including curriculum vitae, relevant certificates and an
outline of scientific background should be sent to:

Dr.-Ing. B. Bouaifi
TU Clausthal
Im Schleeke 112
D-38640 Goslar

Tel: (05321) 3470-0
Fax: (05321) 3470-40
WWW homepage:


From: Wayne Joubert <>
Date: Tue, 6 Aug 1996 13:03:04 -0600 (MDT)
Subject: Postdoctoral Positions at Los Alamos National Laboratory

Post Doctoral Fellowship Positions
Los Alamos National Laboratory

The Scientific Computing Group (CIC-19) at Los Alamos National Laboratory is
currently seeking highly motivated Ph.D. graduates to participate in
the Post Doctoral Fellowship program. Graduates with experience in any
or all of the following categories are encouraged to apply:

Numerical Linear Algebra
Parallel Computer Programming

Experience with Fortran 90, HPF, PVM and MPI on parallel machines such as
the Cray T3D, Connection Machine CM-5, IBM SP-2, SGI Power Challenge and
workstation clusters is desirable. Experience with iterative linear
solver methods such as conjugate gradient methods and incomplete Cholesky
preconditioners is also desirable.

A Ph.D. completed within the last three years or soon to be completed is
required. Appointments have a duration of two years, with the possibility
of extension for a third year.

To receive fullest consideration for all fellowship possibilities, we must
receive your application as soon as possible. Interested individuals are
encouraged to send a copy of their resume by e-mail, fax, or U.S. mail to:

Wayne Joubert
Los Alamos National Laboratory
Group CIC-19, MS B256
Los Alamos, NM 87545

FAX: 505-667-1126

Los Alamos National Laboratory is an equal-opportunity employer.


From: Tamar Schlick <schlick@acf4.NYU.EDU>
Date: Fri, 9 Aug 1996 17:48:26 -0400 (EDT)
Subject: Postdoctoral Position at New York University

Postdoctoral Position
in Computational Mathematics/Structural Biology
at New York University

Applications are invited for a postdoctoral
position involving molecular and Brownian
dynamics simulations of DNA supercoiling
at New York University.

Ph.D. recipients in various areas are welcome,
but a solid background in computing and applied
mathematics, along with strong knowledge in
biophysics or chemistry, is highly desired.

Applications in writing, including full curriculum
vitae, statement of research interests, and three
reference letters should be sent to: Tamar Schlick,
Department of Chemistry and the Courant Institute
of Mathematical Sciences, New York University,
251 Mercer Street, New York, NY 10012 USA.

For further information, contact T.S. at the
email address , by phone
(212-998-3116) or fax (212-995-4152).
You are also invited to check out the group's
home page on:
for related information.


From: Richard Brualdi <>
Date: Mon, 5 Aug 1996 11:03:56 -0500 (CDT)
Subject: Contents, Linear Algebra and its Applications

Contents, Volume 245

Jin-Hsien Wang (Kaohsiung, Taiwan, Republic of China)
Product of Invertible Operators of Quadratic Type 1
Stewart W. Neufeld (Winnipeg, Canada)
A Diameter Bound on the Exponent of a Primitive Directed Graph 27
Ra#l Andrade and Alicia Labra (Santiago, Chile)
On a Class of Baric Algebras 49
I. Gohberg (Tel-Aviv, Israel), M. A. Kaashoek, and J. Kos (Amsterdam, The
The Asymptotic Behavior of the Singular Values of Matrix Powers and
Applications 55
Naum Krupnik (Ramat-Gan, Israel) and Steffen Roch (Leipzig, Germany)
On Invertibility Symbols in Banach Algebras 77
J. Rudolph (Stuttgart, Germany)
Duality in Time-Varying Linear Systems: A Module Theoretic Approach 83
J. W. Hoffmann (Kaiserslautern, Germany) and P. A. Fuhrmann (Beer Sheva, Israel)
On Balanced Realizations of Bounded Real and Positive Real Functions 107
Zhi-hao Cao, Jin-jun Xie (Shanghai, People's Republic of China), and Ren-Cang
Li (Berkeley, California)
A Sharp Version of Kahan's Theorem on Clustered Eigenvalues 147
Kirill A. Kopotun (Edmonton, Alberta, Canada)
A Note on the Convexity of the Sum of Subpermanents 157
Keith G. Woodgate (London, United Kingdom)
Least-Squares Solution of F = PG Over Positive Semidefinite Symmetric P 171
Sambhavi Lakshminarayanan (New York, New York) and R. Chandrasekaran
(Richardson, Texas)
Tessellation and g-Tessellation of Circulants, Q6, and Qt6 191
Ettore Fornasini and Maria Elena Valcher (Padova, Italy)
On the Spectral and Combinatorial Structure of 2D Positive Systems 223
Patrick J. Rabier and Werner C. Rheinboldt (Pittsburgh, Pennsylvania)
Classical and Generalized Solutions of Time-Dependent Linear Differential-
Algebraic Equations 259
M. Koppinen (Turku, Finland)
Three Automorphism Theorems for Triangular Matrix Algebras 295
Patrick J. Rabier and Werner C. Rheinboldt (Pittsburgh, Pennsylvania)
Time-Dependent Linear DAEs With Discontinuous Inputs 305
Nicolau C. Saldanha and Carlos Tomei (Rio de Janeiro, Brazil)
The Accumulated Distribution of Quadratic Forms on the Sphere 335
J. A. Dias de Silva (Lisboa, Portual)
New Conditions for Equality of Decomposable Symmetrized Tensors 353
Author Index 373


From: SIAM <>
Date: Mon, 05 Aug 96 12:12:20 EST
Subject: Contents, SIAM Review

SIAM Review
September 1996, Volume 38, Number 3

On Projection Algorithms for Solving Convex Feasibility Problems
Heinz H. Bauschke and Jonathan M. Borwein

Conjugate Gradient Methods for Toeplitz Systems
Raymond H. Chan and Michael K. Ng

Modelling the Stem Curve of a Palm in a Strong Wind
C. Philipsen, S. Markvorsen, and W. Kliem

Erratum and Reformulation: On the Stem Curve of a Tall Palm in a Strong Wind
Donald F. Winter

Catastrophe Theory Implications for Rightsizing when Planning Interim Solutions
for Progressing from a Partial Mainframe to Client-Server Distributed Databases:
3D Previewing of Possible Problems
Barry S. Thornton and W. T. Hung

The Use of Linear Programming in the Construction of Extremal Solutions to
Linear Inverse Problems
Stephen P. Huestis

The Matrix Exponential
I. E. Leonard


Limit Theorems of Probability Theory: Sequences of Independent Random Variables
(Valentin V. Petrov), Wlodzimierz Bryc

Parametric Statistical Theory (Johann Pfanzagl), Morris L. Eaton

Observational Studies (Paul R. Rosenbaum), Daniel F. Heitjan

Moving Finite Elements (M. J. Baines), Weizhang Huang

Optimization and Dynamical Systems (Uwe Helmke and John B. Moore), Yuji Kodama

Asymptotics of Nonlinearities and Operator Equations (Alexander M.
Krasnosel'skii), Jean Mawhin

Math into LATEX: An Introduction to LATEX and AMS-LATEX (George Grotzes),
Kenneth R. Meyer

An Introduction to Infinite-Dimensional Linear System Theory (R. F. Curtain and
H. J. Zwart), K. A. Morris

Modeling, Analysis and Control of Dynamic Elastic Multi-Link Structures (J. E.
Lagnese, G. Leugeuring, and E. J. P. G. Schmidt), Serge Nicaise

Introduction to Perturbation Methods (M. H. Holmes), Robert E. O'Malley, Jr.

Global Dynamics, Phase Space Transport, Orbits Homoclinic to Resonances, and
Applications (Stephen Wiggins), Kenneth J. Palmer

Perturbation Methods in the Computer Age (David C. Wilcox), R. H. Rand

The Technique of Pseudodifferential Operators (H. O. Cordes), Michael Taylor

Multigrid Methods for Finite Elements (V. V. Shaidurov), Jinchao Xu

Multigrid Methods for Process Simulation (W. Joppich and S. Mijalkovic), Jinchao

Mathematical Go Chilling Gets the Last Point (Elwyn Berlekamp and David Wolfe),
Zhiping You and James A. Yorke





From: SIAM <>
Date: Wed, 07 Aug 96 10:40:21 EST
Subject: Contents, SIAM Numerical Analysis

SIAM Journal on Numerical Analysis
OCTOBER 1996, Volume 33, Number 5

Fast Algorithms for Polynomial Interpolation, Integration, and Differentiation
A. Dutt, M. Gu, and V. Rokhlin

The Runge-Kutta Theory in a Nutshell
Peter Albrecht

A First-Order Exactly Incompressible Finite Element for Axisymmetric Fluid Flow
Barry Bernstein, Kathleen A. Feigl, and Elwood T. Olsen

Two-Grid Discretization Techniques for Linear and Nonlinear PDEs
Jinchao Xu

On the Superconvergence of Galerkin Methods for Hyperbolic IBVP
David Gottlieb, Bertil Gustafsson, Pelle Olsson, and Bo Strand

Convergence of a Boundary Integral Method for Water Waves
J. Thomas Beale, Thomas Y. Hou, and John Lowengrub

Analysis and Application of Fourier-Gegenbauer Method to Stiff Differential
L. Vozovoi, M. Israeli, and A. Averbuch

Piecewise Linear Petrov-Galerkin Error Estimates for the Box Method
Thomas Kerkhoven

The Fourier-Finite-Element Method for Poisson's Equation in Axisymmetric Domains
with Edges
Bernd Heinrich

The Approximate Solution of Defining Equations for Generalized Turning Points
A. Griewank and G. W. Reddien

Sequential Regularization Methods for Higher Index DAEs with Constraint
Singularities: The Linear Index-2 Case
Uri M. Ascher and Ping Lin

A New Class of Discretization Methods for the Solution of Linear
Differential-Algebraic Equations with Variable Coefficients
Peter Kunkel and Volker Mehrmann

Extra Smoothness Requirements for Galerkin Methods for the Wave Equation
Mark C. Haase

Parallel Treatment of a Class of Differential-Algebraic Systems
J. Bahi, E. Griepentrog, and J. C. Miellou

Uniqueness and Nonuniqueness for the Approximation of Quasilinear Elliptic
N. Andre and M. Chipot

Quadrature Error Bounds with Applications to Lattice Rules
Fred J. Hickernell

A Multiple-Exchange Algorithm for Complex Chebyshev Approximation by Polynomials
on the Unit Circle
Ching-Yih Tseng

High-Order Positivity-Preserving Kinetic Schemes for the Compressible Euler
J. L. Estivalezes and P. Villedieu

Qualitative Analysis of Newton's Flow
Vladimir Janovsky and Viktor Seige


End of NA Digest