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
|
! { dg-do run }
! Tests the fix for PR35740, where the trick of interchanging the descriptor
! dimensions to implement TRANSPOSE did not work if it is an argument of
! an elemental function - eg. CONJG. The fix forces a library call for such
! cases. During the diagnosis of the PR, it was found that the scalarizer was
! completely thrown if the argument of TRANSPOSE was a non-variable
! expression; eg a + c below. This is also fixed by the library call.
!
! Contributed by Dominik Muth <dominik.muth@gmx.de>
!
program main
implicit none
complex, dimension(2,2) :: a,b,c,d
a(1,1) = (1.,1.)
a(2,1) = (2.,2.)
a(1,2) = (3.,3.)
a(2,2) = (4.,4.)
!
b = a
b = conjg(transpose(b))
d = a
d = transpose(conjg(d))
if (any (b /= d)) STOP 1
!
d = matmul (b, a )
if (any (d /= matmul (transpose(conjg(a)), a))) STOP 2
if (any (d /= matmul (conjg(transpose(a)), a))) STOP 3
!
c = (0.0,1.0)
b = conjg(transpose(a + c))
d = transpose(conjg(a + c))
if (any (b /= d)) STOP 4
!
d = matmul (b, a + c)
if (any (d /= matmul (transpose(conjg(a + c)), a + c))) STOP 5
if (any (d /= matmul (conjg(transpose(a + c)), a + c))) STOP 6
END program main
|
