@@ -959,6 +959,44 @@ static const char *medium_primes[10] = {
959
959
"8aOicQ5gIbFCarFUgSgzh40LpuZ0jjK1u48/YT+C0h1dAQ8CIEgZjHZT+5/7cCRGmJlo+XCp7S41MSQ2ZNRSJh2texRYtvAXBAZfR8A8twl316P"
960
960
};
961
961
962
+ const mp_digit prime_tab [] = {
963
+ 0x0002 , 0x0003 , 0x0005 , 0x0007 , 0x000B , 0x000D , 0x0011 , 0x0013 ,
964
+ 0x0017 , 0x001D , 0x001F , 0x0025 , 0x0029 , 0x002B , 0x002F , 0x0035 ,
965
+ 0x003B , 0x003D , 0x0043 , 0x0047 , 0x0049 , 0x004F , 0x0053 , 0x0059 ,
966
+ 0x0061 , 0x0065 , 0x0067 , 0x006B , 0x006D , 0x0071 , 0x007F , 0x0083 ,
967
+ 0x0089 , 0x008B , 0x0095 , 0x0097 , 0x009D , 0x00A3 , 0x00A7 , 0x00AD ,
968
+ 0x00B3 , 0x00B5 , 0x00BF , 0x00C1 , 0x00C5 , 0x00C7 , 0x00D3 , 0x00DF ,
969
+ 0x00E3 , 0x00E5 , 0x00E9 , 0x00EF , 0x00F1 , 0x00FB , 0x0101 , 0x0107 ,
970
+ 0x010D , 0x010F , 0x0115 , 0x0119 , 0x011B , 0x0125 , 0x0133 , 0x0137 ,
971
+
972
+ 0x0139 , 0x013D , 0x014B , 0x0151 , 0x015B , 0x015D , 0x0161 , 0x0167 ,
973
+ 0x016F , 0x0175 , 0x017B , 0x017F , 0x0185 , 0x018D , 0x0191 , 0x0199 ,
974
+ 0x01A3 , 0x01A5 , 0x01AF , 0x01B1 , 0x01B7 , 0x01BB , 0x01C1 , 0x01C9 ,
975
+ 0x01CD , 0x01CF , 0x01D3 , 0x01DF , 0x01E7 , 0x01EB , 0x01F3 , 0x01F7 ,
976
+ 0x01FD , 0x0209 , 0x020B , 0x021D , 0x0223 , 0x022D , 0x0233 , 0x0239 ,
977
+ 0x023B , 0x0241 , 0x024B , 0x0251 , 0x0257 , 0x0259 , 0x025F , 0x0265 ,
978
+ 0x0269 , 0x026B , 0x0277 , 0x0281 , 0x0283 , 0x0287 , 0x028D , 0x0293 ,
979
+ 0x0295 , 0x02A1 , 0x02A5 , 0x02AB , 0x02B3 , 0x02BD , 0x02C5 , 0x02CF ,
980
+
981
+ 0x02D7 , 0x02DD , 0x02E3 , 0x02E7 , 0x02EF , 0x02F5 , 0x02F9 , 0x0301 ,
982
+ 0x0305 , 0x0313 , 0x031D , 0x0329 , 0x032B , 0x0335 , 0x0337 , 0x033B ,
983
+ 0x033D , 0x0347 , 0x0355 , 0x0359 , 0x035B , 0x035F , 0x036D , 0x0371 ,
984
+ 0x0373 , 0x0377 , 0x038B , 0x038F , 0x0397 , 0x03A1 , 0x03A9 , 0x03AD ,
985
+ 0x03B3 , 0x03B9 , 0x03C7 , 0x03CB , 0x03D1 , 0x03D7 , 0x03DF , 0x03E5 ,
986
+ 0x03F1 , 0x03F5 , 0x03FB , 0x03FD , 0x0407 , 0x0409 , 0x040F , 0x0419 ,
987
+ 0x041B , 0x0425 , 0x0427 , 0x042D , 0x043F , 0x0443 , 0x0445 , 0x0449 ,
988
+ 0x044F , 0x0455 , 0x045D , 0x0463 , 0x0469 , 0x047F , 0x0481 , 0x048B ,
989
+
990
+ 0x0493 , 0x049D , 0x04A3 , 0x04A9 , 0x04B1 , 0x04BD , 0x04C1 , 0x04C7 ,
991
+ 0x04CD , 0x04CF , 0x04D5 , 0x04E1 , 0x04EB , 0x04FD , 0x04FF , 0x0503 ,
992
+ 0x0509 , 0x050B , 0x0511 , 0x0515 , 0x0517 , 0x051B , 0x0527 , 0x0529 ,
993
+ 0x052F , 0x0551 , 0x0557 , 0x055D , 0x0565 , 0x0577 , 0x0581 , 0x058F ,
994
+ 0x0593 , 0x0595 , 0x0599 , 0x059F , 0x05A7 , 0x05AB , 0x05AD , 0x05B3 ,
995
+ 0x05BF , 0x05C9 , 0x05CB , 0x05CF , 0x05D1 , 0x05D5 , 0x05DB , 0x05E7 ,
996
+ 0x05F3 , 0x05FB , 0x0607 , 0x060D , 0x0611 , 0x0617 , 0x061F , 0x0623 ,
997
+ 0x062B , 0x062F , 0x063D , 0x0641 , 0x0647 , 0x0649 , 0x064D , 0x0653
998
+ };
999
+
962
1000
#define ARR_LENGTH (a ) ((int)(sizeof((a))/sizeof((a)[0])))
963
1001
964
1002
static int test_mp_prime_miller_rabin (void )
@@ -1026,7 +1064,7 @@ static int test_mp_prime_miller_rabin(void)
1026
1064
DO (mp_read_radix (& a , SPSP_2_100_LARGE [i ], 64 ));
1027
1065
for (j = 2u ; j <= (mp_digit )mp_prime_rabin_miller_trials (mp_count_bits (& a )); j ++ ) {
1028
1066
result = false;
1029
- mp_set (& b , (mp_digit )s_mp_prime_tab [j ]);
1067
+ mp_set (& b , (mp_digit )prime_tab [j ]);
1030
1068
DO (mp_prime_miller_rabin (& a , & b , & result ));
1031
1069
}
1032
1070
/* These numbers are not big enough for the heuristics to work */
0 commit comments