slarmm()

real function slarmm	(	real	anorm,
		real	bnorm,
		real	cnorm )

SLARMM

Purpose:

!>
!> SLARMM returns a factor s in (0, 1] such that the linear updates
!>
!>    (s * C) - A * (s * B)  and  (s * C) - (s * A) * B
!>
!> cannot overflow, where A, B, and C are matrices of conforming
!> dimensions.
!>
!> This is an auxiliary routine so there is no argument checking.
!> 

Parameters

[in]	ANORM	!> ANORM is REAL !> The infinity norm of A. ANORM >= 0. !> The number of rows of the matrix A. M >= 0. !>
[in]	BNORM	!> BNORM is REAL !> The infinity norm of B. BNORM >= 0. !>
[in]	CNORM	!> CNORM is REAL !> The infinity norm of C. CNORM >= 0. !> References: C. C. Kjelgaard Mikkelsen and L. Karlsson, Blocked Algorithms for Robust Solution of Triangular Linear Systems. In: International Conference on Parallel Processing and Applied Mathematics, pages 68–78. Springer, 2017.

Definition at line 60 of file slarmm.f.

      IMPLICIT NONE
*     .. Scalar Arguments ..
      REAL               ANORM, BNORM, CNORM
*     .. Parameters ..
      REAL               ONE, HALF, FOUR
      parameter( one = 1.0e0, half = 0.5e+0, four = 4.0e+0 )
*     ..
*     .. Local Scalars ..
      REAL               BIGNUM, SMLNUM
*     ..
*     .. External Functions ..
      REAL               SLAMCH
      EXTERNAL           slamch
*     ..
*     .. Executable Statements ..
*
*
*     Determine machine dependent parameters to control overflow.
*
      smlnum = slamch( 'Safe minimum' ) / slamch( 'Precision' )
      bignum = ( one / smlnum ) / four
*
*     Compute a scale factor.
*
      slarmm = one
      IF( bnorm .LE. one ) THEN
         IF( anorm * bnorm .GT. bignum - cnorm ) THEN
            slarmm = half
         END IF
      ELSE
         IF( anorm .GT. (bignum - cnorm) / bnorm ) THEN
            slarmm = half / bnorm
         END IF
      END IF
      RETURN
*
*     ==== End of SLARMM ====
*

Here is the call graph for this function:

Here is the caller graph for this function: