| 0,0 → 1,67 |
|---|
| #! /usr/bin/env python |
| """ fmath.prime_test |
| |
| Test for prime number routines |
| |
| Copyright © (c) 2002 by Paul A. Lambert |
| Read LICENSE.txt for license information. |
| """ |
| import unittest |
| from fmath.prime import rabin_miller |
| from string import split,atoi |
| |
| class primeTests(unittest.TestCase): |
| """ Tests for the rabin_miller probabilistic primality """ |
| def testKnowPrimes(self): |
| """ Test a few known primes """ |
| for prime in knownPrimes: |
| for testNum in (2,3,5,7,9,13,17,19,21,27,33,101,205): |
| assert( rabin_miller(prime, testNum) ), 'Prime: %d should pass rabin miller test using %d!'%(prime,testNum) |
| def testKnowComposites(self): |
| """ Test a few known primes """ |
| for nonPrime in knownComposites: |
| witness = 1 |
| for testNum in (2,3,5,7,9,13,17,19,21,27,33,101,205): |
| if not rabin_miller(nonPrime, testNum): |
| witness = 0 |
| #print witness, nonPrime, testNum |
| assert( not witness ), 'NonPrime: %d should Fail rabin miller tests!'%nonPrime |
| def test50K(self): |
| """ Test all composites below the first 50000 primes |
| print out the maxium witness number required """ |
| f=open('primes_1st_50k.txt','r') |
| s=f.read() |
| f.close() |
| ps=split(s) |
| primes_1st_50k = [atoi(i) for i in ps] |
| maxWitness = 2 |
| lastPrime = 2 |
| for prime in primes_1st_50k[1:]: |
| for composite in range(lastPrime+2,prime,2): |
| for witness in primes_1st_50k: |
| if not rabin_miller(composite, witness): |
| if witness > maxWitness: |
| maxWitness = witness |
| print 'composite, maxWitness',composite,maxWitness |
| lastPrime = prime |
| break |
| #assert(witness!=primes_1st_50K[-1]),'MR test failed!!!!!!!' |
| lastPrime = prime |
| print 'last tested',prime |
| |
| knownPrimes = [611693, 611707, 611729, 611753, 611791, |
| 611801, 611803, 611827, 611833, 611837, |
| 611839, 611873, 611879, 611887, 611903, |
| 611921, 611927, 611939, 611951, 611953] |
| |
| knownComposites = [611695, 611709, 611723, 611757, 611793, |
| 611805, 611807, 611829, 611835, 611800, |
| 611841, 611875, 611881, 611889, 611905, |
| 611923, 611925, 611941, 611947, 611937] |
| |
| # Make this test module runnable from the command prompt |
| if __name__ == "__main__": |
| unittest.main() |
| |
| |
| |