5 ! Test cases provided by Dominique d'Humieres <dominiq@lps.ens.fr>
6 ! and Michael Richmond <michael.a.richmond@nasa.gov>.
9 intrinsic :: dot_product
, sqrt
14 real, dimension(:), intent(in
) :: r
16 len
= sqrt(dot_product(r
,r
))
21 INTEGER, PARAMETER :: trantb(1,1) = RESHAPE((/1,2/), shape
=(/1,1/))
22 next_state
= trantb(1, 1)
23 END FUNCTION next_state
25 end module vector_calculus