Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=8), | intent(in), | dimension(:) | :: | f | ||
real(kind=8), | intent(in) | :: | dh | |||
integer, | intent(in), | optional | :: | order |
function derivative2(f,dh,order) result(df) real(8),dimension(:),intent(in) :: f real(8),intent(in) :: dh integer,intent(in),optional :: order real(8),dimension(size(f)) :: df integer :: i,L,order_ L=size(f) order_=4;if(present(order))order_=order if(L < order_ + 1) stop "derivative2: L < order+1." select case(order_) case(2) df = derivF2_n222(f,dh) case(4) df = derivF2_n444(f,dh) case(6) df = derivF2_n666(f,dh) case default stop "derivative2: order is not 2,4 or 6" end select end function derivative2