pblas.F90

Go to the documentation of this file.

!! Copyright (C) 2011 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 pblas_oct_m
  implicit none
 
  private
 
  public ::            &
    pblas_gemm,        &
    pblas_herk,        &
    pblas_trsm,        &
    pblas_tran
 
 
  interface pblas_gemm
    subroutine pdgemm(transa, transb, m, n, k, alpha, &
      a, ia, ja, desca, b, ib, jb, descb, &
      beta, c, ic, jc, descc)
      use, intrinsic :: iso_fortran_env
      implicit none
 
      character(len=1) :: transa, transb
      integer          :: m, n, k, ia, ja, ib, jb, ic, jc
      real(real64)          :: alpha, beta
      real(real64)          :: a, b, c
      integer          :: desca, descb, descc
    end subroutine pdgemm
 
    subroutine pzgemm(transa, transb, m, n, k, alpha, &
      a, ia, ja, desca, b, ib, jb, descb, &
      beta, c, ic, jc, descc)
      use, intrinsic :: iso_fortran_env
      implicit none
 
      character(len=1) :: transa, transb
      integer          :: m, n, k, ia, ja, ib, jb, ic, jc
      complex(real64)       :: alpha, beta
      complex(real64)       :: a, b, c
      integer          :: desca, descb, descc
    end subroutine pzgemm
  end interface
 
  interface pblas_herk
    subroutine pdsyrk(uplo, trans, n, k, alpha, a, ia, ja, desca, beta, c, ic, jc, descc)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: uplo
      character(1), intent(in)    :: trans
      integer,      intent(in)    :: n
      integer,      intent(in)    :: k
      real(real64),      intent(in)    :: alpha
      real(real64),      intent(in)    :: a
      integer,      intent(in)    :: ia
      integer,      intent(in)    :: ja
      integer,      intent(in)    :: desca
      real(real64),      intent(in)    :: beta
      real(real64),      intent(inout) :: c
      integer,      intent(in)    :: ic
      integer,      intent(in)    :: jc
      integer,      intent(in)    :: descc
    end subroutine pdsyrk
 
    subroutine pzherk(uplo, trans, n, k, alpha, a, ia, ja, desca, beta, c, ic, jc, descc)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: uplo
      character(1), intent(in)    :: trans
      integer,      intent(in)    :: n
      integer,      intent(in)    :: k
      complex(real64),   intent(in)    :: alpha
      complex(real64),   intent(in)    :: a
      integer,      intent(in)    :: ia
      integer,      intent(in)    :: ja
      integer,      intent(in)    :: desca
      complex(real64),   intent(in)    :: beta
      complex(real64),   intent(inout) :: c
      integer,      intent(in)    :: ic
      integer,      intent(in)    :: jc
      integer,      intent(in)    :: descc
    end subroutine pzherk
 
  end interface pblas_herk
 
  interface pblas_trsm
    subroutine pdtrsm(side, uplo, transa, diag, m, n, alpha, a, ia, ja, desca, b, ib, jb, descb)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: side
      character(1), intent(in)    :: uplo
      character(1), intent(in)    :: transa
      character(1), intent(in)    :: diag
      integer,      intent(in)    :: m
      integer,      intent(in)    :: n
      real(real64),      intent(in)    :: alpha
      real(real64),      intent(in)    :: a
      integer,      intent(in)    :: ia
      integer,      intent(in)    :: ja
      integer,      intent(in)    :: desca
      real(real64),      intent(inout) :: b
      integer,      intent(in)    :: ib
      integer,      intent(in)    :: jb
      integer,      intent(in)    :: descb
    end subroutine pdtrsm
 
    subroutine pztrsm(side, uplo, transa, diag, m, n, alpha, a, ia, ja, desca, b, ib, jb, descb)
      use, intrinsic :: iso_fortran_env
      implicit none
      character(1), intent(in)    :: side
      character(1), intent(in)    :: uplo
      character(1), intent(in)    :: transa
      character(1), intent(in)    :: diag
      integer,      intent(in)    :: m
      integer,      intent(in)    :: n
      complex(real64),   intent(in)    :: alpha
      complex(real64),   intent(in)    :: a
      integer,      intent(in)    :: ia
      integer,      intent(in)    :: ja
      integer,      intent(in)    :: desca
      complex(real64),   intent(inout) :: b
      integer,      intent(in)    :: ib
      integer,      intent(in)    :: jb
      integer,      intent(in)    :: descb
    end subroutine pztrsm
 
  end interface pblas_trsm
 
  !  PDTRAN  transposes a matrix
  !
  !     sub( C ) := beta*sub( C ) + alpha*op( sub( A ) )
  !
  !  where
  !
  !     sub( C ) denotes C(IC:IC+M-1,JC:JC+N-1),
  !
  !     sub( A ) denotes A(IA:IA+N-1,JA:JA+M-1), and, op( X ) = X**T.
  !
  !  Thus, op( sub( A ) ) denotes A(IA:IA+N-1,JA:JA+M-1)**T.
  !
  !  Beta is a scalar, sub( C ) is an m by n submatrix, and sub( A ) is an
  !  n by m submatrix.
  interface pblas_tran
    subroutine pdtran(M, N, ALPHA, A, IA, JA, DESCA, BETA, C, IC, JC, DESCC)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer :: M, N, IA, JA, DESCA, IC, JC, DESCC
      real(real64) :: ALPHA, A, BETA, C
    end subroutine
 
    ! conjugated transpose
    subroutine pztranc(M, N, ALPHA, A, IA, JA, DESCA, BETA, C, IC, JC, DESCC)
      use, intrinsic :: iso_fortran_env
      implicit none
      integer :: M, N, IA, JA, DESCA, IC, JC, DESCC
      complex(real64) :: ALPHA, A, BETA, C
    end subroutine
  end interface pblas_tran
 
end module pblas_oct_m
 
!! Local Variables:
!! mode: f90
!! coding: utf-8
!! End: