!> If scale * sqrt( sumsq ) > tbig then
!> we now require: scl >= sqrt( TINY*EPS ) / sbig on entry,
!> and if scale * sqrt( sumsq ) < tsml then
!> we now require: scl <= sqrt( HUGE ) / ssml on entry,
!> where
!> tbig -- upper threshold for values whose square is representable;
!> sbig -- scaling constant for big numbers; \see la_constants.f90
!> tsml -- lower threshold for values whose square is representable;
!> ssml -- scaling constant for small numbers; \see la_constants.f90
!> and
!> TINY*EPS -- tiniest representable number;
!> HUGE -- biggest representable number.

• Commit e8a6fae
• #570: Revert "Merge pull request #290"
• Branch master

Weslley S. Pereira