NA Digest Sunday, May 20, 2007 Volume 07 : Issue 21

Today's Editor:
Tamara G. Kolda
Sandia National Labs

Submissions for NA Digest:

Mail to

Information via email about NA-NET:

Mail to


From: Nick Trefethen <>
Date: Fri, 18 May 2007 19:38:46 +0100
Subject: Nick Higham elected a Fellow of the Royal Society

Congratulations to Nick Higham of the University of
Manchester, who is one of forty-four new Fellows of
the Royal Society elected May 17. The full list
of new fellows can be found at,
and also includes mathematicians Edwin Perkins,
Terence Tao, and Trevor Wooley.

Nick Trefethen, Oxford U.


From: "Martin Fuchs" <>
Date: Fri, 18 May 2007 16:34:52 +0200 (CEST)
Subject: Design optimization test problems

I am currently working on a tool involving design optimization with
discrete choice variables, of the form:

min f(z) s.t. z in S1xS2

where S1 is a finite subset of R^n and S2 is an m-dimensional box.

I would be very grateful for reference to test problems, preferably coded
in Matlab.

Martin Fuchs
Univ. of Vienna


From: Steven Ashby <>
Date: Sun, 20 May 2007 15:24:54 -0400
Subject: Fernbach Award Nominations - Deadline 30 June 2007

On behalf of the IEEE Computer Society and SC07, I invite you to nominate
outstanding candidates in the field of supercomputing for the Sidney
Fernbach Award. This award recognizes innovative uses of high performance
computing in problem solving. It was established in 1992 in memory of Sidney
Fernbach, one of the pioneers in the development and application of high
performance computers for the solution of large computational problems. A
certificate and $2,000 are awarded for innovative approaches and outstanding
contributions in the application of high performance computers.

The deadline for nominations is 30 JUNE 2007 (one month earlier than in
previous years). The nomination form, plus the required endorsements by
several other persons in the field, must be submitted to the IEEE Computer
Society web site ( Additional information can be
found at under the Technical Program tab;
also see the flyer Should you
have any questions, feel free to contact me via email at
We look forward to receiving your nominations.

Steven Ashby
IEEE Computer Society Fernbach Award Chair
SC07 Awards Chair


From: Bill McLean <>
Date: Mon, 14 May 2007 03:56:38 -0400
Subject: subject: f95 version of gaussq.f

The fortran 95 module gaussquad.f95 can compute the integration
points and weights for any of the classical Gauss quadrature rules,
including the Radau and Lobatto variants. Included is a low-level
routine that can handle any Gauss rule given the coefficients in the 3-term
recurrence relation for the associated orthogonal polynomials.
You can download the source distribution from

My f95 code is a goto-less version of the f77 routines in netlib's
gaussq.f, and should produce nearly identical results. The new package
adds an optional C language interface that relies on a couple of
the interoperability features of the f2003 standard, plus a few test
programs and examples in fortran, C and python.


From: William Mitchell <>
Date: Mon, 14 May 2007 08:36:42 -0400
Subject: PHAML 1.0 available

I am pleased to announce the availability of Version 1.0 of PHAML.

PHAML (Parallel Hierarchical Adaptive MultiLevel) is a Fortran 90
module for the solution of second order linear elliptic partial
differential equations of the form

-( p(x,y) * u ) -( q(x,y) * u ) + r(x,y) * u = f(x,y)
x x y y

with Dirichlet, Neumann (natural), Periodic, or Mixed boundary

where p>0, q>0, r, and f are functions of x and y.

It also solves elliptic eigenvalue problems where f is lambda * u and
the boundary conditions are homogeneous, and systems of elliptic PDEs.
Other types of differential equations (e.g. parabolic, nonlinear) can
be solved using PHAML as the core elliptic solver in an iteration.

PHAML features:

- low and high order finite elements on triangle grids
- h-, p-, and hp-adaptive mesh refinement based on newest node
- hierarchical basis multigrid linear system solver
- refinement-tree based partitioning method for load balancing
- arbitrary 2D connected, bounded domains, including curved
boundaries and holes
- Fortran 90 implementation
- message passing parallelism
- extensive visualization
- optional use of other popular packages for alternate methods

PHAML has been tested on clusters and distributed-memory parallel
computers with many Unix-like operating systems. The distribution
contains examples demonstrating its use for solving elliptic
equations, parabolic equations, nonlinear equations, systems of
elliptic equations, eigenvalue problems, and general domains.

For more information on the PHAML project, and to download the
software, please visit


From: Tim Davis <>
Date: Mon, 14 May 2007 11:20:35 -0400
Subject: KLU: sparse LU factorization for circuit simulation matrices

I would like to announce the release of Version 1.0 of KLU, a
mild-mannered sparse LU factorization algorithm targeted for matrices
arising in circuit simulation. It can be 100's of times faster than
x=A\b in MATLAB for these problems. KLU (the beta version) is used in
Xyce (developed at Sandia National Laboratories) and several commercial
circuit simulators.

See for details.

KLU is a left-looking method that does not exploit supernodes, in
contrast to its Supernodal left-looking cousin, SuperLU. Thus KLU = a
clark Kent LU. Version 1.0 adds a stable API, new features to reduce
the work in computing a suitable block triangular form, support for
64-bit integers, a very simple MATLAB interface, and two minor bug
fixes. SuperLU and other BLAS-based methods are faster than KLU for
matrices arising from 2D/3D geometries (for example), but less so for
circuit matrices which tend to be extremely sparse if ordered properly
and thus gain no benefit from the use of the BLAS.

Tim Davis, University of Florida


From: Peter Forsyth <>
Date: Mon, 14 May 2007 09:25:12 -0400
Subject: Computational Methods in Finance Conference: July 26-27, 2007

Computational Methods in Finance
July 26-27, 2007
University of Waterloo, Waterloo Ontario, Canada

This conference will focus on algorithmic and computational issues
applied to topics in finance. Areas covered include Monte
Carlo methods, numerical solution of partial differential
equations, portfolio optimization and stochastic optimization.

The Conference program is available at

More information, including a registration form,
is available at the Conference Website:


From: Peter Benner <>
Date: Fri, 18 May 2007 20:31:46 +0200
Subject: 1-Day Workshop on Matrix Equations, TU Chemnitz, Jun 2007


1-Day Workshop on Matrix Equations

will be held at Chemnitz University of Technology on June 15. The
workshop is devoted to latest developments in theory and numerical
algorithms for algebraic and differential matrix equations of
Sylvester, Lyapunov, Bernoulli, Riccati type and applications thereof.

The workshop is organized by

Peter Benner (TU Chemnitz)
Heike Fassbender (TU Braunschweig)
Lars Grasedyck (MPI for Mathematics in the Sciences, Leipzig)
Daniel Kressner (TU Berlin)

Talks will be 30 or 45 minutes, including discussion, and will start
10am. The meeting will end with a joint dinner.

Interested participants should contact Peter Benner at directly.


From: "Connie Young" <>
Date: Mon, 14 May 2007 11:27:03 -0400
Subject: SIAM Conference on Parallel Processing for Scientific Computing, Mar 2008

Conference Name: SIAM Conference on Parallel Processing for Scientific
Location: Atlanta, Georgia
Dates: March 12-14, 2008

The Call for Presentations for this conference is now available. Please
visit for more information.

Contact SIAM Conference Department at
with any questions about the conference.


From: Xiaobing Feng <>
Date: Fri, 18 May 2007 17:50:44 -0400
Subject: Computational Methods in Geosciences, Xian, China, July 5-7, 2007

Final Announcement of

on the occasion of Professor Jim Douglas, Jr.'s 80th birthday

July 5-7, 2007, Xi'an Jiaotong University, Xian, Shannxi, China

Workshop Website:

Organizing Committee:

Zhangxin Chen (Southern Methodist University)
Xiaobing Feng (University of Tennessee)
Kaitai Li (Xian Jiaotong University)
Dongwoo Sheen (Seoul National University)
Lihe Wang (University of Iowa)

For other inquires contact

Professor Xiaobing Feng ( <>
Professor Kaitai Li (


From: Luc GIRAUD <>
Date: Tue, 15 May 2007 06:09:04 -0400
Subject: Preconditioning 2007 in Toulouse, Jul 2007


International Conference Centre, Meteo France,
Toulouse, France

July 9-12, 2007

See the conference web site

Please note that the DEADLINE for REGISTRATION is June 15 and that
DEADLINE for the possible HOTEL REGISTRATIONS are also approaching.

A preliminary program is available on the Web.

The conference has been endorsed by the SIAM/SIAGLA and SMAI/GAMNI.

Academic and Industrial financial support for the conference is being
provided by CEA , CERFACS , EADS, ENSEEIHT, INRIA , IRIT , The Boeing
Company and Total.

The Preconditioning 2007 Conference Chairs,

Esmond G. Ng (Lawrence Berkeley National Laboratory)
Yousef Saad (University of Minnesota)
Wei-Pai Tang (The Boeing Company)


From: Kate Watt <>
Date: Tue, 15 May 2007 05:59:23 -0400
Subject: Nonlinearity High-Profile Articles of 2006

I am delighted to announce the publication of a special collection of
articles that highlight the breadth of high-quality research published in
Nonlinearity during 2006. Articles in this collection are freely available
until 31 December 2007 at:

These 'High-Profile Articles' consist of the top-downloaded papers published
in 2006, many of which have already been highly cited. Whether you are
reading these articles for the first time or from renewed interest, we hope
you will find them enjoyable and stimulating.

Information on how to read, write for or subscribe to Nonlinearity can be
found at or you can email the journal
publishing team at


From: Adam Phillips <>
Date: Fri, 11 May 2007 15:04:22 +0100
Subject: Inverse Problems Editorial Board highlights 2006

We are delighted to announce the Inverse Problems Editorial Board
highlights 2006:

The Editorial Board of Inverse Problems have selected articles published
in 2006 to be highlighted on the journal web page and we invite you to
view this taster of the high quality content published last year. This is
intended not as a list of the 'best' articles, but as an interesting and
stimulating reading list. Articles were selected for many reasons, some
contain outstanding research and breakthroughs, some may have an
especially clear exposition and are beautifully presented, others are
instructive, containing results and tools useful to many readers. We are
pleased to make these articles freely available and very much hope that
you will enjoy reading them.


From: Ronald Cools <>
Date: Fri, 18 May 2007 04:30:16 -0400
Subject: Faculty Position in Mathematical Engineering at Comp. Sci., K.U.Leuven

The Department of Computer Science, Faculty of Engineering of the K.U.Leuven,
Belgium announces a vacancy for a full time tenured or tenure-track position
in the field of mathematical engineering.

Candidates should have substantial internationally recognized research
expertise in one or more areas of mathematical engineering and must be able
to bridge the gap between fundamental numerical techniques and engineering
applications. The candidate has expertise in one or more of the following
(1) the design and analysis of fundamental numerical techniques,
(2) the development of scientific software,
(3) numerical modeling, simulation and analysis of complex systems.
He/she should also have teaching experience in these areas. The candidate is
expected to contribute to the new master program on mathematical engineering
at K.U.Leuven and/or to related programs.

Applicants should have demonstrated the ability to supervise research
projects and to attract external research funding. His/her activities
will be integrated in the Numerical Analysis and Applied Mathematics
section of the Department of Computer Science. Information about the
section and its research groups can be found via the departments website .

The candidates must have a PhD in Sciences or Applied Sciences. Initially
candidates can teach in English. They must commit themselves to learning
Dutch, so that they are able to teach in Dutch after a trial period.

The appointment will be made in the rank of assistant professor or higher,
depending on the qualifications of the candidate. The K.U.Leuven is an equal
opportunity employer.

The appointment will start in October 2008. Details about the application
procedure can be found at .
Inquiries about the position should be directed to Pierre Verbaeten, head
of the department of Computer Science (
The evaluation of candidates will begin on October 1, 2007.


From: Brenda Willoughby <>
Date: Fri, 18 May 2007 08:35:31 -0400
Subject: University Lecturer, Mathematical Institute, Oxford University

Mathematical, Physical and Life Sciences Division
Mathematical Institute
in association with The Queen's College

The Mathematical Institute proposes to appoint a University Lecturer in the
Analysis and Computation of Nonlinear Partial Differential Equations with
effect from 1 October 2007 or as soon as possible thereafter. The
successful candidate will be offered a non-tutorial fellowship by The
Queen's College, under arrangements described in the further particulars.
The salary will be on a scale up to £51,095 per annum.

The successful candidate will have a strong record of research of an
international standard of excellence in the analysis and computation of
nonlinear partial differential equations, and will be expected to contribute
effectively to teaching and administration.

Further particulars, containing details of the application procedure and of
the duties, may be obtained from the Administrative Assistant (Vacancies),
Mathematical Institute, 24–29 St Giles', Oxford OX1 3LB (email or by visiting The
closing date for applications is 18 June 2007. Please quote reference
number BK/07/014.

The University is an Equal Opportunities Employer.


From: Neville Ford <>
Date: Wed, 16 May 2007 20:35:04 +0100
Subject: Lecturer/Senior Lecturer Position at University of Chester, UK


£25,334 - £39,033 per annum

The Department of Mathematics at the University of Chester is seeking to appoint
a Lecturer/Senior Lecturer to teach at all levels within the undergraduate and
postgraduate Mathematics curricula and to undertake research. This post
offers the successful candidate an excellent opportunity to contribute
to all aspects of the work of our compact and successful department. You
must possess a PhD in Mathematics. Experience of teaching University
Mathematics and successful publication of research results will be an

Quote reference number: HRMS/221

Candidates may download further details and an application form from our
web site at
HRM Services, University of Chester, Parkgate Road, Chester CH1 4BJ, UK
quoting the relevant reference number.

Closing date: Monday 4th June 2007
Working towards equality of opportunity


From: Pierre Ramet <>
Date: Tue, 15 May 2007 11:33:18 -0400
Subject: Postdoctoral position in MHD simulation at Bordeaux University

The context of this work falls within the scope of ASTER project that intends
to develop and implement methods to improve the simulation of MHD
instabilities that are needed to evaluate mechanisms to control the energy
losses observed in the standard tokamak operating scenario (ITER).

The applications we study are complex in that they require simulations at a
very large scale; several tens or hundreds teraflops of computation using
several terabyte of data are often needed. To carry out these more and more
accurate full-scale simulations without increasing the number of unknowns in a
uniform way, a technique consists in handling a finer grid where the solution
varies abruptly and a coarser grid at other places accordingly to quite
specific criteria.
This approach is broadly used when explicit numerical schemes are considered.
Here, the characteristic time scales may be very different so that an implicit
scheme is needed.

The goal of this work is first to make a state of the art of the mesh
refinement methods based on numerical implicit schemes in particular for
irregular grids, then to study within a parallel framework the principal
difficulties encountered in these methods such as:
management of the hierarchy of the grids which themselves are distributed,
criteria of "refinement/unrefinement", management of the versatility of the
parallel solvers. The effectiveness of such a method will therefore require
the use of dynamic load balancing methods over the set of processors to
compensate the variability of the grid during the simulation.

Please send your CV, a cover letter describing your research interests, and
the names and email addresses of at least 2 references to:
O. Coulaud and/or P. Ramet, INRIA, Univ. Bordeaux 1, 33405 Talence, FRANCE


From: Frederic Nataf <>
Date: Mon, 14 May 2007 17:55:40 +0200
Subject: PhD Position at French Petroleum Institute / Paris VI University

PhD Position at French Petroleum Institute / Paris VI University

Parallel Linear Solvers for Large Sparse Systems arising from
Petroleum Industry

The French Petroleum Institute invites applications for a PhD
position starting in October 2007.
It will support the research on parallel linear solvers for large
sparse systems arising from petroleum industry.

Desirable qualifications:
- strong background in numerical linear algebra and/or
numerical methods for PDEs
- good programming skills

Duration : three years
Location : French Petroleum Institute and University Paris VI, Paris
Salary (gross wage includes medical care) : first year : 1900 Euros,
second year : 2050 Euros, third year : 2200 euros

Send your application (CV + references) to

Frederic Nataf
Jacques Louis Lions Laboratory, University Paris VI


Date: Wed, 16 May 2007 10:29:27 +0200
Subject: PhD position in XFEM-simulation for fluids / Aachen

PhD position in XFEM-simulation for fluids

Chair for Computational Analysis of Technical Systems
Aachen University, Germany

A PhD position is available at Aachen University from July 1st,
2007. The position is part of the project “Numerical methods for
discontinuities in continuum mechanics” which is supported by the
prestigious Emmy-Noether-program of the German Research Society

The applicant will work in the field of computational fluid dynamics
and concentrate on applications where discontinuities are present in
the field variables (e.g. two-fluid flows, flows involving shocks and
small boundary layers). The extended finite element method (XFEM) is,
in general, highly suited for applications involving discontinuities
and has been intensively studied in solid mechanics. The development
and application of the XFEM for fluid mechanics is the aim of this

A thorough knowledge of the classical finite element method (FEM) is
required, and experience in computational fluid dynamics is
preferred. Programming experience is a must, preferably in C.

Please send your application to: Thomas-Peter Fries, Chair for CATS,
Aachen University, Steinbachstr. 53 B, 52074 Aachen, Germany. Your
application by email is welcome: Deadline
for applications is June 8th, 2007.


From: Johannes Kraus <>
Date: Wed, 16 May 2007 03:52:25 -0400
Subject: Ph.D. Position at RICAM, Austria

In the framework of a three-year research project on "Algebraic Multigrid and
Multilevel Methods for Vector-Field Problems", funded by the Austrian Science
Foundation, a Ph.D. position is available at the Johann Radon Institute for
Computational and Applied Mathematics (RICAM) starting with July 1, 2007.

For this challenging position we seek a highly motivated candidate.
A strong background in applied mathematics, in particular, in numerical
methods for partial differential equations and in numerical linear algebra,
and excellent programming skills are required. Familiarity with finite
element methods and (algebraic) multigrid methods is desirable. Some
knowledge of electromagnetics would be advantageous.

Informal enquiries can be made to Dr. Johannes Kraus at:

A complete application including academic records (copies of relevant
documents, a statement about scientific interest and achievements related to
the master thesis) should be sent to
Dr. Johannes K. Kraus
Johann Radon Institute
Austrian Academy of Sciences
Altenbergerstraße 69
A-4040 Linz Austria

Alternatively, the application material can be sent by email to

Closing date for applications is May 27, 2007.
The Austrian Science Foundation is an equal opportunity employer.


From: Romas Baronas <>
Date: Fri, 18 May 2007 14:55:10 +0300
Subject: Contents, Nonlinear Analysis: Modelling and Control

Nonlinear Analysis: Modelling and Control, ISSN 1392-5113,
Volume 12, Number 2, 2007

A free on-line edition is available at:


The Nehari Manifold for p-Laplacian Equation with Dirichlet Boundary Condition,
Pages 143-155.
G.A. Afrouzi, S. Mahdavi, Z. Naghizadeh.

A Study of Steady Buoyancy-Driven Dissipative Micropolar Free
Convection Heat and Mass Transfer in a Darcian Porous Regime with Chemical
Reaction, Pages 157-180.
O. Anwar Bég, R. Bhargava, S. Rawat, H.S. Takhar, Tasweer A. Bég.

Process Distribution in the Network Systems, Pages 181-189.
J. Augutis, E. Uspuras, R. Krikstolaitis, V. Matuzas.

Conjugate Effects of Heat and Mass Transfer on Natural Convection Flow Across an
Isothermal Horizontal Circular Cylinder With Chemical Reaction, Pages 191-201.
Md.A. Hye, Md.M. Molla, M.A.H. Khan.

Natural Convection in a Square Cavity in the Presence of Heated Plate, Pages
P. Kandaswamy, J. Lee, A.K. Abdul Hakeem.

Some Advanced Fiber-Optic Amplitude Modulated Reflection Displacement and
Refractive Index Sensors, Pages 213-225.
V. Kleiza, J. Verkelis.

A Nonlinear Dynamical Model to Study the Removal of Gaseous and Particulate
Pollutants in a Rain System, Pages 227-243.
R. Naresh, S. Sundar.

Computational Modeling of the Amperometric Bioanalytical System for Lipase
Activity Assay: a Time-Dependent Response, Pages 245-251.
M. Puida, F. Ivanauskas, I. Ignatjev, G. Valinèius, V. Razumas.

Nonlinear Spectra for Parameter Dependent Ordinary Differential Equations, Pages
F. Sadyrbaev, A. Gritsans.

On Similarity Solutions and Interface Reactions for a Vector-Valued Stefan
Problem, Pages 269-288.
F.J. Vermolen.

For a paper submission, please refer to

Dr. Romas Baronas, Journal Secretary,
Nonlinear Analysis: Modelling and Control


From: Liz Martin <>
Date: Fri, 11 May 2007 10:59:16 +0100
Subject: Contents, Nonlinearity, volume 20, issue 6, June 2007


Volume 20, Issue 6, June 2007

Individual articles are free for 30 days following their publication on
the web. This issue is available at: URL:

Pages: 1313--1534


A KAM theorem for Hamiltonian networks with long ranged couplings
J Geng and Y Yi

Critical Fujita exponents for a class of quasilinear equations with
homogeneous Neumann boundary data
C Wang, S Zheng and Z Wang

The amplitude equations for the first instability of electro-convection in
nematic liquid crystals in the case of two unbounded space directions
G Schneider and H Uecker

Resonant eigenstates for a quantized chaotic system
S Nonnenmacher and M Rubin

Some semi-linear parabolic equations with nonlinear flux boundary
P Grindrod and E Tuma

Statistical properties of nonlinear shell models of turbulence from linear
advection models: rigorous results
R Benzi, B Levant, I Procaccia and E S Titi

Cnoidal wave solutions to Boussinesq systems
H Chen, M Chen and N V Nguyen

Topology and separation of self-similar fractals in the plane
C Bandt and H Rao

Uniqueness for the Navier--Stokes problem: remarks on a theorem of
Jean--Yves Chemin
P G Lemari\'e-Rieusset

A Borel--Cantelli lemma for intermittent interval maps
S Gou\"ezel

Numerical continuation of normally hyperbolic invariant manifolds
H W Broer, A Hagen and G Vegter

End of NA Digest