NA Digest Monday, June 15, 1998 Volume 98 : Issue 22

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

Submissions for NA Digest:

Mail to

Information about NA-NET:

Mail to

URL for the World Wide Web:

From: Victoria Hernandez and Jorge Estrada <>
Date: Wed, 10 Jun 1998 11:16:50 -0500 (CUBA)
Subject: Euclidean Distance from a Point to a Conic


We are researchers at Institute of Cybernetics, Mathematics and Physics,
(ICIMAF) from Havana, Cuba. We have been working out a new algorithm
to compute the euclidean distance from a point to an arbitrary conic.
Now it is ready and the algorithm works fairly well even when the
points are not near to the conic.
Therefore we would like to know if anyone in the NA community can
figure out (real) situations where it is necessary to compute the
euclidean distance from a point to a curve, but you may not assume
that the point is very close to the curve (unlike the usual fitting
problems). Thank you in advance for your kind attention!
Sincerely yours,
Victoria Hernandez and Jorge Estrada


From: Amie Casper <>
Date: Wed, 10 Jun 1998 10:34:43 -0400
Subject: New Book, The Geometry Toolbox for Graphics and Modeling

The Geometry Toolbox for Graphics and Modeling
Gerald Farin, Dianne Hansford

Now available from A K Peters, Ltd.:

288 pages; Hardcover; ISBN 1-56881-074-1
List price: $48.00

The Geometry Toolbox takes a novel and particularly visual approach to
teaching the basic concepts of two- and three-dimensional geometry. It
explains the geometry essential for today's computer-modeling,
computer-graphics, and animation systems. While the basic theory is
completely covered, the emphasis of the book is not on abstract proofs but
rather on examples and algorithms. The Geometry Toolbox is the ideal text
for professionals who want to get acquainted with the latest geometric
tools. The two chapters on basic curves and surfaces are an ideal stepping
stone into the world of graphics and modeling. It is also a unique
textbook for a modern introduction to linear algebra and matrix theory.

A K Peters, Ltd., 63 South Ave., Natick, MA 01760
Ph: 508/655-9933, Fax: 508/655-5847


From: Alan Hindmarsh <>
Date: Fri, 12 Jun 1998 10:23:12 -0700 (PDT)
Subject: New Book, Incompressible Flow and the Finite Element Method

New Book Announcement:
Incompressible Flow and the Finite Element Method, by

Philip Gresho (UC Lawrence Livermore National Laboratory) and
Robert Sani (University of Colorado), in collaboration with
Michael Engelman (Fluid Dynamics International),
1998, John Wiley & Sons, 1040 pages, $295.

This comprehensive reference work deals with the important details in
the application of the finite element method to incompressible flows.
It addresses the theoretical background and the detailed development
of appropriate numerical methods applied to the solution of a wide
range of incompressible flows. A lucid explanation of the FEM is
combined with extensive coverage of the prototypical fluid mechanics
equation: the advection-diffusion equation. For both this equation
and the equations of principal interest -- the Navier-Stokes equations
-- detailed discussion of both the continuous and discrete equations
is presented, as well as explanations of how to properly march the
time-dependent equations using smart implicit methods. Presentations
of both ODE and DAE solution methods have benefitted from the
significant help of Alan Hindmarsh (LLNL). Boundary and initial
conditions, so important in applications, are carefully described and
discussed, including well-posedness. The important role played by the
pressure, so confusing in the past, is carefully explained.

This book explains and emphasizes consistency in six areas:
* consistent mass matrix
* consistent pressure Poisson equation
* consistent penalty methods
* consistent normal direction
* consistent heat flux
* consistent forces

Fully indexed and referenced, this book is an essential refrence tool
for all researchers, students, and applied scientists in incompressible
fluid mechanics.


From: Willi Schoenauer <>
Date: Mon, 15 Jun 98 11:13:40 +0100
Subject: Internet Book Scientific Supercomputing

The first (hand-written) Internet book on Scientific Supercomputing
now has been finished:

Scientific Supercomputing:

Architecture and Use of Shared and Distributed Memory
Parallel Computers


1. Introduction
2. Prototypes with their Bottlenecks
3. Arithmetic Operations and Memory Bandwidth
4. The CRAY T90
5. The Fujitsu VPP300 (and VPP700)
6. The IBM RS/6000 SP
7. The CRAY T3E and DEC Alpha Processors
8. Performance Analysis
9. Basic Considerations Concerning Data Structures
10. Fortran, Autovectorization and Autoparallelization, Programming Models
11. Recurrences
12. Matrix Operations
13. Systems of Linear Equations on Vector Processors
14. Remarks to the Parallelization of the Gauss Algorithm
15. Tridiagonal Linear Systems
16. Iterative Solution of Linear Systems of Equations
17. Explaining the Gap Between Theoretical Peak Performance and Real
Performance for Supercomputer Architectures
18. References

The book is available under the URL

In the Addendum further architectures and algorithms will be
presented (open-ended).

Prof. Dr. Willi Schoenauer - Universitaet Karlsruhe -
Rechenzentrum - Zirkel 2 - 76128 Karlsruhe -
Voice: (+49) (0)721-6083754 - Fax (+49) (0)721-32550 -


From: Roger Haagmans <>
Date: Mon, 8 Jun 1998 09:59:33 +0100
Subject: School on Wavelets in the Geosciences


Supported by:
International Association of Geodesy
Netherlands Geodetic Commission
Sub-Faculty of Geodetic Engineering DUT
Delft Institute for Earth-Oriented Space Research

October, 4 - 9, 1998
At the Delft Institute for Earth-Oriented Space Research
Faculty of Civil Engineering and Geosciences
Delft University of Technology
Thijsseweg 11
Delft, The Netherlands

Latest information
This text contain the latest information on the School on Wavelets in the
Geosciences. The final team of lecturers consists of Willi Freeden, Matthias
Holschneider and Wim Sweldens. More detailed information on the contents of
the lectures as provided by the lecturers is available and registration
information is extended. The first registrations show an interesting
diversity in scientific disciplines among the participants: so don't
hesitate to register but beware of the registration deadlines and the
limited number of places.

Objectives of the school

The basic objective of the School is provide the necessary information to
understand the potential and limitations of the application of wavelets in
the geosciences. This includes:
* the mathematical representation in one and more dimensions like on the
* the properties as compared to Fourier techniques
* the signal representation and analysis ability
* the use of operators in terms of wavelets
* gaining experiences with wavelets using examples from geosciences in
computer exercises


The course will last for six days and contains three major subjects. Every
subject will be covered in two days time. All topics will be supported by
practical exercises on the computer with examples from geodynamics,
topography representation, gravity field modelling etc.
The lectures and subjects are:
1. Dr.Matthias Holschneider, Laboratoire de G=E9omagn=E9tisme, Institut de
Physique du Globe de Paris, France;

compression and filtering with wavelets.
2. Dr. Wim Sweldens, Mathematical Sciences Research Centre, Lucent
Technologies Bell Laboratories, Murray Hill, NJ, USA

3. Prof. Dr. Willi Freeden, Geomathematics Group. Department of
Mathematics, University of Kaiserslautern, Germany


Scientific Committee
* Prof. Dr. Willi Freeden, Geomathematics Group, Department of
Mathematics, University of Kaiserslautern, Germany.
* Prof. Dr. Roland Klees, Delft Institute for Earth-Oriented Space
Research (DEOS), Faculty of Civil Engineering and Geosciences, Delft
University of Technology, The Netherlands.

Local organising committee and school secretariat address:
Dipl.-Math. Michael Bayer (Geomathematics Group), Dr. Ir.Martin van Gelderen,
Ir. Roger Haagmans and ms. Wil Coops-Luijten (DEOS)
Secretariat IAG School on Wavelets in the Geosciences
c/o Wil Coops-Luijten
DEOS, Faculty of Civil Engineering and Geosciences,
Delft University of Technology
Thijsseweg 11, NL-2629 JA Delft, The Netherlands
Telephone: +31 15 2783289
Fax: +31 15 2783711

More detailed information can be obtained from:


From: Jennifer Mueller <>
Date: Thu, 11 Jun 1998 16:22:28 -0400
Subject: Workshop on Nonlinear and Improperly Posed Problems

Workshop deadline extended

The deadline for abstracts and registration for the

International Workshop on Nonlinear and Improperly Posed Problems
August 15-13, 1998
International Hotel Sapanaca
Kocaeli - Turkey

has been extended to June 30, 1998. For details please see the website


From: Peggy Sicilano <>
Date: Fri, 12 Jun 1998 10:07:27 -0400 (EDT)
Subject: IMACS Conference on Monte Carlo Methods

June, 7-11, 1999, Varna BULGARIA

The Second IMACS Seminar on Monte Carlo methods will be held in Varna,
Bulgaria, International Conference Hall "Joliot Curie" June 7-11, 1999. The
purpose of the seminar is to provide a forum for the presentation of recent
advances in the analysis, implementation and applications of Monte Carlo
simulation techniques and, in particular to stimulate the exchange of
information between specialists in these areas.

The topics should cover both theoretical developments: (random numbers,
quasi-MC methods, numerical MC methods statistical analysis, variance
reduction, perturbation techniques, MC error analysis)
and application fields: (particle transport, reliability analysis,
quantum mechanics, statistical physics, stimulation of random processes
and fields).

Calendar (for submission and selection of contributions):
Deadline for submission of abstracts (one page stating clearly the originality and importance of the results): September 1998
Notification of acceptance (on basis of the abstract): November 1998
Deadline for submission of full papers: January 1999
Notification of acceptance of full papers: April 1999

For further information:
or send email to:


From: Kenneth Appel <>
Date: Mon, 8 Jun 1998 07:44:21 -0700
Subject: Position at University of New Hampshire

The Department of Mathematics of the University of New Hampshire invites
applications for a temporary position (1-3 years).

Preference will be given to candidates whose research areas
include control theory, dynamical systems, or numerical analysis.

Candidates are expected to have a PhD in Mathematics by August 1998 and
demonstrated research ability. Strong commitment to teaching is expected.
The teaching load is expected to be two courses per semester, at least one a
large lecture.
The positions will begin in the fall semester of 1998. Review of applications
will begin on June 20, 1998, but applications will be considered until the
positions are filled.

Resume, e-mail address if possible, and three letters of recommendation should
be sent to Kenneth Appel, Chair, Department of Mathematics, University of
New Hampshire, Durham, N.H., 03824 ( The University of
New Hampshire is an EO/AA employer.


From: Ake Bjorck <>
Date: Fri, 12 Jun 1998 09:46:49 +0200 (MET DST)
Subject: Contents, BIT

CONTENTS BIT Volume 38, No. 3 (September 1998)
ISSN 0006-3835

For information to contributers and about subscriptions see
Available on the Internet as from 1998---contact
Swets & Zeitlinger at

Half-explicit Runge-Kutta methods with explicit stages for
differential-algbraic systems of index 2
M. Arnold, pp. 415--438

Pseudo-symplectic Runge-Kutta methods
A. Aubry and P. Chartier, pp. 439--461

On the variable stepsize approximations of a Cauchy problem for
the evolution equation
N. Yu. Bakaev, pp. 462--485

On the sensitivity of the LU factorization
X.-W. Chang and C. C. Paige, pp. 486--501

Relative perturbation results for eigenvalues and eigenvectors
of diagonalisable matrices
S. C. Eisenstat and I. C. F. Ipsen, pp. 502--509

Relaxed functional iteration techniques for the numerical solution
of M/G/1 type Markov chains
P. Favati and B. Meini, pp. 510--526

Minimum residual methods for augmented systems
B. Fisher, A. Ramage, D. Silvester, and A. J. Wathen, pp. 527--543

A framework for polynomial preconditioners based on fast transforms I:
S. Holmgren and K. Otto, pp. 544--559

A total least squares method for Toeplitz systems of equations
J. Kamm and J. G. Nagy, pp. 560--582

Computing zeros of analytic mappings: A logarithmic residue approach
P. Kravanja, R. Cools, and A. Haegemans, pp. 583--596

A numerical method for the Benjamin--Ono equation
V. Thom\'ee and A. S. Vasudeva Murthy, pp. 597--611


A remark on ``a note on constructing a symmetric matrix with specified
diagonal entries and eigenvalues''
K. D. Ikramov, p. 612

A stability result for general linear methods with characteristic
function having real poles only
G.-Y. Psihoyios and J. R. Cash, pp. 613--618


From: Hans Schneider <>
Date: Sun, 7 Jun 1998 18:05:50 -0500 (CDT)
Subject: Contents, Linear Algebra and its Applications

Journal : Linear Algebra and Its Applications
Volume issue : 277/1-3
Year : 1998
Anticipated Publication Date: 05-JUN-98

pp. 1-9
Convergence of inhomogeneous products of matrices and coefficients of ergodicity
UG Rothblum, J Hartfiel

pp. 11-31
Completely positive matrices with a book-graph
F Barioli

pp. 33-39
Local inversion of matrices with sparse inverses
CR Johnson, M Lundquist

pp. 41-48
On the matrix equation A^l+A^l+k=J_n
Y Wu, Q Li

pp. 49-56
Falsity of wang's conjecture on stars
CS Karuppan Chetty

pp. 57-61
Matrix characterization of MDS linear codes over modules
XD Dong, CB Soh

pp. 63-81
A balanced canonical form for discrete-time minimal systems using characteristic
J Hoffmann

pp. 83-95
Matrices of zeros and ones with the maximum jump number
BO Cheng

pp. 97-125
On lie gradings II
M Havlicek, M Havlicek

pp. 127-134
An index theorem for the product of linear relations
RW Cross

pp. 135-142
A note on the hyperbolic singular value decomposition
BC Levy

pp. 143-148
The image of the adjoint mapping
DW Robinson

pp. 149-185
Primes in several classes of the positive matrices
G Picci, JH Van Schuppen

pp. 187-198
Eigenvalue location for nonnegative and Z-matrices
SM Fallat

pp. 199-236
The truncated hamburger matrix moment problems in the nondegenerate and
degenerate cases, and matrix continued fractions
GN Chen

pp. 237-238
On a determinant result of I. Olkin
M Marcus

pp. 239-251
Matrix sandwich problems
MC Golumbic

pp. 253-269
On the number of invariant polynomails of the product of matrices with
prescribed similarity classes
YL Zhang

pp. 271-289
Multimatriods IV. Chain-group representions
A Bouchet

pp. 291-298
Norms of sampling operators
P Zizler

pp. 299-311
Reverse order law for reflexive generalized inverses of products of matrices
AR De Pierro, Musheng Wei

pp. 313-336
Time-varying discrete Riccati equation in terms of Ben Artzi - Gohberg dichotomy
V Ionescu

pp. 337-356
On matrices satisfying a maximum principle with respect to a cone
MR Weber

pp. 357-360
Review of"Nonnegative Matrices and Applications" by R.B. Bapat and T.E.S.
S Kirkland


From: Mehrmann Volker <volker.mehrmann@Mathematik.TU-Chemnitz.DE>
Date: Wed, 10 Jun 1998 11:27:40 +0200 (METDST)
Subject: Contents, Linear Algebra and its Applications


Journal : LAA Linear Algebra and Its Applications
Volume/issue : 275-276
Year : 1998
Journal no. : 7738
ISSN : 0024-3795
Cover Date : May 15, 1998

pp. 1-1
[B. Shader, S. Kirkland, V. Mehrmann, G. Michler]

pp. 3-18
A max version of the Perron-Frobenius theorem
R.B. Bapat

pp. 19-47
The symplectic eigenvalue problem, the butterfly form, the SR algorithm
and the
Lanczos method
P. Benner, H. Fassbender

pp. 49-75
Compactly supported solutions of two-scale difference equations
L. Berg

pp. 77-106
Some comparison theorems for weak nonnegative splittings of bounded
J.-J. Climent, C. Perea

pp. 107-120
Practical task-oriented parallelism for Gaussian elimination in distributed
G. Cooperman

pp. 121-140
Extension of C*algebras and Moore-Penrose stability of sequences
V.D. Didenko, B. Silbermann

pp. 141-160
The role of formalism in the teaching of the theory of vectors spaces
J.-L. Dorier

pp. 161-177
Differentiable structure of the set of controllable (A,B)t-invariant
J. Ferrer, F. Puerta, X. Puerta

pp. 179-187
Some characterizations of symmetric inverse M-matrices
M. Fiedler

pp. 189-200
Spectral properties of rational matrix functions with nonnegative
K.-H. Foerster

pp. 201-224
Matrix rank based conditions for reachability/controllability of discrete
repetitive processes
K. Galkowski, E. Rogers

pp. 225-248
Representations of Toeplitz-plus-Hankel matrices using trigonometric
transformations with application to fast matrix-vector multiplication
G. Heinig, K. Rost

pp. 249-259
On cones and stability
D. Hershkowitz

pp. 261-279
Modifying the inertia of matrices arising in optimization
N.J. Higham

pp. 281-286
The Gauss-Huard algorithm and LU factorization
W. Hoffmann

pp. 287-314
On structured perturbation of Hermitian matrices
T. Hu, L. Qiu

pp. 315-325
An algorithm for nilpotent completions of partial Jordan matrices
C. Jordan, J.R. Torregrosa

pp. 327-347
Uniform convergence of a fast algorithm for cauchy singular intergral
P. Junghanss, U. Luther

pp. 349-357
Extreme nonnegative matrices
T.J. Laffey

pp. 359-379
Modal coupling in linear control systems using robust eigenstructure
N.K. Nichols, D.M. Littleboy

pp. 381-400
A Modified Block Newton iteration for approximating an invariant
subspace of a
symmetric matrix
R. Loesche, H. Schwetlick, G. Timmermann

pp. 401-415
A projection method for computing the minimum eigenvalue of a
symmetric positive definite Toeplitz matrix
W. Mackens, H. Voss

pp. 417-431
Spectral sensitivity of products of bidiagonals
B.N. Parlett

pp. 433-450
Fast and stable algorithms for discrete spherical Fourier transforms
D. Potts, G. Steidl, M. Tasche

pp. 451-470
Parallel codes for computing the numerical rank
G. Quintana-Orti, E.S. Quintana-Orti

pp. 471-493
Exceptional modules are three modules
C.M. Ringel

pp. 495-507
Digraph based determination of Jordan block size structure of singular
K. Roebenack, K.J. Reinschke

pp. 509-529
Restricted invariant factor assignment under state-feedback
A. Roca, I. Zaballa

pp. 531-536
A note on singular values of Cauchy-Toeplitz matrices
S. Roch, B. Silbermann

pp. 537-549
Inversion formulas and fast algorithms for Lowner-Vandermonde
K. Rost, Z. Vavrin

pp. 551-562
Matrix cyclization over complex polynomials
W. Schmale

pp. 563-577
On matrices for which norm bounds are attained
H. Schneider, H.F. Weinberger

pp. 579-594
Fiedler matrices and their factorization
J.L. Stuart, J. Weaver

pp. 595-615
Sensitivity analysis of the discrete-time algebraic Riccati equation
J.-G. Sun

pp. 617-626
On normal matrices of zeros and ones with fixed row sum
B.-Y. Wang, F. Zhang

pp. 627-629
[Conference report]]
[V. Mehrmann, H. Schneider]
pp. 631-632

Author Index Vol. 275-276


From: Maxim Larin <>
Date: Wed, 10 Jun 98 10:26:36 +0700
Subject: Contents, Siberian Journal of Numerical Mathematics

Siberian Journal of Numerical Mathematics

Editor-in-Chief: A.S.Alekseev
Institute of Comp. Math. and Mathem. Geophysics,
Russian Academy of Sciences, Siberian Branch,
Lavrentjeva 6, Novosibirsk 630090, RUSSIA

Managing Editor: V.A.Vasilenko
Institute of Comp. Math. and Mathem. Geophysics,
Russian Academy of Sciences, Siberian Branch,
Lavrentjeva 6, Novosibirsk 630090, RUSSIA

Scientific secretary: V.P.Il'in
Institute of Comp. Math. and Mathem. Geophysics,
Russian Academy of Sciences, Siberian Branch,
Lavrentjeva 6, Novosibirsk 630090, RUSSIA

Aims and Scope

The Publishing House of Siberian Division of Russian Academy of Sciences
from 1998 began to issue the new quarterly Siberian Journal on Numerical
Mathematics (SibJNM). The editorial board consists of the well-known
specialists in numerical mathematics and its applications from Siberia
and from all over Russia. The aim of this action is to demonstrate, in
concentrated form, to the Russian and International Mathematical Community
the latest most important investigations of Siberian numerical
mathematicians in various scientific and engineering fields. Indeed,
the cooperation with the numerical mathematicians from the other regions
of Russia and from other countries is welcome. Our Journal will be
distributed in Russia, in the Western and Eastern Europe, in the United
States and in Japan.

Main SibJNM Topics

- Numerical Mathematics
- Mathematical and Numerical Modeling
- Applied Informatics
- Automatization of Scientific and Applied Studies

Editorial Office

Editorial Office of SibJNM,
prospect Lavrentieva, 6,
630090, Novosibirsk, Russia.
Tel: (007) 3832-396545
Fax: (007) 3832-343783

Contents of Volume 1, 1998

Number 1, pp.1-98

A.S. Alekseev
On the scientific and educational potential of Siberia in the field of
computational and applied mathematics (foreword of the
Editor-in-Chief), pp.3-4

S.K. Godunov and V.M. Gordienko
The Krylov space and the Kalman equation (in Russian), pp.5-10

A.N. Gorban'
Generalized approximation theorem and computational capabilities of
neural networks (in Russian), pp.11-24

A.N. Konovalov
Conjugate-factorized models in mathematical physics
problems (in Russian), pp.25-58

L.D. Menikhes and V.P. Tanana
The finite-dimensional approximation for the Lavrent'ev method
(in Russian), pp.59-66

G.A. Mikhailov
New Monte Carlo methods for solving boundary value problems (and related
topics), pp.67-76

V.V. Ostapenko
Approximation of Hugoniot's conditions by explicit conservative
difference schemes for non-stationar shock waves
(in Russian), pp.77-87

S.B. Sorokin
Step-by-step inversion method for elasticity problems, pp.89-98

Number 2, pp.99-196

E.G. Bykova and V.V. Shaidurov
A nonuniform difference scheme with fourth order of accuracy
in a domain with smooth boundary (in Russian), pp.99-118

A.V. Voytishek
On the permissible class of interpolations for discrete-stochastic
procedures of global estimation of functions (in Russian),

V.M. Malbackov
On the statistical properties of the hydrodynamic models based
on solutions to the Boussinesq equations (in Russian),

R.V. Namm
On characterization of a limit point in the iterative
prox-regularization method (in Russian), pp.143-152

V.A. Rukavishnikov and A.Yu. Bespalov
On the h-p version of the finite element method
for one-dimensional boundary value problem with singularity
of a solution, pp.153-170

L.I. Skurin
An iterative-marching method for solving problems of fluid and
gas mechanics (in Russian), pp.171-182

T.A. Sushkevich, S.A. Strelkov, A.K. Kulikov, and
S.V. Maksakova
A model of the polarized radiation transfer in a planar layer
with interface of two media (in Russian), pp.183-194

Number 3 (forthcoming), pp.197-2??

I.A. Blatov
On the incomplete factorization for
fast Fourier transform method for the discrete Poisson equation
in a curvelinear boundary domain (in Russian).

A cascadic multigrid algorithm in finite
element method for the three-dimensional Dirichlet problem
(in Russian).

V.A.Debelov, A.M.Matsokin, and S.A.Upol'nikov
Plane subdivision and set operations (in Russian).

A.I. Zadorin
Numerical solution of the equation with a
small parameter and the point source on the infinite interval
(in Russian).

B.G. Mikhailenko and O.N. Soboleva
Absorbing boundary conditions for the elastic theory equations
(in Russian).

V.F. Raputa, A.I. Krylova, and G.A. Platov
An inverse problem of estimating the total emission
for the nonstationary boundary layer of the atmosphere (in Russian).

G.I. Shishkin
Grid approximations of singularly perturbed systems of parabolic
convection-diffusion equations with counterflow.

An additional information you may obtain at


End of NA Digest