# sspt01.f - Man Page

TESTING/LIN/sspt01.f

## Synopsis

### Functions/Subroutines

subroutine sspt01 (uplo, n, a, afac, ipiv, c, ldc, rwork, resid)
SSPT01

## Function/Subroutine Documentation

### subroutine sspt01 (character uplo, integer n, real, dimension( * ) a, real, dimension( * ) afac, integer, dimension( * ) ipiv, real, dimension( ldc, * ) c, integer ldc, real, dimension( * ) rwork, real resid)

SSPT01

Purpose:

``` SSPT01 reconstructs a symmetric indefinite packed matrix A from its
block L*D*L' or U*D*U' factorization and computes the residual
norm( C - A ) / ( N * norm(A) * EPS ),
where C is the reconstructed matrix and EPS is the machine epsilon.```
Parameters

UPLO

```          UPLO is CHARACTER*1
Specifies whether the upper or lower triangular part of the
symmetric matrix A is stored:
= 'U':  Upper triangular
= 'L':  Lower triangular```

N

```          N is INTEGER
The number of rows and columns of the matrix A.  N >= 0.```

A

```          A is REAL array, dimension (N*(N+1)/2)
The original symmetric matrix A, stored as a packed
triangular matrix.```

AFAC

```          AFAC is REAL array, dimension (N*(N+1)/2)
The factored form of the matrix A, stored as a packed
triangular matrix.  AFAC contains the block diagonal matrix D
and the multipliers used to obtain the factor L or U from the
block L*D*L' or U*D*U' factorization as computed by SSPTRF.```

IPIV

```          IPIV is INTEGER array, dimension (N)
The pivot indices from SSPTRF.```

C

`          C is REAL array, dimension (LDC,N)`

LDC

```          LDC is INTEGER
The leading dimension of the array C.  LDC >= max(1,N).```

RWORK

`          RWORK is REAL array, dimension (N)`

RESID

```          RESID is REAL
If UPLO = 'L', norm(L*D*L' - A) / ( N * norm(A) * EPS )
If UPLO = 'U', norm(U*D*U' - A) / ( N * norm(A) * EPS )```
Author

Univ. of Tennessee

Univ. of California Berkeley

NAG Ltd.

Definition at line 109 of file sspt01.f.

## Author

Generated automatically by Doxygen for LAPACK from the source code.

## Referenced By

The man page sspt01(3) is an alias of sspt01.f(3).

Tue Nov 28 2023 12:08:43 Version 3.12.0 LAPACK