NA Digest Saturday, May 15, 2004 Volume 04 : Issue 20

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

Submissions for NA Digest:

Mail to

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


From: Nils Wagner <>
Date: Fri, 14 May 2004 13:05:45 +0200
Subject: Iterative Refinement of Spectral Projectors

Dear experts,

Assuming that a crude initial approximation \tilde{P}_\Gamma(A) to the
spectral projector, which is defined by the matrix integral

P_\Gamma(A) = \frac{1}{2\pi i} \oint_Gamma (z I - A)^{-1} dz (1)

is available (by quadrature for example), how can we refine the
approximations ?

The spectral projector is idempotent

P^2_\Gamma(A)-P_\Gamma(A) = 0 (2)

and the matrix multiplication is commutative

P_\Gamma(A) A - A P_\Gamma(A) = 0 (3)

Does anyone know of any references for this problem ?

Thanks in advance.

Nils Wagner


From: Steven G. Johnson <>
Date: Fri, 14 May 2004 21:10:36 -0400 (EDT)
Subject: Harminv 1.0, Extracting Frequencies Better Than the FFT

ANNOUNCE: Harminv 1.0 -- extracting frequencies better than the FFT

I'm pleased to announce the availability of Harminv 1.0, a free program
and C library for harmonic inversion: decomposing a time-series into a sum
of sinusoids, including exponentially decaying sinusoids.

Harminv is an implementation of the "filter diagonalization method" (FDM)
of Mandelshtam & Taylor (see URL for references), which maps the harmonic
inversion problem onto a small eigen-problem (size proportional to the
number of sinusoids).

This method has been widely employed in physics since its inception in
1997, and is often able to obtain much more robust and accurate solutions
for the frequencies, etcetera, than e.g. direct least-squared fitting of
the data or its FFT spectrum (which are typically very ill-conditioned

I've been happily using this code for a few years now in my own research,
and I finally got around to releasing it. I hope that others find it

Steven G. Johnson


From: Daniel Stewart <>
Date: Tue, 11 May 2004 11:44:30 +0100
Subject: New Book, Orthogonal Polynomials: Computation and Approximation

New from Oxford University Press

Orthogonal Polynomials: Computation and Approximation

Walter Gautschi, Purdue University, USA

* Breaks new ground in the wide field of orthogonal polynomials and their
* Eminently practical, yet based on solid theory
* First available collection of relevant Matlab codes
* Contains ready-to-use computer programs
* Enables many applications of orthogonal polynomials without comprehensive
knowledge of underlying theory
* Provides ample references to literature

This is the first book on constructive methods for, and applications of
orthogonal polynomials. The book begins with a concise introduction to the
theory of polynomials orthogonal on the real line (or a portion thereof),
relative to a positive measure of integration. Topics which are particularly
relevant to computation are emphasized. The second chapter develops
computational methods for generating the coefficients in the basic
three-term recurrence relation. The methods are of two kinds: moment-based
methods and discretization methods. The former are provided with a detailed
sensitivity analysis. Other topics addressed concern Cauchy integrals of
orthogonal polynomials and their computation, a new discussion of
modification algorithms, and the generation of Sobolev orthogonal
polynomials. The final chapter deals with selected applications: the
numerical evaluation of integrals, especially by Gauss-type quadrature
methods, polynomial least squares approximation, moment-preserving spline
approximation, and the summation of slowly convergent series. Detailed
historic and bibliographic notes are appended to each chapter. The book will
be of interest not only to mathematicians and numerical analysts, but also
to a wide clientele of scientists and engineers who perceive a need for
applying orthogonal polynomials.

Contents: Basic Theory; Orthogonal polynomials; Properties of orthogonal
polynomials; Three-term recurrence relation; Quadrature rules; Classical
orthogonal polynomials; Kernal polynomials; Sobolev orthogonal polynomials;
Orthogonal polynomials on the semicircle; Notes to chapter 1; Computational
Methods; Moment-based methods; Discretization methods; Computing Cauchy
integrals of orthogonal polynomials; Modification algorithms; Computing
Sobolev orthogonal polynomials; Notes to chapter 2; Applications;
Quadrature; Least squares approximation; Moment-preserving spline
approximation; Slowly convergent series; Notes to chapter 3.


From: John C. Nash <>
Date: Thu, 13 May 2004 14:43:35 -0400 (EDT)
Subject: Nash and Walker-Smith Book and Software

Some time ago our book

an Integrated System in BASIC

(originally published by Marcel Dekker, New York, 1987) went out of print.
We took over its distribution via Nash Information Services Inc., and
added more text and software as

Example Problems and Software Extensions
(to accompany Nonlinear Parameter Estimation:
an Integrated System in BASIC)

There are still sporadic enquiries about it, but rather than bother with
arranging shipping and payment we are making all the text and software
available for download. You can find a link from

or go directly to the download site

We retain the copyright to all the materials, and request the courtesy of
proper acknowledgement or reference if you make use of them. We expect,
particularly, that some of the example problems and some of the
bibliographic references will prove helpful. Do let us know. FYI, the
publisher is officially

Nash Information Services Inc.
#903, 373 Laurier Ave. E., Ottawa, Ontario, K1N 8X6, Canada
(613) 236 6108

John C. Nash and Mary Walker-Smith


From: Suely Oliveira <>
Date: Wed, 12 May 2004 13:58:10 -0500
Subject: ETNA Special Issue on Combinatorial Scientific Computing

As annouced previously, ETNA will be having a special issue on
Combinatorial Scientific Computing. This issue was planned for the
CSC workshop (following the SIAM Parallel Processing Conference in
San Francisco) last February but is open to any author who may
have a contribution in a related area. For more information about the
special issue refer to:

Please note that the submission deadline was changed from May 31st to
August 9th.

Suely Oliveira, Associate Professor


From: Vasiliy Govorukhin <>
Date: Wed, 12 May 2004 16:02:33 +0400
Subject: Conference in Russia on Mathematical Hydrodynamics

International conference "Mathematical hydrodynamics: models and methods"
dedicated to the 70th Anniversary of Professor V.Yudovich

4-8 October 2004
Rostov on Don

The International conference "Mathematical hydrodynamics: models and methods"
devoted to the 70th Anniversary of professor V.Yudovich will be held in
Rostov-on-Don, Russia from 4th to 8th October of 2004. It will be organised
by Rostov State University .

The main goal of conference is Mathematical Hydrodynamics in a broad context.
The scope of the conference includes but not limited the following topics:

Mathematical theory of fluid dynamics: solvability and uniqueness.
Analytical dynamics and geometric-dufferential methods in hydrodynamics.
Stability of flows for ideal and viscous fluids.
Convective flows.
Asymptotics in hydrodynamics. Vibrodynamics. Parametric resonance.
Boundary layers.
Spectral theory in the stability problems of hydrodynamics.
Qualitative and numerical methods.
Bifurcation Analysis. Transitions. Systems with symmetry and cosymmetry.
Computer experiment methods and results.

Prof. Belokon' A.V. Rostov State University
Acad. Matishov G.G. RAS, Rostov on Don
Prof. Yudovich V.I. Rostov State University

Progam Committee
Prof. Andreev V.K. Krasnoyarsk State University
Acad. Babeshko V.A. Krasnodar State University
Prof. Belokon' A.V. Rostov State University
Prof. Vladimirov V.A. Hall University, UK
Prof. Gertsenshtein S.Ja. Moscow State University
Prof. Iooss G. Nice University, France
Prof. Zaslavsky G.M. New York University, USA
Prof. Karasozen B. Middle East Technical University, Turkey
Prof. Lyubimov D.V. Perm State University
Prof. Nasedkin A.V. Rostov State University
Acad. Plotnikov P.I. Mathematics Institute of RAS, Novosibirsk
Prof. Polezhaev V.I. Istitute of Appl. Mech. RAS, Noscow
Acad. Pukhnachev V.V. Hydrodynamics Institute of RAS, Novosibirsk
Prof. Sedenko V.I. Rostov Economics Academy
Prof. Simonenko I.B. Rostov State University
Prof. Yudovich V.I. Rostov State University

Preliminary list of lectures
Yudovich V.I. Vibrodynamics and vibrogeometry in continuum media
Pukhnachev V.V. Group approach in hydrodynamics: problems and perspectives
Polezhaev V.I. Computer experiment in turbulent convection
Vladimirov V.A. Vorticity effects and spray generation: theory and experiment
Gertsenshtein S.Ja. Nonlinear interaction of waves and stochastic fluid flows
Lyubimov D.V. Thermobibrational convection in compressible fluids
Plotnikov P.I., Iooss G., Toland G. Old problems and new achievments in surface waves theory

Organization Committee:
Dr. Govorukhin V.N. Rostov State University
Dr. Kurakin L.G. Rostov State University
Dr. Morgulis A.B. Rostov State University
Dr. Norkin M.V. Rostov State University
Dr. Tsybulin V.G. Rostov State University
Dr. Tsyvenkova O.A. Rostov State University


Secretary: S. V. Revina

Rostov State University
Mechanical mathematical department
Zorge St., 5,
Rostov on Don, 344090, Russia.


(8632) 221312



From: Tony Drummond <>
Date: May 10, 2004 7:52:52 PM EDT
Subject: Workshop in Berkeley on High End Computer Simulations

Fifth Workshop on the DOE
Advanced Computational Software Collection

Enabling Technologies For High End Computer Simulations
Place: Lawrence Berkeley National Laboratory
Date: August 24-27, 2004

The DOE Advanced CompuTational Software Collection (ACTS Collection, comprises a set of tools aim to simplify the
solution of common and important computational problems and have
substantially benefited a wide range of scientific and industrial
applications. These benefits are accounted not only for running
efficiently in high performing computing environments but also
realizing computation that would not have been possible otherwise.

Despite these successes, there is still a need for a greater
infrastructure to each out academia and industry through a
dissemination and instruction on the state-of-the-art
tools for high performance computing environments and
simultaneously provide an umbrella for tool developers to
receive the feedback from these communities. This
workshop is part of an approach to build such an infrastructure.

The four-day workshop will present an introduction to the
ACTS Collection for application scientists whose research
demand includes either large amounts of computation, the use of
robust numerical algorithms,
or combinations of these. The workshop will include a range
of tutorials on the tools (currently available in the collection and
some deliverables from the DOE SciDAC ISICs), discussion
sessions aimed to solve specific computational needs by the
participants, and hands-on practices using the NERSC's
state-of-the-art computers. We are planning to organize
parallel sessions and group the tutorials by topics, as follows:

+ Direct and Iterative Methods for the solution of linear and
non-linear systems of equations
+ PDE's and Multi-level Methods
+ Numerical Optimization
+ Structured and Unstructured meshes (Generation, Manipulation
and Computation)
+ Development of High Performance Computing applications
+ Performance monitoring and tuning
+ Grid computing
+ Software Interoperability

Who Should attend this workshop?

This workshop is opened to computational scientists from
industry and academia. Registration fees are fully sponsored by
the DOE's office of Science. In addition, DOE will sponsor travel
expenses for a limited number of graduate students and postdoctoral
fellows. This support includes round-trip transportation to and from
Berkeley, local transportation in Berkeley, lodging, breakfasts,
lunches and workshop materials.

The format of this workshop is tailored to the computational
problems and needs of its participants and functionality available
in the ACTS Collection. To effectively meet our goals, the number
of participants is limited and people interested in attending
must fill out an on-line application form no later than June 11, 2004.
In addition, graduate students and postdoctoral fellows applying
for travel support must provide a letter of support from their
advisor or supervisor. This letter must also be submitted on-line
and no later than the June 11, 2004 deadline.

For more information on the workshop, please contact Tony Drummond at
(510) 486-7624 or Osni Marques at (510) 486-5290.

Important Dates to Remember:

Proposal submission deadline: June 25, 2004
Proposal review completed and invitations sent: July 2, 2004
Attendee confirmation of participation deadline: July 16, 2004
Travel arrangements completed for all: July 30, 2004
Workshop Dates: August 24-27, 2004

Relevant Links:

ACTS Collection information:
On-line application form:
On-line support letter submission:


From: Volker Reichelt <>
Date: Wed, 12 May 2004 17:05:23 +0200 (CEST)
Subject: Faculty Position at RWTH Aachen

Full Professor Mathematics
RWTH Aachen
Faculty of Mathematics, Computer Science and Natural Sciences

We are seeking qualified applicants for teaching and research in the area
of Mathematics and Scientific Computing. The starting date is 01.10.2005.
Candidates with a distinguished expertise in Mathematical Modeling,
Numerical Simulation and Optimization, or Large Scale Scientific Computation
are particularly welcome. Scientific priorities could range from the design
and realization of large scale applications to theoretical foundations
in Applied Analysis. Key requirements are the ability and willingness to
conduct interdisciplinary collaborations with colleagues from Science and
Engineering in the newly founded "Centre for Computational Engineering
Science (CCES)". The candidate is expected to teach students in
Computational Engineering Science (CES) and in Mathematics.

The candidate should have a university degree, Ph.D., Habilitation
(particular qualification to teach at professorial level) or equivalent
academic achievements, and didactical skills.

Please send your written application to the

Dean of the Faculty of Mathematics, Computer Science and Natural Sciences
Prof. Dr. A. Krieg
RWTH Aachen
Templergraben 55
D-52056 Aachen

The deadline for applications is 31.07.2004.

Employment equity is one of our aims.

We welcome applications from persons with disabilities.

See also


From: Chandrika Kamath <>
Date: Tue, 11 May 2004 13:44:21 -0700
Subject: Postdoctoral Position at Lawrence Livermore National Laboratory

The Sapphire data mining project in the Center for Applied Scientific
Computing at the Lawrence Livermore National Laboratory is looking for
post-doctoral candidates in advanced image processing. The research focus
will include robust techniques for the identification of objects and
extraction of features from images. In addition, we are interested in
feature extraction in massive three-dimensional multi-variable data from
time-varying semi-structured and unstructured meshes that result from
computer simulations.

More details on the position are at under post-doc jobs, job#002295.

Information on the Sapphire project is available at

US citizens only.

Chandrika Kamath
Sapphire project lead


From: Fran Moshiri <>
Date: Tue, 11 May 2004 17:39:28 -0500
Subject: Postdoctoral Position at Rice University

Pfeiffer-VIGRE Instructorship
Department of Computational and Applied Mathematics
Rice University

The Department of Computational and Applied Mathematics (
invites applications for a postdoctoral Pfeiffer-VIGRE Instructorship.
Pfeiffer-VIGRE Instructors are appointed for two academic years, with
possible extension for a third year. The term of appointment may begin at
any time after July 1, 2004. Candidates should have received a PhD within
18 months preceding the time of appointment and have demonstrated potential
for excellence in both research and teaching. Candidates must be US citizens
or permanent residents.

The Pfeiffer-VIGRE Instructorships honor Professor Emeritus Paul Pfeiffer, a
founding member of the department, and are a component of the Rice VIGRE
program ( Pfeiffer-VIGRE Instructors teach one lecture
course each semester, and work with a team of faculty and graduate and
undergraduate students in one of the VIGRE research areas. The research
area for this appointment encompasses control, design, and inverse problems
and their interaction with computational harmonic analysis ("wavelets").
Candidates for this position will be expected to document interest and
experience in one or more of these fields.

Rice University is a private research university with a long tradition of
excellence in undergraduate and graduate science and engineering education.
The Department of Computational and Applied Mathematics hosts research
programs in optimization, numerical linear algebra, control and inverse
problems, and partial differential equations. Interdisciplinary work is a
fundamental aspect of the Department's program.

Applicants should send a letter of application, current vita, and
descriptions of research plans and teaching experience to

Staffing Committee
Department of Computational and Applied Mathematics
Rice University
6100 Main Street, MS-134
Houston Texas 77005-1892

At least three letters of reference are required; the candidate should have
these sent directly to the above address.

Rice University is an Affirmative Action, Equal Opportunity institution.


From: Peter V. Coveney <>
Date: Fri, 14 May 2004 11:53:50 +0100
Subject: Studentship at University College London

Centre for Computational Science, University College London
& Silicon Graphics Ltd

EPSRC Engineering Doctorate PhD Degree

High performance simulation, visualisation and grid computing

Applications are invited for a four-year EPSRC funded Ph.D. studentship
under a newly instituted Engineering Doctorate at University College
London, co-funded by Silicon Graphics Ltd, to perform research on one or
more of the following: (i) scalability of lattice-Boltzmann and/or
molecular dynamics algorithms and associated codes on future SGI computer
architectures; (ii) high-end visualisation capabilities for computational
steering of large scale simulation based on ray-tracing algorithms; (iii)
lightweight middleware to facilitate deployment of computational steering
for complex modelling and simulations.

The successful candidate will work closely alongside SGI scientific staff
worldwide, with access to upcoming SGI technologies, and with Professor


From: Raimondas Ciegis <>
Date: Thu, 13 May 2004 06:48:43 +0200 (WET)
Subject: Contents, Mathematical Modelling and Analysis


The Baltic Journal
on Mathematical Applications,
Numerical Analysis and Differential Equations

ISSN 1392-6292
Electronical edition:

Volume 9, Number 2, 2004


I.B. Bolotin, K.M. Rasulov
The first basic boundary value problem of Riemann's type for bianalytical
functions in a plane with slots
p. 91-98

S. Chakravarty, P.K. Mandal, A. Mandal
Numerical simulation of unsteady two-layered pulsatile blood flow in a
stenosed flexible artery: effect of peripheral layer viscosity
p. 99-114

R. {\v C}iegis, A. Dement'ev, G. {\v S}ilko
Parallel numerical modelling of short laser pulse compression
p. 115-126

T. Cirulis, H. Kalis, O. Lietuvietis
Comparative analysis by means of finite differences and DM methods for
linearized problem of gyrotrons
p. 127-136

M. Dubatovskaya, S. Rogosin
On an exact description of the Schottky groups of symmetries
p. 137-148

A. Misevi{\v c}ius
An improved hybrid optimization algorithm for the quadratic assignment
p. 149-168

I.V. Rybak
Monotone and conservative difference schemes for elliptic equations with
mixed derivatives
p. 169-178


From: Science Direct <>
Date: Thu, 13 May 2004 06:16:01 -0400 (EDT)
Subject: Contents, Journal of Approximation Theory

Journal of Approximation Theory Volume 126, Issue 1, Pages 1-113 (January 2004)

A tribute to Geza Freud, Pages 1-15
H. N. Mhaskar

Direct and converse results for the weighted rational approximation of functions with inner singularities, Pages 16-35
B. Della Vecchia

Short proofs of some theorems in approximation theory, Pages 36-39
Ingvar Mellander

On Blaschke products associated with n-widths, Pages 40-51
L. Baratchart, V. A. Prokhorov and E. B. Saff

Kolmogorov widths between the anisotropic space and the space of functions with mixed smoothness, Pages 52-59
Heping Wang and Yongsheng Sun

Uniqueness of best [phi]-approximation from the set of splines with infinitely many simple knots, Pages 60-113
A. Damas and M. Marano


From: Science Direct <>
Date: Wed, 12 May 2004 12:27:15 -0400 (EDT)
Subject: Contents, Nonlinear Analysis

Nonlinear Analysis Volume 57, Issue 3, Pages 323-484 (May 2004)

Evolution hemivariational inequalities with hysteresis, Pages 323-340
Leszek Gasiski

Existence and uniqueness of positive radial solutions for the Lane-Emden system, Pages 341-348
R.Robert Dalmasso

Universal bounds for global solutions of a forced porous medium equation, Pages 349-362
Michael Winkler

Stability theorems and existence results for periodic solutions of nonlinear impulsive delay differential equations with variable coefficients, Pages 363-399
Yuji Liu and Weigao Ge

Positive, unbounded and monotone solutions of the singular second Painleve equation on the half-line, Pages 401-419
P. K. Palamides and G. N. Galanis

Stability in a diffusive food-chain model with Michaelis-Menten functional response, Pages 421-433
Z.Zhigui Lin and Michael Pedersen

Solvability of nonlocal boundary value problems for ordinary differential equations of higher order, Pages 435-458
Yuji Liu and Weigao Ge

On solutions for a hyperbolic system with differential inclusion and memory source term on the boundary, Pages 459-472
Jong Yeoul Park and S.H.Sun Hye Park

Existence and optimal estimates of solutions for singular nonlinear Dirichlet problems, Pages 473-484
Zhijun Zhang and Jiangang Cheng

Nonlinear Analysis Volume 57, Issue 2, Pages 153-322 (April 2004)

A parabolic free boundary problem with double pinning, Pages 153-172
G. S. Weiss

Random fixed points of K-set- and pseudo-contractive random maps, Pages 173-181
Naseer Shahzad

Nonlocal Cauchy problems governed by compact operator families, Pages 183-189
Jin Liang, James Liu and T.-J.Ti-Jun Xiao

A study on local properties of Fourier series, Pages 191-197
H.Huseyn Bor

Existence of infinitely many weak solutions for a Neumann problem, Pages 199-209
Giovanni Anello

Abstract existence theorems of positive solutions for nonlinear boundary value problems, Pages 211-227
Yongxiang Li

Minimax control for evolutionary variational bilateral problem, Pages 229-252
Qihong Chen

On the existence of eventually positive solutions of fourth-order quasilinear differential equations, Pages 253-263
Manabu Naito and Fentao Wu

An Hopf-Lax formula for a class of measurable Hamilton-Jacobi equations, Pages 265-286
Fabio Camilli

Trajectories in Godel-type space-times under the action of a vector field, Pages 287-309
Elvira Mirenghi and Maria Tucci

Monotone penalty approximation of extremal solutions for quasilinear noncoercive variational inequalities, Pages 311-322
S. Carl and V.K.Vy K. Le


End of NA Digest