# return False if koncnaResitev1 > nStars: return False, False ''' return koncnaResitev1, koncnaResitev0 / resitev1 # divide Lmax scale numHorSec = 48 # on how many parts should we devide L - how many different Lmax should we take noIterationsPerMaxL = 50000 # How many points N per each file/maxL - should be as big as possible ##trenutno 1000, da jih hitreje generira logPoints = np.linspace( 2, 10, numHorSec) # divide on numHorSec equal parts a scale from 2 to 10 - # for each Lmax create a file with points for maxL in logPoints: array = [] array1 = [] type_dist = "loguniform" for _ in range(0, noIterationsPerMaxL): point, nrsettlements = getPoint(maxL, type_dist) if type(point) != type(False): array.append(point) array1.append(nrsettlements) saveData(array, "/inf_" + type_dist + "_" + str(maxL)) #saveData(array1, "/inf_" + type_dist + "_settlements_" + str(maxL)) print("File: inf_" + type_dist + "_" + str(maxL) + ".txt created") print('done')
# x.append(astrophysicsProbability) #from our distribution calculated distribution # y.append(astrophysicsProbability2) #from our distribution return resitev drawnPoints = 0 numHorSec = 48 noIterationsPerMaxN = 20000 logPoints = np.linspace(0, 4, numHorSec) allPoints = noIterationsPerMaxN * numHorSec fixed_n = [1, 10, 100, 1000, 10000] for maxN in logPoints: # for maxN in fixed_n: array = [] for i in range(0, noIterationsPerMaxN): point = getPoint(maxN) if type(point) != type(False): array.append(point) saveData(array, "inf" + str(maxN)) print("File: inf" + str(maxN) + ".txt created. no of points:" + str(len(array))) drawnPoints = drawnPoints + len(array) pointFraction = (drawnPoints * 100) / allPoints print('done') print('Drawn points: ' + str(drawnPoints) + ' Which is: ' + str(pointFraction) + '%')
return resitev # devide Lmax scale numHorSec = 48 # on how many parts should we devide L - how many different Lmax should we take noIterationsPerMaxL = 50000 # How many points N per each file/maxL - should be as big as possible ##trenutno 1000, da jih hitreje generira logPoints = np.linspace(2, 10, numHorSec) # devide on numHorSec equal parts a scale from 2 to 10 - # x = [] # y=[] # for each Lmax create a file with points for maxL in logPoints: array = [] for _ in range(0, noIterationsPerMaxL): point = getPoint(maxL) if type(point) != type(False): array.append(point) saveData(array, "/inf" + str(maxL)) print("File: inf" + str(maxL) + ".txt created") ''' #for checking fLife x = np.array(x) y = np.array(y) plt.hist(x,np.linspace(-3,3,100),facecolor='blue',alpha=0.75) #plt.show() plt.hist(y,np.linspace(-3,3,100),facecolor='orange',alpha=0.5) plt.show() ''' print('done')
# return getPoint(maxL) return False ''' if resitev < -6: return False return 10**resitev # divide Lmax scale drawnPoints = 0 numHorSec = 20 # on how many parts should we divide L - how many different Lmax should we take noIterationsPerMaxL = 5000000 # How many points N per each file/maxL - should be as big as possible ##trenutno 1000, da jih hitreje generira # logPoints = np.linspace(2, 10, numHorSec) # divide on numHorSec equal parts a scale from 2 to 10 - # allPoints = noIterationsPerMaxL * numHorSec # for each Lmax create a file with points # for maxL in logPoints: array = [] for _ in range(0, noIterationsPerMaxL): point = getPoint() if type(point) != type(False): array.append(point) saveData(array, "/N_Sandberg_no_cut") # print("File: inf" + str(maxL) + ".txt created. Size: " + str(len(array))) # drawnPoints = drawnPoints + len(array) # pointFraction = (drawnPoints * 100) / allPoints print('done') # print('Drawn points: ' + str(drawnPoints) + ' Which is: ' + str(pointFraction) + '%')