[6] | 1 | module m_elstack |
---|
| 2 | |
---|
| 3 | use m_buffer |
---|
| 4 | |
---|
| 5 | private |
---|
| 6 | |
---|
| 7 | ! |
---|
| 8 | ! Simple stack to keep track of which elements have appeared so far |
---|
| 9 | ! |
---|
| 10 | integer, parameter, private :: STACK_SIZE = 40 |
---|
| 11 | |
---|
| 12 | type, public :: elstack_t |
---|
| 13 | private |
---|
| 14 | integer :: n_items |
---|
| 15 | type(buffer_t), dimension(STACK_SIZE) :: data |
---|
| 16 | end type elstack_t |
---|
| 17 | |
---|
| 18 | public :: push_elstack, pop_elstack, reset_elstack, print_elstack |
---|
| 19 | public :: init_elstack |
---|
| 20 | public :: get_top_elstack, is_empty, get_elstack_signature |
---|
| 21 | |
---|
| 22 | interface is_empty |
---|
| 23 | module procedure is_empty_elstack |
---|
| 24 | end interface |
---|
| 25 | private :: is_empty_elstack |
---|
| 26 | |
---|
| 27 | CONTAINS |
---|
| 28 | |
---|
| 29 | !----------------------------------------------------------------- |
---|
| 30 | subroutine init_elstack(elstack) |
---|
| 31 | type(elstack_t), intent(inout) :: elstack |
---|
| 32 | |
---|
| 33 | integer :: i |
---|
| 34 | |
---|
| 35 | elstack%n_items = 0 |
---|
| 36 | do i = 1, STACK_SIZE ! to avoid "undefined status" |
---|
| 37 | call init_buffer(elstack%data(i)) |
---|
| 38 | enddo |
---|
| 39 | end subroutine init_elstack |
---|
| 40 | |
---|
| 41 | !----------------------------------------------------------------- |
---|
| 42 | subroutine reset_elstack(elstack) |
---|
| 43 | type(elstack_t), intent(inout) :: elstack |
---|
| 44 | |
---|
| 45 | integer :: i |
---|
| 46 | |
---|
| 47 | elstack%n_items = 0 |
---|
| 48 | do i = 1, STACK_SIZE |
---|
| 49 | call reset_buffer(elstack%data(i)) |
---|
| 50 | enddo |
---|
| 51 | end subroutine reset_elstack |
---|
| 52 | |
---|
| 53 | !----------------------------------------------------------------- |
---|
| 54 | function is_empty_elstack(elstack) result(answer) |
---|
| 55 | type(elstack_t), intent(in) :: elstack |
---|
| 56 | logical :: answer |
---|
| 57 | |
---|
| 58 | answer = (elstack%n_items == 0) |
---|
| 59 | end function is_empty_elstack |
---|
| 60 | |
---|
| 61 | !----------------------------------------------------------------- |
---|
| 62 | subroutine push_elstack(item,elstack) |
---|
| 63 | type(buffer_t), intent(in) :: item |
---|
| 64 | type(elstack_t), intent(inout) :: elstack |
---|
| 65 | |
---|
| 66 | integer :: n |
---|
| 67 | |
---|
| 68 | n = elstack%n_items |
---|
| 69 | if (n == STACK_SIZE) then |
---|
| 70 | stop "*Element stack full" |
---|
| 71 | endif |
---|
| 72 | n = n + 1 |
---|
| 73 | elstack%data(n) = item |
---|
| 74 | elstack%n_items = n |
---|
| 75 | |
---|
| 76 | end subroutine push_elstack |
---|
| 77 | |
---|
| 78 | !----------------------------------------------------------------- |
---|
| 79 | subroutine pop_elstack(elstack,item) |
---|
| 80 | type(elstack_t), intent(inout) :: elstack |
---|
| 81 | type(buffer_t), intent(out) :: item |
---|
| 82 | |
---|
| 83 | ! |
---|
| 84 | ! We assume the elstack is not empty... (the user has called is_empty first) |
---|
| 85 | ! |
---|
| 86 | integer :: n |
---|
| 87 | |
---|
| 88 | n = elstack%n_items |
---|
| 89 | if (n == 0) then |
---|
| 90 | stop "*********Element stack empty" |
---|
| 91 | endif |
---|
| 92 | item = elstack%data(n) |
---|
| 93 | elstack%n_items = n - 1 |
---|
| 94 | |
---|
| 95 | end subroutine pop_elstack |
---|
| 96 | |
---|
| 97 | !----------------------------------------------------------------- |
---|
| 98 | subroutine get_top_elstack(elstack,item) |
---|
| 99 | ! |
---|
| 100 | ! Get the top element of the stack, *without popping it*. |
---|
| 101 | ! |
---|
| 102 | type(elstack_t), intent(in) :: elstack |
---|
| 103 | type(buffer_t), intent(out) :: item |
---|
| 104 | |
---|
| 105 | ! |
---|
| 106 | ! We assume the elstack is not empty... (the user has called is_empty first) |
---|
| 107 | ! |
---|
| 108 | integer :: n |
---|
| 109 | |
---|
| 110 | n = elstack%n_items |
---|
| 111 | if (n == 0) then |
---|
| 112 | stop "*********Element stack empty" |
---|
| 113 | endif |
---|
| 114 | item = elstack%data(n) |
---|
| 115 | |
---|
| 116 | end subroutine get_top_elstack |
---|
| 117 | |
---|
| 118 | !----------------------------------------------------------------- |
---|
| 119 | subroutine print_elstack(elstack,unit) |
---|
| 120 | type(elstack_t), intent(in) :: elstack |
---|
| 121 | integer, intent(in) :: unit |
---|
| 122 | integer :: i |
---|
| 123 | |
---|
| 124 | do i = elstack%n_items, 1, -1 |
---|
| 125 | write(unit=unit,fmt=*) str(elstack%data(i)) |
---|
| 126 | enddo |
---|
| 127 | |
---|
| 128 | end subroutine print_elstack |
---|
| 129 | |
---|
| 130 | !------------------------------------------------------------- |
---|
| 131 | subroutine get_elstack_signature(elstack,string) |
---|
| 132 | type(elstack_t), intent(in) :: elstack |
---|
| 133 | character(len=*), intent(out) :: string |
---|
| 134 | integer :: i, length, j |
---|
| 135 | |
---|
| 136 | string = "" |
---|
| 137 | j = 0 |
---|
| 138 | do i = 1, elstack%n_items |
---|
| 139 | length = len(elstack%data(i)) |
---|
| 140 | string(j+1:j+1) = "/" |
---|
| 141 | j = j+1 |
---|
| 142 | string(j+1:j+length) = str(elstack%data(i)) |
---|
| 143 | j = j + length |
---|
| 144 | enddo |
---|
| 145 | |
---|
| 146 | end subroutine get_elstack_signature |
---|
| 147 | |
---|
| 148 | end module m_elstack |
---|
| 149 | |
---|
| 150 | |
---|
| 151 | |
---|
| 152 | |
---|
| 153 | |
---|