2 ! Testcase from PR 25396: User defined operators returning arrays.
7 interface operator(.cross
.)
13 ! Cross product between two 3d vectors.
14 pure
function cross(a
, b
)
15 real, dimension(3), intent(in
) :: a
,b
16 real, dimension(3) :: cross
18 cross
= (/ a(2) * b(3) - a(3) * b(2), &
19 a(3) * b(1) - a(1) * b(3), &
20 a(1) * b(2) - a(2) * b(1) /)
32 a
= dot_product (t(:,1), t(:,2) .cross
. t(:,3))