41
42
43
44
45
46
47
48 INTEGER I(2), J(2), K(2)
49
50
51
52
53
54 INTEGER IPOW15, IPOW16, IPOW30
55 parameter( ipow15=2**15, ipow16=2**16, ipow30=2**30 )
56
57
58 INTEGER KT, LT
59
60
61 INTRINSIC mod
62
63
64
65 kt = k(1)*j(1)
66 IF( kt.LT.0 ) kt = (kt+ipow30) + ipow30
67 i(1) = mod(kt,ipow16)
68 lt = k(1)*j(2) + k(2)*j(1)
69 IF( lt.LT.0 ) lt = (lt+ipow30) + ipow30
70 kt = kt/ipow16 + lt
71 IF( kt.LT.0 ) kt = (kt+ipow30) + ipow30
72 i(2) = mod( kt, ipow15 )
73
74 RETURN
75
76
77