1!  { dg-do run }
2!  Tests the fix for PR35740, where the trick of interchanging the descriptor
3!  dimensions to implement TRANSPOSE did not work if it is an argument of
4!  an elemental function - eg. CONJG.  The fix forces a library call for such
5!  cases.  During the diagnosis of the PR, it was found that the scalarizer was
6!  completely thrown if the argument of TRANSPOSE was a non-variable
7!  expression; eg a + c below.  This is also fixed by the library call.
8!
9!  Contributed by Dominik Muth <dominik.muth@gmx.de>
10!
11program main
12  implicit none
13  complex, dimension(2,2) :: a,b,c,d
14  a(1,1) = (1.,1.)
15  a(2,1) = (2.,2.)
16  a(1,2) = (3.,3.)
17  a(2,2) = (4.,4.)
18!
19  b = a
20  b = conjg(transpose(b))
21  d = a
22  d = transpose(conjg(d))
23  if (any (b /= d)) call abort ()
24!
25  d = matmul (b,  a )
26  if (any (d /= matmul (transpose(conjg(a)), a))) call abort ()
27  if (any (d /= matmul (conjg(transpose(a)), a))) call abort ()
28!
29  c = (0.0,1.0)
30  b = conjg(transpose(a + c))
31  d = transpose(conjg(a + c))
32  if (any (b /= d)) call abort ()
33!
34  d = matmul (b,  a + c)
35  if (any (d /= matmul (transpose(conjg(a + c)), a + c))) call abort ()
36  if (any (d /= matmul (conjg(transpose(a + c)), a + c))) call abort ()
37 END program main
38