vecfor_R16P-doctest-23.f90 Source File

program volatile_doctest
use vecfor_R16P
 type(vector_R16P) :: pt(0:3)

 pt(0) = -1 * ey_R16P
 pt(1) = 1 * ex_R16P
 pt(2) = 1 * ey_R16P
 pt(3) = -1 * ex_R16P
 print "(L1)", is_concyclic_R16P(pt(0), pt1=pt(1), pt2=pt(2), pt3=pt(3))
endprogram volatile_doctest