NA Digest Sunday, November 16, 2003 Volume 03 : Issue 46

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

Submissions for NA Digest:

Mail to

Information via e-mail about NA-NET: Mail to


From: John G. Lewis <>
Date: Fri, 14 Nov 2003 15:03:01 -0800
Subject: Change of Address for John Lewis

After 25 years with Boeing, I've decided to try some new high-speed
challenges by joining the Cascade development team at Cray Inc. My new
addresses are:

Dr. John G. Lewis
Cray Inc.
411 First Avenue South, Suite 600
Seattle WA 98104-2860

Telephone: 206-701-2273
Fax: 206-701-2500
Email: jglewis "at"


From: Ben Fang <>
Date: Thu, 13 Nov 2003 18:10:31 -0800
Subject: Update of Our Multi-Dimension Interpolation

In March, we reported our simple, powerful interpolation method for multi-dimensional
space at . The method can be applied to arbitrary
dimensions (1D, 2D, 3D, ...10D,....).

There are newly added info at the same site, which include new sample cases, FAQ, and
Investor's Forum. In Investor's Forum, technical proposal to government agency (NSF) is
included. We are pursuing to apply our method to spherical and cylindrical coordinates to
facilitate 3-D surface reconstructions. Feedbacks are most welcome. Critical comments
from the computing community are appreciated.

Dr. Ben K.K. Fang


From: Nils Wagner <>
Date: Fri, 14 Nov 2003 10:39:24 +0100
Subject: From NASTRAN to MATLAB

I am looking for freely available MATLAB routines for reading NASTRAN
files (*.o4 files (binary output))
As far as I know, there is only a commercial tool for this purpose
(Structural Dynamics Toolbox).

Any hints would be appreciated.
Thanks in advance,

Nils Wagner


From: Francesca Mazzia <>
Date: Wed, 12 Nov 2003 12:59:18 +0100 (CET)
Subject: New Release of the Test Set for IVP Solvers

Dear colleague,

Hereby we announce the release 2.2 of the Test Set for IVP
Solvers, in which a number of Initial Value Problems from several
application fields has been brought together.

The test set consists of a descriptive part and a software part.
The first part describes test problems and reports on the behavior
of a few state-of-the-art solvers when applied to these problems.
The software serves as a platform to test the performance of a
solver on a particular test problem.

With respect to the release 2.1, we have updated all the solvers,
two more solvers (MEBDFI and GAMD) are supported, and five well
known stiff problems have been added. Moreover ODE, DAE and IDE
problems have been gathered into three different lists and link to
separated pages reporting a short description, work precision
diagrams and plots of the solutions. A short description of each
solver together with its related references, has also been added
in the descriptive part of the Test Set.

It is possible to submit new test problems.

The Test Set can be found at

Comments, remarks and suggestions are welcome.

Francesca Mazzia and Felice Iavernaro, Dipartimento di Matematica,
Universita' di Bari, Italy.


From: Oleg Burdakov <>
Date: Wed, 12 Nov 2003 19:10:14 +0100 (MET)
Subject: Subscriptions to Optimization Methods and Software

Dear Colleague,

We would like to inform you that Taylor & Francis, the new publisher of
the journal "Optimization Methods and Software" (OMS), is now offering
individual electronic subscriptions of OMS at an annual price of 100
Dollars. For details see the leaflet:
This price refers to all individual subscribers, not only to the AMS or
SIAM members as mentioned in the leaflet. Online sample copy is
available on the journal home page:

Best regards,
Oleg Burdakov and Andreas Griewank,
Co-Editors of OMS


From: Bernard Philippe <>
Date: Mon, 10 Nov 2003 09:19:51 +0100 (CET)
Subject: HPC Course in Africa

HPC COURSE IN AFRICA (Some scholarships will be granted)
First semester of the UNESCO chair at LAMSIN

Tunis, 2nd of February - 30th of April 2004

Object of the program

This set of courses will describe the state-of-the-art for programming
procedures for large problems in linear algebra.
Some of the courses are devoted to the definition of currently used
methods. Special attention is given to the race for efficiency which
opposes iterative and direct methods for solving large linear systems.
Other courses focus on the parallelization of the methods and on their
implementation on networks of workstations or on multiprocessors.
By analysing all the steps of the numerical simulation, from modelling to
computer implementation, it is shown, using a particular application -
underground water flow - how to build efficient software which can exploit
local or exterior computing resources.

Patrick Amestoy, ENSEEIHT, Toulouse
Jocelyne Erhel, INRIA, Rennes
Jerome Jaffre, INRIA, Rocquencourt
Mohamed Jemni, University of Tunis
Hugues Leroy, INRIA Rennes
Gerard Meurant, CEA/DIF, Bruyeres le Chatel
Bernard Philippe, ENIT/LAMSIN, Tunis
Jean Roberts, INRIA, Rocquencourt
Ahmed Sameh, Purdue University

Nejla Hariga, ENAT, Tunis
Dany Mezher, ESIB, Beyrouth
Ridha Touihri, IPEIM, Monastir

Organizer : Bernard Philippe, ENIT-LAMSIN, Tunis / INRIA-Rennes


This program is open to all researchers and engineers who would like to
acquire the skills necessary for programming numerical models in a way
that best uses the capacities of the computer. It is especially
well-suited to PhD students.
There is no registration fee for the course. Some fellowships for
reimbursing living and travel expenses will be granted to PhD students
from African countries who are sponsored by their universities.

Course organization

The courses will be given in French and/or in English, depending on the
composition of the audience.
Each course is one-week long. It is followed up by lab exercises which are
organized by members of LAMSIN. Lab sequences will conclude with the
installation and use of up-to-date software.


The necessary background in mathematics is the bachelor's level. In
computer science, the attendees should know how to program in at least one
of the languages: C, FORTRAN, Matlab, or SCILAB..


For a full program, see

For registration and/or for a scholarship application, contact

Bernard Philippe

phone: +216 71 874 700 ext : 555

email :


From: Johanna Hunt <>
Date: Mon, 10 Nov 2003 13:33:27 -0000
Subject: Workshop at Hertfordshire on Automatic Differentiation

Joint University of Hertfordshire/Cranfield University
Automatic Differentiation Workshop
Thursday 27th & Friday 28th November 2003

to be held at:
Room LD454
College Lane Campus
University of Hertfordshire
Herts AL10 9AB
England, Europe

Sponsored by the Engineering and Physical Sciences Research Council

This 2-day workshop is the latest in a series that provides a forum for the
presentation of theoretical developments in and applications of Automatic
Differentiation (AD) and adjoint methods.

Updated information can be found at:
and follow link to "Workshops on Automatic Differentiation"


Christian Bischof (RWTH-AACHEN)
'Sensitivities for a Single Drop Simulation'
George Corliss (Marquette University)
'Applying Global Optimization in Structural Engineering'
Paul Cusdin (Queen's University Belfast)
'A Study in AD for Computational Fluid Dynamics'
Chris Hill (MIT) (TBC)
Julia Sternberg (Dresden University)
'A-Revolve: An adaptive memory- and run-time-reduced procedure
for calculating adjoints; with an application to the instationary
Navier-Stokes system '
Mohamed Tadjouddine (Cranfield University)
'Elimination AD Applied to Jacobian Assembly for an Implicit
Compressible CFD solver.'
Bruce Christianson (University of Hertfordshire)
'Random thoughts on modelling uncertainties'
Andreas Griewank (Dresden University of Technology)
'Jacobians and Hessians are Scarcely Matrices!'
Koichi Kubota (Chuo University / University of Hertfordshire)
'Detection of numerically unstable computations with automatic
Uwe Naumann (Argonne National Laboratory)
'The 'No Free Refill' Conjecture'
John Pryce (Cranfield University)
'Why the method can't be fooled by hidden cancellation in code
defining the DAE'
Jean Utke (Argonne)
'Why is the DAG falling apart'


From: Jose Castillo <>
Date: Tue, 11 Nov 2003 17:12:23 -0800
Subject: Panamerican Advanced Studies Institute

PanAmerican Advance Science Institute (PASI II) and Workshop in
Applied and Computational Mathematics (Panam V)

Tegucigalpa, Honduras, Central America
June 13-26, 2004

This NSF-DOE sponsored event will consist of a number of short courses and
it is directed to graduate students, post-docs and faculty of the Americas.

There will be 40 fellowships available for students attending the event: 20
for students from USA Universities and 20 for students of Universities in
the rest of America. These fellowships will cover travel and local expenses
for two weeks, since the grantees are expected to attend the second week

Please visit the event WEB site
at for
details and pre-requisites.

Lecturers that have agreed to give short courses so far are:

Bertil Gustafson, Uppsala University, Sweden
Teresa Larsen, San Diego State University, USA
Godela Scherer, Universidad Simon Bolivar, Venezuela and University of
Reading, UK
Maria C. Pereyra, University of New Mexico, Albuquerque, USA
Antonio Redondo, Los Alamos National Laboratory, New Mexico, USA
Sergio Rojas, Universidad de Los Andes, Merida, Venezuela
Horst Simon, Lawrence Berkeley Lab.
Details on the second week Workshop can be found in the event WEB site:


From: Willi Schoenauer <>
Date: Wed, 12 Nov 2003 18:06:12 +0100
Subject: Block Lecture at Karlsruhe on Scientific Supercomputing

In the week 16-20 Feb. 2004 there is a block lecture

"Scientific Supercomputing: Architecture and Use of Shared and Distributed Memory Parallel
with exercises on different supercomputers at the computer center of the University of
Karlsruhe in Germany. There is no fee. Details are given in

The contents of the lecture and how to obtain the manuscript (paperback, 309 pages)
can be seen via

Willi Schoenauer
Rechenzentrum der Universitaet Karlsruhe
D-76128 Karlsruhe, Germany


From: Bette Byrne <>
Date: Fri, 14 Nov 2003 10:51:06 +0000
Subject: Conference in Oxford on Numerical Methods for Fluid Dynamics

(in association with ERCOFTAC)


The 8th ICFD International Conference on Numerical Methods for Fluid Dynamics organised
by the ICFD (Institute for Computational Fluid Dynamics - in
association with ERCOFTAC(the European Research Community On Flow, Turbulence And
Combustion - will be held at St Catherine's College, Oxford,
from 29 March - 1 April 2004.

The following have agreed to give invited talks: Remi Abgrall (Bordeaux), Jerry Brackbill
(Los Alamos), Tom Hou (Caltech), Tony Hutton (QinetiQ Ltd), Roland Keunings (CESAME, Louvain),
Mike King (BP), Anthony Patera (MIT), Piotr Smolarkiewicz (NCAR Boulder), John Trangenstein
(Duke), Alessandro Veneziani (Milan), Andy Wathen (Oxford). In addition to invited lectures
the programme will include contributed talks of twenty minutes and poster sessions.


In association with ERCOFTAC (the European Research Community On Flow, Turbulence And
Combustion - we are able to offer a number of scholarships for
attendance at the ICFD Conference by students from the European Community. These
will be awarded on a first come, first served basis. If you are interested in further
information on the scholarships please email Bette Byrne, the Conference Secretary on


Two page abstracts for contributed papers should be submitted by 25 NOVEMBER 2003.
All accepted papers will be refereed for inclusion in the Proceedings which will be published
in a Special Issue of the International Journal for Numerical Methods in Fluids


A feature of the meeting will be the fourth award of "The Bill Morton Prize" for a paper on
CFD by a young research worker. The Prize paper, which should be suitable for a 30 minute
talk, will be chosen from eligible contributed papers submitted for the Conference and
marked "for the Bill Morton Prize". The shortlist of papers will be presented by the authors
at a special session at the Conference. Joint papers must be accompanied by a statement
from collaborators confirming that the author made a substantial contribution.

Candidates should be under the age of 31 on 1 April 2004. Each candidate should include a
statement of his/her age on that date and should indicate that he/she would be available to
present his/her paper at the Conference. The Prize will be presented at the Conference dinner.

Further details can be found at

Abstracts should be addressed to:

Mrs B Byrne
Institute for Computational Fluid Dynamics
Wolfson Building, Parks Road
Oxford OX1 3QD, UK.
Phone: 44-(0) 1865-273883
Fax: 44-(0) 1865-273839


From: Michele Benzi <>
Date: Fri, 14 Nov 2003 11:35:31 -0500 (EST)
Subject: Conference at Emory on Preconditioning



DATE: MAY 19-21, 2005 [right before the Householder Symposium XVI]

PLACE: Emory University, Atlanta, GA

This conference is a sequel to three successful "Preconditioning" meetings,
the first held in Minneapolis in 1999, the second in Tahoe City in 2001, and
the third in Napa in 2003. The theme of the conference is preconditioning
methods, in a broad sense, with some emphasis on industrial problems and

A forthcoming announcement will give details on paper submission,
important dates, program committee members, and invited speakers.

The organizers,

Michele Benzi, Emory University
Esmond Ng, Lawrence Berkeley National Laboratory
Yousef Saad, University of Minnesota
Wei-Pai Tang, The Boeing Company


From: Eilish Hathaway <>
Date: Fri, 14 Nov 2003 11:37:22 -0800
Subject: IPAM Undergraduate Summer Program

Undergraduate Summer Program: We are very excited to announce the fourth
year of an innovative summer program for undergraduates at UCLA's
Institute for Pure and Applied Mathematics. We are recruiting talented
national and international undergraduates and would be grateful for your
assistance in bringing this program to their attention. Our program
called Research in Industrial Projects for Students ("RIPS") creates
teams of 3-5 undergraduates paired with faculty mentors and industry
liasons whose goal is to solve real-world industrial problems. Sponsors
included companies such as Pixar and national laboratories such as Los
Alamos. Last summer we completed our 3rd program and as in the previous
two years, it was very successful and the undergraduates who
participated were enthusiastic about their experience. More information,
including comments from students and information on prior projects can
be found on our website at

Questions can be addressed to

You can view and/or print out a poster for the program (and we'd be
grateful if you could post it in an area visible to undergraduates) at
this url


From: Jorge More' <>
Date: Thu, 13 Nov 2003 12:18:40 -0600
Subject: Wilkinson Fellowship in Scientific Computing


The Mathematics and Computer Science Division of Argonne National
Laboratory invites outstanding candidates to apply for the Wilkinson
Fellowship in Scientific Computing. The appointment is for one year
and may be renewed for a second year.

This fellowship was created in memory of Dr. James Hardy Wilkinson,
F.R.S., who had a close association with the Mathematics and Computer
Science Division as a consultant and guiding spirit for the EISPACK
and LINPACK projects. The Wilkinson Fellowship is intended to
encourage scientists actively engaged in state-of-the-art research in
scientific computing. Candidates must have received a recent Ph.D.
prior to the beginning of the appointment. The benefits of the
appointment include a highly competitive salary, moving expenses, and
a generous professional travel allowance.

The appointment will be in the Mathematics and Computer Science
Division, which has strong programs in scientific computing, software
tools, and computational mathematics. Of special interest are
algorithms and software for linear algebra, optimization, differential
equations, computational differentiation, and unstructured mesh
computations; software tools for parallel computing; and numerical
methods for computational science problems. For details about
current projects, see

Internationally recognized for innovative research in high-performance
computing, the Mathematics and Computer Science Division supports an
excellent computational environment that includes access to a
scientific visualization laboratory and a distributed computing
laboratory. For further information, see

Argonne is located in the southwestern Chicago suburbs, offering the
advantages of affordable housing, good schools, and easy access to the
cultural attractions of the city.

Resumes should be sent by e-mail to, referencing
position mcs-305721. Please include a statement of research goals and
the names and addresses of three references. The closing date for
applications is January 12, 2004. The applications will be reviewed by
a selection committee and a candidate announced in April 2004.

Argonne is an affirmative action/equal opportunity employer.


From: Kees Vuik <>
Date: Mon, 10 Nov 2003 15:40:24 +0100
Subject: Postdoctoral Position at TU Delft

Parallel solvers for the 2D and 3D acoustic and electromagnetic wave equations

Physical Chemistry and Molecular Thermodynamics/Numerical Mathematics,
Julianalaan 136, Delft 2628 BL, The Netherlands

Function description:
Postdoctoral Position: Rigorous modeling of 3D wave propagation in
inhomogeneous media for geophysical and optical problems

The research in this project aims at finding efficient parallel solvers for the
2D and 3D acoustic and electromagnetic wave equations. Both iterative and
direct solvers are studied. Within the numerical mathematics group new
preconditioners are being developed to reduce the number of iterations required
for an iterative solution of the 3D Helmholtz equation. The time needed to
perform a single iteration is still very large however. The task of the
postdoctoral fellow involves the parallellisation of the iterative solution of
the 3D Helmholtz equation using domain decomposition techniques. The proposed
computer architecture for parallellisation is a Linux cluster. The work will
be carried out in close collaboration with researchers at Shell/SIEP in
Rijswijk and Philips in Eindhoven.

The research is part of a joint project of the sections Physical Chemistry and
Molecular Thermodynamics ( of the Department of
Chemical Engineering (Prof. Dr. S.W. de Leeuw), Numerical Mathematics
( at the Department of Applied
Mathematics (Dr. C. Vuik), Shell/SIEP in Rijswijk (Dr. R. Plessix) and Philips
Natlab in Eindhoven (Prof. H.P. Urbach).

We are looking for an enthusiastic and motivated candidate, who recently
obtained (or soon will obtain) a PhD. Degree in Numerical Mathematics,
Computational Science, Computational Physics or related discipline. He/she has
an interest in Computational Modelling and Seismics. Good communicative skills
and fluency in English are essential.

Working conditions:
Enrollment and salary are according to the Employment Agreement Dutch
Universities. Enrollment will start as soon as possible. The postdoctoral fellow
will be appointed for a period of two years.

How to apply:
Letters of application (by email or ordinary mail) including a CV to:
Mrs. D.T. Landstra
Julianalaan 136
2628 BL Delft
The Netherlands

Additional information:
Prof. Dr. S.W. de Leeuw Dr. C. Vuik
phone: +31 (0)15-2785037 +31 (0)15-2785530


From: Chris Tisdell <>
Date: Wed, 12 Nov 2003 13:28:26 +1100
Subject: PhD Scholarship at UNSW, Sydney in


The School of Mathematics, UNSW is pleased to announce a PhD scholarship for
work on the project "Dynamic Equations on Time Scales".

Dynamic equations on time scales (measure chains) is a new and revolutionary
theory that models phenomena whose behaviour can change smoothly at one
level, but irregularly at another. This project will investigate these types
of equations and look at ideas such as existence, uniqueness and behaviour of
solutions. The theory is expected to apply to such "stop-start" behaviour as:
insect population models, stock markets, crop harvests and combustion.

The project would be ideal for students who have a background and interest in
differential equations, difference equations and nonlinear analysis.

The project is funded through the Australian Research Council's Discovery
Grants and funding is available for up to three years.

Scholarships carry a stipend of $23,294 AUD per annum for a maximum of three
years. (It does not cover tuition fees, currently $23,040 AUD per year. Fee
waivers will be available for PhD scholars who are Australian citizens or
Australian residents.)

Applicants should be highly motivated, and have a clear idea of their goals
and ambitions. They should show a commitment to fundamental research in

Successful candidates will work under the supervision of Dr Chris Tisdell.

To apply, please send one hard copy of each of (i) a full curriculum vitae,
(ii) a covering letter and (iii) the names and contact details of two
referees (e-mail addresses are essential) by normal mail to:

Dr Chris Tisdell
School of Mathematics
UNSW Sydney 2052

(Do not ask referees to send their references directly to me. I will contact

The deadline for applications is Wednesday, 24th December 2003.

Further information can be obtained from:

Dr Chris Tisdell
School of Mathematics
UNSW Sydney 2052


From: Science Direct <>
Date: Wed, 12 Nov 2003 06:39:11 -0500 (EST)
Subject: Contents, Journal of Approximation Theory

Journal of Approximation Theory Volume 124, Issue 2, Pages 139-282 (October 2003)

On the values of continued fractions: q-series, Pages 139-153
Tapani Matala-aho

Convergence and asymptotic behavior of Jackson integrals associated with irreducible
reduced root systems, Pages 154-180
Masahiko Ito

Sequence spaces of spline functions on subsets and l[infin]-spaces, Pages 181-193
Peter Wingren

Finding the projection of a point onto the intersection of convex sets via projections
onto half-spaces, Pages 194-218
Lev M. Bregman, Yair Censor, Simeon Reich and Yael Zepkowitz-Malachi

Strong converse inequalities for Baskakov operators, Pages 219-231
Shunsheng Guo and Qiulan Qi

Approximation of the rth differential operator by means of linear shape preserving
operators of finite rank, Pages 232-241
S. P. Sidorov

Multi-node higher order expansions of a function, Pages 242-253
Han Xuli

Normal indices in Nikishin systems, Pages 254-262
A. Branquinho, J. Bustamante, A. Foulquie Moreno and G. Lopez Lagomasino

Determining radii of meromorphy via orthogonal polynomials on the unit circle, Pages 263-281
D. Barrios Rolania, G. Lopez Lagomasino and E. B. Saff


From: Communications in Math Sciences <>
Date: Mon, 10 Nov 2003 10:02:16 -0600 (CST)
Subject: Contents, Communications in Mathematical Sciences

Communications in Mathematical Sciences (CMS), Vol 1, No. 4, 2003
* CMS is online at *

Stanley Osher and Richard Tsai
Review Article: Level Set Methods and Their Applications in Image Science
Laurent Boudin, Laurent Desvillettes and Renaud Motte
A modeling of compressible droplets in a fluid
Fanghua Lin and Tai-Chia Lin
Multiple time scale dynamics in coupled Ginzburg-Landau equations
Jingxue Yin, Peidong Lei and Zhouqun Wu
Uniqueness of BV entropy solutions for high dimensional quasilinear
parabolic equations with arbitrary degeneracy
Ioannis G. Kevrekidis, C. William Gear, James M. Hyman, Panagiotis G.
Kevrekidis, Olof Runborg, and Constantinos Theodoropoulos
Equation-Free, Coarse-Grained Multiscale Computation: Enabling
Microscopic Simulators to Perform System-Level Tasks
Li Chen, Ling Hsiao and Yong Li
Strong Solution to a Kind of Cross Diffusion Parabolic System
Philippe G. LeFloch and Mai Duc Thanh
The Riemann Problem for Fluid Flows in a Nozzle with Discontinous
Weizhu Bao, Norbert J. Mauser and H.P. Stimming
Effective One Particle Quantum Dynamics of Electrons: a Numerical
Study of the Schrodinger-Poission-Xalpha Model
Jack Xin and Yingyong Qi
A PDE based two level model of the masking property of human ear


From: Joe Traub <>
Date: Thu, 13 Nov 2003 20:44:49 -0700 (MST)
Subject: Contents, Journal of Complexity

Volume 19, Number 6, December 2003
2003 Best Paper Award Committee
Guest Editors' Preface

Five- and six-dimensional lattice rules generated by structured matrices
Ronald Cools, Hilde Govaert

Worst case complexity of multivariate Feynman-Kac path integration
Marek Kwas, Younming Li

I-binomial scrambling of digital net and sequences
Shu Tezuka, Henri Faure

Where does smoothness count the most for Fredholm equations of the second
kind with noisy information?
Arthur G. Werschulz

Regular Articles

A recursive algorithm for the infinity-norm fixed point problem
Spencer Shellman, K. Sikorski


End of NA Digest