NA Digest Tuesday, February 22, 2000 Volume 00 : Issue 08

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: Petter Abrahamsen <>
Date: Mon, 14 Feb 2000 09:10:21 +0100
Subject: LAPACK++ and Cholesky decomposition

We are currently developing some C/C++ programs that include solving
symmetric linear positive definite (SPD) systems. In addition, we
need the Cholesky decomposition of the associated SPD matrices for
further use in some matrix expressions. Speed is quite critical since
we will solve many different fairly large equation systems based on
the same Cholesky decomposition.

We have previously used LAPLACK++ for solving linear equation systems
and thought this would be an appropriate package. However, we have run
into major problems with our installation. Partly we have problems
figuring out which methods are available for Cholesky decomposition from
the somewhat erroneous and lacking documentation. Also, we have problems
compiling some of the code and even linking it because some header files
seem to lack the associated .cc files necessary to make an archive. (Is
this a local problem or a problem with the LAPACK++ package?).
So, is there anybody who have successfully used LAPACK++ for Cholesky
decomposition of SPD matrices and then used the solution in further
matrix operations including solving linear equation systems?
Alternatively, have any of you suggestions to an alternative free
library for use in C++ or C? (Preferably C++).

Sincerely yours
Petter Abrahamsen
Norwegian Computing Center
P.O.Box 114 Blindern
N-0314 Oslo, Norway


From: Constructive Approximation <>
Date: Wed, 16 Feb 2000 16:41:52 -0500 (EST)
Subject: Selected Papers of J.L. Walsh

We are pleased to announce the availability of "Joseph L. Walsh
Selected Papers", published by Springer-Verlag (ISBN 0-387-98782-7),
with commentaries by Q.I. Rahman, F. Schipp, P. Gauthier, D. Gaier,
T.J. Rivlin, E.B. Saff, and W. Schempp.


From: Craig Loehl <>
Date: Thu, 17 Feb 2000 11:32:20 EST
Subject: Global Optimization for Mathematica

Global Optimization 3.0 Released

Loehle Enterprises announces the release of v. 3.0 of Global Optimization for
Mathematica. Version 3.0 features 4 functions including a generalized
hill-climbing algorithm that accepts any mix of real, integer, and discrete
variables and is very robust to local minima, a 0-1 integer solver using the
Interchange Method, and an allocation solver for nonlinear investment
problems. It is ideal for engineering design problems, finance, nonlinear
model fitting, and other applications. For more information contact
630-579-ll90, or


From: Sharon Filipowski <>
Date: Fri, 18 Feb 2000 08:17:02 -0800
Subject: Optimization Prize for Young Researchers

CALL FOR NOMINATIONS: Optimization Prize for Young Researchers

PRINCIPAL GUIDELINE. The Optimization Prize for Young Researchers,
established in 1998 and administered by the Optimization Section (OS)
within the Institute for Operations Research and Management Science
(INFORMS), is awarded annually at the INFORMS Fall National Meeting
to one (or more) young researchers for the most outstanding paper in
optimization that is submitted to or published in a refereed professional
journal. The prize serves as an esteemed recognition of promising
colleagues who are at the beginning of their academic or industrial career.

DESCRIPTION OF THE AWARD. The Optimization award includes a cash
amount of US$1,000 and a citation certificate. The award winners will be
invited to give a one-hour lecture of the winning paper at the INFORMS
Fall National Meeting in the year the award is made. It is expected that
the winners will be responsible for the travel expenses to present the paper
at the INFORMS meeting.

ELIGIBILITY. The authors and paper must satisfy the following three
conditions to be eligible for the prize:

(a) the paper must either be published in a refereed professional journal
no more than three years before the closing date of nomination, or be
submitted to and received by a refereed professional journal no more than
three years before the closing date of nomination;

(b) all authors must have been awarded their terminal degree within five
years of the closing date of nomination;

(c) the topic of the paper must belong to the field of optimization in its
broadest sense.

NOMINATION. Nominations should be sent before July 15, 2000 to
Robert J. Vanderbei
Dept. of Operations Research and Financial Engineering
Princeton University
Princeton, NJ 08544
Nominations should be accompanied by a supporting letter.


From: Jud Wolfskill <wolfskil@MIT.EDU>
Date: Fri, 18 Feb 2000 14:14:07 -0400
Subject: Two New Books, Using MPI and Using MPI-2

The following is a book which readers of this list might find of
interest. For more information please visit

Using MPI
Portable Parallel Programming with the Message Passing Interface
second edition
William Gropp, Ewing Lusk, and Anthony Skjellum

Using MPI-2
Advanced Features of the Message Passing Interface
William Gropp, Ewing Lusk, and Rajeev Thakur

The Message Passing Interface (MPI) specification is widely used for
solving significant scientific and engineering problems on parallel
computers. There exist more than a dozen implementations on computer
platforms ranging from IBM SP-2 supercomputers to clusters of PCs
running Windows NT or Linux ("Beowulf" machines). The initial MPI
Standard document, MPI-1, was recently updated by the MPI Forum. The
new version, MPI-2, contains both significant enhancements to the
existing MPI core and new features.

Using MPI is a completely up-to-date version of the authors' 1994
introduction to the core functions of MPI. It adds material on the new
C++ and Fortran 90 bindings for MPI throughout the book. It contains
greater discussion of datatype extents, the most frequently
misunderstood feature of MPI-1, as well as material on the new
extensions to basic MPI functionality added by the MPI-2 Forum in the
area of MPI datatypes and collective operations.

Using MPI-2 covers the new extensions to basic MPI. These include
parallel I/O, remote memory access operations, and dynamic process
management. The volume also includes material on tuning MPI
applications for high performance on modern MPI implementations.

William Gropp and Ewing Lusk are Senior Computer Scientists in the
Mathematics and Computer Science Division at Argonne National
Laboratory. Anthony Skjellum is Associate Professor of Computer Science
and Director of the High Performance Computing Laboratory at
Mississippi State University. Rajeev Thakur is an Assistant Computer
Scientist in the Mathematics and Computer Science Division at Argonne
National Laboratory.


From: Gian Antonio D'Addetta <>
Date: Tue, 15 Feb 2000 13:05:29 +0100 (MET)
Subject: Symposium on Continuous and Discontinuous Modelling of Materials

International Symposium on Continuous and Discontinuous Modelling of Cohesive
Frictional Materials in Stuttgart, Germany from 27th to 28th of April.

Dear Colleagues,

We would like to invite interested researchers working in the field of
computational material modelling of geomaterials, experimental investigations of
geomaterials or related fields to participate in the symposium "Continuous and
Discontinuous Modelling of Cohesive Frictional Materials".

This symposium focuses on various aspects of modelling techniques for
cohesive frictional materials and is planned as an informative seminarlike
meeting. Following the preceding workshops of the research group
"Modelling of Cohesive Frictional Materials" this meeting is projected as a
subject-spanning symposium with international audience. Invited
national and international experts in targeted research areas will enlight
current developments and problems in the numerical modelling of
cohesive frictional materials and provide a deeper understanding of the
micro and macro description of geomaterials. Lectures from the members
of the research group "Modelling of Cohesive Frictional Materials" will
give an overview on the research topics of this group. Advances reached
within the first research period of two years will be presented as well as
the cooperation between the various institutes.

We would very much appreciate if you could directly approach potential
participants and inform them about this symposium.

For further information on invited guests and a preliminary list of lectures
of the symposium please visit the CDM2000-homepage at:

Thank you for your effort.

Gian Antonio D'Addetta
Institute of Structural Mechanics - University of Stuttgart - Germany


From: Sabrina Nordt <>
Date: Tue, 15 Feb 2000 16:45:13 +0000
Subject: Workshop on Numerical Methods for Free Boundary Problems

VIII. Workshop on
Numerical methods for free boundary problems
Berlin (Germany), October 5 --7, 2000

organized by
Freie Universitaet Berlin, Institut fuer Mathematik I
Weierstrass Institute for Applied Analysis and Stochastics
(WIAS), Berlin

Scientific committee: G. Dziuk (Freiburg)
C.M. Elliott (Brighton)
R.H. Nochetto (College Park)
C. Verdi (Milano)

Organizing committee: R. Kornhuber (FU-Berlin),
J. Sprekels (WIAS, Berlin)

The meeting is intended to stimulate the interplay of
mathematical modelling, numerical analysis, numerical
computation and real world problems in the field of
free boundary problems, such as phase transformations,
superconductivity or fluid interfaces.

The workshop continues the series of meetings on
numerical and computational aspects of free boundary
problems started with the ESF/FBP scientific program.
Preliminary workshops took place at Jyvaeskylae (1993),
Pont-a'-Mousson (1994), Pavia and Freiburg (1995),
Lamoura (1996), Ittingen (1997), Braga (1998), and Freiburg (1999).

A tentative list of invited speakers contains:

R.F. Almgren, E. Baensch, J.W. Barret, J. Blowey, F. Bornemann,
J.-M. Vanden-Broeck, P. Colli-Franzone, K. Deckelnick, G. Dziuk,
C.M. Elliott, W. Engl, D. Hoemberg, K.-H. Hoffmann, M. Kimura,
O. Klein, R. Klein, R. Krause, J. Lowengrub, R.H. Nochetto,
M. Paolini, L. Prigozhin, J. Rappaz, G. Savare', A. Schmidt,
G.J. Schmitz, J.A. Sethian, M. Shelley, V. Styles, T. Tiihonen,
A.-K. Tornberg, C. Verdi

The conference fee of 75 DM includes coffee/tea and conference dinner.

The workshop is sponsored by - Deutsche Forschungsgemeinschaft*
- WIAS Berlin
- TMR*
* grant pending
Contact address:

Sabrina Nordt | Phone: #49 (30) 838-75407
Freie Universitaet Berlin | Fax: #49 (30) 838-54977
Fachbereich Mathematik + Informatik | e-mail:
Institut fuer Mathematik I |
Arnimallee 2-6, |
D-14195 Berlin (Germany) |



From: Stephane Gerbi <>
Date: Wed, 16 Feb 2000 13:34:24 +0100
Subject: Savoisiennes Numerical Analysis Days


to be held 25-26 may 2000 at:
University of Savoie, Campus Scientifique, Le Bourget du Lac, France

Organised by C. Bourdarias, S. Gerbi, M. Gisclon
from the Mathematics Department of the Universty of Savoie, Chambery, France.

List of Speakers:

F. Bouchut (Ecole Normale Superieure, Paris),
P. Fabrie (University of Bordeaux),
I. Faille (Insitut Francais du Petrole, Paris),
T. Gallouet (University of Marseille),
E. Grenier (Ecole Normale Superieure de Lyon),
O.Gues (University of Nice),
R. Herbin (University of Marseille),
F. Nataf (Ecole Polytechnique, Paris),
D. Serre (Ecole Normale Superieure de Lyon),
J.P.Vila (Institut des Sciences Appliquees de Toulouse).

The program of these Days will include 10 invited 1 hour lectures.
Posters presentation by young reseachers is also organised.

There are NO registration fees but the inscription is recommended.

For more information send e-mail to:
or consult the Web site:

Kindly inform your colleagues interested in this field.

Stephane Gerbi
Laboratoire de Mathematiques
Universite de Savoie
Campus Scientifique
75376 Le Bourget du Lac Cedex
tel: +33 (0) 4 79 75 87 27
fax: +33 (0) 4 79 75 81 42


From: Simon Chandler-Wilde <>
Date: Fri, 18 Feb 2000 23:56:59 -0000
Subject: Workshop on Computational Acoustics and Electromagnetics

An Interdisciplinary Workshop on
14-15 September 2000, University of Bath, UK

Forming part of the Second International Conference on
Boundary Integral Methods: Theory and Applications,
12-16 September 2000

held under the auspices of the Institute of Mathematics
and its Applications and Co-Sponsored by the
Institution of Electrical Engineers and the Institute of Acoustics.


This interdisciplinary workshop, supported by the UK Engineering
and Physical Sciences Research Council,
aims to connect engineers, mathematicians, and industrialists
working on innovative boundary element algorithms and their
application in computational acoustics and electromagnetics,
and to explore the opportunities for increased collaboration in this field.

The following have agreed to give invited talks at the workshop:
W.C. Chew (University of Illinois), A. Kirsch (University of Karlsruhe),
S.P. Walker (Imperial College, London), L.Wrobel (Brunel University, London).

Workshop Organising and Scientific Committee:
Simon Chandler-Wilde (Brunel University, Chair), Ke Chen (University of
Liverpool), Penny Davies (University of Strathclyde), Ivan Graham (University
of Bath), Mireille Levy (Rutherford Appleton Laboratory).

The larger conference includes additional invited talks from:

C. Constanda (Strathclyde University), W. Hackbusch (Max-Planck-Institute,
Leipzig), T. Hou (California Institute of Technology), A. Rathsfeld
(Weierstrass Institute, Berlin), S.A. Sauter (University of Zurich),
I.H. Sloan (University of New South Wales, Sydney), W.L. Wendland
(University of Stuttgart).

Future announcements about the workshop and further details/
application forms can be found at the web site

Future details concerning the larger conference can be found at:


From: Gene Golub <golub@sccm.Stanford.EDU>
Date: Tue, 15 Feb 100 13:55:30 PST
Subject: Forsythe Fellowship at Stanford

Scientific Computing and Computational Mathematics Program
2 Year Acting Assistant Professorship

The School of Engineering at Stanford University has established
Forsythe Fellows in the Scientific Computing and Computational
Mathematics (SCCM) Program. These positions have a duration of two
years and are jointly funded through the School of Engineering and the
research support of faculty members in the SCCM Program. An optional
third year of funding may be available, contingent upon financial

The SCCM program is an interdisciplinary graduate program awarding MS
and PhD degrees. The successful candidate will be expected to carry
out a program of research in scientific computation. The current areas
of research of members of the core faculty include solution of pde's,
linear algebra, optimization and fluid mechanics. He/she will teach
two courses within the SCCM program each academic year and will also
help in the administration of the weekly SCCM seminar. Typically one
of these courses may be at an advanced level and in the field of
interest of the successful candidate while the other will be an
undergraduate, or beginning graduate, level course in computational
mathematics or numerical analysis. It is desirable for the candidate
to have some knowledge of parallel computation.

For further information about the SCCM Program see the URL:

The starting date for this Fellowship is Sept 15, 2000 or possibly earlier.

Please mail applications to:

Forsythe Fellow Committee
Computer Science Dept.
Gates 2B
Stanford University
Stanford CA 94305
Attn. G Golub

The applicants should have three references sent to the Committee.
Stanford University is an equal opportunity employer.
Closing Date: April 1, 2000.


From: Charles Romine <>
Date: Sun, 13 Feb 2000 14:51:09 -0500
Subject: Householder Fellowship at Oak Ridge National Laboratory

Fellowship in Scientific Computing

The Computer Science & Mathematics (CSM) Division at the Oak Ridge
National Laboratory (ORNL) invites outstanding candidates to apply
for the 2000 Alston S. Householder Fellowship in Scientific Computing.
The Fellowship honors Dr. Alston S. Householder, founding director
of the Mathematics Division (now CSM) at ORNL and recognizes his
seminal research contributions to the fields of numerical analysis
and scientific computing.

Funding for the Householder Fellowship comes from the Computational
Mathematics Project, which is supported by the Office of Mathematical,
Information, and Computational Sciences of the U.S. Department of Energy.

Each Householder Fellowship is a one-year appointment (potentially
renewable for a second year) that provides access to state-of-the-art
computational facilities (high-performance workstations and parallel
architectures), and collaborative research opportunities in active
research programs in advanced scientific computing and computational
sciences. The Fellowship offers a competitive salary, fringe benefits,
and travel opportunities.

The purpose of the Householder Fellowship is to promote innovative
research in scientific computing on advanced computer architectures and
to facilitate technology transfer from the laboratory research
environment to industry and academia through advanced training of new
computational scientists.

For more information, visit our web site at

Position requires a Ph.D. in computer science, mathematics, or statistics
(completed no more than three years prior to the appointment). Must have
a strong background and research interest in large-scale scientific
computing. To apply, send resume, statement of research, and three
letters of recommendation by March 31, 2000, to: Householder Fellowship,
Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, TN 37831-6359.

Finalists for the Fellowship will be invited to visit ORNL to present a
seminar and visit the area. The selection committee's final decision will
be made in May 2000. The selected Fellow must be available to begin the
appointment during CY 2000.

ORNL, a multiprogram research facility managed by Lockheed Martin Energy
Research Corp. for the U.S. Department of Energy, is an equal opportunity
employer committed to building and maintaining a diverse work force.


From: Benedict Tanyi <>
Date: Tue, 15 Feb 2000 10:25:38 +0000
Subject: Position at the Fujitsu European Centre, UK


Financial/Mathematical Modeller:

FECIT (a subsidiary of Fujitsu Ltd., Japan) is a
multidisciplinary research centre devoted to the development
of information technology on the latest high performance
parallel computers (visit for more
information on FECIT's research activities).

A vacancy exits within the Financial Engineering
Research Group at FECIT for a researcher in
Financial Risk Management/Modelling. Our scope is
fairly broad, but experience in the following areas
are most desirable:

1. Portfolio Optimization Theory
2. Derivative Pricing (mathematics of)
3. Artificial Intelligence, particularly Neural Networks,
Genetic Programming, etc., with applications in Finance.
4. Good Object-Oriented Programming skills are also desirable.

Practical experience in one or more of the above areas and
clear mathematical analytical skills are preferred, but not
necessary. Applications are encouraged from recent Ph.D.s or
suitable qualified graduates.

It is also essential that the applicants have:
- the ability to do independent research;
- the ability and desire to work in teams of individuals
with diverse backgrounds;
- enthusiasm for working on applications;
- good communication skills.

The salary for this position will be in the range of 22,000 to 30,000 pounds,
depending on qualification. The benefits package includes 25 days annual
holiday, a pension scheme, private health insurance, death in services benefits,
subsidized Health Club fees and subsidized canteen facilities.

Interested applicants should send their CV's to Mrs. Edna Davis, at the
address below. Her E-mail is Further details about the
position can be obtained from Dr. Benedict Tanyi,

Mrs Edna Davis
Fujitsu European Centre for Information Technology
2 Longwalk Road
Stockley Park
Middlesex UB11 1AB
United Kingdom


From: Jean-Davis Benamou <>
Date: Thu, 17 Feb 2000 17:30:26 +0100
Subject: Postdoctoral Position at INRIA, France

We invite applications either at the post-doc level or higher to join our
team at INRIA-Rocquencourt, France (near Paris) to work on the development
of a scientific software library dedicated to Geometrical Optics and
related Hamilton-Jacobi problems. The person chosen for this post will
participate in the object oriented architecture design of the library
and will be in charge of coding the classes.

He will also need to interact with people interested in using or contributing
to the library. The position is for a period of one to two years (possibly
longer) starting as soon as possible. The salary will depend on experience.

Candidates should have some background in numerical methods for PDE's and
ODE's and experience in coding such methods. The candidate should either
have knowledge of OO programming and C++ or be willing to train himself
in such coding techniques. The possibilty for the candidate to split his
time between the development of the library and his own research provided
it is related to G.O. can be discussed.

Interested people should contact


From: Loyce Adams <>
Date: Wed, 16 Feb 2000 10:11:15 -0800 (PST)
Subject: Postdoctoral Position at University of Washington

VIGRE Postdoctoral Position
University of Washington
Department of Applied Mathematics

Applications are invited for one Postdoctoral Fellow starting
September 16, 2000. Applicants are required to have earned a
Ph.D. degree between March 15, 1999 and the start date. This
position is part of an NSF VIGRE program in Applied and
Computational Mathematical Sciences (ACMS) joint between the
Applied Mathematics, Mathematics, and Statistics departments.
The current opening is in the Department of Applied Mathematics,
but the recipient will have opportunities for research collaboration
within these mathematical sciences units and other applied departments
on campus. More details may be found at


From: Hans Schneider <>
Date: Sun, 13 Feb 2000 21:41:09 -0600 (CST)
Subject: Contents, Linear Algebra and Its Applications

Journal: Linear Algebra and Its Applications
ISSN : 0024-3795
Volume : 306
Issue : 1-3
Date : 22-Feb-2000

pp 1-13
Equivalence constants for matrix norms: a problem of Goldberg
A Tonge

pp 15-24
Classes of Schur D-stable matrices
R Fleming

pp 25-31
A representation theorem for algebras with commuting involutions
M Cabrera

pp 33-44
Sparsity of orthogonal matrices with restrictions
GS Cheon, BL Shader

pp 45-57
On sums of three square-zero matrices
K Takahashi, PEIYUAN Wu

pp 59-86
The Nevanlinna-Pick interpolation problems and power moment problems
for matrix-valued functions III: The infinitely many data case

pp 87-102
Generalized totally positive matrices
M Fiedler

pp 103-121
On almost regular tournament matrices
C Eschenbach, JR Weaver

pp 123-130
Dual graphs and knot invariants
M Lien, W Watkins

pp 131-143
Decomposing a matrix into circulant and diagonal factors
M Schmid, R Steinwandt

pp 145-154
Centrogonal matrices
O Krafft

pp 155-163
Spectral clustering properties of block multilevel Hankel matrices
D Fasino

pp 165-182
On condensed forms for partially commuting matrices
YUA Alpin, L Elsner

pp 183-188
A tree whose complement is not eigensharp
VL Watts

pp 189-202
Growth in Gaussian elimination for weighing matrices, W(n,n-1)
CHRISTOS Koukouvinos

pp 203-209
Semidefiniteness without real symmetry
CR Johnson, RB Reams

pp 211-211
Author index


From: Iain Duff <>
Date: Fri, 18 Feb 2000 09:57:34 GMT
Subject: Contents, IMA Journal of Numerical Analysis


Contents of Volume 20, Number 1.

E.P. Stephan and T. Tran
Domain decomposition algorithms for indefinite weakly singular integral
equations: the $h$ and $p$ versions.
pp 1-24

M. Reifenberg and J-P. Berrut
Numerical solution of boundary integral equations by means of
attenuation factors.
pp 25-46

E.F. Toro and S.J. Billett
Centred TVD schemes for hyperbolic conservation laws.
pp 47-79

Z. Chen, T. Shih and X. Yue
Numerical methods for Stefan problems with prescribed convection
and nonlinear flux.
pp 81-98

P.W. Hemker, G.I. Shishkin and L.P. Shishkina
$\epsilon$-uniform schemes with high-order time-accuracy for
parabolic singular perturbation problems.
pp 99-121

Wang Xinghua.
Convergence of Newton's method and uniqueness of the solution of
equations in Banach space.
pp 123-134

J. Schropp
One-step and multistep procedures for constrained minimization problems.
pp 135-152

Huang Chengming
Dissipativity of Runge-Kutta methods for dynamical systems with delays.
pp 153-166


End of NA Digest