NA Digest Friday, June 27, 2003 Volume 03 : Issue 26

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

Submissions for NA Digest:

Mail to

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


From: Arieh Iserles <>
Date: Mon, 23 Jun 2003 09:43:37 +0100
Subject: Leslie Fox Prize

The 11th Leslie Fox Prize Competition was held at the Isaac Newton
Institute in Cambridge on Friday, 20th June.

The finalists presented their talks in the following (randomly chosen)

"Monotone numerical schemes for nonlinear elliptic partial differential

Melvin LEOK
"Foundations of computational geometric mechanics"

"A posteriori error estimation for finite element discretization of
parameter identification problems"

Tatjana STYKEL
"Balanced truncation model reduction for semidiscretized Stokes equation"

"Adaptive mollifiers for high resolution recovery of piecewise
smooth data from its spectral information"

"A parallel multilevel partition of unity method for elliptic partial
differential equations"

All talks were of an exceedingly high quality, thereby setting the
adjudicating committee, consisting of Arieh Iserles (chair), Alastair
Spence and Endre Suli, an unenviable task. After substantive
deliberation, the adjudicators awarded Leslie Fox Prizes to all
finalists, with the First Prize being awarded to Jared Tanner.

Congratulations to all the finalists.


From: Charles Wampler <>
Date: Tue, 24 Jun 2003 14:35:38 -0400
Subject: A Singular Eigenvalue Problem

The following singular eigenvalue problem arises in the solution of a
problem in geometry. I have never seen this type of problem in standard
linear algebra texts, but perhaps there is some existing theory to handle it.

We have NxN matrices A,B,C,D with known entries. We wish to find the
eigenvalues x and eigenvectors v that satisfy
lim_{q->0} [ A + qB + x(C + qD)] v = 0.
If A and C are general, then this is just a generalized eigenvalue problem
with q=0. In the case at hand, the eigenvalue problem is singular at q=0,
because A and C have a common null space. However, for general q, the
problem is nonsingular, so the limit values must exist.

We could consider the characteristic equation, gathering terms by powers
of q as
det[ A + qB + x(C + qD)] = p_0(x) + q*p_1(x) + ... + q^N*p_N(x).
In the case at hand, p_0(x) is identically zero. The eigenvalues we seek
are the roots of the first p_k(x) that is not trivial. This is surely a bad
way to approach the problem numerically.

I would appreciate any ideas on how to solve this problem in a stable and
efficient manner, by say, converting it to a nonsingular eigenvalue problem.
The matrices A,B,C,D are all size N=40.

Charles Wampler
General Motors R&D Center


From: Nils Wagner <>
Date: Wed, 25 Jun 2003 11:12:53 +0200
Subject: Large-scale Generalized Eigenvalue Problems

I am looking for technical reports illustrating the practical experience
in the numerical solution of large-scale (n>10^5) generalized eigenvalue
problems of the form

A x = \lambda B x

with real matrices A, B, where A is asymmetric and B is SPD. Suitable
solvers for this problem class are of particular interest.

Any hints would be appreciated.
Thanks in advance,

Nils Wagner


From: Leslie Greengard <>
Date: Mon, 23 Jun 2003 21:16:09 -0400 (EDT)
Subject: Fast Multipole Codes

MadMax Optics Inc. is pleased to announce the release of FMMPART3D,
a Fast Multipole-based particle code for 3D Coulombic/gravitational
interactions. This release is a FORTRAN callable library routine
which evaluates all pairwise interactions due to N monopole and/or
dipole sources at N locations in O(N) time. The user can set the
accuracy to 3, 6, 9, or 12 digits.

A second routine, FMMPARTTARG3D, evaluates the potential and its gradient
induced by the sources at an arbitrary set of target locations.

The library and user's manual are available for free trial at
The trial license expires on May 1, 2004.

The 2D FMM Toolbox for MatLab (TM) is also available at this website.


From: Andre Weideman <>
Date: Fri, 27 Jun 2003 08:50:22 +0200 (GMT)
Subject: Revised Version of DMSUITE

An updated version of DMSUITE (A MATLAB Differentiation Matrix
Suite) by Satish Reddy and myself is now available. Larry Shampine
informed us that several of the functions failed to execute
under MATLAB Release 13, and these corrections, along with a few
smaller ones, have now been made.

The suite is freely available for download from the MathWorks Web Site:
please go to
If you find the suite useful, please complete the review on that page.

DMSUITE also has a Web Site (and hopefully soon a guest book) at
This Web Site replaces the old one at, which
will soon cease to exist. Please update bookmarks accordingly.

For those readers unfamiliar with DMSUITE, here is an abstract
from the accompanying paper, which was published in ACM TOMS,
Vol. 26, pp. 465--519 (2000):

Abstract: A software suite consisting of seventeen MATLAB functions
for solving differential equations by the spectral collocation
(a.k.a. pseudospectral) method is presented. It includes functions for
computing derivatives of arbitrary order corresponding to Chebyshev,
Hermite, Laguerre, Fourier, and sinc interpolants. Auxiliary functions
are included for incorporating boundary conditions, performing
interpolation using barycentric formulas, and computing roots of
orthogonal polynomials. It is demonstrated how to use the package
for solving eigenvalue, boundary value, and initial value problems
arising in the fields of special functions, quantum mechanics, nonlinear
waves, and hydrodynamic stability.

JAC (Andre') Weideman
Department of Applied Mathematics
University of Stellenbosch
Stellenbosch 7600
South Africa


From: Aliona Dreglea <>
Date: Tue, 24 Jun 2003 14:25:52 +0100
Subject: Change of Address for D.N. Sidorov

new current address:
Senior Research Fellow
Department of Applied Mathematics
Energy Systems Institute (SB)
Russian Academy of Sciences
130 Lermontov Str, Irkutsk 33 Russia
office phone: +7 3952 428 440
fax: +7 3952 426 796
e-mail :

former address:
Department of Electronic and Electrical
Engineering, Trinity College Dublin,
Dublin 2, Ireland


From: Ya-xiang Yuan <>
Date: Mon, 23 Jun 2003 10:09:15 +0800
Subject: Conference in China Numerical Linear Algebra and Optimization

International Conference on Numerical Linear Algebra and Optimization
(GuiLin, China, Oct 7-10, 2003)

An International Conference on Numerical Linear Algebra and Optimization of
Computational Mathematics and Scientific/Engineering Computing, Academy of
Mathematics and Systems Sciences, Chinese Academy of Sciences, and Guangxi
University. Its sponsors include Chinese Natural Science Foundation, Chinese
Mathematical Society, Chinese Society of Industrial and Applied Mathematics.
Invited lectures on recent advances of numerical optimization and numerical
linear algebra will be given.

INVITED SPEAKERS include(tentative list):
M. Fukushima (Kyoto Univ, Japan)
G. Golub (Stanford Univ, USA)
B.S. He (Nanjing Univ. China)
T. Luo (McMaster Univ., Canada)
M.J.D. Powell (Univ of Cambridge, UK)
M. Raydon (Univ Central de Venezuela, Venezuela)
E.W. Sachs (Univ Trier, Germany)
T. Terkelay (McMaster Univ., Canada)
L.Q. Qi (Hong Kong Polytechnic)
R. S. Varga (Kent St Univ., USA)
G. Yu (Univ of Texas, USA)
J.Z. Zhang (City Univ, Hong Kong)
S.L. Zhang (Tokyo Univ., Japan)

Limited number of short (about 20 minutes) contributed talks will be accepted
for presentation. Papers on theoretical, computational and practical aspects
of numerical optimization and numerical algebra are welcome.

For more details, please visit
or contact:

Prof. Ya-xiang Yuan and Dr. Yuhong Dai
Institute of Computational Mathematics and Scientific Computing
Chinese Academy of Sciences, P.O. Box 2719, Beijing 100080, China
Tel: +86-10-255-9001, +86-10-6254-5820 FAX: +86-10-254-2485


From: T. Terlaky <>
Date: Tue, 24 Jun 2003 10:37:00 -0400
Subject: McMaster Optimization Conference

3rd Annual McMaster Optimization Conference:
Theory and Applications
(MOPTA 03)
July 30 - August 1, 2003, McMaster University
Hamilton, Ontario, Canada

The 3rd annual McMaster Optimization Conference (MOPTA 03) will be
held at the campus of McMaster University. It will be hosted by the
Advanced Optimization Lab at the Department of Computing and Software
and it is co-sponsored by the Fields Institute and MITACS.

The conference aims to bring together a diverse group of people from
both discrete and continuous optimization, working on both theoretical
and applied aspects. We aim to bring together researchers from both the
theoretical and applied communities who do not usually get the chance
to interact in the framework of a medium-scale event.

Preliminary Program will be posted on the MOPTA WEB site before June 30, 2003.

REGISTRATION -- via the conference WEB site:
Deadline for early registration: June 30, 2003.

On behalf of the Organizing Committee
Tam=E1s Terlaky, (Chair, McMaster University)

Further information is available at

Distinguished researchers will give one-hour long invited talks on topics
of wide interest. Confirmed invited speakers include:

Laurent El Ghaoui, University of California, Berkeley, CA
"A Guided Tour of Robust Optimization"

Lisa K. Fleischer, Carnegie Mellon University, Pittsburg, PA
Approximately Optimal Control of Fluid Networks

Anthony Vannelli, University of Waterloo, Waterloo, ON

Masakazu Kojima, Tokyo Institute of Technology, Tokyo, Japan -
"Parallel Implementation of Primal-dual
Interior-point Methods for Semidefinite Programs"

George Nemhauser, Georgia Institute of Technology, Atlanta, GA
"A Branch-and-Cut Approach to Nonconvex Combinatorial Optimization"

Stratos Pistikopoulos, Imperial College, London, UK <bio_pistikopoulos.htm>
A guided tour to multi-parametric mixed integer and continuous programming"

Margaret H. Wright, Courant Institute, New York University, NY
"Direct Search Methods: What's the Fuss?"


From: Heike Fassbender <>
Date: Wed, 25 Jun 2003 09:40:59 +0200
Subject: GAMM Workshop on Applied and Numerical Linear Algebra

The 3. GAMM Workshop on Applied and Numerical Linear Algebra
will be held in Braunschweig, Germany, September 12.-13. 2003.

The special emphasis of the workshop is on 'Applications in
Medicine and Biology'. But all other aspects of applied and
numerical linear algebra are most welcome.

Invited speakers (confirmed):
Sabine Van Huffel (Katholieke Universiteit Leuven, Belgium)
Christof Schuette (Freie Universitaet Berlin, Germany)
Bernd Fischer (Medizinische Universitaet Luebeck, Germany)

The workshop will consist of three invited talks and 20 - 25
contributed talks of 25 minutes. The deadline for submission
of abstracts is July, 15.

More information on deadlines and registration can be found on the

Heike Fassbender


From: Giuseppe Ali <>
Date: Wed, 25 Jun 2003 11:57:51 +0200
Subject: Workshop in Naples on Optimization in Electromagnetism

First Announcement

MACSI-net Working Groups 02 and 13
Coupled Problems/Model Reduction
Electromagnetics in Telecommunication

Workshop on Optimization and Coupled Problems in Electromagnetism

September 22-23, 2003
Complesso Monte S. Angelo
Universit=E0 degli Studi Federico II
Naples, Italy

For information go to Events in the MACSI-net web page:
or go directly to the web page:


From: Simon Chandler-Wilde <>
Date: Fri, 27 Jun 2003 10:43:42 +0100 (BST)
Subject: Faculty Position at the University of Reading

The University of Reading
Department of Mathematics
Lecturer 'B' in Applied Mathematics/Numerical Analysis - Ref. AC0327
#26,270 to #33,679 p.a.

We are seeking an outstanding candidate to take up a Lectureship in
Applied Mathematics/Numerical Analysis, available from 1st October 2003.
The post has a particular emphasis on research in the first four years,
with correspondingly reduced teaching and administrative duties.

Applied Mathematics at Reading has an excellent research record in the
mathematics of waves, numerical analysis, reaction-diffusion theory,
computational and theoretical fluid dynamics, applied analysis and control
theory. It achieved a grade 5 in the last Research Assessment Exercise. It
has strong research links with the grade 5* Department of Meteorology
(within the new School of Mathematics, Meteorology and Physics) together
with the Met Office and several industrial partners. Candidates for the
Lectureship should have a growing research reputation in an area of
applied mathematics related to one or more of the current research
activities of the Department. Some priority may be placed on the
appointment of a candidate whose research overlaps with that of Dr Simon
Chandler-Wilde, who takes up the post of Professor of Applied Mathematics
at Reading from 1st September 2003.

The successful candidate will be expected to make a significant
contribution to the department's research output, to supervise
postgraduate students and postdocs, and to undertake teaching and
administrative duties, with a reduced (50%) teaching and administrative
load in the first four years. The appointment may be subject to a
probationary period in the first instance. Informal enquiries may be made
to Simon Chandler-Wilde (tel 0118 378 8991 - phone in the mornings if
possible), e-mail, web For further information about the department

Application form and further particulars available from the Personnel
Office, The University of Reading, Whiteknights, PO Box 217, Reading, RG6
6AH, telephone 0118 378 6771 (answerphone). Email
giving full name and address. Application Forms and Further particulars
are also available from Closing date for
applications:25 July 2003. Please quote Reference AC0327.


From: Arvid Naess <>
Date: Fri, 27 Jun 2003 13:47:28 +0200
Subject: PhD Position at Norwegian University of Science and Technology


Potential candidates for a PhD position at the Centre of Excellence for
Ships and Ocean Structures at the Norwegian University of Science and
Technology are cordially invited to submit their expression of interest.
Further information will then be provided (after July 20). The centre web

The ideal candidate is one with a good background in numerical analysis,
hydrodynamics or structural analysis and with a desire to participate in the
development of methods for response analysis of structures subjected to the
ocean environment. This would also involve the use of probabilistic methods.

With kind regards,
Arvid Naess


From: Tracy Paroz <>
Date: Fri, 27 Jun 2003 08:38:08 +1000
Subject: Software Engineer Position in Australia

Software Engineer - Grand Challenge Computational Science

Software Engineering, Parallel Computing, Computational Science Software
and Visualisation, Software Integration

Australian Computational Earth Systems Simulator (ACcESS) - UQ St Lucia

The Australian Computational Earth Systems Simulator (ACcESS) is a new
Major National Research Facility that is developing computational
software for simulating the dynamics of the solid earth, and
establishing a thematic parallel supercomputer with TFlops performance.

We are seeking three or more highly motivated candidates for the
position of a Software Engineer to work together with the team
developing the ACcESS software system. The successful candidates will be
assigned to software projects that will address grand challenge
scientific and computational problems in global to micro-scale earth
dynamics. These range from mantle convection and plate-mantle dynamics,
to earthquake physics and complexity, to mineralisation processes, to
rock micro-mechanics underlying earthquake and mining processes.

Applicants must have a strong background in advanced software
engineering, and a demonstrated ability in high performance computing
and/or scientific/computational software development and scientific
visualisation. The successful applicant will be based at the facility
headquarters at the University of Queensland.

Applicants should possess a degree in a relevant discipline, advanced
programming skills in C/C++ and object oriented software design and
implementation. Experience in high performance/parallel computing and
software development for discrete simulation methods - particularly
finite elements and/or discrete elements - is highly desirable.

This is an initial three year fixed term full-time appointment at
Academic Level B or C depending on level of experience and
qualifications. Three or more positions are available.

The remuneration package will be in the range of $65,249 - $77,483
(Level B) or $79,930 to $92,164 (Level C) per annum, which includes
employer superannuation contributions of 17%.

Obtain the position description and selection criteria here or contact
Ms Leonie Chen (07) 3365 7923 or Telephone
Professor Peter Mora, (07) 3365 7418, to discuss the position. Send
applications to the Human Resource Officer, Faculty of Engineering,
Physical Sciences and Architecture at the address below, or email

Closing date for applications:11 July 2003 Reference Number: 3006890


From: Ed Saff <>
Date: Tue, 24 Jun 2003 17:02:32 -0500 (CDT)
Subject: Contents, Constructive Approximation

Constructive Approximation Contents: vol. 19, #3, 2003.

"Markov-Berstein-Type Inequality for Trigonometric Polynomials
with Respect to Doubling Weights on [-w,w]" by T. Erdelyi, pp. 329-338.

"Tangential Markov Inequalities on Transcendental Curves" by L.P. Bos,
A. Brudnyi, N. Levenberg, and V. Totik, pp. 339-354.

"Characterizations of K-Functionals Built from Franctional Powers of
Infinitesimal Generators of Semigroups" by W. Trebels and U. Westphal,
pp. 503-528.

"Self-Affine Sets and Graph-Directed Systems" by X.-G. He, K.-S. Lau,
and H. Rao, pp. 373-397.

"An Extended Faber Operator" by J. Muller,
pp. 399-410.

"Zero Distribution of Bergman Orthogonal Polynomials for certain Planar
Domains" by A.L. Levin, E.B. Saff, and N.S. Stylianopoulos, pp. 411-436.

"(sqrt)3-Subdivision Schemes: Maximal Sum Rule Orders" by Q.T. Jiang, P.
Oswald, and S.D. Riemenschneider, pp. 437-463.

"Constrained Aproximation in Banach Spaces" by M. Smith, pp. 465-476.


From: Vladik Kreinovich <>
Date: Tue, 24 Jun 2003 19:01:42 -0600 (MDT)
Subject: Contents, Reliable Computing

Reliable Computing
Volume 9, issue 5, 2003

Special Issue
Proceedings of the Validated Computing 2002
May 23-25, 2002, Toronto, Canada
Guest Editor: R. Baker Kearfott

Interval Global Optimization in Solvent Design
Luke E. K. Achenie, Manish Sinha

Extended Interval Power Function
Walter Kraemer, Juergen Wolff von Gudenberg

Are There Easy-to-Check Necessary and Sufficient
Conditions for Straightforward Interval Computations
To Be Exact?
Vladik Kreinovich, Luc Longpre, James J. Buckley

Some Computer Assisted Proofs for Solutions of
the Heat Convection Problems
Mitsuhiro T. Nakao, Yoshitaka Watanabe,
Nobito Yamamoto, Takaaki Nishida

Reliable Computation of Frequency Response Plots
for Nonrational Transfer Functions to
Prescribed Accuracy
Paluri S. V. Nataraj, Jayesh J. Barve

Accelerated Enclosure Methods for Ordinary Free
Boundary Problems
Uwe Schaefer


Interval-Related Talks at the 2003 IEEE International
Conference on Fuzzy Systems


From: Oleg Burdakov <>
Date: Thu, 26 Jun 2003 20:19:00 +0200 (MEST)
Subject: Contents, Optimization Methods and Software

Table of Contents
Optimization Methods and Software (OMS), Volume 18, Number 3 (June, 2003)

Yuhong Dai, John D. Lamb and Wenbin Liu
Novel supervisor-searcher cooperation algorithms for
minimisation problems with strong noise

Abdel-Rahman Hedar and Masao Fukushima
Minimizing multimodal functions by simplex coding genetic algorithms

Quin Ni
A globally convergent method of moving asymptotes with
trust region technique

Michael Hinze and Thomas Slawig
Adjoint gradients compared to gradients from algorithmic
differentiation in instantaneous control of the Navier-Stokes

M. Kocvara and M. Stingl
PENNON: A code for convex nonlinear and semidefinite programming

Janos D. Pinter
Globally optimized calibration of nonlinear models:
techniques, software and applications

Forthcoming papers and complete table of contents for the journal OMS:

Online sample copy and Instructions for authors :


End of NA Digest