| Directory: | ./ |
|---|---|
| File: | misc/arth_m.f90 |
| Date: | 2022-01-11 19:19:34 |
| Exec | Total | Coverage | |
|---|---|---|---|
| Lines: | 12 | 28 | 42.9% |
| Branches: | 12 | 36 | 33.3% |
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | MODULE arth_m | ||
| 2 | |||
| 3 | IMPLICIT NONE | ||
| 4 | |||
| 5 | INTEGER, PARAMETER, private:: NPAR_ARTH=16, NPAR2_ARTH=8 | ||
| 6 | |||
| 7 | INTERFACE arth | ||
| 8 | ! Returns an arithmetic progression, given a first term "first", an | ||
| 9 | ! increment and a number of terms "n" (including "first"). | ||
| 10 | |||
| 11 | MODULE PROCEDURE arth_r, arth_i | ||
| 12 | ! The difference between the procedures is the kind and type of | ||
| 13 | ! arguments first and increment and of function result. | ||
| 14 | END INTERFACE | ||
| 15 | |||
| 16 | private arth_r, arth_i | ||
| 17 | |||
| 18 | CONTAINS | ||
| 19 | |||
| 20 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | pure FUNCTION arth_r(first,increment,n) |
| 21 | |||
| 22 | REAL, INTENT(IN) :: first,increment | ||
| 23 | INTEGER, INTENT(IN) :: n | ||
| 24 | REAL arth_r(n) | ||
| 25 | |||
| 26 | ! Local: | ||
| 27 | INTEGER :: k,k2 | ||
| 28 | REAL :: temp | ||
| 29 | |||
| 30 | !--------------------------------------- | ||
| 31 | |||
| 32 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | if (n > 0) arth_r(1)=first |
| 33 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | if (n <= NPAR_ARTH) then |
| 34 | ✗ | do k=2,n | |
| 35 | ✗ | arth_r(k)=arth_r(k-1)+increment | |
| 36 | end do | ||
| 37 | else | ||
| 38 |
2/2✓ Branch 0 taken 7 times.
✓ Branch 1 taken 1 times.
|
8 | do k=2,NPAR2_ARTH |
| 39 | 8 | arth_r(k)=arth_r(k-1)+increment | |
| 40 | end do | ||
| 41 | 1 | temp=increment*NPAR2_ARTH | |
| 42 | k=NPAR2_ARTH | ||
| 43 | 2 | do | |
| 44 |
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 2 times.
|
3 | if (k >= n) exit |
| 45 | 2 | k2=k+k | |
| 46 |
5/6✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 24 times.
✓ Branch 4 taken 24 times.
✓ Branch 5 taken 2 times.
|
52 | arth_r(k+1:min(k2,n)) = temp + arth_r(1:min(k,n-k)) |
| 47 | 2 | temp=temp+temp | |
| 48 | 1 | k=k2 | |
| 49 | end do | ||
| 50 | end if | ||
| 51 | |||
| 52 | END FUNCTION arth_r | ||
| 53 | |||
| 54 | !************************************* | ||
| 55 | |||
| 56 | ✗ | pure FUNCTION arth_i(first,increment,n) | |
| 57 | |||
| 58 | INTEGER, INTENT(IN) :: first,increment,n | ||
| 59 | INTEGER arth_i(n) | ||
| 60 | |||
| 61 | ! Local: | ||
| 62 | INTEGER :: k,k2,temp | ||
| 63 | |||
| 64 | !--------------------------------------- | ||
| 65 | |||
| 66 | ✗ | if (n > 0) arth_i(1)=first | |
| 67 | ✗ | if (n <= NPAR_ARTH) then | |
| 68 | ✗ | do k=2,n | |
| 69 | ✗ | arth_i(k)=arth_i(k-1)+increment | |
| 70 | end do | ||
| 71 | else | ||
| 72 | ✗ | do k=2,NPAR2_ARTH | |
| 73 | ✗ | arth_i(k)=arth_i(k-1)+increment | |
| 74 | end do | ||
| 75 | ✗ | temp=increment*NPAR2_ARTH | |
| 76 | k=NPAR2_ARTH | ||
| 77 | ✗ | do | |
| 78 | ✗ | if (k >= n) exit | |
| 79 | ✗ | k2=k+k | |
| 80 | ✗ | arth_i(k+1:min(k2,n))=temp+arth_i(1:min(k,n-k)) | |
| 81 | ✗ | temp=temp+temp | |
| 82 | ✗ | k=k2 | |
| 83 | end do | ||
| 84 | end if | ||
| 85 | |||
| 86 | END FUNCTION arth_i | ||
| 87 | |||
| 88 | END MODULE arth_m | ||
| 89 |