1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
|
program implied_do_loops6
real, allocatable :: Q(:)
real :: R(5)
integer :: rank
allocate(Q(4))
Q = 5.0
Q([1, 3]) = -4.0
print *, Q
R = [abs(Q), 0.0_4]
print *, R
if( any(R /= [4.0, 5.0, 4.0, 5.0, 0.0]) ) error stop
rank = maxval([0_4, trueloc(Q > 0)])
print *, rank
if( rank /= 4 ) error stop
contains
function trueloc(x) result(loc)
logical, intent(in) :: x(:)
integer(4), allocatable :: loc(:)
integer :: i, j
allocate(loc(count(x)))
j = 1
do i = 1, size(x)
if( x(i) ) then
loc(j) = i
j = j + 1
end if
end do
end function trueloc
end program
|