logspace Function

public function logspace(start, stop, num, base) result(array)

Arguments

Type IntentOptional Attributes Name
real(kind=8) :: start
real(kind=8) :: stop
integer :: num
real(kind=8), optional :: base

Return Value real(kind=8), (num)


Calls

proc~~logspace~~CallsGraph proc~logspace logspace proc~linspace~2 linspace proc~logspace->proc~linspace~2

Source Code

  function logspace(start,stop,num,base) result(array)
    real(8)          :: start,stop,array(num)
    integer          :: num,i
    ! logical,optional :: iend
    ! logical          :: endpoint_
    real(8),optional :: base
    real(8)          :: base_
    real(8)          :: A,B
    base_= 10.d0;if(present(base))base_=base
    if(num<0)stop "logspace: N<0, abort."
    A=start;if(start==0.d0)A=1.d-12
    B=stop;if(stop==0.d0)B=1.d-12
    A=log(A)/log(base_) ; B=log(B)/log(base_)
    array = linspace(A,B,num=num,iend=.true.)
    array = base_**array
  end function logspace