lacau/Projet_Boites
2
0
Fork 0
Code Issues Pull requests Releases Wiki Activity

chore: clean up "zero padding"

Browse source
This commit is contained in:
Paul ALNET 2023-06-04 08:11:26 +02:00
parent 7bee845a97
commit 0cdc13b869
1 changed files with 4 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
Probas.py
View file

@ -88,11 +88,13 @@ def stats_NFBP_iter(R, N):
HSumVariance[n] += sim["H"][n]**2 HSumVariance[n] += sim["H"][n]**2
T=sim['T'] T=sim['T']
V=sim['V'] V=sim['V']
for i in range(N): # ensure that T, V have the same length as Sum_T, Sum_V
for i in range(N - sim['i']):
T.append(0) T.append(0)
V.append(0) V.append(0)
Sum_T=[x+y for x,y in zip(Sum_T,T)] Sum_T=[x+y for x,y in zip(Sum_T,T)]
Sum_V=[x+y for x,y in zip(Sum_V,V)] Sum_V=[x+y for x,y in zip(Sum_V,V)]
Sum_T=[x/R for x in Sum_T] Sum_T=[x/R for x in Sum_T]
Sum_V=[round(x/R,2) for x in Sum_V] Sum_V=[round(x/R,2) for x in Sum_V]
#print(Sum_V) #print(Sum_V)
@ -102,7 +104,7 @@ def stats_NFBP_iter(R, N):
# TODO clarify line below # TODO clarify line below
print(" {} * {} iterations of T".format(R,N),'\n') print(" {} * {} iterations of T".format(R,N),'\n')
for n in range(N): for n in range(min(N, 10)):
Hn = HSum[n]/R # moyenne Hn = HSum[n]/R # moyenne
HVariance = sqrt(HSumVariance[n]/(R-1) - Hn**2) # Variance HVariance = sqrt(HSumVariance[n]/(R-1) - Hn**2) # Variance
print("Index of bin containing the {}th item (H_{}) : {} (variance {})".format(n, n, Hn, HVariance)) print("Index of bin containing the {}th item (H_{}) : {} (variance {})".format(n, n, Hn, HVariance))