line3n.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 shape_line3n
  integer, parameter, private :: kreal = kind(0.0d0)
 
contains
  subroutine shapefunc_line3n(lcoord,func)
    real(kind=kreal), intent(in) :: lcoord(1)
    real(kind=kreal) :: func(3)
    func(1) =-0.5d0*(1.d0-lcoord(1))*lcoord(1)
    func(2) = 0.5d0*(1.d0+lcoord(1))*lcoord(1)
    func(3) = 1.0d0-lcoord(1)*lcoord(1)
  end subroutine
 
  subroutine shapederiv_line3n(lcoord,func)
    real(kind=kreal), intent(in) :: lcoord(1)
    real(kind=kreal) :: func(3,1)
    func(1,1) = lcoord(1)-0.5d0
    func(2,1) = lcoord(1)+0.5d0
    func(3,1) =-2.d0*lcoord(1)
  end subroutine
 
  subroutine shape2ndderiv_line3n(func)
    real(kind=kreal) :: func(3,1,1)
    func(1,1,1) = 1.d0
    func(2,1,1) = 1.d0
    func(3,1,1) = -2.d0
  end subroutine
 
end module