lumatrix.f90

Path: libsrc/lumatrix/lumatrix.f90
Last Update: Mon Aug 19 17:49:31 +0900 2013

This file provides the following module.

Required files

libsrc/lumatrix/lumatrix.f90  

Methods

ludecomp21   ludecomp32   lusolve211   lusolve212   lusolve322   lusolve323  

Included Modules

dc_message

Public Instance methods

ludecomp21( alu, kp )
Subroutine :
alu(:,:) :real(8), intent(inout)
: 入力／ＬＵ行列
kp(size(alu,1)) :integer, intent(out)
: ピボット

ALU(NDIM,NDIM), KP(NDIM) NDIM x NDIM の行列を LU 分解. ＬＵ行列は 入力行列に上書きされる.

[Source]

subroutine ludecomp21(alu,kp)
  !
  ! ALU(NDIM,NDIM), KP(NDIM)
  ! NDIM x NDIM の行列を LU 分解.
  ! ＬＵ行列は 入力行列に上書きされる.
  !
  use dc_message

  real(8), intent(inout) :: alu(:,:)                  ! 入力／ＬＵ行列
  integer, intent(out)   :: kp(size(alu,1))           ! ピボット

  if ( size(alu,1) > size(alu,2) ) then
    call MessageNotify('E','ludecomp', 'The third dimension is less than the second')
  elseif( size(alu,1) < size(alu,2) ) then
    call MessageNotify('W','ludecomp', 'The third dimension is grater than the second')
  endif

  !" 行列のＬＵ分解（部分ピボット選択）
  call LUMAK1( alu, kp, size(alu,1) )

end subroutine ludecomp21
ludecomp32( alu, kp )
Subroutine :
alu(:,:,:) :real(8), intent(inout)
: 入力／ＬＵ行列
kp(size(alu,1),size(alu,2)) :integer, intent(out)
: ピボット

ALU(JDIM,NDIM,NDIM), KP(JDIM,NDIM) NDIM x NDIM の行列 JDIM 個を一度に LU 分解. ＬＵ行列は 入力行列に上書きされる.

[Source]

subroutine ludecomp32(alu,kp)
  !
  ! ALU(JDIM,NDIM,NDIM), KP(JDIM,NDIM)
  ! NDIM x NDIM の行列 JDIM 個を一度に LU 分解.
  ! ＬＵ行列は 入力行列に上書きされる.
  !
  use dc_message

  real(8), intent(inout) :: alu(:,:,:)                  ! 入力／ＬＵ行列
  integer, intent(out)   :: kp(size(alu,1),size(alu,2)) ! ピボット

  if ( size(alu,2) > size(alu,3) ) then
    call MessageNotify('E','ludecomp', 'The third dimension is less than the second')
  elseif( size(alu,2) < size(alu,3) ) then
    call MessageNotify('W','ludecomp', 'The third dimension is grater than the second')
  endif

  !" 行列のＬＵ分解（部分ピボット選択）
  call LUMAKE( alu, kp, size(alu,1), size(alu,2) )

end subroutine ludecomp32
lusolve211( alu, kp, b ) result(lusolve211)
Function :
lusolve211(size(b)) :real(8)
: 解
alu(:,:) :real(8), intent(in)
: 入力／ＬＵ行列
kp(:) :integer, intent(in)
: ピボット
b(:) :real(8), intent(in)
: 右辺ベクトル

ALU(NDIM,NDIM), KP(NDIM), B(NDIM) NDIM x NDIM 型行列の連立方程式 A X = B を 1 個の B に対して計算する.

[Source]

function lusolve211(alu,kp,b)
  !
  ! ALU(NDIM,NDIM), KP(NDIM), B(NDIM)
  ! NDIM x NDIM 型行列の連立方程式
  ! A X = B を 1 個の B に対して計算する.
  !
  use dc_message

  real(8), intent(in)  :: alu(:,:)              ! 入力／ＬＵ行列
  integer, intent(in)  :: kp(:)                 ! ピボット
  real(8), intent(in)  :: b(:)                  ! 右辺ベクトル

  real(8) :: lusolve211(size(b))                   ! 解

  lusolve211 = b
  call LUSOL2( lusolve211, alu , kp, 1, size(b) )

end function lusolve211
lusolve212( alu, kp, b ) result(lusolve212)
Function :
lusolve212(size(b,1),size(b,2)) :real(8)
: 解
alu(:,:) :real(8), intent(in)
: 入力／ＬＵ行列
kp(:) :integer, intent(in)
: ピボット
b(:,:) :real(8), intent(in)
: 右辺ベクトル

ALU(NDIM,NDIM), KP(NDIM), B(JDIM,NDIM) NDIM x NDIM 型行列の連立方程式 A X = B を JDIM 個の B に対して計算する.

[Source]

function lusolve212(alu,kp,b)
  !
  ! ALU(NDIM,NDIM), KP(NDIM), B(JDIM,NDIM)
  ! NDIM x NDIM 型行列の連立方程式
  ! A X = B を JDIM 個の B に対して計算する.
  !
  use dc_message

  real(8), intent(in)  :: alu(:,:)              ! 入力／ＬＵ行列
  integer, intent(in)  :: kp(:)                 ! ピボット
  real(8), intent(in)  :: b(:,:)                ! 右辺ベクトル

  real(8) :: lusolve212(size(b,1),size(b,2))       ! 解

  lusolve212 = b
  call LUSOL2( lusolve212, alu , kp, size(b,1), size(b,2) )

end function lusolve212
lusolve322( alu, kp, b ) result(lusolve322)
Function :
lusolve322(size(b,1),size(b,2)) :real(8)
: 解
alu(:,:,:) :real(8), intent(in)
: 入力／ＬＵ行列
kp(:,:) :integer, intent(in)
: ピボット
b(:,:) :real(8), intent(in)
: 右辺ベクトル

ALU(JDIM,NDIM,NDIM), KP(JDIM,NDIM), B(JDIM,NDIM) NDIM x NDIM 型行列を JDIM 個並べた連立方程式 A X = B をひとつの B の並びに対して計算する.

[Source]

function lusolve322(alu,kp,b)
  !
  ! ALU(JDIM,NDIM,NDIM), KP(JDIM,NDIM), B(JDIM,NDIM)
  ! NDIM x NDIM 型行列を JDIM 個並べた連立方程式
  ! A X = B をひとつの B の並びに対して計算する.
  !
  use dc_message

  real(8), intent(in)  :: alu(:,:,:)                   ! 入力／ＬＵ行列
  integer, intent(in)  :: kp(:,:)                      ! ピボット
  real(8), intent(in)  :: b(:,:)                       ! 右辺ベクトル

  real(8) :: lusolve322(size(b,1),size(b,2))             ! 解

  lusolve322 = b
  call LUSOLV( lusolve322, alu , kp, 1, size(b,1), size(b,2) )

end function lusolve322
lusolve323( alu, kp, b ) result(lusolve323)
Function :
lusolve323(size(b,1),size(b,2),size(b,3)) :real(8)
: 解
alu(:,:,:) :real(8), intent(in)
: 入力／ＬＵ行列
kp(:,:) :integer, intent(in)
: ピボット
b(:,:,:) :real(8), intent(in)
: 右辺ベクトル

ALU(JDIM,NDIM,NDIM), KP(JDIM,NDIM), B(IDIM,JDIM,NDIM) NDIM x NDIM 型行列を JDIM 個並べた連立方程式 A X = B を IDIM 個の B に対して計算する.

[Source]

function lusolve323(alu,kp,b)
  !
  ! ALU(JDIM,NDIM,NDIM), KP(JDIM,NDIM), B(IDIM,JDIM,NDIM)
  ! NDIM x NDIM 型行列を JDIM 個並べた連立方程式
  ! A X = B を IDIM 個の B に対して計算する.
  !
  use dc_message

  real(8), intent(in)  :: alu(:,:,:)                   ! 入力／ＬＵ行列
  integer, intent(in)  :: kp(:,:)                      ! ピボット
  real(8), intent(in)  :: b(:,:,:)                     ! 右辺ベクトル

  real(8) :: lusolve323(size(b,1),size(b,2),size(b,3)) ! 解

  lusolve323 = b
  call LUSOLV( lusolve323, alu , kp, size(b,1), size(b,2), size(b,3) )

end function lusolve323