NA Digest Sunday, May 2, 2004 Volume 04 : Issue 18

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

Submissions for NA Digest:

Mail to

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


From: Joe Grcar <>
Date: Fri, 30 Apr 2004 16:51:14 -0700
Subject: Phil Colella Elected to NAS

Congratulations to Phil Colella on being elected
to the National Academy of Science. Phil has made
fundamental contributions to algorithms and software
for solving all manner of partial differential equations.

Joseph Grcar
Lawrence Berkeley National Laboratory


From: Kaj Madsen <>
Date: Sat, 1 May 2004 19:24:23 +0200 (METDST)
Subject: Cleve Moler Honorary Doctor

Cleve Moler Honorary Doctor
In the presence of Her Majesty, Queen Margrethe of Denmark,
Cleve Moler has been appointed Honorary Doctor (doctor technices,
honoris causa) at the Technical University of Denmark.
The degree was given at a ceremony on April 30.

It is motivated by Cleve's outstanding contributions in
research and software development, and especially for his
ground-breaking work with the software environment Matlab, which
has revolutionized the art of numerical computing, for the benefit
of scientists and engineers in a huge number of subjects throughout
the world.

Kaj Madsen
Informatics and Mathematical Modelling
Technical University of Denmark


From: Jane Cullum <>
Date: Thu, 29 Apr 2004 17:00:02 -0600 (MDT)
Subject: Change of Address for Jane Cullum

As of April 30, 2004 my new mailing address is

P.O. Box 626
Hawley, PA 18428-0626

and although my LANL mail should still be forwarded to me,
please use my AOL address

Jane Cullum


From: Reinhard Nabben <nabben@math.TU-Berlin.DE>
Date: Tue, 27 Apr 2004 21:58:30 +0200 (CEST)
Subject: Change of Address for Reinhard Nabben

My new address is:

Prof. Dr. Reinhard Nabben
Technische Universitaet Berlin
Institut fuer Mathematik, MA 3-3
Strasse des 17. Juni 136
D-10623 Berlin

Tel.: 49-(0)30-314-29291

Office: Math. Building Room 372


From: Joel Malard <>
Date: Sat, 01 May 2004 15:31:15 -0700
Subject: ACRE, Parallel Covariance Component Estimation Code

A couple of people have asked recently for a copy of the parallel
(restricted/residual) maximum (REML/ML) gradient algorithms from the paper:

J.M. Malard, "Parallel Restricted Maximum Likelihood Estimation
for Linear models with a Dense Exogenous Matrix", Parallel
Computing, 28, pp343-53, 2002.

The code has been upgraded to PETSc 2.2.0 and TAO 1.6 and is available
by sending an email at This software solves
covariance component estimation problems for linear models were the
residual vector comes from a normal distribution. The Cholesky
factorization of the covariance matrix is a sparse matrix. The problem
must be framed in the following form, which underlines that REML and ML
estimation can be viewed as a next step in complexity after solving
linear least squares.

Given a dense matrix A and a response vector b, find the
Best Linear Unbiased Estimator x and an upper triangular
matrix L such that Ax=b+e and the matrix LL' (L times
the transpose of L) is equal to the expected value of xx'.
Both L and A are assumed full rank.

It is customary in statistical modeling to split the matrix A into fixed
and random effects. The two formulations are equivalent but no script is
provided to do the conversion.

The purpose of this project was to demonstrate that linear estimation
algorithms such as REML can scale to a few hundred processors on a
distributed memory platform. The dll webpage contains some additional information.

If anyone needs an implementation of the REML Hessian matrix using the
forward differentiation mode, it has existed in the past, send an email
to the above address. My current priority is to allow for a singular
matrix L. Comments, bugs reports and suggestions are welcome at the same
email address.

With best regards,

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


From: David Bailey <>
Date: Sat, 1 May 2004 17:10:18 -0700
Subject: New Book on Experimental Mathematics

I am pleased to announce that the second volume of our work experimental
mathematics is now out, available from the publisher or from
The second volume gives numerous case studies of experimental mathematics in
action, ranging from sequences, series, products, integrals, Fourier series,
zeta functions, partitions, primes and polynomials. Some advanced numerical
techniques are also presented. The reference info is:

Authors: Jonathan M. Borwein, David H. Bailey and Roland Girgensohn

Title: "Experimentation in Mathematics: Computational Paths to Discovery"

Publisher: A. K. Peters (, 2004

If you would like a better taste of what these books are like, a 75-page
condensed version is available FREE at:

David H Bailey


From: Alessia Ritorto <>
Date: Wed, 28 Apr 2004 16:50:31 +0000
Subject: New Book, Arithmetic of Infinity

Sergeyev Ya.D., Arithmetic of infinity
Edizioni Orizzonti Meridionali, 2003, 112 pp., ISBN 88-89064-01-3

The book presents a new type of arithmetic that allows us to execute
arithmetical operations with infinite numbers in the same manner as we
are used to do with finite ones. The problem of infinity is considered in a
coherent way different from those proposed by Georg Cantor, Abraham Robinson,
and John Conway. However, the new approach does not contradict Cantor, it
complements his theory. The new viewpoint gives detailed answers to many
questions and paradoxes regarding infinite and infinitesimal quantities.
Particularly, simple applications of the new approach to limit theory,
measure theory, and set theory are given. The book opens new interesting
perspectives for the theory of computations too.

The book is mainly addressed to mathematicians, physicists, computer scientists,
and students. However, it is written in such a way that any person having a
high school education and who is interested in the foundations of these
sciences will be able to understand it easily.


From: Gene Golub <>
Date: Thu, 29 Apr 2004 10:10:18 -0700 (PDT)
Subject: Presentations from Two Recent Conferences


You can find the presentations from two recent conferences at the
following website.

Solution Methods for Saddle Point Systems in Computational Mechanics.

BIRS workshop on Model Reduction Problems and Matrix Methods.



From: Peter V. Coveney <>
Date: Fri, 30 Apr 2004 16:39:16 +0100
Subject: 2nd Annual RealityGrid Workshop

The 2nd Annual RealityGrid Workshop will take place
at the Royal Society of London on 15-16 June 2004.

The meeting will feature a mixture of talks by investigators within the
project and invited external
national and international speakers who are working in related fields. The
format will allow ample time for informal discussions
and the viewing of various demonstrations.

Details are now available at:

from where you can navigate to the programme and registration pages.

Please be sure to register if you wish to attend. Last year, we had to turn
some people away as the event proved to be very popular.

Peter Coveney


From: Vasco Brattka <>
Date: Fri, 30 Apr 2004 10:10:10 +0200 (CEST)
Subject: Conference in Germany on Real Numbers and Computers

*** Extended deadline ***

November 15-17, 2004
Schloss Dagstuhl, Germany

Web site:

Submission of manuscripts: May 15, 2004
Notification of acceptance: July 15, 2004
Revised Papers required by: September 1, 2004
Conference: November 15-17, 2004


The aim of the symposia on "Real Numbers and Computers" is to bring
together specialists from various research areas, all concerned with
problems related to computations based on real numbers. These
computations may use any number system implemented by a software
package or in hardware, including floating and fixed point, serial,
on line, continued fractions, exact, multiple precision, interval
and stochastic arithmetic.

Results are sought on both applied and fundamental questions.
Important topics discussed during these conferences include but are
not limited to:

* Foundation and properties of number systems
* Computability and complexity
* Formal aspects and automatic proof checking
* Links with number theory and automata theory
* Basic arithmetic operations
* Implementation of the standard and special functions
* Engineering of floating and fixed point algorithms
* Symbolic manipulation of numbers
* Accuracy and reliability for applications and industry
* Robust geometric algorithms and exact geometric computation
* Hardware design support and implementations

The conference will feature invited lectures and contributed talks.
Original research results and insightful analyses of current concerns
are solicited for submission. Survey and tutorial articles may be
suitable for submission if clearly identified as such.

Invited lectures:

Benno Fuchssteiner, Paderborn, Germany
Simon Plouffe, Montreal, Canada
Stefan Schirra, Magdeburg, Germany


Authors are invited to submit papers (Postscript or PDF) no later than May 15,
2004 by email to

The proceedings will be distributed at the conference.

A special issue on Real Numbers and Computers of the journal
Theoretical Informatics and Applications
is to follow the conference. Authors of the conference will be strongly
encouraged to propose expanded work for review and inclusion.

Steering Committee:

Jean-Claude Bajard, Montpellier, France
Jean-Marie Chesneaux, Paris, France
Marc Daumas, Lyon, France
Christiane Frougny, Paris, France
Peter Kornerup, Odense, Denmark (Chair)
Dominique Michelucci, Dijon, France
Jean-Michel Muller, Lyon, France


From: Willy Doerfler <>
Date: Mon, 26 Apr 2004 13:58:26 +0200
Subject: HP-FEM Summer School in Karlsruhe

Theory - Applications - Programming
20.09 - 24.09.2004, Karlsruhe, Germany

Ainsworth (Glasgow), Carstensen (Berlin),
Doerfler (Karlsruhe), Heuveline (Heidelberg),
Melenk (Reading), Monk (Delaware), Wieners (Karlsruhe)

Topics: We present the theoretical basis of the
hp-finite element method like interpolation
and error estimates, quadrature rules, solvers and
preconditioners. Applications in electrodynamics and
continuum mechanics will be considered. In the practical
part we explain some datastructure concepts and present
programming examples.
The audience is expected to have basic knowledge in theory
of elliptic partial differential equations and standard finite
element methods including some programming experience.

G. Alefeld, W. Doerfler, A. Rieder, C. Wieners (Karlsruhe)
supported by
o Centre for Scientific Computing and Mathematical Modelling (IWRMM)
in Karlsruhe
o Scientific Computing in Baden-Wuerttemberg (WIR)

More information is provided under

Contact: W. Doerfler, Inst. Appl. Math. 2, D-78128 Karlsruhe


From: Science Direct <>
Date: Tue, 27 Apr 2004 10:15:19 -0400 (EDT)
Subject: Contents, Linear Algebra and its Applications

Linear Algebra and its Applications Volume 381, (1 April 2004)

System theoretic based characterisation and computation of the least common multiple of a set of polynomials, Pages 1-23
Nicos Karcanias and Marilena Mitrouli

Conjugacy invariants of , Pages 25-35
B. Foreman

Inertia theorems for pairs of matrices, Pages 37-52
Cristina Ferreira and Fernando C. Silva

A property concerning the Hadamard powers of inverse M-matrices, Pages 53-60
Shencan Chen

A limit theorem for sets of stochastic matrices, Pages 61-76
Anne Condon and Michael Saks

On the powers of a vectorial matroid - an equivalence relation, Pages 77-95
Rosario Fernandes

A decreasing sequence of upper bounds on the largest Laplacian eigenvalue of a graph, Pages 97-116
Oscar Rojo and Hector Rojo

The k-edge connected subgraph problem I: Polytopes and critical extreme points, Pages 117-139
M. Didi Biha and A. R. Mahjoub

A least squares approach to reduce stable discrete linear systems preserving their stability, Pages 141-163
Sven Feldmann and Patrick Lang

Positive extension problems for a class of structured matrices, Pages 165-195
Vladimir Bolotnikov and Paul A. Smith

Spectral shorted matrices, Pages 197-217
Jorge Antezana, Gustavo Corach and Demetrio Stojanoff

Parametric extensions of Shannon inequality and its reverse one in Hilbert space operators, Pages 219-235
Takayuki Furuta

Some facets of the polytope of even permutation matrices, Pages 237-244
Jeffrey Hood and David Perkinson

Eventually nonnegative matrices are similar to seminonnegative matrices, Pages 245-258
Sarah Carnochan Naqvi and Judith J. McDonald

Numerical ranges of unbounded operators arising in quantum physics, Pages 259-279
N. Bebiano, R. Lemos and J. da Providencia


From: Science Direct <>
Date: Tue, 27 Apr 2004 10:15:19 -0400 (EDT)
Subject: Contents, Journal of Approximation Theory

Orthonormal polynomials for generalized Freud-type weights and higher-order Hermite-Fejer interpolation polynomials, Pages 1-38
T. Kasuga and R. Sakai

(C,[alpha]) summability of Walsh-Kaczmarz-Fourier series, Pages 39-60
Peter Simon

Multivariate tight affine frames with a small number of generators, Pages 61-73
Maura Salvatori and P.M.Paolo M. Soardi

Function-valued Korovkin systems without quasiconcavity and set-valued Korovkin systems without convexity, Pages 74-82
Pedro Teran

Weierstrass' theorem with weights, Pages 83-107
Ana Portilla, Yamilet Quintana, J.M.Jose M. Rodriguez and Eva Touris

Kergin interpolation in Banach spaces, Pages 108-123
Lars Filipsson

Where do homogeneous polynomials on l1n attain their norm?, Pages 124-133
David Perez-Garcia and Ignacio Villanueva

Erratum to \"Multivariate approximating averages\" [J. Approx. Theory 125 (2003) 85-105], Page 134
E. Belinsky, F. Dai and Z. Ditzian


End of NA Digest