hecmw_allocate.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2019 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
 
module hecmw_allocate
contains
 
  subroutine hecmw_allocate_matrix(hecMAT, mat, nBlock)
    use hecmw_util
 
    implicit none
 
    type (hecmwST_matrix):: hecMAT, mat
    integer:: nBlock
    integer:: ierr
 
    allocate( mat%AL( nblock*hecmat%NPL ), stat=ierr )
    if ( ierr /= 0 ) then
      write(0,*) 'Allocation error: Not enough memory for matrix!'
      stop
    end if
    allocate( mat%AU( nblock*hecmat%NPU ), stat=ierr )
    if ( ierr /= 0 ) then
      write(0,*) 'Allocation error: Not enough memory for matrix!'
      stop
    end if
    allocate( mat%D ( nblock*hecmat%NP  ), stat=ierr )
    if ( ierr /= 0 ) then
      write(0,*) 'Allocation error: Not enough memory for matrix!'
      stop
    end if
 
    mat%AL = 0.0d0
    mat%AU = 0.0d0
    mat%D  = 0.0d0
 
  end subroutine hecmw_allocate_matrix
 
  subroutine hecmw_allocate_vector_i(vector, size)
    use hecmw_util
    implicit none
 
    integer(kind=kint), pointer:: vector(:)
    integer(kind=kint):: size
    integer(kind=kint):: ierr
 
    allocate( vector( size ), stat=ierr )
    if ( ierr /= 0 ) then
      write(0,*) 'Allocation error: Not enough memory for integer array!'
      stop
    end if
 
    vector = 0
  end subroutine hecmw_allocate_vector_i
 
  subroutine hecmw_allocate_vector_r(vector, size)
    use hecmw_util
    implicit none
 
    real(kind=kreal), pointer:: vector(:)
    integer(kind=kint):: size
    integer(kind=kint):: ierr
 
    allocate( vector( size ), stat=ierr )
    if ( ierr /= 0 ) then
      write(0,*) 'Allocation error: Not enough memory for real array!'
      stop
    end if
 
    vector = 0.0d0
  end subroutine hecmw_allocate_vector_r
 
end module hecmw_allocate