summaryrefslogtreecommitdiff
path: root/gcc/testsuite/gfortran.dg/transpose_conjg_1.f90
blob: 3b28827b38ab2f5aca9ae73cac1e4a37b7fcdedc (plain)
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)) call abort ()
!
  d = matmul (b,  a )
  if (any (d /= matmul (transpose(conjg(a)), a))) call abort ()
  if (any (d /= matmul (conjg(transpose(a)), a))) call abort ()
!
  c = (0.0,1.0)
  b = conjg(transpose(a + c))
  d = transpose(conjg(a + c))
  if (any (b /= d)) call abort ()
!
  d = matmul (b,  a + c)
  if (any (d /= matmul (transpose(conjg(a + c)), a + c))) call abort ()
  if (any (d /= matmul (conjg(transpose(a + c)), a + c))) call abort ()
 END program main