NA Digest Sunday, May 13, 1990 Volume 90 : Issue 19

Today's Editor: Cleve Moler

Today's Topics:


From: The Editor
Date: Sun May 13 20:42:07 PDT 1990
Subject: Re: Fortran, C++ or Something Else?

Editor's Note:
Baker Kearfott sent in a collection of responses to his query last
week about programming languages for teaching numerical analysis
and mathematical software. Judging from the number of respones,
there is quite a bit of interest in this topic. But it was not
clear that all the authors of responses intended their remarks for
general consumption. So, let's hold off another week and invite
anyone interested in responding for the NA News to send their comments
to Kerfott,, or to directly to
The responses in Baker's collection were from David Goldberg, Layne
Watson, George Corliss, Louis Howell, Bill Behrman, Paul Frederickson,
David Gay, David Hough, Cleve Moler and "Mike".


From: Nelson H. F. Beebe <>
Date: Wed 9 May 90 14:12:33-MDT
Subject: Machine-readable Bibliography Data Bases

With the help of the GNU awk utility, I recently converted
ACM Transactions on Mathematical Software (TOMS) algorithm
620 (References and Keywords for Collected Algorithms of the
ACM) to BibTeX format, with corrections, additions of
accents in names, and math mode sequences. This covers only
the 137 algorithms published in TOMS (one is from TOPLAS)
between June 1975 and June 1985 (I added the last year myself,
but ran out of steam in keeping it up-to-date); ordinary
papers in TOMS are NOT covered.

I would like to extend this work to cover the full contents
of TOMS, and other journals in numerical analysis and
computer science. Since this kind of work should only need
to be done once, the collections would naturally be made
FREELY AVAILABLE via standard electronic channels, including

I have software that can convert from UNIX bib/refer or
Bibix formats to BibTeX format. Tib is fairly close to
bib/refer, and Scribe format was the basis of the BibTeX
format design. Thus, collections in other formats can very
likely be massaged, largely automatically, into a consistent
BibTeX form, and with the help of a modest amount of awk
programming, back into other formats as required.

So, do any of you have bibliography data base files covering
journals in the cited areas that you would be willing to
share with the community? And perhaps those of you who are
journal editors can lean on the publishers to make, and keep
up-to-date, machine-readable copies of journal contents
which can be made freely available electronically.

Please correspond with me directly, and I'll summarize
results to this list. I will be out of the country for the
rest of May, so you won't get a response back from me until

--Nelson Bebee
Center for Scientific Computing
South Physics
University of Utah
Salt Lake City, UT 84112
(801) 581-5254


From: Nasser M. Khraishi <khraishi@isl.Stanford.EDU>
Date: Thu, 10 May 90 22:53:29 PDT
Subject: Character Arguments in the BLAS.

I have been lately trying to write a program in C that calls the
Fortran level 1, 2, and 3 BLAS that I obtained from netlib.

After reading the reports that go with the level 2 and 3 BLAS, I am
not convinced that the character argument is the best way of passing
information about the job to be performed by a subroutine.
My main objection is due to the fact that different "other language"
compilers, such as C, handle character strings passed to and from Fortran

Thus, the main issue of portability of BLAS is jeopardized by
the practice of passing the character strings. Actually, from reading
the reports, it is obvious that the authors ran into this issue even for
Fortran compilers, since case of characters is a factor for some compilers.

In view of this, and lacking an agreement on how different languages
should interact, I do think that going back to the more portable
practice of passing integers corresponding to JOB information
(as is currently done in LINPACK) is the better way to go.

Nasser M. Khraishi
Information Systems Lab.
Stanford University.


From: Aaron Fogelson <>
Date: Fri, 11 May 1990 15:28:03 MDT
Subject: Eigenvalues of Elliptic Systems

Dear Colleagues,

Does anyone know of software to determine the eigenvalues LAMDA of
a pair of coupled elliptic equations of the form:


where U = (u1,u2), the Laplacian is in two variables, D is a diagonal matrix,
and K is not a diagonal matrix. The first component of U, u1, satisfies
homogeneous Dirichlet boundary conditions while the second component, u2,
satisfies homogeneous Neumann conditions.

I am asking this question on behalf of my colleague Hans Othmer, so
please address replies to

Thank you.

Aaron Fogelson


From: Gene Golub <>
Date: Fri, 11 May 1990 17:49:07 PDT
Subject: Reports from Stanford

The following reports are available on a limited basis. Please send
your request to

Gene Golub

- Manuscript NA-90-04 April 1990
Quadratically constrained least squares and quadratic problems, by
Gene H. Golub and Urs Von Matt.

- Manuscript NA-90-05 May 1990
Jacobi matrices for sums of weight functions, by Sylvan Elhay, Gene
H. Golub and Jaroslav Kautsky.

- Manuscript NA-90-06 May 1990
The nonsymmetric Lanczos algorithm and controllability, by Daniel
Boley and Gene Golub.


From: Dianne Oleary <oleary@cs.UMD.EDU>
Date: Thu, 10 May 90 10:50:04 -0400
Subject: Monroe Martin Prize in Applied Mathematics


UNIVERSITY OF MARYLAND, College Park is pleased to announce the
IN APPLIED MATHEMATICS (including numerical analysis) by a young
research worker. Candidates must be residents of North America
and not more than 35 years of age at the filing deadline. Submitted
papers must be by a single author and have been published, or
accepted for publication, in the open literature subsequent to
January 31, 1985. The work must not have been performed in
connection with the completion of requirements for an academic

Applications from qualified candidates, or nominations, are
solicited for the Monroe Martin Prize. Entries should include a
copy of the paper or contribution, with a covering letter, and be
submitted on or before July 31, 1990, to


The award will be announced by November 1, 1990. The recipient
will be asked to present his or her work at the Monroe Martin Lecture
at the University of Maryland in December, 1990, and will be awarded
a prize of $1000 plus travel expenses.

The Monroe Martin Prize was established to commemorate the
achievements of Professor Monroe Martin, former Director of the
Institute for Fluid Dynamics and Applied Mathematics, and Chair
of the Mathematics Department at the University of Maryland.
Previous prize winners are Neil Berger, Marshall Slemrod, and
Jonathan Goodman.


From: Karl Knapp <>
Date: Mon, 7 May 90 13:01:55 CDT
Subject: Positions at Numerical Algorithms Group

The Numerical Algorithms Group, Inc (NAG) in suburban Chicago is
seeking one or more Numerical Analysts to join its technical staff. NAG
develops and markets the NAG Fortran Library, the world's most comprehensive
Mathematics and Statistics subroutine library. NAG also develops and
markets a number of other numerical software products: statistics
packages (GLIM & GENSTAT), the NAG Ada Library, a Graphics Library,
optimization packages, and other specialty libraries.

NAG is seeking candidates with an M.S. in Numerical Analysis or
related field with broad interests. The successful candidate(s) will
be responsible for supporting customers (primarily researchers), porting
products to new hardware/software environments, and developing new
products. Positions generally require occasional travel to the parent
office in Oxford, England. NAG has a network of Sun workstations, though
much of our work is done on external state-of-the-art systems, typically
vector and/or parallel computers. Previous UNIX experience is a big plus.
Starting salary: $35K - $48K. NAG Inc provides an excellent work
environment and has very good benefits: 401K plan, five weeks of
vacation, family medical and dental coverage.

NAG Inc is also seeking a Systems Manager to administer its internal
network of Sun workstations. Please contact NAG for further details.

Please send your resume, in confidence, to:

Numerical Algorithms Group, Inc
attn: Technical Manager
1400 Opus Place, Suite 200
Downers Grove, IL 60515-5702

phone: (708) 971-2337
FAX: (708) 971-2706


From: Jon Kettenring <>
Date: Tue, 8 May 90 17:21:09 -0400
Subject: Interface'91 Preliminary Notice


On the eve of Interface'90, planning is underway for the 1991 meeting, in
Seattle on Apr 21-24, 1991. The general purpose of these interface meetings
is to explore topics of interest to both the computer science and statistics
communities, topics such as visualization of data, parallel computing, neural
networks, large data problems, data base design, programming paradigms,
computer-human interfaces, numerical computations, algorithms, etc.

The purpose of this message is to solicit program ideas. Please send your
suggestions by e-mail to: (or bellcore!iface91).

If you'll be at Interface'90 in East Lansing, MI, on May 16-19 this year
and have program ideas for next year, please look me up.

(We're sending this notice several ways. Sorry if you get more than one copy!)

Jon Kettenring, Program Chair for Interface'91
Bellcore, 445 South Street,
Morristown, NJ, 07960. Email: or bellcore!jon


From: Milo Dorr <>
Date: Mon, 7 May 90 11:12:24 PDT
Subject: Advances in Mathematics, Computations, and Reactor Physics


American Nuclear Society (Mathematics & Computation Division
and Reactor Physics Division) International Topical Meeting:


April 28 - May 1, 1991
Green Tree Marriott
Pittsburgh, PA

Program Chairmain: Schedule:

I. K. Abu-Shumays Summaries Due Aug. 20, 1990
RT-Mathematics, 34F Full Papers Due Jan. 15, 1991
Bettis Atomic Power Laboratory Meeting Starts Apr. 28, 1991
P. O. Box 79
West Mifflin, PA 15122-0079 BANQUET - GUEST PROGRAMS - TOURS

AUTHORS: Submit 1000 word summaries to the program chairman.

Papers are solicited on all topics of interest to the membership
of the Mathematics & Computation Division and the Reactor
Physics Division of the American Nuclear Society, including
the following topics:

Advances in Supercomputers
Networks and Communications
Distributed Computing
Parallel Processing and Vectorization
Computing Practices/Software Engineering
Graphics and Animation
High Speed Computing Applied to Reactor Operation
Real-time Systems/Simulation Models
Artificial Intelligence and Expert Systems
Optimization of Nuclear Systems
Robotics and Manufacturing Applications
Reactor Design and Operation
Neural Networks

. . .

Advances in Computational Methods
Solution of Linear and Nonlinear Systems
Numerical Optimization Schemes
Deterministic and Stochastic Transport Methods
Numerical Methods for Radiative, Neutral, and
Charged Particle Transport
Diffusion Theory and Kinetics
Numerical Methods in Fluid Dynamics,
Structural Mechanics, and Thermal Hydraulics
Computational and Mathematical Problems of
Reactor Dynamics and Safety Analysis
Computational Aspects of Power Distribution,
Control, and Optimization
Methods for Probabilistic Risk Assessment

. . .

Reactor Physics, Theory, and Analysis Methods
Advanced Reactor Concepts/Designs
Space Reactor and Research Reactor Physics
Physics and Safety Aspects of Core Design
Reactor Design Validation and Operating Experience
Physics of Reactor Design
Physics of Criticality Safety
Point and Space-Time Core Models for Transient Analysis
Nuclear Plant Analyzers, Data Bases & Code Systems
Improvement and Validation of Plant Simulation Codes
Instrumentation, Control, and Monitoring
Fuel Cycle and Fuel Management
Shielding and Core Management
Nuclear Data
Integral Experiments - Measurements and Analysis
Reactor Physics Calculations Versus Measurements
Sensitivity Theory and Data Adjustment Techniques for
Reactor Analysis
Physics and Computational Aspects of Fusion

. . .

Paper and Poster Sessions on Methods, Algorithms,
and Computer Performance Benchmarks
Poster Session on Computer Programs which offer new
or improved methods or options (only code abstracts
required, such as published in NS&E)

. . .

Best Benchmarks and Programs Paper Awards
Best Student Paper Awards

The new Pittsburgh and its environs offer outstanding technical,
cultural, and scenic attractions. The city and the conference
sponsors cordially invite you to contribute to our 1991
Topical, and bring your family to participate in all other

I. K. ABU-SHUMAYS (PROGRAM CHAIRMAN, 412/476-6469, FTS 721-6469)


From: Douglas N. Arnold <>
Date: Sat, 12 May 90 22:12:18 EDT
Subject: Course on MODULEF Finite Element Library

Pennsylvania State University
July 16-20, 1989

Club MODULEF was created in 1974 to bring together researchers
from academics, government, and industry interested in developing
a finite element program library with the following features:

o access to source code, enabling users to add and modify

o modular organization

o rapid incorporation of current theoretical developments

At present the club has about 250 institutional members from around the
world. The MODULEF library presently contains over 300,000 lines of
FORTRAN code and is being expanded continually. Though based in
France, at the Institut National de Recherche en Informatique et en
Automatique (I.N.R.I.A.), there is presently a major effort to make all
the documentation available in English.

Courses in the use of MODULEF have been offered regularly in France.
Because of the increasing American membership in Club MODULEF, such a
course is being offered for the first time in the United States. The
aim of the course is to introduce new users to the MODULEF library and
teach them both to use existing modules and to incorporate new modules.


I. Using existing modules
o A sample physical problem
- variational formulation
- approximation by finite elements
- steps in the implementation

o The MODULEF library
- organization of the library
- principle modules
- interactive use
- writing a calling program

o Examples and hands-on use

II. Creating new modules
We will explain how to create new modules to
be used with existing modules in order to
expand the capability of the library


The course will be held at the Pennsylvania State University in
University Park, Pennsylvania. It will run from Monday to Friday,
July 16-20, 1990.


Tuition is $200 for full-time academic participants and $500 for
industrial participants. Tuition includes a copy of "MODULEF: A
Modular Library of Finite Elements", a 350 page tutorial and reference
PARTICIPANTS TO REGISTER. Tuition must be paid at the time of


DOUGLAS ARNOLD. Professor Arnold is a leading researcher in numerical
analysis and differential equations at the Pennsylvania State
University and is a user and contributor to the MODULEF library.

MICHEL BERNADOU. Professor Bernadou is a leading researcher in
numerical analysis and computational shell theory and the director of
the MODULEF project at I.N.R.I.A.

PAUL LOUIS GEORGE. Dr. George is a full time researcher with the
Modulef project at I.N.R.I.A. He works on mesh generation techniques
and is the principal author of the extensive mesh generation facilities
in the MODULEF code.

MARINIA VIDRASCU. Dr. Vidrascu has been a full-time researcher with
the MODULEF project at I.N.R.I.A. since its inception and is the
largest single contributor to the design and coding of the MODULEF
library. She works in nonlinear elasticity.

Questions about the course may be directed to

Douglas Arnold,
Department of Mathematics
Pennsylvania State University
University Park, PA 16802

telephone: (814) 865-0246
electronic mail:

For information on membership in Club MODULEF and obtaining the Modulef
software library contact

I.N.R.I.A electronic mail:
Domaine de Voluceau, B.P. 105
78153 Le Chesnay Cedex, FRANCE
Telephone : (1) 39635432


From: ICIAM '91 <>
Date: Mon, 7 May 90 09:38 EDT
Subject: Second International Conference on Industrial and Applied Mathematics

The Second International Conference on Industrial and Applied Mathematics

ICIAM Update

o More than 1200 applied and computational mathematicians, engineers and
scientists worldwide have already sent in requests for abstract and
minisymposium proposal forms for ICIAM 91.

o More than 35% of the requests are from outside the U.S.A.

o More than 60 representatives throughout the world are distributing ICIAM 91

o Proposals and abstracts already received encompass exciting and recent
advances in applied and computational mathematics.

o It's not too early!

ICIAM 91 is truly anticipated as an international forum for applied and
computational mathematicians and users of mathematics. Be a part of it....
organize a minisymposium....contribute a paper....come to Washington D.C. If
you still haven't told us of your interests, send in the attached form and
make sure you stay on the ICIAM 91 mailing list.

ICIAM 91--July 8-12, 1991, Washington, D.C., USA

Message to SIAM Members:

ICIAM 91 is an important international event that will occur every four years.
It is meant to spotlight the applied and computational community and focus
worldwide attention on mathematics professionals. There will be no SIAM Annual
Meeting in 1991. Make plans instead to come to ICIAM 91.

Invited Presentations*

Dynamic Energy Minimization and Phase Transformations in Solids
John M. Ball
Department of Mathematics
Heriot-Watt University, Scotland

Intermediate Asymptotics in Micromechanics
Grigorii I. Barenblatt
Institute of Oceanography
USSR Academy of Sciences, USSR

Mathematical and Computational Models of Vision
Michael Brady
Department of Engineering Science
University of Oxford, Great Britain

Aspects of Constrained Optimization
Andrew R. Conn
Department of Combinatorics and Optimization
University of Waterloo, Canada

Modulation Equations Arising in the
Mechanics of Continuous Media
Wiktor Eckhaus
Mathematisch Instituut
Rijksuniversiteit te Utrecht, The Netherlands

Some Industrial Problems Solved in the Framework of the
European Consortium of Mathematics in Industry
Antonio Fasano
Istituto Mathematico "U. Dini"
Universita Degli Studi (Firenze), Italy

Numerical Solution of Linear Algebraic Equations
Arising from Convection-Diffusion Equations
Gene H. Golub
Department of Computer Science
Stanford University, USA

Discrete Mathematics in Manufacturing
Martin Groetschel
Institut fuer Mathematik
Universitaet Augsburg, West Germany

Computational Methods for Chemically Reacting Flow
Thomas J.R. Hughes
Division of Applied Mechanics
Stanford University, USA

Interior Point Methods in Optimization
Narendra K. Karmarkar
Mathematical Sciences Research Center
AT&T Bell Laboratories, USA

Forced Nonlinear Surface Waves
Klaus W. Kirchgassner
Institut fuer Mathematik
Universitaet Stuttgart, West Germany

Optimal Control and Viscosity Solutions
Pierre-Louis Lions
Universite de Paris IX, France

Wavelets and Applications to Numerical Analysis
Yves F. Meyer
Universite Paris Dauphine, France

Dynamics of Patterns, Waves and Interfaces from the
Reaction-Diffusion Aspect
Masayasu Mimura
Department of Mathematics
Hiroshima University, Japan

Generation of Biological Pattern and the Formulation of
Morphogenetic Flows
James D. Murray
Center for Mathematical Biology, University of Oxford, Great Britain, and
Department of Applied Mathematics, University of Washington, USA

Radar Architectures--From Microwave Processing to Computational Power
Gabriel Ruget
Defense and Control Systems Division
Thomson CSF, France

Advances in Strongly Polynomial Linear Programming Algorithms
Eva Tardos
School of Operations Research and Industrial Engineering
Cornell University, USA

Applications of Massively Parallel Computing
David J. Wallace
Department of Physics
University of Edinburgh, Scotland

Computational Modeling of Flows in Porous Media
Mary F. Wheeler
Department of Mathematics
University of Houston, USA

Hierarchical Finite Elements and Related Preconditioners for
Elliptic Partial Differential Equations
Harry Yserentant
Fachbereich Mathematik
Universitaet Dortmund, West Germany

(*Titles of presentations tentative.)

Conference Topics

Applied and Computational Mathematics, Computer Science, Applied Probability
and Statistics, Scientific Computing, and Applications in Engineering and the
Biological, Chemical, and Physical Sciences.

Minisymposia Proposals

The ICIAM 91 Program Committee invites you, as a potential organizer, to
submit a proposal for a minisymposium.

A minisymposium is a session of 3-5 speakers focusing on a single topic.
Sessions may be scheduled for two or three hours. The organizer of a
minisymposium invites the speakers and decides on the topics they are to
address. If you would like to organize a minisymposium, you must submit a
proposal on an ICIAM minisymposium form. Acceptance of proposals will be by
decision of the ICIAM 91 Program Committee.

Deadline date for minisymposium proposals: July 31, 1990

Contributed Presentations - Posters and Lectures

Participate in ICIAM 91 by submitting a paper, which you may present in
lecture or poster format. The ICIAM 91 Program Committee is encouraging
contributors to present their papers in poster form to increase communication
among participants, foster the development of international friendships, and
reduce the need for large numbers of parallel sessions.

Authors will have approximately 15 minutes for contributed presentations
(lecture format), with an additional 5 minutes for questions. Alternatively,
they may elect the poster format that encourages interactive discussions with
individuals interested in their work using flip charts and other visual aids.

If you desire to present a paper (lecture or poster format), you must submit a
summary not exceeding 100 words on an ICIAM 91 contributed paper/poster
presentation form or facsimile. You may also submit an abstract via e-mail.
Macros are available in LaTeX or TeX. To receive macros via e-mail contact Papers will be reviewed by the program committee.
Everyone who submits a paper will be notified by mail regarding acceptance.

Deadline date for submission of contributed presentation forms: September 30,

Registration Information and Announcements

To obtain future announcements, the ICIAM 91 program, and registration
information, please contact SIAM.

If you have questions about ICIAM 91, write to:

ICIAM 91 Conference Manager
c/o SIAM
3600 University City Science Center
Philadelphia, PA 19104-2688, USA

FAX (TELECOPY): (215) 386-7999
PHONE: (215) 382-9800

Exhibitor Information

Booths and combined exhibits will be available to companies wishing to display
their products and/or services. To ensure that you receive an exhibitor kit,
please contact SIAM.


The conference will take place in Washington, D.C., at the Sheraton
Washington, where SIAM is holding a block of rooms for those wishing to stay
at the conferene site. The Sheraton is a modern air-conditioned hotel with
2,000 guest rooms, 25 meeting rooms, an exercise room, a large outdoor
swimming pool, and many restaurants and lounge areas. It is located on the
city's "Metro" subway system which provides easy inexpensive access to most
points of interest.

The hotel, on 12 acres of a wooded park-like setting, is within walking
distance of the National Zoo, a few blocks from the National Cathedral,
Embassy Row, historic Georgetown and Rock Creek Park. The Sheraton Washington
boasts a multi-lingual staff fluent in a total of 20 languages, foreign
currency exchange services, telex facilities, and electrical adapters. A U.S.
post office is on the premises.

For those seeking more economical accommodations, there will be a limited
number of dormitory rooms available on a first come first served basis at
George Washington University, just a few Metro stops away from the Sheraton

Program Committee

Robert E. O'Malley Jr., Chair
Robert Azencott
Franco Brezzi
Carlo Cercignani
John E. Dennis, Jr.
Peter J. Deuflhard
Iain S. Duff
J.C.R. Hunt
P.A. Raviart
Ivar Stakgold
Adriaan H.P. Van der Burgh
Wolfgang L. Wendland

Standing Committee

Roger Temam, Chair
Vinicio Boffi
C. William Gear, ICIAM 91 President
Douglas S. Jones
Patrick Lascaux
Robert M. M. Mattheij
James McKenna
Jean-Claude Nedelec
Helmut Neunzert
Ronald A. Scriven
Alberto Tesei
Wolfgang Walter

Sponsoring Societies

Gesellschaft fur Angewandte Mathematik und Mechanik (GAMM)
Institute of Mathematics and Its Applications (IMA)
Societa Italiana di Matematica Applicata e Industriale (SIMAI)
Societe de Mathematiques Appliquees et Industrielles (SMAI)
Society for Industrial and Applied Mathematics (SIAM)

With the cooperation of

Institut National de Recherche en Informatique et en Automatique (INRIA)

Hosted by SIAM


End of NA Digest