NA Digest Saturday, September 20, 2003 Volume 03 : Issue 38

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

Submissions for NA Digest:

Mail to

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


From: Joel Malard <>
Date: Mon, 15 Sep 2003 11:38:37 -0700
Subject: Monte Carlo Speedup

It is being told that the best way to achieve perfect speedup of Monte Carlo
simulations is to simply do as many runs as there are processors. This sounds
like a half truth because if the two runs use the same random numbers there is
no speedup at all.

Are there known results that quantify speedup as a function of statistical
independence among random number streams?

Would you know of any real-life example where incorrect conclusions were drawn
from Monte Carlo because of invalid assumptions about the sequence of random

With best regards,

Joel M. Malard, Ph.D.
Senior Research Scientist II
Pacific Northwest National Laboratory
Battelle Boulevard, PO Box 999
Mail Stop K1-85
Richland, WA 99352


From: Nils Wagner <>
Date: Thu, 18 Sep 2003 11:22:38 +0200
Subject: Spectral Projectors for Matrix Polynomials

Dear NA readers,

The spectral projector for a linear matrix polynomial (i.e.,
generalized) is

P_\Gamma(A,B) = \frac{1}{2\pi i} \oint\limits_\Gamma (\lambda B-A)^{-1}
B d\lambda

for A x = \lambda B x. Is it possible to find a similar expression for
quadratic matrix polynomials (A + \lambda B + \lambda^2 C)x = 0 ?

It is well known, that a matrix polynomial can be reduced to a (large)
generalized problem, but
is there any other possibility with respect to the representation of the
spectral projector ?

Does anyone know of any references for this problem ?

Nils Wagner


From: Yves Renard <>
Date: Mon, 15 Sep 2003 18:18:10 +0200
Subject: GMM++, A Generic C++ Matrix Library

GMM++ 1.5 Available for Download

This generic matrix library is built as an interface for already existing
class of vectors and matrices. It provides also predefined type of dense,
skyline and sparse matrices. It is largely inspired from MTL/ITL but we hope
a little bit simpler to use. matrix x matrix multiplication works for any
type of matrices (even with mixed types).

An important extension is the possibility to access to sub-matrices for read
or write operations with any kind of interfaced matrix.

Linear solvers has been imported from ITL (CG, BICGSTAB, GMRES, QMR, dense
LU, ...) with classical preconditionners (Incomplete cholesky, ILU, ILUT ...)

The performance is comparable with the one of MTL, slightly better on some
solvers we optimized.

Inpout / output with standard formats is also provided.

This library is available under LGPL licenseand may be downloaded at:

A short documentation exists.

Yves Renard (
Dept de Mathematiques, INSA de Toulouse
Complexe Scientifique de Rangueil
31077 Toulouse Cedex, FRANCE


From: Marielba Rojas <>
Date: Fri, 19 Sep 2003 16:33:07 +0200
Subject: MATLAB Software for Large-Scale Trust-Regions Subproblems

Dear Colleagues,

We are pleased to announce the availability of our Matlab software for
Large-Scale Trust-Region Subproblems.

The software implements the LSTRS method for large-scale quadratic problems
with a quadratic constraint.

To obtain the software please go to:

Marielba Rojas
Sandra Santos
Dan Sorensen


From: Nadja English <>
Date: Fri, 19 Sep 2003 10:40:12 -0400
Subject: New Book, A Numerical Library in JAVA

New From CRC Press
A Numerical Library in JAVA for Scientists and Engineers by Hang T. Lau

This groundbreaking handbook represents the only inexpensive numerical library
available for the use of JAVA in scientific computation. It provides a library
of Java-based procedures that perform fundamental mathematical computations
critical to researchers. It includes procedural descriptions of linear algebra,
the numerical solutions of ordinary and partial differential equations,
optimization, parameter estimation, and special functions of mathematical
physics. The handbook also enables the user to experiment with the many JAVA
procedures provided in order to observe computational success and simulate

ISBN 158484304; Catalog Number C4304; 1,088 pp.; Pub Date 8/27/03; $99.95/


From: Abhijit Gosavi <>
Date: Sat, 13 Sep 2003 18:31:24 -0400
Subject: New Book, Reinforcement Learning and Simulation Optimization

Colleagues and Friends:

Kluwer Academic has recently published a new book on simulation-based
optimization and reinforcement learning.
More information about the book can be found at the following website.

Here is a brief description of the book's contents.


This book provides a unified treatment to static and dynamic models of
simulation-based optimization. It aims at providing an introduction to
exciting developments in the theory of simulation-based optimization in an
accessible format. The book divides the field of simulation-based
optimization into two parts: static and dynamic optimization. The methods
discussed for dynamic optimization are primarily Reinforcement Learning and
Learning Automata, while those for static optimization are Simultaneous
Perturbation and Meta-heuristics.


From: Russ Merris <>
Date: Fri, 19 Sep 2003 15:43:03 -0700
Subject: New Book on Combinatorics

New Book

Russell Merris, Combinatorics, Second Edition, Wiley-Interscience, 2003
ISBN 0-471-26296-X

Table of Contents (Asterisks indicate optional sections.)
Chapter One: The Mathematics of Choice
1.1 The Fundamental Counting Principle
1.2 Pascal's Triangle
*1.3 Elementary Probability
*1.4 Error Correcting Codes
1.5 Combinatorial Identities
1.6 Four Ways to Choose
1.7 The Binomial and Multinomial Theorems
1.8 Partitions
1.9 Elementary Symmetric Functions
*1.10 Combinatorial Algorithms

Chapter Two: The Combinatorics of Finite Functions
2.1 Stirling Numbers of the Second Kind
2.2 Bells, Balls, and Urns
2.3 The Principle of Inclusion and Exclusion
2.4 Disjoint Cycles
2.5 Stirling Numbers of the First Kind

Chapter Three: Polya's Theory of Enumeration
3.1 Function Composition
3.2 Permutation Groups
3.3 Burnside's Lemma
3.4 Symmetry Groups
3.5 Color Patterns
3.6 Polya's Theorem
3.7 The Cycle Index Polynomial

Chapter Four: Generating Functions
4.1 Difference Sequences
4.2 Ordinary Generating Functions
4.3 Applications of Generating Functions
4.4 Exponential Generating Functions
4.5 Recursive Techniques

Chapter Five: Enumeration in Graphs
5.1 The Pigeonhole Principle
*5.2 Edge Colorings and Ramsey Theory
5.3 Chromatic Polynomials
*5.4 Planar Graphs
5.5 Matching Polynomials
5.6 Oriented Graphs
5.7 Graphic Partitions

Chapter Six: Codes and Designs
6.1 Linear Codes
6.2 Decoding Algorithms
6.3 Latin Squares
6.4 Balanced Incomplete Block Designs

A1 Symmetric Polynomials
A2 Sorting Algorithms
A3 Matrix Theory
Hints and Answers to Selected Odd Numbered Exercises

Highlights of Second Edition enhancements: Smoothed and polished
exposition with a sharpened focus on key ideas; expanded discussion of
linear condes; new optional section on algorithms; greatly expanded
hints and answers section; and many new exercises and examples.


From: Joanna Littleton <>
Date: Thu, 18 Sep 2003 09:51:26 -0400
Subject: SIAM W. T. and Idalia Reid Prize


SIAM will present the W. T. and Idalia Reid Prize at the SIAM Annual Meeting
scheduled to be held July 12-16, 2004, in Portland, Oregon.

The award will be given for research in, or other contributions to, the broadly
defined areas of differential equations and control theory. The prize may be
given either for a single notable achievement or for a collection of such
achievements. Committee Chair John Burns wishes to stress the breadth of the
eligible fields. He welcomes nominations in all areas of differential
equations (ODEs, PDEs, etc.), both analytical and numerical, as well as in
control theory.


The prize is awarded to any member of the scientific community who meets the
general guidelines of the prize description above.

Description of Award

The award consists of an engraved medal and a $10,000 cash prize.


A letter of nomination, including a description of achievement(s), should be
sent to the address below. Nominations must be received in the SIAM office by
January 30, 2004.

Professor John A. Burns
Chair, Reid Prize Selection Committee
c/o Joanna Littleton
3600 University City Science Center
Philadelphia, PA 19104-2688
Telephone: 215-382-9800 ext. 303
Fax: 215-386-7999

Selection Committee

Members of the selection committee are John A. Burns (Virginia Polytechnic
Institute and State University), Chair; H. Thomas Banks (North Carolina State
University); Max Gunzburger (Florida State University); James P. Keener
(University of Utah); Alan J. Laub (University of California, Davis).


From: Bette Byrne <>
Date: Thu, 18 Sep 2003 10:06:46 +0100
Subject: Conference at Oxford on Computational Fluid Dynamics


The 8th ICFD International Conference on Numerical Methods for Fluid Dynamics
organised by the ICFD (Institute for Computational Fluid Dynamics) will be
held at St Catherine's College, Oxford, from 29 March - 1 April 2004.

The following have agreed to give invited talks: Remi Abgrall (Bordeaux),
Jerry Brackbill (Los Alamos), Tom Hou (Caltech), Tony Hutton (QinetiQ Ltd),
Andy Keane (Southampton), Roland Keunings (CESAME, Louvain), Mike King (BP),
Anthony Patera (MIT), Piotr Smolarkiewicz (NCAR Boulder), John Trangenstein
(Duke), Alessandro Veneziani (Milan), Andy Wathen (Oxford). In addition to
invited lectures the programme will include contributed talks of twenty
minutes and poster sessions.


Two page abstracts for contributed papers should be submitted by 25 NOVEMBER
2003, stating a preference for oral or poster presentation. NOTIFICATION OF
ACCEPTANCE WILL BE GIVEN BY 21 JANUARY 2004. All accepted papers (including
those for posters) will be required at the meeting for distribution immediately
after the Conference on CD. All accepted papers will be refereed for inclusion
in the Proceedings which will be published, as at the last Conference, in a
Special Issue of the International Journal for Numerical Methods in Fluids.


A feature of the meeting will be the fourth award of "The Bill Morton Prize"
for a paper on CFD by a young research worker. The Prize paper, which should
be suitable for a 30 minute talk, will be chosen from eligible contributed
papers submitted for the Conference and marked "for the Bill Morton Prize".
The shortlist of papers will be presented by the authors at a special session
at the Conference. Joint papers must be accompanied by a statement from
collaborators confirming that the author made a substantial contribution.

Candidates should be under the age of 31 on 1 April 2004. Each candidate
should include a statement of his/her age on that date and should indicate
that he/she would be available to present his/her paper at the Conference.

The Prize will be presented at the Conference dinner.

Further details can be found at

Abstracts should be addressed to:

Mrs B Byrne
Institute for Computational Fluid Dynamics
Wolfson Building, Parks Road
Oxford OX1 3QD, UK.
Phone: 44-(0) 1865-273883
Fax: 44-(0) 1865-273839


From: Sizwe Mabizela <>
Date: Mon, 15 Sep 2003 15:25:04 +0200
Subject: Conference in Cape Town on Topology and its Applications

Dear NA-netter

I wish to draw your attention to the announcement of the "19th Summer
Conference on Topology and its Applications" which will be held at the
University of Cape Town from 5 July to 9 July 2004. I particularly wish
to draw your attention to a special session on Topology in Analysis and
Geometry which will be held as part of this conference. Full details are
obtainable on the following website:
This is an excellent opportunity for your to visit the magnificent
southern tip of the African continent!

Sizwe Mabizela.


From: J. Levesley <>
Date: Tue, 16 Sep 2003 09:19:39 +0100 (BST)
Subject: Will Light Memorial Conference

First Announcement of the

Will Light Memorial Conference
University of Leicester, UK,
December 18-19, 2003

18th: Jeremy Levesley (Leicester, UK)
19th: Rick Beatson (Canterbury, NZ)
Rob Brownlee (Leicester, UK)
Manfred von Golitschek (Wurzburg, Germany)
Tim Goodman (Dundee, UK)
George Phillips (St Andrews, UK)
Robert Schaback (Gottingen, Germany)

For further information see the conference website at

The conference is sponsored by the London Mathematical Society and
the Centre for Mathematical Modelling at the University of Leicester.


From: Knut Andreas Lie <>
Date: Mon, 15 Sep 2003 10:10:28 +0200
Subject: Winter School in Norway on Adaptive Methods for PDEs

Dr. Holms Hotel, Geilo, Norway, March 7-12, 2004

Dr. Mats G. Larson (Chalmers University of Technology, Sweden):
Introduction to adaptive methods, in particular:
* finite element methods
* adaptivity
* error estimation

Prof. Rolf Rannacher (University of Heidelberg, Germany):
Goal-oriented error estimation and adaptive methods and its
practical use for numerical simulation and optimization in
fluid dynamics, structural mechanics, chemistry and astrophysics.


Knut-Andreas Lie (SINTEF)
Aslak Tveito (Simula Research Laboratory)

Research Council of Norway through the BeMatA project

Please send your registration as soon as possible to
Ms. Mariann Fjeldstad, SINTEF Applied Mathematics,
P.O. Box 124 Blindern, N-0314 Oslo, Norway.
Phone: (+47) 22 06 76 54
Fax: (+47) 22 06 73 50

Please supply name and contact information (address, phone, fax, email).
Deadline for registration is January 15 2004.


From: Catherine Larkins <>
Date: Mon, 15 Sep 2003 09:50:20 -0600
Subject: Faculty Positions at University of Colorado, Boulder

- Assistant, Associate or full Professor - We invite applications for
a tenured or tenure-track position, to begin in August 2004 or
sooner. Applicants at all levels will be considered, but PhD
required. Strong preference will be given to candidates with
demonstrated research expertise in applied analysis, especially as
analysis intersects the other research areas of the department.
These areas of research include computational mathematics, dynamical
systems, nonlinear waves, analysis of differential equations,
physical applied mathematics, statistics and applied probability, and
mathematical biology. Further information about the department can
be found at: http// Excellence in teaching is also
expected. Applicants should send a letter of application, a current
curriculum vitae, a statement of research interests, and an AMS cover
sheet (see http// ) to Chair of Analysis
Committee, Department of Applied Mathematics, University of Colorado,
Boulder, CO, 80309-0526. Three letters of recommendation should
also be sent to this address. This is a continuation of an existing
search, which will end when the position is filled. The University
also be sent to this address. This is a continuation of an existing
search, which will end when the position is filled. The University
of Colorado at Boulder is committed to diversity and equality in
education and employment.

- Instructor - We invite applications for a non-tenure-track
instructorship, to begin in August, 2004. This is initially a
one-year position, typically renewed for up to two additional years
based on satisfactory performance in teaching and research.
Preference will be given to candidates with research interests in
applied analysis or statistics/probability, but strong candidates in
other areas of applied mathematics will also be considered. Areas of
research expertise in the department include computational
mathematics, dynamical systems, nonlinear waves, analysis of
differential equations, physical applied mathematics, statistics and
applied probability, and mathematical biology. Further information
about the department can be found at: http//
Applicants should send a letter of application, a current curriculum
vitae, a statement of research interests, an AMS cover sheet (see
http// ) to Chair of Instructor Committee,
Department of Applied Mathematics, University of Colorado, Boulder,
CO 80309-0526. Three letters of recommendation should also be sent
to this address. One of these letters should address the applicant's
teaching qualifications. This is a continuation of an existing
search, which will end when the position is filled. The University of
Colorado at Boulder is an equal opportunity/nondiscrimination


From: Luca Formaggia <>
Date: Mon, 15 Sep 2003 19:27:21 +0200
Subject: PhD Position at Politecnico di Milano, Italy

PhD position open at Politecnico di Milano (Italy)

Please note: this positions is open only to Citizens of EU countries or
associated states (ITALY EXCLUDED). Please, do not send your application
if you are not eligible, if in doubt please consult

The position is within the frame of the EU funded project HaeMOdel

Short description: Position leading to a PhD in Applied Mathematics (3
years). The applicant should hold a degree in Mathematics or Engineering
or Physics or related subjects and have a good attitude towards
numerical methods and scientific computation. He/she will be involved in
the research activity of the HaeMOdel group, aimed at the development of
efficent tools for the numerical simulation of the human cardiovascular
system and will participate to the training activities of the project.
Skill in numerical techniques and programming is welcomed.

Applicants should write (possibly by e-mail) before October the 7th 2003 to:

Alessandro Veneziani Mox Mathematics Department Via Bonardi,9 20133

Please send a detailed CV and provide also the name of one person that
might be contacted for references.

In accordance to Italian law, the position is subject to the succesful
participation to an examination that will be held at Politecnico di
Milano, Milan, Italy (expected date in November 2003); aimed to
assessing the candidate capabilities. The examination may be given in
Italian or in English.


From: Igor Boglaev <>
Date: Thu, 18 Sep 2003 09:31:41 +1200
Subject: Postdoctoral Position at Massey University

A position of Postdoctoral Fellow in scientific computation is available
at Institute of Fundamental Sciences,
Massey University (Palmerston North, New Zealand), in the group of Igor

The project concerns with monotone algorithms for solving nonlinear
singularly perturbed problems and involves high performance algorithms.
An applicant must have a good background in numerical partial
differential equations and parallel scientific computing.
An applicant should have a PhD in a related field.

The appointment will be for two years. Interested candidates should send
a CV to the address below, and should also arrange
to have three letters of recommendation sent.

Please send all correspondence to:
Human Resources Section
(attn Lynette Hensman)
Massey University, Private Bag 11 222,
Palmerston North, New Zealand

For more information, please visit hhtp://
or contact Igor Boglaev at


From: Joke Blom <>
Date: Fri, 19 Sep 2003 21:32:10 +0200 (CEST)
Subject: PhD Position at CWI, Amsterdam

CWI is an internationally renowned research institute in mathematics and
computer science, located in Amsterdam, The Netherlands. The focus is on
fundamental research problems, derived from societal needs.


The research group "Nonlinear PDEs: Analysis and Scientific Computing"
at CWI has a 4-year position for a

PhD Student in Numerical Mathematics/Scientific Computing

on the project "Mesoscale simulation paradigms in the Silicon Cell".
The research topic is the development of numerical methods and software for
an adaptive grid method for reaction-diffusion PDE systems with irregular,
dynamically changing, interfaces. The aim of the project is to compare this
continuum approach with a particle-based method (Lattice Boltzmann, developed
by a second PhD student) for problems with large spatial and temporal scales
such as arise, e.g., in the study of biological cells.


Candidates should have an MSc in (numerical) mathematics, scientific
computing or computational science, (computational) physics, computational
fluid dynamics, or a comparable expertise. Candidates are also expected to
have an active interest in life science applications.


Please send your application to the
Personnel Department CWI,
P.O. Box 94079,
1090 GB Amsterdam (NL),
with mentioning of "PhD Vacancy MAS1 - PDELB".
Your application should provide a concise description of your mathematical
education, your CV, examination marks and if possible a copy of your
M.Sc thesis. Students who are about to finish their M.Sc thesis are also
encouraged to apply. Consideration of applications will start immediately and
will continue until the position is filled.


More information about the project can be found at or contact
Joke Blom, e-mail:, tel.: +31-20-5924263, or
Prof. dr. J.G. Verwer,, +31-20-5924095.


From: Science Direct <>
Date: Mon, 15 Sep 2003 12:28:49 -0400 (EDT)
Subject: Contents, Linear Algebra and its Applications

Linear Algebra and its Applications Volume 373, Pages 1-332 (1 November 2003)
Combinatorial Matrix Theory Conference (Pohang, 2002)
Pohang, South Korea, 14 January - 17 January 2002
Edited by S.-G. Hwang, A.R. Kraeuter, B.L. Shader and J.-Y. Shao


Special Issue on the Combinatorial Matrix Theory Conference, Pages 1-3
Suk-Geun HwangArnold R. KrauterBryan L. ShaderJia-Yu Shao

A note on limit points for algebraic connectivity, Pages 5-11
Steve Kirkland

Matrix completion problems for pairs of related classes of matrices, Pages 13-29
Leslie Hogben

On graphs with algebraic connectivity equal to minimum edge density, Pages 31-50
Shaun M. Fallat, Steve Kirkland and Sukanta Pati

Linear preservers and diagonal hypergraphs, Pages 51-65
Richard A. Brualdi

A note on k-primitive directed graphs, Pages 67-74
LeRoy B. Beasley and Steve Kirkland

The linear algebra of the k-Fibonacci matrix, Pages 75-87
Gwang-Yeon Lee and Jin-Soo Kim

Matrices determined by a linear recurrence relation among entries, Pages 89-99
Gi-Sang Cheon, Suk-Geun Hwang, Seog-Hoon Rim and Seok-Zun Song

Linear criteria for lifting automorphisms of elementary abelian regular
coverings, Pages 101-119
Shao-Fei Du, Jin Ho Kwak and Ming-Yao Xu

An approach to solving Ak=J-I, Pages 121-142
Yaokun Wu and Qiao Li

Basic matrices, Pages 143-151
Miroslav Fiedler

A lower bound on the maximum permanent in [Lambda]nk, Pages 153-167
Ian M. Wanless

Generalized exponents of boolean matrices, Pages 169-182
Bolian Liu

On k-hypertournament matrices, Pages 183-195
Youngmee Koh and Sangwook Ree

Extremes of permanents of (0,1)-matrices, Pages 197-210
Seok-Zun Song, Suk-Geun Hwang, Seog-Hoon Rim and Gi-Sang Cheon

Sparse orthogonal matrices, Pages 211-222
Gi-Sang Cheon, Suk-Geun Hwang, Seog-Hoon Rim, Bryan L. Shader and Seok-Zun Song

Number of nonzero entries of S2NS matrices and matrices with signed generalized
inverses, Pages 223-239
Jia-Yu Shao, Jin-Ling He and Hai-Ying Shan

Which graphs are determined by their spectrum?, Pages 241-272
Edwin R. van Dam and Willem H. Haemers

Asymptotic enumeration of 0-1 matrices with equal row sums and equal column
sums, Pages 273-287
Brendan D. McKay and Xiaoji Wang

Factorizations of matrices over semirings, Pages 289-296
Han Hyuk Cho and Suh-Ryung Kim

The maximal determinant and subdeterminants of +/-1 matrices, Pages 297-310
Jennifer Seberry, Tianbing Xia, Christos Koukouvinos and Marilena Mitrouli

Inverse eigenvalue problems and lists of multiplicities of eigenvalues for
matrices whose graph is a tree: the case of generalized stars and double
generalized stars, Pages 311-330
Charles R. Johnson, Antonio Leal Duarte and Carlos M. Saiago


End of NA Digest