NA Digest Monday, January 31, 2011 Volume 11 : Issue 05

Today's Editor:
Daniel M. Dunlavy
Sandia National Labs

Submissions for NA Digest:

Mail to

Information via email about NA-NET:

Mail to


From: Pavel Solin <>
Date: Sun, 30 Jan 2011 18:50:06 -0800
Subject: Runge Kutta Brain Teasers

Difficulty level 1:
Can you tell which of the following names does not
belong to the rest?


Difficulty level 2:
Find an embedded implicit Runge-Kutta method.

I would be very grateful for any help. My search encompasses
embedded fully implicit as well as diagonally implicit, or ESDIRK
methods. So far, it appears that virtually all existing embedded
RK methods are explicit, see for example the Wikipedia page:


Pavel Solin

Pavel Solin
University of Nevada, Reno
Home page:
FEMTEC 2011:


From: Andrea Piacentini <>
Date: Wed, 26 Jan 2011 10:38:54 -0500
Subject: CERFACS' PALM coupler now open source (LGPL 3.0)

PALM is a dynamic parallel coupler developed at CERFACS.
Originally designed for data assimilation, PALM is now used for a
variety of applications in many fields of multi-physics coupling
ranging from oceanography to hydraulics, from hydrology to
agronomy, from aeronautics to space engineering and so on. PALM
based applications run on every size of computer, from laptop to
massively parallel or vector supercomputers. More than 250
scientists and engineers have been trained to PALM usage, and
every week new research laboratories, industries or computing
service companies join the PALM users' community. After more
than ten years of developments and thanks to the positive and
constructive users' feedback, PALM is now ready to become open
source under the version 3.0 of the LGPL license. For this reason
it has been renamed O-PALM, as a shorthand for Open-PALM.
You can find further detailed informations on O-PALM, including
contacts, on the dedicated web site

You can download the most recent O-PALM version (currently 4.0.0)
The release includes the PALM library, the PrePALM graphic user
interface and the documentation (including a tutorial).

You do not need to sign any agreement, but you are kindly
requested to register as a user. This will contribute to make the
O-PALM users' community getting more and more structured and to
share resources and hints. CERFACS regularly provides training
sessions with a frequency that will suit the users' needs. The
dates, conditions and fares for the training sessions are
constantly updated on the PALM web site:
CERFACS shareholders and CNRS institutions benefit a specific
training fare.


From: Dan Baruth <>
Date: Thu, 27 Jan 2011 19:42:47 -0500
Subject: Guaranteed Extended Precision - Timing Benchmark

Hi all,

As my quest for function evaluation time benchmarks (NA-Digest,
V. 10, # 51) was not successful, I have published one for the
Bessel function of the first kind, for guaranteed precision of 64
bits (i.e. 80-bit floating point). See details at\XP_Bessel.htm. Also, upon
request (for testing purposes) the Bessel function is available
for free.

Your comments and suggestions are welcome.

D. Baruth


From: David Silvester <>
Date: Wed, 26 Jan 2011 18:01:06 -0500
Subject: IFISS Version 3.1 Software Release

We are happy to announce the release of our

Incompressible Flow & Iterative Solver Software (IFISS) ver 3.1

The IFISS software toolbox may be used to generate typical linear
systems arising from finite element discretizations of steady and
unsteady diffusion, convection-diffusion, Stokes flow and
Navier-Stokes flow problems. It has built in geometric and
algebraic multigrid methods, state-of-the-art Krylov subspace
solvers, and it includes a variety of appropriate preconditioning

New features in Version 3.1 include
- a posteriori error estimation for Q2 approximation of diffusion
- a posteriori error estimation for Q2-Q1 and Q2-P1 (Navier-)
Stokes approximation,
- a proof of concept "optimal" Stokes solver (est_minres) with a
built in stopping criterion,
- fully implicit self-adaptive time stepping for stabiiized Q1-P0
and Q1-Q1 approximation of the Navier-Stokes equations.

IFISS is freely available and can be run with Matlab and Octave.
For download, user guide, and other information see or

David Silvester, Howard Elman, Alison Ramage and many IFISS


From: George Anastassiou <>
Date: Tue, 25 Jan 2011 06:08:31 -0500
Subject: New Book, Intelligent Mathematics: Computational Analysis

George A. Anastassiou of U. Memphis, USA, published a new
monograph, titled:

=93Intelligent Mathematics: Computational Analysis=94

Publisher: Springer, Series: Intelligent Systems Reference
Library, Vol. 5, 1st Edition., 2011, XVIII, 802 p. 1 illus.,


From: "Joppich, Wolfgang" <>
Date: Tue, 25 Jan 2011 04:13:38 -0500
Subject: Multigrid Course, Germany, Apr, 2011

There will be a MULTIGRID COURSE, Introduction to standard
methods, at the University of Applied Sciences Bonn-Rhein-Sieg,
from Monday 18.04.2011 to Wednesday 20.04.2011. For more
information contact or
look at

Please, pass this information to colleagues which might be
interested in such a course.


From: Pamela Bye <>
Date: Fri, 28 Jan 2011 10:50:04 -0500
Subject: IMA Math Education of Engineers, UK, Apr 2011

Mathematical Education of Engineers (MEE =96 2011)
Monday 11 April 2011, Loughborough University

This conference is the seventh MEE to be held at Loughborough for
the exchange of ideas, practices and pedagogy in the mathematical
education of engineers. This conference aims to reflect the
progress and experiences of initiatives within the teaching of
mathematics to engineers in recent years and to debate areas of
known concern. This meeting will provide an opportunity for
participants to network with specialist colleagues and learn
together from current practice. A one-day format is led by two
main contemporary themes:
i) How might we best teach Maths to Engineers? Interactive
discussion meeting including keynote presentations from
experienced colleagues: Dr Les Mustoe (Mathematics) and
Professor Graham Sander (Engineering)
ii) What Maths should we teach Engineers? Keynote talk =96
=91Requirements for Mathematical Education for Engineers from
an Industrial Point of View=92 by Dr Dirk Labuhn, Airbus
Operations GmbH, Bremen, Germany
These keynote events will be supported by parallel sessions
covering case studies; pedagogic outcomes from practitioners,
enhancements from the use of technology; experiences and
reflective studies on enhancing interaction and different
learning methods. Contributions from accepted speakers will
include curriculum development, improving understanding of
engineering mathematics and meeting the skills expectations of
graduate engineers.

Further details:


From: Ronald Boisvert <>
Date: Wed, 26 Jan 2011 11:45:51 -0500
Subject: Special Functions in the 21st Century, Apr 2011

The International Conference on Special Functions in the 21st
Century: Theory and Applications will be held in Washington DC
on April 6-8, 2011.

The conference will provide a forum for the exchange of
expertise, experience and insights among researchers in special
functions. Participants will include contributors to the
recently released NIST Digital Library of Mathematical
Functions (DLMF). It will also provide an opportunity for DLMF
users to interact with its creators and to explore potential
areas of fruitful future developments. The conference will also
celebrate the career of Professor Frank W. J. Olver, Mathematics
Editor of the DLMF.

Travel support for students and postdoctoral or early career
researchers is available.

Organizing Committee:
Daniel Lozier, National Institute of Standards and Technology
Adri Olde Daalhuis, University of Edinburgh
Nico Temme, CWI, Amsterdam
Roderick Wong, City University of Hong Kong

For details, see


From: Lisa Mc Goldrick <>
Date: Mon, 31 Jan 2011 09:33:40 -0500
Subject: 3rd Intl Conf on Num Meth for Finance, Ireland, Jun 2011

3rd International Conference on Numerical Methods for Finance
June 8th - 10th 2011, Kemmy Business School, Limerick, Ireland

The aim of this conference is to attract leading researchers,
both practitioners and academics, to discuss new and relevant
numerical methods for the solution of practical computational
problems in finance. Discussion of new products and new models
requiring new numerical methods for their solution is also
welcome. The conference is focused on topics of current interest,
and the programme is organised around the talks given by
distinguished invited speakers. Individual contributions may be
submitted to either an oral or a poster session. Participants are
encouraged to propose minisymposia consisting of a number of
papers on a given topic. Offers to organise panel discussions are
also welcome. Formal refereed proceedings will be published after
the conference.

Scientific Programme Committee
Rama Cont (Ecole Polytechnique, Palaiseau)
David Edelman (Smurfit Business School, Dublin)
Desmond Higham (University of Strathclyde, Glasgow)
Olivier Pironneau (University of Paris, Jussieu)
Wil Schilders (Technical University, Eindhoven)

You are now invited to submit an abstract (4 pages maximum) of an
oral or poster contribution to The
deadline for receipt of abstracts is 18th May 2011, but early
submission is strongly encouraged.

For further information, email or


From: Ilse Ipsen <>
Date: Mon, 24 Jan 2011 11:00:18 -0500
Subject: Graduate Student Modeling Workshop, Jul 2011

The 17th Industrial Mathematical & Statistical Modeling (IMSM)
Workshop for Graduate Students will take place at North Carolina
State University, 7-15 July 2011. The workshop is sponsored by
the Statistical and Applied Mathematical Science Institute
(SAMSI) and the Center for Research in Scientific
Computation (CRSC).

The IMSM workshop exposes graduate students in mathematics,
engineering, and statistics to exciting real-world problems from
industry and government. Besides giving students experience in
the team approach to problem solving, the IMSM workshop can help
them to decide what kind of professional career they want.

Local expenses and travel expenses will be covered for students
at US institutions. The application deadline is 15 April
2011. Information is available at
and questions can be directed to

With best regards,
Ilse Ipsen, Pierre Gremaud, and Ralph Smith


From: Paolo Bientinesi <>
Date: Wed, 26 Jan 2011 09:39:49 +0100
Subject: Corrected URL for Aachen Conf (AC.CES), Jul 2011

In the last issue of the Digest I announced the Aachen Conference
on Computational Engineering Science (AC.CES) July 13-15, 2011,
RWTH Aachen. Unfortunately there was a typo in the website

The correct link is

I take the opportunity to provide the list of confirmed speakers:

* Lorenz T. Biegler (Carnegie Mellon University)
* Peter G. Binev (University of South Carolina)
* Kai-Uwe Bletzinger (TU Munich)
* Albert Cohen (Universit=E9 Paris 6)
* Christodoulos A. Floudas (Princeton University)
* Omar Ghattas (University of Texas at Austin)
* Max Gunzburger (Florida State University)
* Matthias Heinkenschloss (Rice University)
* Barbara Kaltenbacher (University of Graz)
* Rudibert King (TU Berlin)
* Sandro Macchietto (Imperial College London)
* Anthony T. Patera (MIT)
* Rolf Rannacher (Heidelberg University)
* Stefan Ulbrich (TU Darmstadt)
* Karen Willcox (MIT)
* Enrique Zuazua (Basque Center for Applied Mathematics)


From: "marco lapegna" <>
Date: Tue, 25 Jan 2011 12:24:29 +0100
Subject: Hierarchical Parallelism in HPC Systems, Poland, Sep 2011

INVITATION AND CALL FOR PAPER: Workshop on Models, Algorithms and
Methodologies for Hierarchical Parallelism in new HPC Systems
Torun, Poland, September 11-14, 2011

In conjunction with PPAM 2011

From an architectural point of view, a High Performance Computing
system can be described by means of a hierarchical multi-level
structure: at the highest level there are several systems
connected among them by geographic networks (System level), an
intermediate level is composed by the nodes in a single system
communicating among them by means of dedicated fast networks or
high performance switches (Node level), at the lowest level,
finally, there are several computing elements, computing cores as
well as graphic accelerators, sharing resources in a single
CPU (Core level). These architecture levels have very different
features and they require different algorithmic development
methodologies. For such a reason, the development of algorithms
and scientific software for these system implies a suitable
combination of several methodologies to deal with the different
kinds of parallelism corresponding to each architectural
level. The general aim is then the development of hybrid and
hierarchical algorithms, able to be aware of the underlying

Submission of Papers: April 30, 2011
Notification of Acceptance: June 15, 2011
Conference: September 11-14, 2011
Camera-Ready Papers: Oct 31, 2011

Workshop webpage:


From: Josef Dick <>
Date: Wed, 26 Jan 2011 18:38:36 -0500
Subject: MCQMC 2012 Conference, Sydney, Feb 2012

The 10th International Conference on Monte Carlo and Quasi-Monte
Carlo Methods in Scientific Computing (MCQMC 2012) will be held
in Sydney, Feb 13-17, 2012.

The web site is now live. To
receive further announcements please go the web site, click on
the "mailing list" tab, and sign up.

MCQMC is a biennial conference devoted to Monte Carlo and
quasi-Monte Carlo methods and their interactions and
applications. (In brief, quasi-Monte Carlo methods replace the
random choices that characterize the Monte Carlo method by well
chosen deterministic choices.) For more information, click on
the "Background" tab on the web site.

The plenary speakers are Pierre Del Moral, Mike Giles, Fred
J. Hickernell, Aicke Hinrichs, Michael Lacey, Kerrie Mengersen,
Andreas Neuenkirch, Art B. Owen, Leszek Plaskota, Eckhard
Platen. See the web site for more information.

The web site includes a call for special sessions. Proposals for
special sessions should contain
* a short description of the theme and scope of the session
* the name(s) of the organizer(s)
* the names of four speakers (preferably from different
institutions) who have agreed to participate
The deadline for submitting proposals is July 1, 2011. There
will be a later call for contributed talks.

Ian H. Sloan, Frances Y. Kuo, Josef Dick, Gareth Peters
University of New South Wales,


From: Wim Michiels <>
Date: Wed, 26 Jan 2011 14:11:37 -0500
Subject: Postdoc Position, K.U. Leuven

Postdoc position: Numerical simulation of complementarity
systems, with applications to power electronic circuits

Power electronics are widely present in applications where one
form of electrical energy is converted into another. Examples
include power supplies in consumer electronics, industrial
electric motor drives, electro-heating, lighting and
energy-efficient interfaces between renewable energy resources
and the distribution grid.

Power electronic converters make use of semiconductor components
operating at a high switching frequency to reach the desired
dynamics at a slower time scale. Modeling and simulation methods
are indispensable for the design and analysis of these
converters. Because of the large difference in time scales,
non-smooth methods are needed to make a simulation possible in a
reasonable amount of time. Recent research results indicate the
potential of complementarity systems, a subclass of hybrid
systems, in modeling and simulation power electronic systems.

In the context sketched above we have a vacancy for a junior
postdoc (less than 5 years experience after PhD) whose task is to
develop numerical methods for time-integration (simulation) and
control of switched complementarity systems. The emphasis lies on
generic methods and software, whose application is not
exclusively limited to power electronic systems (other
applications include mechanical systems and dynamic

The postdoc position result from a joint project between the
Scientific Computing Research Group (htt://
and the Division Electa ( of
the K.U. Leuven. The requirements are a strong background in
computational mathematics, in particular numerical
time-integration (simulation), a strong interest in programming
and strong collaboration skills. The position will be initially
offered for one year.

Application instructions can be found at


From: Yukihiko Nakata <>
Date: Thu, 27 Jan 2011 18:48:14 +0100
Subject: PhD, Structured Population Dynamics, BCAM

BCAM, Basque Center for Applied Mathematics, is looking for one
Ph. D student position to join the research line Mathematical
Biology & Molecular Simulation for the following project funded
by Spanish Ministry of Science and Innovation (MICINN).

Reference code: MTM2010-18318
Principal Investigator: Philipp Getto

1. Description of the project

The project is devoted to the computation and analysis of regions
in parameter planes that characterize stability and instability
of equilibria for so called structured population models. In this
context, in particular, the principle of linearized stability and
Hopf bifurcations for infinite dimensional dynamical systems play
a role. Our fields of application include population dynamics in
cell biology, ecology and epidemiology.

2. General information

For information on the project, please check:

For information on the research line, please check:

For an overview of activities of BCAM, please check:

BCAM is in the process of relocating to the center of Bilbao, a
modern city with a rich cultural life as well as with many
possibilities for outdoor activities.

3. Required skill for candidates

Applicants should have completed a bachelor degree in mathematics
or related sciences or be currently enrolled in the last year.

4. Application process

The project is funded by the Spanish Ministry of Science and
Innovation (MICINN) and interested candidates should apply via
the website of MICINN

2.Click =93AYUDAS Y CONVOCATORIAS=94 on the right column.

3.Click =93LIA de Recursos Humanos=94 in =93Plan Nacional de I+D+i
2008-2011=94 on the left column.

4.Click =93Formaci=F3n de Personal Investigador (FPI-MICINN)=94 on the
center column.

The call is expected to be open in the first half of February. We
therefore suggest candidates to send us their contact details as
soon as possible to (Philipp Getto) with a copy
to (Yukihiko Nakata) so that we can provide
further information on the application process once the call is


From: Colin Cotter <>
Date: Mon, 24 Jan 2011 10:42:58 -0500
Subject: PhD, Stochastic Subgrid Models, Imperial College London

Funding is available for a PhD studentship supervised by Dr Pavel
Berloff (Department of Mathematics) and Dr Colin Cotter
(Department of Aeronautics) at Imperial College London on the
subject of "Stochastic models of unresolved eddies for ocean and
climate models", funded by the Grantham Institute for Climate

Project description: This project will investigate a new approach
to one of the Grand Challenges of 21st century science, namely
how to accurately model the effects of ocean eddies that are too
small to be represented on grids of ocean models. Due to
nonlinear interactions, even though these eddies are small, their
effects are felt on the global scale; without them, models fail
to predict ocean circulation patterns and their response to
global warming correctly. Recently, a very promising direction
has been to model subgrid eddies stochastically, i.e. as random
processes. In this project, we will use a computational dynamical
systems algorithm called the "Ensemble Kalman Filter", commonly
used in inverse problems, data assimilation and machine learning,
to construct such models; this innovative approach will lead to
new understanding of subgrid eddies and the potential of more
accurate ocean and climate predictions.

For more details see:

Students wishing to apply should contact Dr Colin Cotter
( before 13th February 2011.
International students are welcome and elligible for the
bursary, but will need to contribute the difference between
EU and international fees.


From: Emma Avery <>
Date: Mon, 31 Jan 2011 09:45:59 +0000
Subject: Contents, Inverse Problems, 27(2)


Volume 27, Issue 2, February 2011

Article numbers: 025001--025012

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

Reconstruction of river bed topography from free surface data
using a direct numerical approach in one-dimensional shallow
water flow. A F Gessese, M Sellier, E Van Houten and G Smart

Inversion formulas for the broken-ray Radon transform. Lucia
Florescu, Vadim A Markel and John C Schotland

Stability results for backward parabolic equations with
time-dependent coefficients. Dinh Nho H\`ao and Nguyen Van Duc

Detection and imaging in strongly backscattering randomly layered
media. R Alonso, L Borcea, G Papanicolaou and C Tsogka

Inverse problems for pseudo-Jacobi matrices: existence and
uniqueness results. N Bebiano and J da Provid\^encia

Sharp converse results for the regularization error using
distance functions. Jens Flemming, Bernd Hofmann and Peter

Identification of moving pointwise sources in an
advection--dispersion--reaction equation. M Andrle, F Ben
Belgacem and A El Badia

Source localization in electromyography using the inverse
potential problem. Kees van den Doel, Uri M Ascher and Dinesh K

State space regularization in the nonstationary inverse problem
for diffuse optical tomography. P Hiltunen, S S\"arkk\"a, I
Nissil\"a, A Lajunen and J Lampinen

Tensor completion and low-n-rank tensor recovery via convex
optimization. Silvia Gandy, Benjamin Recht and Isao Yamada

Identification of surface impedance of thin dielectric objects
from far-field data. Noam Zeev

Reconstruction of a function from its spherical (circular) means
with the centers lying on the surface of certain polygons and
polyhedra. Leonid Kunyansky


From: Marcin Paprzycki <>
Date: Sat, 29 Jan 2011 08:16:03 -0500
Subject: Contents, SCPE, 11(4)

Scalable Computing: Practice and Experience
Volume 11, no. 4 (December 2010)

Special Issue

Network Management in Distributed Systems
Editors: Jemal H. Abawajy, Mukaddim Pathan, Al-Sakib Khan Pathan
and Mustafizur Rahman


Special Issue Papers:

* Introduction to the Special Issue: Network Management in
Distributed Systems. Jemal H. Abawajy, Mukaddim Pathan, Al-Sakib
Khan Pathan and Mustafizur Rahman, page iii

* Implementation and Evaluation of a NAT-Gateway for the General
Internet Signaling Transport Protocol. Roland Bless and Martin
R=F6hricht, pages 329=96343

* An Efficient and User Privacy-Preserving Routing Protocol for
Wireless Mesh Networks. Jaydip Sen, pages 345=96358

* ERCTP: End-to-End Reliable and Congestion Aware Transport Layer
Protocol for Heterogeneous WSN. Atif Sharif, Vidyasagar
M. Potdar and A. J. D Rathnayaka, pages 359=96371

* Implementing Self-Protection in Distributed Grid Environment.
Inderpreet Chopra and Maninder Singh, pages 373=96383

Research Papers:

* Reputation-based Majority Voting for Malicious Grid Resources
Tolerance. Ahmed Bendahmane, Mohammad Essaaidi, Ahmed El
Moussaoui and Ali Younes

* A Graph-Based Approach to Context Matching. Andrei Olaru and
Adina Magda Florea, pages 393=96399

Book Reviews:

* CUDA by Example: An Introduction to General-Purpose GPU
Programming. Reviewed by Jie Cheng, pages 401=96402

End of NA Digest