| ! Really test where inside forall with temporary |
| program evil_where |
| implicit none |
| type t |
| logical valid |
| integer :: s |
| integer, dimension(:), pointer :: p |
| end type |
| type (t), dimension (5) :: v |
| integer i |
| |
| allocate (v(1)%p(2)) |
| allocate (v(2)%p(8)) |
| v(3)%p => NULL() |
| allocate (v(4)%p(8)) |
| allocate (v(5)%p(2)) |
| |
| v(:)%valid = (/.true., .true., .false., .true., .true./) |
| v(:)%s = (/1, 8, 999, 6, 2/) |
| v(1)%p(:) = (/9, 10/) |
| v(2)%p(:) = (/1, 2, 3, 4, 5, 6, 7, 8/) |
| v(4)%p(:) = (/13, 14, 15, 16, 17, 18, 19, 20/) |
| v(5)%p(:) = (/11, 12/) |
| |
| forall (i=1:5,v(i)%valid) |
| where (v(i)%p(1:v(i)%s).gt.4) |
| v(i)%p(1:v(i)%s) = v(6-i)%p(1:v(i)%s) |
| end where |
| end forall |
| |
| if (any(v(1)%p(:) .ne. (/11, 10/))) call abort |
| if (any(v(2)%p(:) .ne. (/1, 2, 3, 4, 17, 18, 19, 20/))) call abort |
| if (any(v(4)%p(:) .ne. (/1, 2, 3, 4, 5, 6, 19, 20/))) call abort |
| if (any(v(5)%p(:) .ne. (/9, 10/))) call abort |
| |
| v(1)%p(:) = (/9, 10/) |
| v(2)%p(:) = (/1, 2, 3, 4, 5, 6, 7, 8/) |
| v(4)%p(:) = (/13, 14, 15, 16, 17, 18, 19, 20/) |
| v(5)%p(:) = (/11, 12/) |
| |
| forall (i=1:5,v(i)%valid) |
| where (v(i)%p(1:v(i)%s).le.4) |
| v(i)%p(1:v(i)%s) = v(6-i)%p(1:v(i)%s) |
| end where |
| end forall |
| |
| if (any(v(1)%p(:) .ne. (/9, 10/))) call abort |
| if (any(v(2)%p(:) .ne. (/13, 14, 15, 16, 5, 6, 7, 8/))) call abort |
| if (any(v(4)%p(:) .ne. (/13, 14, 15, 16, 17, 18, 19, 20/))) call abort |
| if (any(v(5)%p(:) .ne. (/11, 12/))) call abort |
| |
| ! I should really free the memory I've allocated. |
| end program |
