view gcc/testsuite/gfortran.dg/transpose_conjg_1.f90 @ 158:494b0b89df80 default tip

...
author Shinji KONO <kono@ie.u-ryukyu.ac.jp>
date Mon, 25 May 2020 18:13:55 +0900
parents 84e7813d76e9
children
line wrap: on
line source

!  { 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