Testing Notebook#
import numpy as np
import EnBiD_ananke
EnBiD_ananke.__version__
'0.2.0.post2'
Creating mock coordinates
np.random.seed(0)
points = np.random.randn(100000, 3)
Test enbid
name='demo'
rho = EnBiD_ananke.enbid(points, name=name)
Dimensions = 3
Normalization constant of Kernel type 3: 0.596831
Reading ASCII format file: to_enbid
Read 300000 records
Allocated 3.8147 MByte for particle storage.
Type = 1 Particles = 100000
Total particles = 100000
Allocated 22.1252 MByte for Binary-Tree and Entropy bins.
Scaling Co-ordinates as x[i]=x[i]/h[i] with h[i]->
1 1 1
Starting to Build Tree .......
Particle Type = 1 First node number = 1
Last node number = 199999
Total number of nodes = 199999
Treebuild time = 0.195773 s
Density Calculation. Smoothing .....
Evaluated = 1 % Time Left = 1.178880 s of 1.190800 s Par no = 15331 Density = 2.822111e+03
Evaluated = 2 % Time Left = 1.164522 s of 1.188300 s Par no = 77799 Density = 1.442662e+03
Evaluated = 3 % Time Left = 1.112352 s of 1.146767 s Par no = 20716 Density = 1.347107e+03
Evaluated = 4 % Time Left = 1.096189 s of 1.141875 s Par no = 36598 Density = 1.518448e+03
Evaluated = 5 % Time Left = 1.091159 s of 1.148600 s Par no = 54783 Density = 1.418503e+03
Evaluated = 6 % Time Left = 1.106572 s of 1.177217 s Par no = 18231 Density = 5.910374e+03
Evaluated = 7 % Time Left = 1.093124 s of 1.175414 s Par no = 77863 Density = 3.802778e+03
Evaluated = 8 % Time Left = 1.081034 s of 1.175050 s Par no = 29694 Density = 3.112394e+03
Evaluated = 9 % Time Left = 1.072393 s of 1.178467 s Par no = 49839 Density = 6.421483e+02
Evaluated = 10 % Time Left = 1.058073 s of 1.175650 s Par no = 63429 Density = 4.052199e+03
Evaluated = 11 % Time Left = 1.043643 s of 1.172645 s Par no = 34096 Density = 2.971922e+03
Evaluated = 12 % Time Left = 1.031143 s of 1.171767 s Par no = 73880 Density = 1.911730e+03
Evaluated = 13 % Time Left = 1.015834 s of 1.167638 s Par no = 59488 Density = 3.942927e+03
Evaluated = 14 % Time Left = 1.016858 s of 1.182407 s Par no = 17747 Density = 4.903259e+03
Evaluated = 15 % Time Left = 1.005680 s of 1.183167 s Par no = 50409 Density = 7.975875e+02
Evaluated = 16 % Time Left = 0.993015 s of 1.182175 s Par no = 66709 Density = 3.939722e+03
Evaluated = 17 % Time Left = 0.982928 s of 1.184265 s Par no = 58854 Density = 7.299768e+02
Evaluated = 18 % Time Left = 0.978344 s of 1.193117 s Par no = 6149 Density = 3.353193e+03
Evaluated = 19 % Time Left = 0.975416 s of 1.204232 s Par no = 13280 Density = 3.527213e+03
Evaluated = 20 % Time Left = 0.967216 s of 1.209035 s Par no = 77854 Density = 7.278303e+03
Evaluated = 21 % Time Left = 0.954330 s of 1.208029 s Par no = 84128 Density = 4.000489e+03
Evaluated = 22 % Time Left = 0.943653 s of 1.209827 s Par no = 12510 Density = 2.366131e+03
Evaluated = 23 % Time Left = 0.930737 s of 1.208765 s Par no = 38046 Density = 2.401039e+03
Evaluated = 24 % Time Left = 0.921387 s of 1.212367 s Par no = 76317 Density = 2.604788e+03
Evaluated = 25 % Time Left = 0.907695 s of 1.210276 s Par no = 55671 Density = 2.979146e+03
Evaluated = 26 % Time Left = 0.898664 s of 1.214427 s Par no = 24840 Density = 3.602203e+03
Evaluated = 27 % Time Left = 0.887568 s of 1.215863 s Par no = 34979 Density = 1.324684e+03
Evaluated = 28 % Time Left = 0.871556 s of 1.210511 s Par no = 17339 Density = 1.698872e+03
Evaluated = 29 % Time Left = 0.860217 s of 1.211590 s Par no = 59059 Density = 3.364317e+03
Evaluated = 30 % Time Left = 0.848971 s of 1.212833 s Par no = 48861 Density = 4.071747e+02
Evaluated = 31 % Time Left = 0.835533 s of 1.210935 s Par no = 34025 Density = 2.219551e+03
Evaluated = 32 % Time Left = 0.824722 s of 1.212844 s Par no = 92378 Density = 2.220241e+03
Evaluated = 33 % Time Left = 0.813035 s of 1.213503 s Par no = 14610 Density = 2.642163e+03
Evaluated = 34 % Time Left = 0.802889 s of 1.216518 s Par no = 79246 Density = 2.082532e+03
Evaluated = 35 % Time Left = 0.793517 s of 1.220814 s Par no = 2398 Density = 3.694222e+03
Evaluated = 36 % Time Left = 0.784734 s of 1.226167 s Par no = 24906 Density = 3.338249e+03
Evaluated = 37 % Time Left = 0.770851 s of 1.223592 s Par no = 64780 Density = 2.756157e+02
Evaluated = 38 % Time Left = 0.757783 s of 1.222250 s Par no = 56355 Density = 5.440238e+03
Evaluated = 39 % Time Left = 0.745475 s of 1.222110 s Par no = 75129 Density = 1.397896e+03
Evaluated = 40 % Time Left = 0.732711 s of 1.221205 s Par no = 37669 Density = 2.296944e+02
Evaluated = 41 % Time Left = 0.719517 s of 1.219541 s Par no = 44245 Density = 1.726854e+02
Evaluated = 42 % Time Left = 0.708527 s of 1.221619 s Par no = 51931 Density = 4.815193e+03
Evaluated = 43 % Time Left = 0.697771 s of 1.224181 s Par no = 43963 Density = 3.715361e+03
Evaluated = 44 % Time Left = 0.685939 s of 1.224914 s Par no = 56586 Density = 4.460856e+03
Evaluated = 45 % Time Left = 0.676918 s of 1.230782 s Par no = 3177 Density = 4.708334e+02
Evaluated = 46 % Time Left = 0.663603 s of 1.228917 s Par no = 16436 Density = 3.666953e+03
Evaluated = 47 % Time Left = 0.652705 s of 1.231543 s Par no = 27577 Density = 1.474241e+03
Evaluated = 48 % Time Left = 0.638816 s of 1.228517 s Par no = 76070 Density = 1.748841e+03
Evaluated = 49 % Time Left = 0.626664 s of 1.228778 s Par no = 5116 Density = 3.177198e+03
Evaluated = 50 % Time Left = 0.614867 s of 1.229758 s Par no = 36593 Density = 3.739479e+03
Evaluated = 51 % Time Left = 0.601677 s of 1.227937 s Par no = 10533 Density = 4.443296e+03
Evaluated = 52 % Time Left = 0.590598 s of 1.230438 s Par no = 64236 Density = 1.314060e+03
Evaluated = 53 % Time Left = 0.578380 s of 1.230623 s Par no = 31107 Density = 4.358351e+03
Evaluated = 54 % Time Left = 0.566096 s of 1.230670 s Par no = 37390 Density = 8.293836e+02
Evaluated = 55 % Time Left = 0.552226 s of 1.227196 s Par no = 44193 Density = 6.911039e+02
Evaluated = 56 % Time Left = 0.540732 s of 1.228964 s Par no = 73564 Density = 1.592207e+03
Evaluated = 57 % Time Left = 0.527950 s of 1.227819 s Par no = 84805 Density = 3.435408e+03
Evaluated = 58 % Time Left = 0.516406 s of 1.229567 s Par no = 70836 Density = 1.086561e+03
Evaluated = 59 % Time Left = 0.503204 s of 1.227358 s Par no = 21538 Density = 9.843017e+02
Evaluated = 60 % Time Left = 0.491015 s of 1.227568 s Par no = 26397 Density = 2.703689e+03
Evaluated = 61 % Time Left = 0.479890 s of 1.230520 s Par no = 59479 Density = 4.639280e+03
Evaluated = 62 % Time Left = 0.468054 s of 1.231753 s Par no = 73940 Density = 6.484902e+03
Evaluated = 63 % Time Left = 0.455467 s of 1.231025 s Par no = 81529 Density = 5.558739e+02
Evaluated = 64 % Time Left = 0.442311 s of 1.228675 s Par no = 83734 Density = 3.109286e+03
Evaluated = 65 % Time Left = 0.429437 s of 1.226997 s Par no = 36656 Density = 1.391414e+03
Evaluated = 66 % Time Left = 0.416495 s of 1.225023 s Par no = 57321 Density = 3.723538e+02
Evaluated = 67 % Time Left = 0.403467 s of 1.222664 s Par no = 35213 Density = 3.780156e+03
Evaluated = 68 % Time Left = 0.392212 s of 1.225701 s Par no = 65080 Density = 1.219758e+03
Evaluated = 69 % Time Left = 0.379816 s of 1.225254 s Par no = 16606 Density = 2.603137e+02
Evaluated = 70 % Time Left = 0.367315 s of 1.224424 s Par no = 87732 Density = 1.023583e+03
Evaluated = 71 % Time Left = 0.355805 s of 1.226956 s Par no = 5640 Density = 5.597383e+02
Evaluated = 72 % Time Left = 0.344073 s of 1.228875 s Par no = 51991 Density = 2.662520e+03
Evaluated = 73 % Time Left = 0.331514 s of 1.227874 s Par no = 26385 Density = 4.305937e+01
Evaluated = 74 % Time Left = 0.319136 s of 1.227493 s Par no = 37310 Density = 3.993851e+03
Evaluated = 75 % Time Left = 0.306509 s of 1.226087 s Par no = 97131 Density = 1.369731e+03
Evaluated = 76 % Time Left = 0.294285 s of 1.226238 s Par no = 3103 Density = 2.436843e+03
Evaluated = 77 % Time Left = 0.281684 s of 1.224768 s Par no = 25522 Density = 1.605624e+03
Evaluated = 78 % Time Left = 0.269137 s of 1.223405 s Par no = 910 Density = 5.151483e+01
Evaluated = 79 % Time Left = 0.256613 s of 1.222023 s Par no = 35125 Density = 2.473162e+01
Evaluated = 80 % Time Left = 0.244173 s of 1.220924 s Par no = 37920 Density = 9.703310e+02
Evaluated = 81 % Time Left = 0.231878 s of 1.220473 s Par no = 24227 Density = 3.733124e+02
Evaluated = 82 % Time Left = 0.219556 s of 1.219821 s Par no = 45909 Density = 2.877809e+03
Evaluated = 83 % Time Left = 0.207273 s of 1.219324 s Par no = 51688 Density = 4.196748e+03
Evaluated = 84 % Time Left = 0.195298 s of 1.220692 s Par no = 95852 Density = 5.509703e+02
Evaluated = 85 % Time Left = 0.183087 s of 1.220664 s Par no = 39450 Density = 1.983139e+02
Evaluated = 86 % Time Left = 0.171183 s of 1.222823 s Par no = 8808 Density = 1.392768e+02
Evaluated = 87 % Time Left = 0.158956 s of 1.222831 s Par no = 57392 Density = 3.091543e+03
Evaluated = 88 % Time Left = 0.146844 s of 1.223806 s Par no = 83393 Density = 3.849348e+03
Evaluated = 89 % Time Left = 0.134889 s of 1.226373 s Par no = 68997 Density = 4.455536e+03
Evaluated = 90 % Time Left = 0.122733 s of 1.227448 s Par no = 79817 Density = 4.377665e+01
Evaluated = 91 % Time Left = 0.110395 s of 1.226743 s Par no = 83472 Density = 2.215261e+03
Evaluated = 92 % Time Left = 0.098346 s of 1.229477 s Par no = 24893 Density = 4.393429e+03
Evaluated = 93 % Time Left = 0.086143 s of 1.230797 s Par no = 64593 Density = 1.252936e+03
Evaluated = 94 % Time Left = 0.073938 s of 1.232510 s Par no = 11536 Density = 2.153480e+03
Evaluated = 95 % Time Left = 0.061726 s of 1.234769 s Par no = 37158 Density = 3.735359e+01
Evaluated = 96 % Time Left = 0.049340 s of 1.233819 s Par no = 43390 Density = 3.406644e+01
Evaluated = 97 % Time Left = 0.036988 s of 1.233329 s Par no = 10529 Density = 3.171351e+03
Evaluated = 98 % Time Left = 0.024696 s of 1.235413 s Par no = 7315 Density = 1.658222e+03
Evaluated = 99 % Time Left = 0.012360 s of 1.237261 s Par no = 32071 Density = 5.807275e+02
Evaluated = 100 % Time Left = -0.000012 s of 1.239249 s Par no = 35331 Density = 1.974778e+01
Total Smoothing Time = 1.23928 s
1 3
writing snapshot file ....
to_enbid_d3n64.est
done with snapshot.
Total Time = 1.92481 s
/home/docs/checkouts/readthedocs.org/user_builds/py-enbid-ananke/envs/v0.2.0.post2/lib/python3.9/site-packages/EnBiD_ananke/__init__.py:237: FutureWarning: The 'delim_whitespace' keyword in pd.read_table is deprecated and will be removed in a future version. Use ``sep='\s+'`` instead
usedvals = pd.read_table(path / CONSTANTS.usedvalues, header=None, delim_whitespace=True,
Show result
print(rho)
[ 654.357 71.6832 4765.37 ... 3993.97 4092.39 415.127 ]