Testing Notebook#


Author: Adrien CR Thob

© 2022 Adrien CR Thob

This notebook is part of the py-EnBiD-ananke project, which is licensed under the GNU General Public License v2.0 or later (GPL-2.0-or-later). The full copyright notice, including terms governing use, modification, and redistribution, is contained in the files LICENSE and COPYRIGHT, which can be found at the root of the source code distribution tree

Updated On: 2026-04-16


import numpy as np
import enbid_ananke
enbid_ananke.__version__
'0.3.4.dev0+15.gc800f25'

Creating mock coordinates

np.random.seed(0)
points = np.random.randn(100000, 3)
velocities = np.random.randn(points.shape[0], 3)
mass = np.exp(np.random.randn(points.shape[0]))

Test enbid

name='demo'
rho = enbid_ananke.enbid(points, velocities=velocities, mass=mass, name=name, caching=False, verbose=True)
Dimensions = 6
Normalization constant of Kernel type 3: 0.774037
Reading GADGET format file: to_enbid 
Allocated 3.8147 MByte for particle storage.
....done with reading.
Type = 1 Particles = 100000
Total particles = 100000
WARNING! particles with multiple mass results might not be accurate
Allocated 33.1879 MByte for Binary-Tree and Entropy bins.

Scaling Co-ordinates as x[i]=x[i]/h[i] with h[i]->
1 1 1 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.2117 s 

Density Calculation. Smoothing .....
Evaluated =   1 % Time Left = 5.268331 s of 5.321600 s Par no = 80305 Density = 7.141296e-01 
Evaluated =   2 % Time Left = 5.344963 s of 5.454100 s Par no = 93926 Density = 2.079772e+01 
Evaluated =   3 % Time Left = 5.632667 s of 5.806933 s Par no = 37988 Density = 2.752346e+02 
Evaluated =   4 % Time Left = 5.777460 s of 6.018250 s Par no = 26785 Density = 1.703962e+02 
Evaluated =   5 % Time Left = 5.870710 s of 6.179760 s Par no = 3855 Density = 2.293737e+02 
Evaluated =   6 % Time Left = 5.876676 s of 6.251850 s Par no = 46956 Density = 2.849157e+02 
Evaluated =   7 % Time Left = 5.911109 s of 6.356100 s Par no = 96759 Density = 1.150033e+02 
Evaluated =   8 % Time Left = 5.911499 s of 6.425612 s Par no = 88921 Density = 1.166282e+02 
Evaluated =   9 % Time Left = 5.870558 s of 6.451233 s Par no = 72202 Density = 1.477496e+02 
Evaluated =  10 % Time Left = 5.853112 s of 6.503530 s Par no = 42648 Density = 9.552992e+01 
Evaluated =  11 % Time Left = 5.815397 s of 6.534227 s Par no = 6388 Density = 1.553956e+02 
Evaluated =  12 % Time Left = 5.714640 s of 6.493983 s Par no = 32000 Density = 2.751964e+02 
Evaluated =  13 % Time Left = 5.696233 s of 6.547469 s Par no = 76572 Density = 2.044967e+02 
Evaluated =  14 % Time Left = 5.619463 s of 6.534336 s Par no = 58036 Density = 8.976174e+00 
Evaluated =  15 % Time Left = 5.552191 s of 6.532067 s Par no = 50258 Density = 1.440317e+02 
Evaluated =  16 % Time Left = 5.502890 s of 6.551137 s Par no = 68992 Density = 1.080991e+02 
Evaluated =  17 % Time Left = 5.421968 s of 6.532571 s Par no = 74604 Density = 4.118754e+01 
Evaluated =  18 % Time Left = 5.367440 s of 6.545739 s Par no = 24512 Density = 8.734037e+01 
Evaluated =  19 % Time Left = 5.322146 s of 6.570632 s Par no = 2467 Density = 1.471868e+01 
Evaluated =  20 % Time Left = 5.260174 s of 6.575300 s Par no = 17203 Density = 1.095730e+02 
Evaluated =  21 % Time Left = 5.207279 s of 6.591576 s Par no = 62546 Density = 1.659647e+02 
Evaluated =  22 % Time Left = 5.119113 s of 6.563050 s Par no = 14503 Density = 9.163900e+01 
Evaluated =  23 % Time Left = 5.050147 s of 6.558717 s Par no = 54143 Density = 1.006519e+02 
Evaluated =  24 % Time Left = 4.998271 s of 6.576758 s Par no = 55616 Density = 9.019440e+01 
Evaluated =  25 % Time Left = 4.941570 s of 6.588848 s Par no = 24138 Density = 9.377963e+01 
Evaluated =  26 % Time Left = 4.880325 s of 6.595123 s Par no = 69109 Density = 1.089671e+02 
Evaluated =  27 % Time Left = 4.809766 s of 6.588811 s Par no = 23852 Density = 1.178417e+02 
Evaluated =  28 % Time Left = 4.723734 s of 6.560832 s Par no = 78410 Density = 1.415703e+02 
Evaluated =  29 % Time Left = 4.651448 s of 6.551428 s Par no = 99211 Density = 7.028931e+01 
Evaluated =  30 % Time Left = 4.599324 s of 6.570557 s Par no = 27607 Density = 1.573432e+02 
Evaluated =  31 % Time Left = 4.546046 s of 6.588568 s Par no = 81375 Density = 2.733477e+01 
Evaluated =  32 % Time Left = 4.498049 s of 6.614875 s Par no = 79457 Density = 1.742180e+01 
Evaluated =  33 % Time Left = 4.428421 s of 6.609682 s Par no = 81012 Density = 5.446134e+01 
Evaluated =  34 % Time Left = 4.346337 s of 6.585459 s Par no = 79989 Density = 5.379813e+01 
Evaluated =  35 % Time Left = 4.275592 s of 6.577934 s Par no = 84097 Density = 1.963045e+02 
Evaluated =  36 % Time Left = 4.226604 s of 6.604172 s Par no = 35932 Density = 5.158055e+01 
Evaluated =  37 % Time Left = 4.177450 s of 6.630978 s Par no = 27555 Density = 6.102217e+01 
Evaluated =  38 % Time Left = 4.130638 s of 6.662426 s Par no = 26268 Density = 6.534772e+01 
Evaluated =  39 % Time Left = 4.073645 s of 6.678215 s Par no = 42879 Density = 1.262151e+02 
Evaluated =  40 % Time Left = 4.013000 s of 6.688445 s Par no = 55482 Density = 3.073591e+01 
Evaluated =  41 % Time Left = 3.957395 s of 6.707563 s Par no = 72476 Density = 7.540024e+01 
Evaluated =  42 % Time Left = 3.894746 s of 6.715195 s Par no = 76397 Density = 5.395208e+01 
Evaluated =  43 % Time Left = 3.824479 s of 6.709730 s Par no = 12939 Density = 3.746408e+01 
Evaluated =  44 % Time Left = 3.754819 s of 6.705155 s Par no = 94789 Density = 2.172225e+01 
Evaluated =  45 % Time Left = 3.682976 s of 6.696442 s Par no = 35983 Density = 2.090884e+01 
Evaluated =  46 % Time Left = 3.594259 s of 6.656159 s Par no = 39379 Density = 1.390368e+01 
Evaluated =  47 % Time Left = 3.525702 s of 6.652394 s Par no = 57737 Density = 6.803721e+01 
Evaluated =  48 % Time Left = 3.458105 s of 6.650329 s Par no = 46337 Density = 5.049719e+01 
Evaluated =  49 % Time Left = 3.393766 s of 6.654573 s Par no = 91051 Density = 7.031203e+00 
Evaluated =  50 % Time Left = 3.327086 s of 6.654306 s Par no = 84718 Density = 2.533191e+01 
Evaluated =  51 % Time Left = 3.257530 s of 6.648157 s Par no = 66270 Density = 6.864194e+01 
Evaluated =  52 % Time Left = 3.198578 s of 6.663842 s Par no = 51768 Density = 8.701906e+01 
Evaluated =  53 % Time Left = 3.126399 s of 6.652055 s Par no = 50751 Density = 2.802341e+00 
Evaluated =  54 % Time Left = 3.055248 s of 6.641989 s Par no = 83825 Density = 7.740121e+01 
Evaluated =  55 % Time Left = 2.987339 s of 6.638678 s Par no = 83533 Density = 1.241632e+02 
Evaluated =  56 % Time Left = 2.908619 s of 6.610648 s Par no = 40668 Density = 1.836257e+01 
Evaluated =  57 % Time Left = 2.832755 s of 6.587954 s Par no = 58938 Density = 1.347231e+02 
Evaluated =  58 % Time Left = 2.759591 s of 6.570612 s Par no = 79239 Density = 4.573686e+01 
Evaluated =  59 % Time Left = 2.697652 s of 6.579800 s Par no = 74985 Density = 2.979457e+01 
Evaluated =  60 % Time Left = 2.631718 s of 6.579460 s Par no = 80940 Density = 4.740730e+00 
Evaluated =  61 % Time Left = 2.559853 s of 6.563893 s Par no = 37556 Density = 3.059815e+01 
Evaluated =  62 % Time Left = 2.491364 s of 6.556394 s Par no = 29279 Density = 1.618676e+01 
Evaluated =  63 % Time Left = 2.427956 s of 6.562221 s Par no = 80126 Density = 1.950052e+01 
Evaluated =  64 % Time Left = 2.362360 s of 6.562292 s Par no = 15070 Density = 5.965748e+01 
Evaluated =  65 % Time Left = 2.292472 s of 6.550108 s Par no = 27233 Density = 6.144010e+01 
Evaluated =  66 % Time Left = 2.228148 s of 6.553570 s Par no = 91836 Density = 2.171463e+00 
Evaluated =  67 % Time Left = 2.162703 s of 6.553845 s Par no = 91573 Density = 2.905033e+00 
Evaluated =  68 % Time Left = 2.096851 s of 6.552865 s Par no = 65169 Density = 5.529618e+00 
Evaluated =  69 % Time Left = 2.030796 s of 6.551167 s Par no = 62593 Density = 3.232186e+00 
Evaluated =  70 % Time Left = 1.966137 s of 6.554010 s Par no = 71165 Density = 3.316961e+01 
Evaluated =  71 % Time Left = 1.901687 s of 6.557766 s Par no = 57553 Density = 9.221316e+00 
Evaluated =  72 % Time Left = 1.837191 s of 6.561629 s Par no = 66318 Density = 5.028376e+00 
Evaluated =  73 % Time Left = 1.772463 s of 6.564919 s Par no = 35690 Density = 5.694636e+01 
Evaluated =  74 % Time Left = 1.708471 s of 6.571293 s Par no = 14894 Density = 4.638654e+01 
Evaluated =  75 % Time Left = 1.642954 s of 6.572077 s Par no = 76643 Density = 3.303287e+00 
Evaluated =  76 % Time Left = 1.579079 s of 6.579768 s Par no = 15021 Density = 1.765642e+01 
Evaluated =  77 % Time Left = 1.515484 s of 6.589345 s Par no = 79810 Density = 3.118361e+01 
Evaluated =  78 % Time Left = 1.450287 s of 6.592515 s Par no = 40455 Density = 4.149385e+01 
Evaluated =  79 % Time Left = 1.385809 s of 6.599403 s Par no = 72504 Density = 1.423325e+01 
Evaluated =  80 % Time Left = 1.321958 s of 6.610121 s Par no = 85345 Density = 1.895958e+01 
Evaluated =  81 % Time Left = 1.258561 s of 6.624354 s Par no = 75587 Density = 1.096318e+02 
Evaluated =  82 % Time Left = 1.194204 s of 6.634835 s Par no = 79939 Density = 1.737358e+01 
Evaluated =  83 % Time Left = 1.131225 s of 6.654654 s Par no = 75612 Density = 4.784624e+01 
Evaluated =  84 % Time Left = 1.066115 s of 6.663633 s Par no = 64827 Density = 3.080563e+00 
Evaluated =  85 % Time Left = 1.001019 s of 6.673902 s Par no = 75103 Density = 4.445605e+00 
Evaluated =  86 % Time Left = 0.934056 s of 6.672307 s Par no = 74984 Density = 1.383137e+01 
Evaluated =  87 % Time Left = 0.867899 s of 6.676657 s Par no = 90670 Density = 5.545300e+00 
Evaluated =  88 % Time Left = 0.802951 s of 6.691813 s Par no = 72460 Density = 7.752909e+01 
Evaluated =  89 % Time Left = 0.738412 s of 6.713445 s Par no = 54642 Density = 1.674362e+00 
Evaluated =  90 % Time Left = 0.672884 s of 6.729513 s Par no = 74375 Density = 3.624767e+00 
Evaluated =  91 % Time Left = 0.606207 s of 6.736386 s Par no = 45222 Density = 2.378857e+02 
Evaluated =  92 % Time Left = 0.539343 s of 6.742627 s Par no = 93733 Density = 3.569487e+00 
Evaluated =  93 % Time Left = 0.472969 s of 6.757667 s Par no = 72029 Density = 3.585799e+00 
Evaluated =  94 % Time Left = 0.406127 s of 6.769910 s Par no = 41294 Density = 3.801584e+00 
Evaluated =  95 % Time Left = 0.339291 s of 6.787175 s Par no = 68481 Density = 7.842824e+00 
Evaluated =  96 % Time Left = 0.271685 s of 6.793830 s Par no = 21797 Density = 1.532215e+00 
Evaluated =  97 % Time Left = 0.203648 s of 6.790522 s Par no = 91023 Density = 5.706758e+00 
Evaluated =  98 % Time Left = 0.135858 s of 6.796315 s Par no = 98577 Density = 4.745541e+00 
Evaluated =  99 % Time Left = 0.067928 s of 6.799553 s Par no = 33226 Density = 6.323213e+00 
Evaluated = 100 % Time Left = -0.000068 s of 6.807822 s Par no = 89304 Density = 1.046962e-02 

Total Smoothing Time = 6.80785 s
1 3

writing snapshot file ....
to_enbid_d6n64.est
done with snapshot.

Total Time = 7.03237 s 

Show result

print(rho)
[  8.581022    6.9317503 141.92      ...  56.267246  266.7389
   8.202949 ]