CLAQR1 man page on Scientific
Printed from http://www.polarhome.com/service/man/?qf=CLAQR1&af=0&tf=2&of=Scientific
CLAQR1(1) LAPACK auxiliary routine (version 3.2) CLAQR1(1)
NAME
SYNOPSIS
SUBROUTINE CLAQR1( N, H, LDH, S1, S2, V )
COMPLEX S1, S2
INTEGER LDH, N
COMPLEX H( LDH, * ), V( * )
COMPLEX ZERO
PARAMETER ( ZERO = ( 0.0e0, 0.0e0 ) )
REAL RZERO
PARAMETER ( RZERO = 0.0e0 )
COMPLEX CDUM, H21S, H31S
REAL S
INTRINSIC ABS, AIMAG, REAL
REAL CABS1
CABS1( CDUM ) = ABS( REAL( CDUM ) ) + ABS( AIMAG( CDUM ) )
IF( N.EQ.2 ) THEN
S = CABS1( H( 1, 1 )-S2 ) + CABS1( H( 2, 1 ) )
IF( S.EQ.RZERO ) THEN
V( 1 ) = ZERO
V( 2 ) = ZERO
ELSE
H21S = H( 2, 1 ) / S
V( 1 ) = H21S*H( 1, 2 ) + ( H( 1, 1 )-S1 )* ( ( H( 1, 1
)-S2 ) / S )
V( 2 ) = H21S*( H( 1, 1 )+H( 2, 2 )-S1-S2 )
END IF
ELSE
S = CABS1( H( 1, 1 )-S2 ) + CABS1( H( 2, 1 ) ) +
CABS1( H( 3, 1 ) )
IF( S.EQ.ZERO ) THEN
V( 1 ) = ZERO
V( 2 ) = ZERO
V( 3 ) = ZERO
ELSE
H21S = H( 2, 1 ) / S
H31S = H( 3, 1 ) / S
V( 1 ) = ( H( 1, 1 )-S1 )*( ( H( 1, 1 )-S2 ) / S ) + H(
1, 2 )*H21S + H( 1, 3 )*H31S
V( 2 ) = H21S*( H( 1, 1 )+H( 2, 2 )-S1-S2 ) + H( 2, 3
)*H31S
V( 3 ) = H31S*( H( 1, 1 )+H( 3, 3 )-S1-S2 ) + H21S*H(
3, 2 )
END IF
END IF
END
PURPOSE
LAPACK auxiliary routine (versioNovember 2008 CLAQR1(1)
[top]
List of man pages available for Scientific
Copyright (c) for man pages and the logo by the respective OS vendor.
For those who want to learn more, the polarhome community provides shell access and support.
[legal]
[privacy]
[GNU]
[policy]
[cookies]
[netiquette]
[sponsors]
[FAQ]
Polarhome, production since 1999.
Member of Polarhome portal.
Based on Fawad Halim's script.
....................................................................
|
Vote for polarhome
|