| 1234567891011121314151617181920212223242526272829303132333435 |
- '''
- generating a test case for sigmas in pallas field.
- the output in csv format: f,Sigma,sigma1,sigma2 respectively for 2-term approximation.
- '''
- from lottery import *
- import numpy as np
- SEP=','
- def calc_sigmas(f, Sigma):
- k=N_TERM
- x = (1-f)
- c = math.log(x)
- neg_c = -1*c
- sigmas = [(int((neg_c/Sigma)**i * (L/fact(i)))) for i in range(1, k+1)]
- return sigmas
- sigmas = calc_sigmas(0.5, 1000)
- assert(len(sigmas)==2)
- with open("pallas_unittests.csv", 'w') as file:
- buf = ''
- for f in np.arange(0.01, 0.99, 10):
- for total_stake in np.arange(100, 1000, 10):
- sigmas = calc_sigmas(f, total_stake)
- line=str(f) + SEP + \
- str(total_stake) + SEP + \
- '{:x}'.format(sigmas[0]) + SEP + \
- '{:x}'.format(sigmas[1]) + '\n'
- buf+=line
- file.write(buf)
|
