f04asf
Scheduled for withdrawal at Mark 28.
Replaced by f07fbf.
f04abf and f04asf have been replaced by f07fbf for performance. The replacement routine is threaded by NAG and may also be threaded in the vendor library (BLAS/LAPACK).
Old Code
CALL f04asf(A,LDA,B,N,C,WK1,WK2,IFAIL)
New Code
CALL f04asf_wrap(a,lda,b,n,c,wk1,wk2,ifail)
 Subroutine f04asf_wrap(a,lda,b,n,c,wk1,wk2,ifail)
  !     .. Use Statements ..
   Use nag_library, Only: dposvx, nag_wp
  !     .. Scalar Arguments ..
   Integer, Intent (In)             :: lda, n
   Integer, Intent (Inout)          :: ifail
  !     .. Array Arguments ..
   Real (Kind=nag_wp), Intent (Inout) :: a(lda,n), b(n)
   Real (Kind=nag_wp), Intent (Out) :: c(n), wk1(1), wk2(1)
  !     .. Local Scalars ..
   Real (Kind=nag_wp)               :: rcond
   Integer                          :: info, ldaf, m
   Character (1)                    :: equed
  !     .. Local Arrays ..
   Real (Kind=nag_wp)               :: s(1), ferr(1), berr(1)
   Real (Kind=nag_wp), Allocatable  :: af(:,:), work(:)
   Integer, Allocatable             :: iwork(:)

   ldaf = n
   m = 1
   Allocate (af(ldaf,n),work(3*n),iwork(n))
  !     The NAG name equivalent of dposvx is f07fbf
   Call dposvx('N','Upper',n,m,a,lda,af,ldaf,equed,s,b,n,  &
       c,n,rcond,ferr,berr,work,iwork,info)
   wk1(1) = rcond
   wk2(1) = berr(1)

   ifail = info
      End Subroutine f04asf_wrap
f04atf
Scheduled for withdrawal at Mark 28.
Replaced by f07abf.
f04aef and f04atf have been replaced by f07abf for performance. The replacement routine is threaded by NAG and may also be threaded in the vendor library (BLAS/LAPACK).
Old Code
Call f04atf(A,LDA,B,N,C,AA,LDAA,WKS1,WKS2,IFAIL)
New Code
CALL f04atf_wrap(a,lda,b,n,c,aa,ldaa,wks1,wks2,ifail)
 Subroutine f04atf_wrap(a,lda,b,n,c,aa,ldaa,wks1,wks2,ifail)
  !     .. Use Statements ..
   Use nag_library, Only: dgesvx, nag_wp
  !     .. Scalar Arguments ..
   Integer, Intent (In)             :: lda, ldaa, n
   Integer, Intent (Inout)          :: ifail
  !     .. Array Arguments ..
   Real (Kind=nag_wp), Intent (Inout) :: a(lda,*), b(n)
   Real (Kind=nag_wp), Intent (Out)   :: c(n), aa(ldaa,n), wks1(1),
  wks2(1)
  !     .. Local Scalars ..
   Real (Kind=nag_wp)               :: rcond
   Integer                          :: info, m
   Character (1)                    :: equed
  !     .. Local Arrays ..
   Real (Kind=nag_wp)               :: cscl(1), rscl(1), ferr(1),
  berr(1)
   Real (Kind=nag_wp), Allocatable  :: work(:)
   Integer, Allocatable             :: ipiv(:), iwork(:)

   m = 1
   Allocate (work(4*n),ipiv(n),iwork(n))
  !     The NAG name equivalent of dgesvx is f07abf
   Call dgesvx('N','N',n,m,a,lda,aa,ldaa,ipiv,equed,rscl,cscl,b,n, &
       c,n,rcond,ferr,berr,work,iwork,info)

   ifail = info
   wks1(1) = rcond
   wks2(1) = berr(1)

      End Subroutine f04atf_wrap
f04eaf