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() 12END 13