gcc/testsuite/gfortran.dg/matmul_8.f03

   1 ! { dg-do run }
   2 ! Transformational intrinsic MATMUL as initialization expression.
   3
   4   REAL, PARAMETER :: PI = 3.141592654, theta = PI/6.0
   5
   6   REAL, PARAMETER :: unity(2,2) = RESHAPE([1.0, 0.0, 0.0, 1.0], [2, 2])
   7   REAL, PARAMETER :: m1(2,2)    = RESHAPE([COS(theta), SIN(theta), -SIN(theta), COS(theta)], [2, 2])
   8   REAL, PARAMETER :: m2(2,2)    = RESHAPE([COS(theta), -SIN(theta), SIN(theta), COS(theta)], [2, 2])
   9   REAL, PARAMETER :: m(2,2)     = MATMUL(m1, m2)
  10
  11   IF (ANY(ABS(m - unity) > EPSILON(0.0))) CALL abort()
  12 END