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
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
! { dg-do run }
!
! Test functionality of pointer class arrays:
! ALLOCATE with source, ASSOCIATED, DEALLOCATE, passing as arguments for
! ELEMENTAL and non-ELEMENTAL procedures, SELECT TYPE and LOWER/UPPER.
!
type :: type1
integer :: i
end type
type, extends(type1) :: type2
real :: r
end type
class(type1), pointer, dimension (:) :: x
allocate(x(2), source = type2(42,42.0))
call display(x, [1], [2], t2 = [type2(42,42.0),type2(42,42.0)])
call display(x, [1], [2], t2 = [type2(111,99.0),type2(111,99.0)])
if (associated (x)) deallocate (x)
allocate(x(1:4), source = [(type2(i,42.0 + float (i)), i = 1, 4)])
call display(x, [1], [4], t2 = [(type2(i,42.0 + float (i)), i = 1, 4)])
call display(x, [1], [4], t2 = [(type2(111,99.0), i = 1, 4)])
if (any (disp (x) .ne. [99.0,99.0,99.0,99.0])) call abort
if (associated (x)) deallocate (x)
allocate(x(1:4), source = type1(42))
call display(x, [1], [4], t1 = [(type1(42), i = 1, 4)])
call display(x, [1], [4], t1 = [type1(42),type1(99),type1(42),type1(42)])
if (any (disp (x) .ne. [0.0,0.0,0.0,0.0])) call abort
if (associated (x)) deallocate (x)
contains
subroutine display(x, lower, upper, t1, t2)
class(type1), pointer, dimension (:) :: x
integer, dimension (:) :: lower, upper
type(type1), optional, dimension(:) :: t1
type(type2), optional, dimension(:) :: t2
select type (x)
type is (type1)
if (present (t1)) then
if (any (x%i .ne. t1%i)) call abort
else
call abort
end if
x(2)%i = 99
type is (type2)
if (present (t2)) then
if (any (x%i .ne. t2%i)) call abort
if (any (x%r .ne. t2%r)) call abort
else
call abort
end if
x%i = 111
x%r = 99.0
end select
call bounds (x, lower, upper)
end subroutine
subroutine bounds (x, lower, upper)
class(type1), pointer, dimension (:) :: x
integer, dimension (:) :: lower, upper
if (any (lower .ne. lbound (x))) call abort
if (any (upper .ne. ubound (x))) call abort
end subroutine
elemental function disp(y) result(ans)
class(type1), intent(in) :: y
real :: ans
select type (y)
type is (type1)
ans = 0.0
type is (type2)
ans = y%r
end select
end function
end
|
