| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | module array_utils | ||
| 2 | implicit none | ||
| 3 | |||
| 4 | private | ||
| 5 | public :: init_derivative_factors | ||
| 6 | |||
| 7 | contains | ||
| 8 | |||
| 9 | !> Initialize factorial-based derivative factors for polynomial derivatives | ||
| 10 | !> derf1(k) = (k-1) | ||
| 11 | !> derf2(k) = (k-1)*(k-2) | ||
| 12 | !> derf3(k) = (k-1)*(k-2)*(k-3) | ||
| 13 | ✗ | pure subroutine init_derivative_factors(ns_max, derf1, derf2, derf3) | |
| 14 | integer, intent(in) :: ns_max | ||
| 15 | double precision, intent(out) :: derf1(ns_max), derf2(ns_max), derf3(ns_max) | ||
| 16 | integer :: k | ||
| 17 | |||
| 18 | ✗ | do k = 1, ns_max | |
| 19 | ✗ | derf1(k) = dble(k - 1) | |
| 20 | ✗ | derf2(k) = dble((k - 1)*(k - 2)) | |
| 21 | ✗ | derf3(k) = dble((k - 1)*(k - 2)*(k - 3)) | |
| 22 | end do | ||
| 23 | |||
| 24 | ✗ | end subroutine init_derivative_factors | |
| 25 | |||
| 26 | end module array_utils | ||
| 27 |