@@ -12,59 +12,67 @@ def sample_seed(seed):
1212 subtraction = abs (seed )
1313 mj = MSEED - subtraction
1414
15- sa0_42 = (mj - 1 ) % MBIG
16- sa0__8 = (1 - sa0_42 ) % MBIG
17- sa0_29 = (sa0_42 - sa0__8 ) % MBIG
18- sa0_50 = (sa0__8 - sa0_29 ) % MBIG
19- sa0_16 = (sa0_29 - sa0_50 ) % MBIG
20- sa0_37 = (sa0_50 - sa0_16 ) % MBIG
21- sa0__3 = (sa0_16 - sa0_37 ) % MBIG
22- sa0_24 = (sa0_37 - sa0__3 ) % MBIG
23- sa0_45 = (sa0__3 - sa0_24 ) % MBIG
24- sa0_11 = (sa0_24 - sa0_45 ) % MBIG
25- sa0_32 = (sa0_45 - sa0_11 ) % MBIG
26- sa0_53 = (sa0_11 - sa0_32 ) % MBIG
27- sa0_19 = (sa0_32 - sa0_53 ) % MBIG
28- sa0_40 = (sa0_53 - sa0_19 ) % MBIG
29- sa0__6 = (sa0_19 - sa0_40 ) % MBIG
30- sa0_27 = (sa0_40 - sa0__6 ) % MBIG
31- sa0_48 = (sa0__6 - sa0_27 ) % MBIG
32- sa0_14 = (sa0_27 - sa0_48 ) % MBIG
33- sa0_35 = (sa0_48 - sa0_14 ) % MBIG
34- sa0__1 = (sa0_14 - sa0_35 ) % MBIG
35- sa0_22 = (sa0_35 - sa0__1 ) % MBIG
36- sa0_43 = (sa0__1 - sa0_22 ) % MBIG
37- sa0__9 = (sa0_22 - sa0_43 ) % MBIG
38- sa0_30 = (sa0_43 - sa0__9 ) % MBIG
39- sa0_51 = (sa0__9 - sa0_30 ) % MBIG
40- sa0_17 = (sa0_30 - sa0_51 ) % MBIG
41- sa0_38 = (sa0_51 - sa0_17 ) % MBIG
42- sa0__4 = (sa0_17 - sa0_38 ) % MBIG
43- sa0_25 = (sa0_38 - sa0__4 ) % MBIG
44- sa0_46 = (sa0__4 - sa0_25 ) % MBIG
45- sa0_12 = (sa0_25 - sa0_46 ) % MBIG
46- sa0_33 = (sa0_46 - sa0_12 ) % MBIG
47- sa0_54 = (sa0_12 - sa0_33 ) % MBIG
48- sa0_20 = (sa0_33 - sa0_54 ) % MBIG
49- sa0_41 = (sa0_54 - sa0_20 ) % MBIG
50- sa0__7 = (sa0_20 - sa0_41 ) % MBIG
51- sa0_28 = (sa0_41 - sa0__7 ) % MBIG
52- sa0_49 = (sa0__7 - sa0_28 ) % MBIG
53- sa0_15 = (sa0_28 - sa0_49 ) % MBIG
54- sa0_36 = (sa0_49 - sa0_15 ) % MBIG
55- sa0__2 = (sa0_15 - sa0_36 ) % MBIG
56- sa0_23 = (sa0_36 - sa0__2 ) % MBIG
57- sa0_44 = (sa0__2 - sa0_23 ) % MBIG
58- sa0_10 = (sa0_23 - sa0_44 ) % MBIG
59- sa0_31 = (sa0_44 - sa0_10 ) % MBIG
60- sa0_52 = (sa0_10 - sa0_31 ) % MBIG
61- sa0_18 = (sa0_31 - sa0_52 ) % MBIG
62- sa0_39 = (sa0_52 - sa0_18 ) % MBIG
63- sa0__5 = (sa0_18 - sa0_39 ) % MBIG
64- sa0_26 = (sa0_39 - sa0__5 ) % MBIG
65- sa0_47 = (sa0__5 - sa0_26 ) % MBIG
66- sa0_13 = (sa0_26 - sa0_47 ) % MBIG
67- sa0_34 = (sa0_47 - sa0_13 ) % MBIG
15+ # note: coefecents are fibinocci numbers
16+ # a = 5
17+ # b = 3
18+ # for i in range(5, 55):
19+ # v = getvarname(i) # sa0_50, sa0_16,..
20+ # aa = (a * ((-1)**i) * -1) % MBIG
21+ # bb = (b * (-1)**i) % MBIG
22+ # print(f"{v} = ({aa} + {bb}*mj) % MBIG")
23+ # c = a
24+ # a += b
25+ # b = c
26+ sa0_50 = (5 + 2147483644 * mj ) % MBIG
27+ sa0_16 = (2147483639 + 5 * mj ) % MBIG
28+ sa0_37 = (13 + 2147483639 * mj ) % MBIG
29+ sa0__3 = (2147483626 + 13 * mj ) % MBIG
30+ sa0_24 = (34 + 2147483626 * mj ) % MBIG
31+ sa0_45 = (2147483592 + 34 * mj ) % MBIG
32+ sa0_11 = (89 + 2147483592 * mj ) % MBIG
33+ sa0_32 = (2147483503 + 89 * mj ) % MBIG
34+ sa0_53 = (233 + 2147483503 * mj ) % MBIG
35+ sa0_19 = (2147483270 + 233 * mj ) % MBIG
36+ sa0_40 = (610 + 2147483270 * mj ) % MBIG
37+ sa0__6 = (2147482660 + 610 * mj ) % MBIG
38+ sa0_27 = (1597 + 2147482660 * mj ) % MBIG
39+ sa0_48 = (2147481063 + 1597 * mj ) % MBIG
40+ sa0_14 = (4181 + 2147481063 * mj ) % MBIG
41+ sa0_35 = (2147476882 + 4181 * mj ) % MBIG
42+ sa0__1 = (10946 + 2147476882 * mj ) % MBIG
43+ sa0_22 = (2147465936 + 10946 * mj ) % MBIG
44+ sa0_43 = (28657 + 2147465936 * mj ) % MBIG
45+ sa0__9 = (2147437279 + 28657 * mj ) % MBIG
46+ sa0_30 = (75025 + 2147437279 * mj ) % MBIG
47+ sa0_51 = (2147362254 + 75025 * mj ) % MBIG
48+ sa0_17 = (196418 + 2147362254 * mj ) % MBIG
49+ sa0_38 = (2147165836 + 196418 * mj ) % MBIG
50+ sa0__4 = (514229 + 2147165836 * mj ) % MBIG
51+ sa0_25 = (2146651607 + 514229 * mj ) % MBIG
52+ sa0_46 = (1346269 + 2146651607 * mj ) % MBIG
53+ sa0_12 = (2145305338 + 1346269 * mj ) % MBIG
54+ sa0_33 = (3524578 + 2145305338 * mj ) % MBIG
55+ sa0_54 = (2141780760 + 3524578 * mj ) % MBIG
56+ sa0_20 = (9227465 + 2141780760 * mj ) % MBIG
57+ sa0_41 = (2132553295 + 9227465 * mj ) % MBIG
58+ sa0__7 = (24157817 + 2132553295 * mj ) % MBIG
59+ sa0_28 = (2108395478 + 24157817 * mj ) % MBIG
60+ sa0_49 = (63245986 + 2108395478 * mj ) % MBIG
61+ sa0_15 = (2045149492 + 63245986 * mj ) % MBIG
62+ sa0_36 = (165580141 + 2045149492 * mj ) % MBIG
63+ sa0__2 = (1879569351 + 165580141 * mj ) % MBIG
64+ sa0_23 = (433494437 + 1879569351 * mj ) % MBIG
65+ sa0_44 = (1446074914 + 433494437 * mj ) % MBIG
66+ sa0_10 = (1134903170 + 1446074914 * mj ) % MBIG
67+ sa0_31 = (311171744 + 1134903170 * mj ) % MBIG
68+ sa0_52 = (823731426 + 311171744 * mj ) % MBIG
69+ sa0_18 = (1634923965 + 823731426 * mj ) % MBIG
70+ sa0_39 = (1336291108 + 1634923965 * mj ) % MBIG
71+ sa0__5 = (298632857 + 1336291108 * mj ) % MBIG
72+ sa0_26 = (1037658251 + 298632857 * mj ) % MBIG
73+ sa0_47 = (1408458253 + 1037658251 * mj ) % MBIG
74+ sa0_13 = (1776683645 + 1408458253 * mj ) % MBIG
75+ sa0_34 = (1779258255 + 1776683645 * mj ) % MBIG
6876
6977
7078 sa1__1 = (sa0__1 - sa0_32 ) % MBIG
@@ -74,18 +82,18 @@ def sample_seed(seed):
7482 sa1__5 = (sa0__5 - sa0_36 ) % MBIG
7583 sa1__6 = (sa0__6 - sa0_37 ) % MBIG
7684 sa1__7 = (sa0__7 - sa0_38 ) % MBIG
77- sa1__8 = (sa0__8 - sa0_39 ) % MBIG
85+ sa1__8 = (( 1 - ( mj - 1 )) - sa0_39 ) % MBIG
7886 sa1__9 = (sa0__9 - sa0_40 ) % MBIG
7987 sa1_10 = (sa0_10 - sa0_41 ) % MBIG
80- sa1_11 = (sa0_11 - sa0_42 ) % MBIG
88+ sa1_11 = (sa0_11 - ( mj - 1 ) ) % MBIG
8189 sa1_12 = (sa0_12 - sa0_43 ) % MBIG
8290 sa1_13 = (sa0_13 - sa0_44 ) % MBIG
8391 sa1_14 = (sa0_14 - sa0_45 ) % MBIG
8492 sa1_15 = (sa0_15 - sa0_46 ) % MBIG
8593 sa1_16 = (sa0_16 - sa0_47 ) % MBIG
8694 sa1_17 = (sa0_17 - sa0_48 ) % MBIG
8795 sa1_18 = (sa0_18 - sa0_49 ) % MBIG
88- sa1_19 = (sa0_19 - sa0_50 ) % MBIG
96+ sa1_19 = (sa0_19 - ( 5 - 3 * mj ) ) % MBIG
8997 sa1_20 = (sa0_20 - sa0_51 ) % MBIG
9098 sa1_21 = (1 - sa0_52 ) % MBIG
9199 sa1_22 = (sa0_22 - sa0_53 ) % MBIG
@@ -95,7 +103,7 @@ def sample_seed(seed):
95103 sa1_26 = (sa0_26 - sa1__2 ) % MBIG
96104 sa1_27 = (sa0_27 - sa1__3 ) % MBIG
97105 sa1_28 = (sa0_28 - sa1__4 ) % MBIG
98- sa1_29 = (sa0_29 - sa1__5 ) % MBIG
106+ sa1_29 = (( 2 * mj - 3 ) - sa1__5 ) % MBIG
99107 sa1_30 = (sa0_30 - sa1__6 ) % MBIG
100108 sa1_31 = (sa0_31 - sa1__7 ) % MBIG
101109 sa1_32 = (sa0_32 - sa1__8 ) % MBIG
@@ -108,15 +116,15 @@ def sample_seed(seed):
108116 sa1_39 = (sa0_39 - sa1_15 ) % MBIG
109117 sa1_40 = (sa0_40 - sa1_16 ) % MBIG
110118 sa1_41 = (sa0_41 - sa1_17 ) % MBIG
111- sa1_42 = (sa0_42 - sa1_18 ) % MBIG
119+ sa1_42 = (( mj - 1 ) - sa1_18 ) % MBIG
112120 sa1_43 = (sa0_43 - sa1_19 ) % MBIG
113121 sa1_44 = (sa0_44 - sa1_20 ) % MBIG
114122 sa1_45 = (sa0_45 - sa1_21 ) % MBIG
115123 sa1_46 = (sa0_46 - sa1_22 ) % MBIG
116124 sa1_47 = (sa0_47 - sa1_23 ) % MBIG
117125 sa1_48 = (sa0_48 - sa1_24 ) % MBIG
118126 sa1_49 = (sa0_49 - sa1_25 ) % MBIG
119- sa1_50 = (sa0_50 - sa1_26 ) % MBIG
127+ sa1_50 = (( 5 - 3 * mj ) - sa1_26 ) % MBIG
120128 sa1_51 = (sa0_51 - sa1_27 ) % MBIG
121129 sa1_52 = (sa0_52 - sa1_28 ) % MBIG
122130 sa1_53 = (sa0_53 - sa1_29 ) % MBIG
0 commit comments