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

Today's Editor:
Cleve Moler
The MathWorks, Inc.
moler@mathworks.com

### Submissions for NA Digest:

Mail to na.digest@na-net.ornl.gov.

Information via e-mail about NA-NET: Mail to na.help@na-net.ornl.gov.

-------------------------------------------------------

From: Arieh Iserles <A.Iserles@damtp.cam.ac.uk>
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)
order:

"Monotone numerical schemes for nonlinear elliptic partial differential
equations"

Melvin LEOK
"Foundations of computational geometric mechanics"

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

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

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

Marc SCHWEITZER
"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 <charles.w.wampler@gm.com>
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 <nwagner@mecha.uni-stuttgart.de>
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.

Nils Wagner

------------------------------

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

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 trial license expires on May 1, 2004.

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

------------------------------

From: Andre Weideman <weideman@dip.sun.ac.za>
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:
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
http://dip.sun.ac.za/~weideman/research/differ.html
This Web Site replaces the old one at ucs.orst.edu, 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

------------------------------

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

Senior Research Fellow
Department of Applied Mathematics
Energy Systems Institute (SB)
130 Lermontov Str, Irkutsk 33 Russia
office phone: +7 3952 428 440
fax: +7 3952 426 796
e-mail : dsidorov@isem.sei.irk.ru

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

------------------------------

From: Ya-xiang Yuan <yyx@lsec.cc.ac.cn>
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)
M.J.D. Powell (Univ of Cambridge, UK)
M. Raydon (Univ Central de Venezuela, Venezuela)
E.W. Sachs (Univ Trier, Germany)
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.

http://lsec.cc.ac.cn/~yyx/nlao2003.html
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
e-mail: yyx@lsec.cc.ac.cn, dyh@lsec.cc.ac.cn

------------------------------

From: T. Terlaky <terlaky@mcmaster.ca>
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
http://www.cas.mcmaster.ca/~mopta

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.

SCOPE
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.

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

REGISTRATION -- via the conference WEB site:
http://www.cas.mcmaster.ca/~mopta/
Deadline for early registration: June 30, 2003.

On behalf of the Organizing Committee
Tam=E1s Terlaky, terlaky@mcmaster.ca (Chair, McMaster University)

Further information is available at
http://www.cas.mcmaster.ca/~mopta/

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>
"ON-LINE OPTIMIZATION VIA OFF-LINE OPTIMIZATION!
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 <h.fassbender@tu-bs.de>
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.

webpage

http://www.tu-bs.de/institute/numerik/GAMM03

Yours,
Heike Fassbender

------------------------------

From: Giuseppe Ali <ali@iam.na.cnr.it>
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:
http://www.macsinet.org
or go directly to the web page: http://www.na.iac.cnr.it/OCPEM2003

------------------------------

From: Simon Chandler-Wilde <Simon.Chandler-Wilde@brunel.ac.uk>
Date: Fri, 27 Jun 2003 10:43:42 +0100 (BST)
Subject: Faculty Position at 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
(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
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 Simon.Chandler-Wilde@brunel.ac.uk, web
www.brunel.ac.uk/~mastsnc). For further information about the department
see www.extra.rdg.ac.uk/Maths/index.asp.

Application form and further particulars available from the Personnel
giving full name and address. Application Forms and Further particulars
are also available from www.reading.ac.uk/Jobs. Closing date for
applications:25 July 2003. Please quote Reference AC0327.

------------------------------

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

PhD POSITION AT CENTRE OF EXCELLENCE AT THE 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
site: www.cesos.ntnu.no

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

arvidn@math.ntnu.no

------------------------------

From: Tracy Paroz <paroz@quakes.uq.edu.au>
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 l.chen@epsa.uq.edu.au. 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
l.chen@epsa.uq.edu.au.

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

------------------------------

From: Ed Saff <esaff@math.vanderbilt.edu>
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.

------------------------------

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
conference
Guest Editor: R. Baker Kearfott

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

Extended Interval Power Function
Walter Kraemer, Juergen Wolff von Gudenberg
339-347

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

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

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

Accelerated Enclosure Methods for Ordinary Free
Boundary Problems
Uwe Schaefer
391-403

Information

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

------------------------------

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

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
247-264

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

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

Michael Hinze and Thomas Slawig
differentiation in instantaneous control of the Navier-Stokes
equations
299-315

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

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