NA Digest Saturday, November 30, 2002 Volume 02 : Issue 47

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

Submissions for NA Digest:

Mail to

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


From: G. W. Stewart <>
Date: Mon, 25 Nov 2002 12:00:30 -0500 (EST)
Subject: Allocatable Arrays in Fortran 95

With the help of my student Roger Lee, I am preparing a package called
Matwrap. It is a Fortran 95 wrapper that links matrix operations and
decompositions to Lapack and the Blas. A preliminary version for real
matrices is near completion. But a final design decision has to be
made before it can be released, and I would like some input.

At the heart of Matwrap is a defined type for a matrix represented
within a rectangular array. The type has a component for the array
containing the matrix. Now by strict Fortran 95 conventions that
array must be a pointer array, which can be allocated and deallocated.
The difficulty with this is that the array must be explicitly
deallocated when it is no longer needed to avoid memory leaks, an
awkward way of doing things at best. There are also problems with
deallocating intermediate values generated by overloaded and new

Now there is an extension of Fortran 95 that allows allocatable arrays
in function results, dummy arguments, and defined types. (For more
see Chapter 13 of Metcalf and Reid's Fortran 90/95 Explained.) Since
such arrays are automatically deallocated when they are no longer
needed, this extension provides a solution to the problems of the last
paragraph. However, if the extension is not widely implemented, it is
not an option.

I would like feedback from Fortran 95 users about whether the
compilers they use support the extension and, just as important, if
the extension works. (The f95 compiler on my Sun workstation supports
the extension, but it does not work properly for returned function
values.) Please be specific about the compiler--source and version.
Thank you for your help.

Pete Stewart


From: Gravvanis George <>
Date: Mon, 25 Nov 2002 20:59:13 +0200
Subject: Mini-symposium in Athens on Preconditioning

Call for Papers
Mini-Symposium Organizer: Dr. G. A. Gravvanis

Within the framework of the Conference
The Sixth Hellenic European Research on Computer Mathematics and its
September 25-27, 2003, Athens University of Economics and Business,
Athens, Greece

Topics of the session: Topics of interest include, but no limited to, the
following: Approximate Inverses, Finite Difference and Finite Element
method, Finite Volume method, Boundary element method, Domain
Decomposition methods, Multigrid Methods, Parallel iterative methods,
Numerical Software Description and Environments, Advanced Numerical
Algorithms and Scientific Computing, Applications in Science and

Prospective authors are invited to submit an extended two-page abstract by
20th February 2003 on topics related to the theme of the Session (PDF,
Postscript, or Microsoft Word documents) to the Session Organizer by
email only (

Recommendations concerning the format of the accepted paper to be published
in the HERCMA Conference Proceedings will be sent to authors. Papers must
not have been previously published or currently submitted for publication
elsewhere. Soon after the Conference, extended forms of selected papers will
be considered for the HERMIS-=EC=F0 special issue publication and other
journals publication.

February 20, 2003: Deadline for the submission of the extended abstract
April 10, 2003: Notification of acceptance
May 30, 2003: Deadline for the submission of the complete paper (camera ready)
September 25-27, 2003: Sixth HERCMA Conference

Dr. George A. Gravvanis
Department of Information and Communication Systems Engineering
University of the Aegean
GR-83200 Karlovasi, Samos, GREECE

For more information about the conference a link has been created at


From: Hoang Duc Minh <>
Date: Tue, 26 Nov 2002 18:23:25 +0100 (CET)
Subject: Conference in Hanoi on Scientific Computing

Modelling, Simulation and Optimization of Complex Processes

March 10-14, 2003
Hanoi, Vietnam

- mathematical modelling
- numerical simulation
- methods for optimization and control
- parallel computing: architectures, algorithms, tools, environment
- symbolic computing
- software development
- applications of scientific computing in
physics, mechanics, chemistry, and biology
environmental and hydrology problems
transport, logistics and site location
communication networks, production scheduling
industrial and commercial problems

Uri Ascher (Vancouver), Roland Bulirsch (Munich),
Zhiming Chen (Beijing), Felix L. Chernousko (Moscow),
Peter Deuflhard (Berlin), Andreas Griewank (Dresden),
Ulrich Langer (Linz), Nguyen Van Lien (Hanoi),
Masayasu Mimura (Hiroshima), Bijan Mohammadi (Montpellier),
Michael Robert Osborne (Canberra), Mauricio G.C. Resende (Florham Park)

The conference will take place at the Institute of Mathematics,
18 Hoang Quoc Viet Street, Cau Giay District, Hanoi, Vietnam

Dr. Phan Thanh An
Institute of Mathematics, P.O. Box 631 - Bo Ho, Hanoi, Vietnam
Phone: +84-4-7563474 (ext.: 212)
Fax: +84-4-7564303

Deadline for registration and submission of abstracts: December 10, 2002
Notification of acceptance: January 10, 2003
Deadline for submission of full papers for the conference proceedings
which will published by Springer: May 15, 2003

Detailed information and online registration can be found at the


From: Slimane Adjerid <>
Date: Wed, 27 Nov 2002 14:20:16 -0500 (EST)
Subject: Conference at RPI Honoring Joe Flaherty


Conference on Adaptive Methods for Partial Differential
Equations and Large-scale Computations
In Honor of the 60th Birthday of Joseph Flaherty
Rensselaer Polytechnic Institute Troy, New York
October 11 and 12, 2003

Conference Topics:

-Efficient adaptive refinement strategies
-Data structures for parallel and adaptive computations
-A posteriori error estimation
-Discontinuous Galerkin (DG) methods
-Grid and mesh generation
-Multiscale applications

Invited Speakers:

-Ivo Babuksa (UT-Austin)
-Randy Bank (UCSD)
-Marsha Berger (NYU)
-Martin Berzins (University of Leeds)
-Rupak Biswas (NASA Ames)
-Bernardo Cockburn (University of Minnesota, Twin Cities)
-Karen Devine (Sandia)
-Herbert Keller (UCSD,UCLA)
-William Mitchell (NIST)
-Tinsley Oden (UT-Austin)
-Robert O'Malley (University of Washington)
-Mark Shephard (Rensselaer)
-Jan Verwer (Netherlands)
-Mary Wheeler (UT-Austin)


Slimane Adjerid (Virginia Tech)
Peter Moore (Southern Methodist University)
James Teresco (Williams College)

Contact Information:


From: Lothar Reichel <>
Date: Fri, 29 Nov 2002 16:42:10 -0500 (EST)
Subject: Conference at Kent State for ETNA

On the occasion of the 10th birthday of the Electronic Journal on
Numerical Analysis (ETNA), the conference

"ETNA: Following the flows of Numerical Analysis"

will be held in Kent, OH, May 29-31, 2003. The meeting, in the tradition of
past Kent numerical analysis conferences, will give participants an opportunity
to present their work and to interact in a friendly and relaxed atmosphere.

A web site for the conference, including a list of participants, is under
construction. It can be found at

Please send e-mail to or if you
are interested in participating in the conference.

Daniela Calvetti and Lothar Reichel

on behalf of the organizing committee.

ETNA is available on the web at as well as on


From: Iain Duff <>
Date: Thu, 28 Nov 2002 09:13:32 GMT
Subject: Numerical Algorithms Topic at EuroPar 2003

Submission Deadline: February 9, 2003

EURO-PAR 2003 Topic 11
Numerical Algorithms and
Scientific Engineering Problems

Klagenfurt, Austria
August 26-29, 2003


Fast and robust parallel algorithms for the basic problems of numerical
mathematics are crucial for solving the problems in computational science and
engineering that arise today.

This workshop will be a forum for the presentation and discussion of new
developments in the field of parallel numerical algorithms, covering all
aspects from basic algorithms, software design and prototyping to efficient
implementation on modern parallel architectures and performance analysis.
Because of its importance for High Performance Scientific Computing
applications, the parallel direct or iterative solution of large systems
of linear and nonlinear equations will be of prime interest.

However, contributions dealing with problems other than numerical linear and
nonlinear algebra, or with general topics concerning parallel numerical
methods are also welcome.

* numerical linear algebra
* large sparse or dense linear systems and eigensystems
* nonlinear systems
* fast transforms (wavelets, FFT)
* discretized partial differential equations


November 8, 2002 First call for papers
February 9, 2003 Final date for submissions
May 1, 2003 Notification of acceptance
July 15, 2003 Early registration deadline


Global Chair

Prof. Iain Duff
Computational Science and Engineering Department
Rutherford Appleton Laboratory
Oxfordshire, UK

Vice Chairs

Prof. Henk van der Vorst
Mathematical Institute
Utrecht University, The Netherlands

Dr. Luc Giraud
Toulouse, France

Local Chair

Prof. Peter Zinterhof
Department of Scientific Computing
Salzburg University, Austria


Authors are requested to submit their papers by using our submission site (please indicate TOPIC 11).
The site will open not later than at the beginning of January 2003.
All accepted papers will be available at the conference in the proceedings
published by Springer-Verlag in the LNCS series.


Euro-Par is the annual series of international conferences dedicated to the
promotion and advancement of all aspects of parallel and distributed computing.
Euro-Par 2003 will represent 19 topics in the categories of hardware, software,
algorithms and applications. The conference normally attracts more than 300
participants coming from universities, research centres and industry.


From: Iain Duff <>
Date: Thu, 28 Nov 2002 15:18:26 GMT
Subject: Alan Curtis Day in Oxford

Alan Curtis Meeting in Oxford. December 5 2002

This meeting is to commemorate the 80th birthday of Alan Curtis (formerly of
AERE Harwell). It will be held on Thursday 5 December in the Lecture Theatre
of Oxford Computing Laboratory, Wolfson Building, Parks Road, Oxford.

If you are unable to attend but want me to pass your greetings to Alan, you are
welcome to email these to me at

Iain Duff (RAL)
Opening remarks

Mike Powell (Cambridge)
Some developments of work with Alan on cubic splines

Kevin Burrage (Queensland)
Stochastic models and simulations for chemically reacting systems

John Reid (RAL)
Sparse matrix research at Harwell and the Rutherford Appleton

Ian Jones (AEA PLC)
Computational fluid dynamics and the role of stiff solvers

Lawrence Daniels (Hyprotech UK Ltd)
Current work with Alan on ODE solvers for HSL

At some point in the day Alan himself will give a short address.


From: Brynjulf Owren <>
Date: Fri, 29 Nov 2002 11:10:12 +0100
Subject: New Talent Invited Speaker at SciCADE'03

SciCADE'03, June 30 - July 4, 2003, Trondheim, Norway.

Call for papers.

We have reserved a slot in the format of the plenary talks
for a
New talent invited speaker

All local expenses will be covered for this speaker, and the registration
fee will be waived as well. Those wishing to be considered for the new
talent award are invited to submit a manuscript to the conference organisers
by April 1, 2003.
Please e-mail manuscripts for the new talent session to

The author of the best paper will be asked to give an invited lecture at
Scicade '03. The criteria for eligibility are

1. candidates must be under 35 as of 1 January 2003; and
2. they should not have completed a PhD more than 4 years ago, that is
before January 1st, 1999.

Jointly authored papers may be submitted by an individual candidate.
In such a case the entry must primarily be the original work of the
candidate, and must be accompanied by a statement from the co-authors
agreeing to the submission and detailing the contribution of the
candidate to the paper.

The entries will be considered by an Adjudicating Committee: its
members are those of the Scientific Committee for SciCADE'03, see

For more information see also the conference home page or send an email to the address below.

The webpages are now open for registration and submission of talks.

For SciCADE'03
Brynjulf Owren


From: Ron Boisvert <>
Date: Fri, 29 Nov 2002 14:33:19 -0500
Subject: Postdoctoral Positions at NIST

The National Research Council Associateship Program is accepting
applications for awards for postdoctoral research to be conducted in
residence at the National Institute of Standards and Technology in
Gaithersburg, Maryland and Boulder, Colorado. Among the topic areas
for such research are applied mechanics, computational biology, materials
science, modeling of communications systems, image analysis, applied
optimization, Monte Carlo methods, data mining, mathematical software,
and quantum computing. For details see
The deadline for applications is February 1. The NRC program at NIST is
restricted to US citizens.


From: Heinz W. Engl <>
Date: Sun, 24 Nov 2002 10:24:21 +0100 (CET)
Subject: Postdoctoral Positions at RICAM, Austria


RICAM is a research institute going into operation on January 1, 2003,
and will be gradually built up to a total of 25 postdoc positions in five

Scientific Computing - Computational Methods for Direct Field Problems
(Prof.Ulrich Langer,

Scientific Computing - Inverse Problems (Prof.Heinz Engl,

Symbolic Computation (Prof.Bruno Buchberger,

Analysis of Partial Differential Equations (Prof.Peter Markowich,

Financial Mathematics (Prof.Gerhard Larcher,;
Prof.Walter Schachermayer,

In addition to ongoing research, which will be partly interdisciplinary
between the areas mentioned, there will be an international visitors'
program in the form of special semesters on application fields.

The institute will be housed on the campus of the Johannes Kepler
Universitaet in Linz, a town of about 200.000 on the Danube, very close to
the Austrian Alps, and half-way between Vienna and Salzburg.

We are looking for PostDocs with a strong interest in one of the fields
above, who are also willing to work in an interdiscplinary environment;
doctorate in mathematics or a closely related field required. In addition,
we are also looking for Senior Postdocs, who could lead a group and would
also be expected to attract additional funding; they should have a
will be English. The positions are initially for up to three years, one
substantial publication record beyond the doctorate. The working language
renewal for three more years is possible depending on achievement.

Please send applications with personal and scientific data, copies of
relevant documents and a statement about scientific interests and
achievements to

Prof. Heinz W. Engl
Director, RICAM
Austrian Academy of Sciences
c/o Kepler University
A-4040 Linz

and a copy to the area leader in whose group you want to work. There is no
deadline, the announcement remains open until all positions are filled; we
expect to be able to fill the first position on March 1, 2003.


From: Jan Hesthaven <>
Date: Mon, 25 Nov 2002 14:25:16 -0500 (EST)
Subject: Postdoctoral Position at Brown University

Postdoctoral Position in Numerical PDE's at Brown University.

A position as a postdoctoral researcher in the Division of Applied
Mathematics at Brown University will be open from summer/fall 2003.

The position is part of an ongoing effort to develop high-
order adaptive methods for the solution of time-dependent
conservation laws on general unstructured grids. The
research involves areas such as space-time error-estimation,
advanced time-integration schemes, spatial hybrid scheme, and
practical strategies for adaptive solutions of time-
dependent problems on parallel computers.

The position is available for two years with the possibility
of extension pending sufficient external funding. The
candidate is expected to have or receive shortly a PhD in
computational science or applied mathematics with a strong
background in numerical solution of partial differential
equations. The ideal applicant will have knowledge and
experience with error-estimation and adaptivity as well as
high performance computing. Thorough experience with one or
more programming languages (C, C++, FORTRAN) is expected.
A background in electromagnetics, gas-dynamics, or plasma-
physics is desired but required.

Interested applicants can get further information about
the position by contacting

Jan S. Hesthaven
Division of Applied Mathematics
Brown University, Box F
Providence, RI 02912
(401) 863-2671

Electronic applications (please use postscript, PDF, or
Word formats) including curriculum vitae, list of publications,
a statement of research interest, and names of 3 references
with e-mail, should be sent to the email given above before
January 15, 2003. However, the search will be ongoing until
the position is filled.

Brown University is an equal opportunity employer and
encourages women and minorities to apply.


From: Vladik Kreinovich <>
Date: Sat, 23 Nov 2002 15:50:37 -0700 (MST)
Subject: Contents, Reliable Computing

Reliable Computing
Volume 9, issue 1, 2003

Mathematical Research

Diagrammatic Analysis of Interval Linear Equations:
Part I: Basic Notions and the One-Dimensional Case
Zenon Kulpa

Controlling the Wrapping Effect in the Solution
of ODEs for Asteroids
Jens Hoefkens, Martin Berz, Kyoko Makino


Taylor Forms - Use and Limits
Arnold Neumaier


COCOS'02: A Workshop on Global Constrained Optimization and Constraint
R. Baker Kearfott


From: M2AN <>
Date: Tue, 26 Nov 2002 11:12:49 +0100
Subject: Contents, Mathematical Modelling and Numerical Analysis

Mathematical Modelling and Numerical Analysis
Vol. 36 No. 5 (September/October 2002)

Olivier Pironneau

*Mathematical and Computational Framework for Reliable Real-Time
Solution of Parametrized Partial Differential Equations
Christophe Prud'homme, Dimitrios V. Rovas, Karen Veroy and Anthony T. Patera

*FER/SubDomain: An Integrated Environment for Finite Element Analysis
using Object-Oriented Approach
Zhi-Qiang Feng and Jean-Michel Cros

*Automatic Differentiation Platform: Design
Christele Faure

*Object oriented design philosophy for scientific computing
R.B. Devloo and Gustavo C. Longhin

*C++ Tools to construct our user-level language
Frederic Hecht

*Distributed Objects for Parallel Numerical Applications
Francoise Baude, Denis Caromel and David Sagnol

*EasyMSG: Tools and techniques for an adaptive overlapping in SPMD programming
Pascal Have

*GO++: A modular Lagrangian/Eulerian software for Hamilton Jacobi
equations of geometric optics type
Jean-David Benamou and Philippe Hoch

*An object-oriented approach to the design of fluid mechanics software
Christophe Calvin, Olga Cueto and Philippe Emonot

*Some new technics regarding the parallelisation of ZeBuLoN, an
object oriented finite element code for structural mechanics
Frederic Feyel

*Concepts-An Object-Oriented Software Package for Partial
Differential Equations
Philipp Frauenfelder and Christian Lage

*New trends in coupled simulations featuring domain decomposition and
Philippe d'Anfray, Laurence Halpern and Juliette Ryan


End of NA Digest