1DLARRA(1)           LAPACK auxiliary routine (version 3.1)           DLARRA(1)
2
3
4

NAME

6       DLARRA - splitting points with threshold SPLTOL
7

SYNOPSIS

9       SUBROUTINE DLARRA( N, D, E, E2, SPLTOL, TNRM, NSPLIT, ISPLIT, INFO )
10
11           IMPLICIT       NONE
12
13           INTEGER        INFO, N, NSPLIT
14
15           DOUBLE         PRECISION SPLTOL, TNRM
16
17           INTEGER        ISPLIT( * )
18
19           DOUBLE         PRECISION D( * ), E( * ), E2( * )
20

PURPOSE

22       Compute  the  splitting  points with threshold SPLTOL.  DLARRA sets any
23       "small" off-diagonal elements to zero.
24
25

ARGUMENTS

27       N       (input) INTEGER
28               The order of the matrix. N > 0.
29
30       D       (input) DOUBLE PRECISION array, dimension (N)
31               On entry, the N diagonal elements of the tridiagonal matrix T.
32
33       E       (input/output) DOUBLE PRECISION array, dimension (N)
34               On entry, the first (N-1) entries contain the subdiagonal  ele‐
35               ments  of  the  tridiagonal matrix T; E(N) need not be set.  On
36               exit, the entries E( ISPLIT( I ) ), 1 <= I <= NSPLIT,  are  set
37               to zero, the other entries of E are untouched.
38
39       E2      (input/output) DOUBLE PRECISION array, dimension (N)
40               On  entry,  the  first (N-1) entries contain the SQUARES of the
41               subdiagonal elements of the tridiagonal matrix  T;  E2(N)  need
42               not  be set.  On exit, the entries E2( ISPLIT( I ) ), 1 <= I <=
43               NSPLIT, have been set to zero
44
45               SPLTOL (input) DOUBLE PRECISION The  threshold  for  splitting.
46               Two criteria can be used:
47               SPLTOL<0 : criterion based on absolute off-diagonal value
48               SPLTOL>0 : criterion that preserves relative accuracy
49
50               TNRM (input) DOUBLE PRECISION The norm of the matrix.
51
52       NSPLIT  (output) INTEGER
53               The number of blocks T splits into. 1 <= NSPLIT <= N.
54
55       ISPLIT  (output) INTEGER array, dimension (N)
56               The  splitting  points,  at which T breaks up into blocks.  The
57               first block consists of rows/columns 1 to ISPLIT(1), the second
58               of  rows/columns  ISPLIT(1)+1  through ISPLIT(2), etc., and the
59               NSPLIT-th consists of rows/columns  ISPLIT(NSPLIT-1)+1  through
60               ISPLIT(NSPLIT)=N.
61
62       INFO    (output) INTEGER
63               = 0:  successful exit
64

FURTHER DETAILS

66       Based on contributions by
67          Beresford Parlett, University of California, Berkeley, USA
68          Jim Demmel, University of California, Berkeley, USA
69          Inderjit Dhillon, University of Texas, Austin, USA
70          Osni Marques, LBNL/NERSC, USA
71          Christof Voemel, University of California, Berkeley, USA
72
73
74
75
76 LAPACK auxiliary routine (versionNo3v.e1m)ber 2006                       DLARRA(1)
Impressum