NA Digest Sunday, November 28, 1993 Volume 93 : Issue 44

Today's Editor:

Cleve Moler
The MathWorks, Inc.

Submissions for NA Digest:

Mail to

Information about NA-NET:

Mail to


From: manolo@titan (Manuel Galan Moreno)
Date: Mon, 22 Nov 93 12:43:13 GMT
Subject: LSQR, Paige and Saunders

I am interested in receiving updated documentation & references about
Paige & Saunders LSQR method for linear system resolution.
Thanks a lot.


From: P.E. Van Camp <>
Date: Wed, 24 Nov 93 12:46:42 +0100
Subject: Computing Eigenvalues of Large, Complex Matrices

I am looking for software to find some of the eigenvalues and -vectors
(the lowest 10-30) of a large (dimension 2000-5000 or more) complex
Hermitian matrix. The matrix is not sparse (altough it contains a number
of zero entries) and usually has some two- or threefold degenerate
I checked Netlib etc. and could'nt find anything useful.

Many thanks,

P.E. Van Camp
University of Antwerp - RUCA
Department of Physics
Research Group TSM
roenenborgerlaan 171
B-2020 Antwerpen, Belgium
Tel. +32/3/2180316
Fax +32/3/2180217
E-mail: vancamp@banruc60.bitnet


From: Nick Trefethen <>
Date: Sat, 27 Nov 93 11:57:35 -0500
Subject: Upstate NA Day

Final schedule

456 Engineering and Theory Center
Cornell University
Ithaca, NY
Thursday, December 9, 1993

On Thursday, December 9, the Cornell Computer Science Department and Cornell
Theory Center will host an informal get-together for numerical analysts.
Graduate students are especially encouraged to attend. For further information
contact Nick Trefethen (, 607-255-4222) or Cindy Robinson
(, 607-255-0985).

9:30 DIRK LAURIE, Potchefstroom University for C.H.E., South Africa
Imperfect periodizing of functions for numerical integration --
must all derivatives vanish?

10:00 CHARLIE VAN LOAN, Department of Computer Science, Cornell University
An eigenvalue problem from adaptive optics

11:00 TOBY DRISCOLL, Center for Applied Mathematics, Cornell University
Interactive Schwarz-Christoffel mapping in Matlab, and an
application in iterative linear algebra

11:30 SAUL TEUKOLSKY, Department of Physics, Cornell University
Numerical Recipes: physicists versus numerical analysts?

2:00 NILOUFER MACKEY, Department of Computer Science, SUNY Buffalo
Convergence of the quaternion-Jacobi method

2:30 STEVE VAVASIS, Department of Computer Science, Cornell University
An accelerated interior point method based on layered least-squares

3:00 THOMAS COLEMAN, Department of Computer Science, Cornell University
Linearly constrained optimization and projected preconditioned
conjugate gradients

4:00 ANNE TREFETHEN, Cornell Theory Center, Cornell University
Parallel Matlab and pipe Poiseuille pseudospectra

4:30 ARIEH ISERLES, University of Cambridge, England
Qualitative analysis of discretized ODEs


From: Ashok Gupta <>
Date: Mon, 22 Nov 93 08:04:01 +0000 (UTC)
Subject: General Purpose Parallel Computing

The British Computer Society
Parallel Processing Specialist Group
General Purpose Parallel Computing

A One Day Open Meeting with Invited and Contributed Papers
22 December 1993, University of Westminster, London, UK

Invited speakers :
Les Valiant, Harvard University
Bill McColl, PRG, University of Oxford, UK
David May, Inmos, UK

A key factor for the growth of parallel computing is the availability of port-
able software. To be portable, software must be written to a model of machine
performance with universal applicability. Software providers must be able to
provide programs whose performance will scale with machine and application size
according to agreed principles. This environment presupposes a model of paral-
lel performance, and one which will perform well for irregular as well as regu-
lar patterns of interaction. Adoption of a common model by machine architects,
algorithm & language designers and programmers is a precondition for general
purpose parallel computing.

Valiant's Bulk Synchronous Parallel (BSP) model provides a bridge between appli-
cation, language design and architecture for parallel computers. BSP is of the
same nature for parallel computing as the Von Neumann model is for sequential
computing. It forms the focus of a project for scalable performance parallel
architectures supporting architecture independent software. The model and its
implications for hardware and software design will be described in invited and
contributed talks.

The PPSG, founded in 1986, exists to foster development of parallel architec-
tures, languages and applications & to disseminate information on parallel pro-
cessing. Membership is completely open; you do not have to be a member of the
British Computer Society. For further information about the group contact ei-
ther of the following :

Chair : Mr. A. Gupta Membership Secretary: Dr. N. Tucker
Philips Research Labs, Crossoak Lane, Paradis Consultants, East Berriow,
Redhill, Surrey, RH1 5HA, UK North Hill, Nr. Launceston, Cornwall, PL15 7NL, UK


From: James F Blowey <>
Date: Mon, 22 Nov 1993 14:28:26 +0000 (GMT)
Subject: Northern England Universities Numerical Analysis Colloquium

The Northern England Universities Numerical Analysis Colloquium will
take place in Van Mildert College, University of Durham on Thursday
January 6 1994.

The invited speakers are:

Arieh Iserles (Cambridge) `Numerical analysis of delay differential
equations with variable delays'


Alistair Watson (Dundee) `Some robust methods for parameter estimation'

There will be time for a few contributed talks. Further information
may be obtained from


From: Yousef Saad <>
Date: Tue, 23 Nov 1993 08:42:23 -0600 (CST)
Subject: Domain-Based Parallelism and Problem Decomposition


On April 25 and 26, 1994, the Minnesota Supercomputer Institute will
host a Workshop on Domain-Based Parallelism and Problem Decomposition
Methods in Computational Science and Engineering, organized in
cooperation with the SIAM Activity Group on Supercomputing and the
Army High Performance Computing Research Center. The workshop will
bring together researchers from many scientific and engineering
disciplines who exploit `problem decomposition techniques' in various
continuous or discrete senses to develop parallel algorithms. The
invited speakers will take a broad perspective suitable for an
interdisciplinary audience of engineers, physical scientists, and
applied mathematicians. There will be no contributed papers.
However, a poster session will be organized for those participants who
wish to present their work. Attendence to the workshop is limited, so
registration is required. Early registration is recommended.

For information or to obtain a registration form, contact one of the
organizers, David Keyes (ICASE) or Youcef Saad (U. of Minnesota, CS)
or Donald Truhlar (U. of Minnesota, Chemistry), or contact the
Symposium Administrator, Michael Olesen, Supercomputer Institute, 1200
Washington Avenue South, Minneapolis, MN 55415, phone):(612) 624-1356,
fax: (612) 624-8861, e-mail:

Invited speakers & lecture topics:

ZLATKO BACIC (New York University): ``Bound States of Strongly
Coupled Multidimensional Molecular Hamiltonians by the Discrete
Variable Representation Approach''.

PETTER BJORSTAD (Institute for Informatics in Bergen, Norway):
``Domain Decomposition and Parallel Computing with Applications to
Oil Reservoir Simulation''.

XIAO-CHUAN CAI (University of Kentucky): ``Domain Decomposition
Methods for Nonsymmetric and Indefinite Systems of Equations''.

CHARBEL FARHAT (University of Colorado at Boulder): ``Tailoring DD
Methods for Fast Parallel Coarse Grid Solution, Multiple Right Hand
Sides, and Adaptive Refinement''.

ROLAND GLOWINSKI (University of Houston): ``Merging Domain
Decomposition and Domain Imbedding Methods for Solving Elliptic and
Time Dependent Problems''.

WILLIAM GROPP (Argonne National Laboratories): ``Parallel Software
for Domain Decomposition Algorithms''.

GEORGES JOLICARD (Universite de Franche-Comte): ``The Wave Operator
Theory as a Tool to Define Active Spaces and Simplified Dynamics in
Large Size Quantum Spaces''.

DANA KNOLL (Idaho National Engineering Laboratory): ``Newton
-Krylov-Schwarz Methods Applied to the Tokamak Edge Plasma Fluid

ANDREW LUMSDAINE (Notre Dame University): ``Domain Decomposition in
Space and Time: Accelerated Waveform Methods for Semiconductor
Device Simulation''.

ALFIO QUARTERONI (Universita degli Studi di Milano in Italy):
``Domain Decomposition Methods for Wave Propagation Problems''.

FRANCOIS-XAVIER ROUX (Office National d'Etudes et de Recherches
Aerospatiales in Chatillon, France): ``Parallel Implementation of a
Domain Decomposition Method for Non-Linear Elasticity Problems''.

DAVID SCHWENKE (NASA Ames Research Center): ``New Ideas for Problem
Decomposition Methods in Quantum Mechanical Reactive Scattering''.

ELLEN STECHEL (Sandia National Laboratories): ``Large-Scale
Electronic Structure Calculations for Very Large Systems''.

TAYFUN TEZDUYAR (University of Minnesota): ``Parallel Finite
Element Computation of Flow Problems''.

V. VENKATAKRISHNAN (Institute of Computer Applications in Science
and Engineering in Virginia): ``Parallel Implicit Methods for
Aerodynamic Applications on Unstructured Grids''.

H.-J. WERNER (Universitat Bielefeld): ``Electron Correlation in

JACOB WHITE (Massachusetts Institute of Technology): ``Aspects of
Multipole and FFT Accelerators for the Rapid Solution of
Three-Dimensional Potential Integral Equations, with Applications''.


From: Alan J. Laub <>
Date: Tue, 23 Nov 93 12:13:47 PST
Subject: Position at Univ. California, Santa Barbara

University of California, Santa Barbara
Electrical and Computer Engineering

Applications are invited for a tenure-track faculty
position in engineering computation, available
effective July 1, 1994. The level of the position
is open. It is expected that the successful
candidate will interact most strongly with faculty
in communications, control, and/or signal and image
processing, but significant interaction with researchers
in computer engineering, computer science, solid state,
and wave electronics are also possible and desirable.
Specific computational areas of interest include
partial differential equations, ordinary differential
equations, optimization, wavelets, and stochastic analysis.
Normally, completion of a doctorate is required at
the time of the appointment. Candidates should have
an established research reputation or outstanding
research potential, a record of or the ability to
attract external research funding, and a strong
commitment to teaching at the undergraduate and
graduate levels. Applicants should send their
resumes and the names and addresses of at least
four professional references to:
Faculty Search Committee
Dept. of Electrical and Computer Engineering
University of California
Santa Barbara, CA 93106-9560
Applications will be received until the position is
filled. UCSB is an Equal Opportunity/Affirmative
Action employer.


From: Bob Russell <>
Date: 23 Nov 93 18:53 -0800
Subject: Position at Simon Fraser University


The Department of Mathematics and Statistics of Simon Fraser University
invites applications for a tenure track position in Applied and
Computational Mathematics at the Assistant Professor level starting
September 1, 1994. We are looking for applicants with an interest in
non-linear applied mathematics; particularly those whose research focuses
on differential equations and their applications. Applicants will be
expected to have completed a PhD degree at the time of appointment and to
have demonstrated a strong teaching and research potential. (This position
is subject to final budgetary approval.)
Applications, including curriculum vitae, descriptive statements on
research plans and teaching activities should be sent, by 15 February 1994,
Dr. Katherine Heinrich, Chair
Department of Mathematics and Statistics
Simon Fraser University
Burnaby, BC V5A 1S6
Please arrange for three letters of reference to be sent directly from the
Simon Fraser University is committed to the principle of equity in
employment and offers equal employment opportunities to all qualified
applicants. In accordance with Canadian immigration requirements, this
advertisement is directed to those who are eligible at the time of the
application for employment in Canada.>


From: Iain Duff <>
Date: Thu, 25 Nov 93 10:43:33 GMT
Subject: Contents: IMA Journal of Numerical Analysis

IMA Journal of Numerical Analysis
Volume 14. Number 1
Jan. 1994

Simpson R B
Testing for effects of asymmetry and instability on preconditioned iterations
of conjugate gradient type.

Qin Sheng
Global error estimates for exponential splitting

Levesley J, Hough D M, and Chandler-Wilde S N
A Chebyshev collocation method for solving Symm's integral equation for
conformal mapping: a partial error analysis.

Mastroianni G and Monegato G
Nystrom interpolants based on the zeros of Legendre polynomials for a
non-compact integral operator equation.

Semper B
Locking in finite-element approximations to long thin extensible beams.

Gastaldi F and Gastaldi L
On a domain decomposition for the transport equation: theory and finite element

Wu W, Spence A, and Cliffe K A
Steady-state/Hopf mode interaction at a symmetry-breaking Takens-Bogdanov point.


End of NA Digest