dens_2dsquare Function

public function dens_2dsquare(x, ts) result(dos)

Arguments

Type IntentOptional Attributes Name
real(kind=8), intent(in) :: x
real(kind=8), intent(in), optional :: ts

Return Value real(kind=8)


Source Code

  function dens_2dsquare(x,ts) result(dos)
    real(8),intent(in)          :: x
    real(8),intent(in),optional :: ts
    real(8)                     :: wband,y,kint,eint,dos
    real(8),parameter           :: pi=acos(-1d0)
    wband=4.d0;if(present(ts))wband=4.d0*ts
    dos=0d0
    if(abs(x)>wband)return
    y  = sqrt(1d0 - (x/wband)**2)
    dos= 2/pi**2/wband*EllipticK(y)
  end function dens_2dsquare