lapack.F90

Go to the documentation of this file.

!! Copyright (C) 2009 X. Andrade
!!
!! This program is free software; you can redistribute it and/or modify
!! it under the terms of the GNU General Public License as published by
!! the Free Software Foundation; either version 2, or (at your option)
!! any later version.
!!
!! This program is distributed in the hope that it will be useful,
!! but WITHOUT ANY WARRANTY; without even the implied warranty of
!! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!! GNU General Public License for more details.
!!
!! You should have received a copy of the GNU General Public License
!! along with this program; if not, write to the Free Software
!! Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
!! 02110-1301, USA.
!!
 
#include "global.h"
 
! -----------------------------------------------------------------------
! -----------------------------------------------------------------------
 
module lapack_oct_m
  use, intrinsic :: iso_fortran_env
  implicit none
 
  public ! only interfaces in this module
 
  interface lapack_potrf
    subroutine dpotrf(uplo, n, a, lda, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: uplo
      integer,      intent(in)    :: n
      real(real64),      intent(inout) :: a
      integer,      intent(in)    :: lda
      integer,      intent(out)   :: info
    end subroutine dpotrf
 
    subroutine zpotrf(uplo, n, a, lda, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: uplo
      integer,      intent(in)    :: n
      complex(real64),   intent(inout) :: a
      integer,      intent(in)    :: lda
      integer,      intent(out)   :: info
    end subroutine zpotrf
  end interface lapack_potrf
 
  interface
    subroutine dpstrf(uplo, n, a, lda, piv, rank, tol, work, info)
      import real64
      implicit none
      character(len=1), intent(in)    :: uplo
      integer,          intent(in)    :: n
      real(real64),     intent(inout) :: a
      integer,          intent(in)    :: lda
      integer,          intent(out)   :: piv
      integer,          intent(out)   :: rank
      real(real64),     intent(in)    :: tol
      real(real64),     intent(inout) :: work
      integer,          intent(out)   :: info
    end subroutine dpstrf
 
    subroutine zpstrf(uplo, n, a, lda, piv, rank, tol, work, info)
      import real64
      implicit none
      character(len=1), intent(in)       :: uplo
      integer,          intent(in)       :: n
      complex(real64),  intent(inout)    :: a
      integer,          intent(in)       :: lda
      integer,          intent(out)      :: piv
      integer,          intent(out)      :: rank
      real(real64),     intent(in)       :: tol
      real(real64),     intent(inout)    :: work
      integer,          intent(out)      :: info
    end subroutine zpstrf
  end interface
 
  interface lapack_sygv
    subroutine dsygv(itype, jobz, uplo, n, a, lda, b, ldb, w, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1),      intent(in)    :: jobz, uplo
      integer,           intent(in)    :: itype, n, lda, ldb, lwork
      real(real64),      intent(inout) :: a, b
      real(real64),      intent(out)   :: w, work
      integer,           intent(out)   :: info
    end subroutine dsygv
  end interface lapack_sygv
 
  interface lapack_sygvd
    subroutine dsygvd(itype, jobz, uplo, n, a, lda, b, ldb, w, work, lwork, iwork, liwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1),      intent(in)    :: jobz, uplo
      integer,           intent(in)    :: itype, n, lda, ldb, lwork, liwork
      real(real64),      intent(inout) :: a, b
      real(real64),      intent(out)   :: w, work
      integer,           intent(out)   :: iwork, info
    end subroutine dsygvd
  end interface lapack_sygvd
 
  interface lapack_hegv
    subroutine zhegv(itype, jobz, uplo, n, a, lda, b, ldb, w, work, lwork, rwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1),      intent(in)    :: jobz, uplo
      integer,           intent(in)    :: n, itype, lda, ldb, lwork
      complex(real64),   intent(inout) :: a, b
      real(real64),      intent(out)   :: w, rwork
      complex(real64),   intent(out)   :: work
      integer,           intent(out)   :: info
    end subroutine zhegv
  end interface lapack_hegv
 
  interface lapack_hegvd
    subroutine zhegvd(itype, jobz, uplo, n, a, lda, b, ldb, w, work, lwork, rwork, lrwork, iwork, liwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1),      intent(in)    :: jobz, uplo
      integer,           intent(in)    :: n, itype, lda, ldb, lwork, lrwork, liwork
      complex(real64),   intent(inout) :: a, b
      real(real64),      intent(out)   :: w, rwork
      complex(real64),   intent(out)   :: work
      integer,           intent(out)   :: iwork, info
    end subroutine zhegvd
  end interface lapack_hegvd
 
  interface
    subroutine dgeev(jobvl, jobvr, n, a, lda, wr, wi, vl, ldvl, vr, ldvr, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: jobvl, jobvr
      integer,      intent(in)    :: n, lda, ldvl, ldvr, lwork
      real(real64),      intent(inout) :: a
      real(real64),      intent(out)   :: wr, wi, vl, vr
      real(real64),      intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine dgeev
 
    subroutine zgeev(jobvl, jobvr, n, a, lda, w, vl, ldvl, vr, ldvr, work, lwork, rwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: jobvl, jobvr
      integer,      intent(in)    :: n, lda, ldvl, ldvr, lwork
      complex(real64),   intent(inout) :: a
      complex(real64),   intent(out)   :: w, vl, vr
      real(real64),      intent(out)   :: rwork
      complex(real64),   intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine zgeev
  end interface
 
  interface lapack_gesvx
    subroutine dgesvx(fact, trans, n, nrhs, a, lda, af, ldaf, ipiv, equed, r, &
      c, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: fact, trans
      integer,      intent(in)    :: n, nrhs, lda, ldaf, ldb, ldx
      real(real64),     intent(inout) :: a, af, r, c, b
      integer,      intent(inout) :: ipiv
      real(real64),     intent(out)   :: x, ferr, berr, work
      real(real64),     intent(out)   :: rcond
      character(1), intent(inout) :: equed
      integer,      intent(out)   :: iwork
      integer,      intent(out)   :: info
    end subroutine dgesvx
 
    subroutine zgesvx (fact, trans, n, nrhs, a, lda, af, ldaf, ipiv, equed, r, &
      c, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: fact, trans
      integer,      intent(in)    :: n, nrhs, lda, ldaf, ldb, ldx
      complex(real64),  intent(inout) :: a, af, b
      real(real64),     intent(inout) :: r, c
      integer,      intent(inout) :: ipiv
      real(real64),     intent(out)   :: ferr, berr
      real(real64),     intent(out)   :: rcond
      complex(real64),  intent(out)   :: x, work
      character(1), intent(inout) :: equed
      real(real64),     intent(out)   :: rwork
      integer,      intent(out)   :: info
    end subroutine zgesvx
  end interface lapack_gesvx
 
  interface lapack_syev
    subroutine dsyev(jobz, uplo, n, a, lda, w, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: jobz, uplo
      integer,      intent(in)    :: n, lda, lwork
      real(real64),      intent(inout) :: a
      real(real64),      intent(out)   :: w, work
      integer,      intent(out)   :: info
    end subroutine dsyev
  end interface lapack_syev
 
  interface lapack_heev
    subroutine zheev(jobz, uplo, n, a, lda, w, work, lwork, rwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: jobz, uplo
      integer,      intent(in)    :: n, lda, lwork
      complex(real64),   intent(inout) :: a
      real(real64),      intent(out)   :: w, rwork
      complex(real64),   intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine zheev
  end interface lapack_heev
 
  interface
    subroutine dsyevx(jobz, range, uplo, n, a, lda, &
      vl, vu, il, iu, abstol, m, w, z, ldz, work, lwork, iwork, ifail, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer,      intent(in)    :: n, lda, il, iu, ldz, lwork
      character(1), intent(in)    :: jobz, range, uplo
      integer,      intent(out)   :: m, iwork, ifail, info
      real(real64),     intent(in)    :: vl, vu, abstol
      real(real64),     intent(inout) :: a
      real(real64),     intent(out)   :: w, z, work
    end subroutine dsyevx
 
    subroutine zheevx(jobz, range, uplo, n, a, lda, &
      vl, vu, il, iu, abstol, m, w, z, ldz, work, lwork, iwork, ifail, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer,      intent(in)    :: n, lda, il, iu, ldz, lwork
      character(1), intent(in)    :: jobz, range, uplo
      integer,      intent(out)   :: m, iwork, ifail, info
      real(real64),     intent(in)    :: vl, vu, abstol
      real(real64),     intent(out)   :: w
      complex(real64),  intent(inout) :: a
      complex(real64),  intent(out)   :: z, work
    end subroutine zheevx
  end interface
 
  interface lapack_geqrf
    subroutine dgeqrf(m, n, a, lda, tau, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer, intent(in)    :: lda, lwork, m, n
      real(real64), intent(inout) :: a
      real(real64), intent(out)   :: tau, work
      integer, intent(out)   :: info
    end subroutine dgeqrf
 
    subroutine zgeqrf(m, n, a, lda, tau, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer, intent(in)       :: lda, lwork, m, n
      complex(real64), intent(inout) :: a
      complex(real64), intent(out)   :: tau, work
      integer, intent(out)      :: info
    end subroutine zgeqrf
  end interface lapack_geqrf
 
  interface lapack_orgqr
    subroutine dorgqr(m, n, k, a, lda, tau, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer, intent(in)    :: k, lda, lwork, m, n
      real(real64), intent(in)    :: tau
      real(real64), intent(inout) :: a
      real(real64), intent(out)   :: work
      integer, intent(out)   :: info
    end subroutine dorgqr
 
    subroutine zungqr(m, n, k, a, lda, tau, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer,    intent(in)    :: k, lda, lwork, m, n
      complex(real64), intent(in)    :: tau
      complex(real64), intent(inout) :: a
      complex(real64), intent(out)   :: work
      integer,    intent(out)   :: info
    end subroutine zungqr
  end interface lapack_orgqr
 
  interface lapack_sygvx
    subroutine dsygvx(itype, jobz, range, uplo, n, a, lda, b, ldb, vl, vu, il, iu, abstol, &
      m, w, z, ldz, work, lwork, iwork, ifail, info)
      use, intrinsic :: iso_fortran_env
      implicit none
 
      integer,             intent(in)    :: itype
      character(len=1),    intent(in)    :: jobz
      character(len=1),    intent(in)    :: range
      character(len=1),    intent(in)    :: uplo
      integer,             intent(in)    :: n
      real(real64),             intent(inout) :: a
      integer,             intent(in)    :: lda
      real(real64),             intent(inout) :: b
      integer,             intent(in)    :: ldb
      real(real64),             intent(in)    :: vl
      real(real64),             intent(in)    :: vu
      integer,             intent(in)    :: il
      integer,             intent(in)    :: iu
      real(real64),             intent(in)    :: abstol
      integer,             intent(out)   :: m
      real(real64),             intent(out)   :: w
      real(real64),             intent(out)   :: z
      integer,             intent(in)    :: ldz
      real(real64),             intent(out)   :: work
      integer,             intent(in)    :: lwork
      integer,             intent(out)   :: iwork
      integer,             intent(out)   :: ifail
      integer,             intent(out)   :: info
    end subroutine dsygvx
  end interface lapack_sygvx
 
  interface lapack_hegvx
    subroutine zhegvx(itype, jobz, range, uplo, n, a, lda, b, ldb, vl, vu, il, iu, abstol, &
      m, w, z, ldz, work, lwork, rwork, iwork, ifail, info)
      use, intrinsic :: iso_fortran_env
      implicit none
 
      integer,             intent(in)    :: itype
      character(len=1),    intent(in)    :: jobz
      character(len=1),    intent(in)    :: range
      character(len=1),    intent(in)    :: uplo
      integer,             intent(in)    :: n
      complex(real64),          intent(inout) :: a
      integer,             intent(in)    :: lda
      complex(real64),          intent(inout) :: b
      integer,             intent(in)    :: ldb
      real(real64),             intent(in)    :: vl
      real(real64),             intent(in)    :: vu
      integer,             intent(in)    :: il
      integer,             intent(in)    :: iu
      real(real64),             intent(in)    :: abstol
      integer,             intent(out)   :: m
      real(real64),             intent(out)   :: w
      complex(real64),          intent(out)   :: z
      integer,             intent(in)    :: ldz
      complex(real64),          intent(out)   :: work
      integer,             intent(in)    :: lwork
      real(real64),             intent(out)   :: rwork
      integer,             intent(out)   :: iwork
      integer,             intent(out)   :: ifail
      integer,             intent(out)   :: info
    end subroutine zhegvx
  end interface lapack_hegvx
 
  interface lapack_gelss
    subroutine dgelss(m, n, nrhs, a, lda, b, ldb, s, rcond, rank, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer, intent(in)    :: m
      integer, intent(in)    :: n
      integer, intent(in)    :: nrhs
      real(real64),intent(inout) :: a
      integer, intent(in)    :: lda
      real(real64),intent(inout) :: b
      integer, intent(in)    :: ldb
      real(real64),intent(out)   :: s
      real(real64),intent(in)    :: rcond
      integer, intent(out)   :: rank
      real(real64),intent(out)   :: work
      integer, intent(in)    :: lwork
      integer, intent(out)   :: info
    end subroutine dgelss
 
    subroutine zgelss(m, n, nrhs, a, lda, b, ldb, s, rcond, rank, work, lwork, rwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer, intent(in)    :: m
      integer, intent(in)    :: n
      integer, intent(in)    :: nrhs
      complex(real64),   intent(inout) :: a
      integer, intent(in)    :: lda
      complex(real64),   intent(inout) :: b
      integer, intent(in)    :: ldb
      real(real64),intent(out)   :: s
      real(real64),intent(in)    :: rcond
      integer, intent(out)   :: rank
      complex(real64),   intent(out)   :: work
      integer, intent(in)    :: lwork
      real(real64),intent(out)   :: rwork
      integer, intent(out)   :: info
    end subroutine zgelss
  end interface lapack_gelss
 
  interface lapack_getrf
    subroutine dgetrf (m, n, a, lda, ipiv, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer,      intent(in)    :: m, n, lda
      real(real64),     intent(inout) :: a
      integer,      intent(out)   :: ipiv
      integer,      intent(out)   :: info
    end subroutine dgetrf
 
    subroutine zgetrf (m, n, a, lda, ipiv, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer,      intent(in)    :: m, n, lda
      complex(real64),  intent(inout) :: a
      integer,      intent(out)   :: ipiv
      integer,      intent(out)   :: info
    end subroutine zgetrf
  end interface lapack_getrf
 
  interface lapack_getri
    subroutine dgetri(n, a, lda, ipiv, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer,      intent(in)    :: n, lda, lwork
      real(real64),     intent(inout) :: a
      integer,      intent(in)    :: ipiv
      real(real64),     intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine dgetri
 
    subroutine zgetri(n, a, lda, ipiv, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer,      intent(in)    :: n, lda, lwork
      complex(real64),  intent(inout) :: a
      integer,      intent(in)    :: ipiv
      complex(real64),  intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine zgetri
  end interface lapack_getri
 
  interface lapack_sytrf
    subroutine dsytrf(uplo, n, a, lda, ipiv, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: uplo
      integer,      intent(in)    :: n, lda, lwork
      real(real64),     intent(inout) :: a
      integer,      intent(out)   :: ipiv
      real(real64),     intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine dsytrf
 
    subroutine zsytrf(uplo, n, a, lda, ipiv, work, lwork, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: uplo
      integer,      intent(in)    :: n, lda, lwork
      complex(real64),  intent(inout) :: a
      integer,      intent(out)   :: ipiv
      complex(real64),  intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine zsytrf
  end interface lapack_sytrf
 
  interface lapack_sytri
    subroutine dsytri (uplo, n, a, lda, ipiv, work, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: uplo
      integer,      intent(in)    :: n, lda
      real(real64),     intent(inout) :: a
      integer,      intent(in)    :: ipiv
      real(real64),     intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine dsytri
 
    subroutine zsytri (uplo, n, a, lda, ipiv, work, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: uplo
      integer,      intent(in)    :: n, lda
      complex(real64),  intent(inout) :: a
      integer,      intent(in)    :: ipiv
      complex(real64),  intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine zsytri
  end interface lapack_sytri
 
  interface lapack_stev
    subroutine dstev (jobz, n, d, e, z, ldz, work, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: jobz
      integer,      intent(in)    :: n
      real(real64), intent(inout) :: d, e
      real(real64), intent(out)   :: z
      integer,      intent(in)    :: ldz
      real(real64), intent(out)   :: work
      integer,      intent(out)   :: info
    end subroutine dstev
  end interface lapack_stev
 
  interface lapack_steqr
    subroutine zsteqr (compz, n, d, e, z, ldz, work, info)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1),    intent(in)    :: compz
      integer,         intent(in)    :: n
      real(real64),    intent(inout) :: d, e
      complex(real64), intent(inout) :: z
      integer,         intent(in)    :: ldz
      real(real64),    intent(out)   :: work
      integer,         intent(out)   :: info
    end subroutine zsteqr
  end interface lapack_steqr
 
 
end module lapack_oct_m
 
!! Local Variables:
!! mode: f90
!! coding: utf-8
!! End: