{ "cells": [ { "cell_type": "markdown", "id": "2f00fe78", "metadata": {}, "source": [ "# Testing Notebook" ] }, { "cell_type": "code", "execution_count": 1, "id": "f130bf6c", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'0.1.3b2.dev2'" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import numpy as np\n", "import EnBiD_ananke\n", "EnBiD_ananke.__version__" ] }, { "cell_type": "markdown", "id": "a8026efc", "metadata": {}, "source": [ "Creating mock coordinates" ] }, { "cell_type": "code", "execution_count": 2, "id": "3a74890b", "metadata": {}, "outputs": [], "source": [ "np.random.seed(0)\n", "points = np.random.randn(100000, 3)" ] }, { "cell_type": "markdown", "id": "d0809afe", "metadata": {}, "source": [ "Test enbid" ] }, { "cell_type": "code", "execution_count": 3, "id": "fa0e3b78", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Dimensions = 3\n", "Normalization constant of Kernel type 3: 0.596831\n", "Reading ASCII format file: to_enbid \n", "Read 300000 records\n", "Allocated 3.8147 MByte for particle storage.\n", "Type = 1 Particles = 100000\n", "Total particles = 100000\n", "Allocated 22.1252 MByte for Binary-Tree and Entropy bins.\n", "\n", "Scaling Co-ordinates as x[i]=x[i]/h[i] with h[i]->\n", "1 1 1 \n", "\n", "Starting to Build Tree .......\n", "Particle Type = 1 First node number = 1\n", " Last node number = 199999\n", "Total number of nodes = 199999\n", "Treebuild time = 0.076733 s \n", "\n", "Density Calculation. Smoothing .....\n", "Evaluated = 1 % Time Left = 0.572214 s of 0.578000 s Par no = 15331 Density = 2.822111e+03 \n", "Evaluated = 2 % Time Left = 0.579272 s of 0.591100 s Par no = 77799 Density = 1.442662e+03 \n", "Evaluated = 3 % Time Left = 0.563241 s of 0.580667 s Par no = 20716 Density = 1.347107e+03 \n", "Evaluated = 4 % Time Left = 0.563514 s of 0.587000 s Par no = 36598 Density = 1.518448e+03 \n", "Evaluated = 5 % Time Left = 0.558556 s of 0.587960 s Par no = 54783 Density = 1.418503e+03 \n", "Evaluated = 6 % Time Left = 0.572360 s of 0.608900 s Par no = 18231 Density = 5.910374e+03 \n", "Evaluated = 7 % Time Left = 0.568862 s of 0.611686 s Par no = 77863 Density = 3.802778e+03 \n", "Evaluated = 8 % Time Left = 0.566173 s of 0.615413 s Par no = 29694 Density = 3.112394e+03 \n", "Evaluated = 9 % Time Left = 0.561990 s of 0.617578 s Par no = 49839 Density = 6.421483e+02 \n", "Evaluated = 10 % Time Left = 0.553242 s of 0.614720 s Par no = 63429 Density = 4.052199e+03 \n", "Evaluated = 11 % Time Left = 0.545200 s of 0.612591 s Par no = 34096 Density = 2.971922e+03 \n", "Evaluated = 12 % Time Left = 0.538136 s of 0.611525 s Par no = 73880 Density = 1.911730e+03 \n", "Evaluated = 13 % Time Left = 0.528124 s of 0.607046 s Par no = 59488 Density = 3.942927e+03 \n", "Evaluated = 14 % Time Left = 0.527598 s of 0.613493 s Par no = 17747 Density = 4.903259e+03 \n", "Evaluated = 15 % Time Left = 0.521208 s of 0.613193 s Par no = 50409 Density = 7.975875e+02 \n", "Evaluated = 16 % Time Left = 0.516814 s of 0.615263 s Par no = 66709 Density = 3.939722e+03 \n", "Evaluated = 17 % Time Left = 0.511088 s of 0.615776 s Par no = 58854 Density = 7.299768e+02 \n", "Evaluated = 18 % Time Left = 0.508380 s of 0.619983 s Par no = 6149 Density = 3.353193e+03 \n", "Evaluated = 19 % Time Left = 0.505165 s of 0.623668 s Par no = 13280 Density = 3.527213e+03 \n", "Evaluated = 20 % Time Left = 0.499162 s of 0.623960 s Par no = 77854 Density = 7.278303e+03 \n", "Evaluated = 21 % Time Left = 0.494805 s of 0.626343 s Par no = 84128 Density = 4.000489e+03 \n", "Evaluated = 22 % Time Left = 0.488238 s of 0.625955 s Par no = 12510 Density = 2.366131e+03 \n", "Evaluated = 23 % Time Left = 0.483473 s of 0.627896 s Par no = 38046 Density = 2.401039e+03 \n", "Evaluated = 24 % Time Left = 0.479183 s of 0.630513 s Par no = 76317 Density = 2.604788e+03 \n", "Evaluated = 25 % Time Left = 0.471687 s of 0.628924 s Par no = 55671 Density = 2.979146e+03 \n", "Evaluated = 26 % Time Left = 0.467093 s of 0.631215 s Par no = 24840 Density = 3.602203e+03 \n", "Evaluated = 27 % Time Left = 0.460910 s of 0.631393 s Par no = 34979 Density = 1.324684e+03 \n", "Evaluated = 28 % Time Left = 0.452611 s of 0.628636 s Par no = 17339 Density = 1.698872e+03 \n", "Evaluated = 29 % Time Left = 0.446447 s of 0.628807 s Par no = 59059 Density = 3.364317e+03 \n", "Evaluated = 30 % Time Left = 0.441915 s of 0.631317 s Par no = 48861 Density = 4.071747e+02 \n", "Evaluated = 31 % Time Left = 0.434970 s of 0.630400 s Par no = 34025 Density = 2.219551e+03 \n", "Evaluated = 32 % Time Left = 0.429881 s of 0.632188 s Par no = 92378 Density = 2.220241e+03 \n", "Evaluated = 33 % Time Left = 0.425383 s of 0.634909 s Par no = 14610 Density = 2.642163e+03 \n", "Evaluated = 34 % Time Left = 0.419759 s of 0.636009 s Par no = 79246 Density = 2.082532e+03 \n", "Evaluated = 35 % Time Left = 0.415245 s of 0.638849 s Par no = 2398 Density = 3.694222e+03 \n", "Evaluated = 36 % Time Left = 0.411053 s of 0.642281 s Par no = 24906 Density = 3.338249e+03 \n", "Evaluated = 37 % Time Left = 0.403818 s of 0.640992 s Par no = 64780 Density = 2.756157e+02 \n", "Evaluated = 38 % Time Left = 0.397368 s of 0.640926 s Par no = 56355 Density = 5.440238e+03 \n", "Evaluated = 39 % Time Left = 0.391881 s of 0.642438 s Par no = 75129 Density = 1.397896e+03 \n", "Evaluated = 40 % Time Left = 0.385560 s of 0.642610 s Par no = 37669 Density = 2.296944e+02 \n", "Evaluated = 41 % Time Left = 0.379621 s of 0.643437 s Par no = 44245 Density = 1.726854e+02 \n", "Evaluated = 42 % Time Left = 0.374100 s of 0.645012 s Par no = 51931 Density = 4.815193e+03 \n", "Evaluated = 43 % Time Left = 0.369008 s of 0.647393 s Par no = 43963 Density = 3.715361e+03 \n", "Evaluated = 44 % Time Left = 0.363408 s of 0.648955 s Par no = 56586 Density = 4.460856e+03 \n", "Evaluated = 45 % Time Left = 0.358502 s of 0.651833 s Par no = 3177 Density = 4.708334e+02 \n", "Evaluated = 46 % Time Left = 0.353145 s of 0.653985 s Par no = 16436 Density = 3.666953e+03 \n", "Evaluated = 47 % Time Left = 0.347449 s of 0.655577 s Par no = 27577 Density = 1.474241e+03 \n", "Evaluated = 48 % Time Left = 0.340208 s of 0.654258 s Par no = 76070 Density = 1.748841e+03 \n", "Evaluated = 49 % Time Left = 0.334020 s of 0.654953 s Par no = 5116 Density = 3.177198e+03 \n", "Evaluated = 50 % Time Left = 0.328331 s of 0.656676 s Par no = 36593 Density = 3.739479e+03 \n", "Evaluated = 51 % Time Left = 0.321242 s of 0.655610 s Par no = 10533 Density = 4.443296e+03 \n", "Evaluated = 52 % Time Left = 0.314676 s of 0.655588 s Par no = 64236 Density = 1.314060e+03 \n", "Evaluated = 53 % Time Left = 0.308304 s of 0.655979 s Par no = 31107 Density = 4.358351e+03 \n", "Evaluated = 54 % Time Left = 0.301848 s of 0.656206 s Par no = 37390 Density = 8.293836e+02 \n", "Evaluated = 55 % Time Left = 0.294812 s of 0.655153 s Par no = 44193 Density = 6.911039e+02 \n", "Evaluated = 56 % Time Left = 0.288907 s of 0.656621 s Par no = 73564 Density = 1.592207e+03 \n", "Evaluated = 57 % Time Left = 0.282201 s of 0.656296 s Par no = 84805 Density = 3.435408e+03 \n", "Evaluated = 58 % Time Left = 0.276374 s of 0.658048 s Par no = 70836 Density = 1.086561e+03 \n", "Evaluated = 59 % Time Left = 0.269181 s of 0.656556 s Par no = 21538 Density = 9.843017e+02 \n", "Evaluated = 60 % Time Left = 0.262623 s of 0.656575 s Par no = 26397 Density = 2.703689e+03 \n", "Evaluated = 61 % Time Left = 0.256569 s of 0.657885 s Par no = 59479 Density = 4.639280e+03 \n", "Evaluated = 62 % Time Left = 0.250463 s of 0.659131 s Par no = 73940 Density = 6.484902e+03 \n", "Evaluated = 63 % Time Left = 0.243657 s of 0.658551 s Par no = 81529 Density = 5.558739e+02 \n", "Evaluated = 64 % Time Left = 0.236579 s of 0.657183 s Par no = 83734 Density = 3.109286e+03 \n", "Evaluated = 65 % Time Left = 0.229474 s of 0.655658 s Par no = 36656 Density = 1.391414e+03 \n", "Evaluated = 66 % Time Left = 0.222403 s of 0.654147 s Par no = 57321 Density = 3.723538e+02 \n", "Evaluated = 67 % Time Left = 0.215533 s of 0.653149 s Par no = 35213 Density = 3.780156e+03 \n", "Evaluated = 68 % Time Left = 0.209652 s of 0.655184 s Par no = 65080 Density = 1.219758e+03 \n", "Evaluated = 69 % Time Left = 0.202872 s of 0.654446 s Par no = 16606 Density = 2.603137e+02 \n", "Evaluated = 70 % Time Left = 0.196184 s of 0.653970 s Par no = 87732 Density = 1.023583e+03 \n", "Evaluated = 71 % Time Left = 0.190398 s of 0.656566 s Par no = 5640 Density = 5.597383e+02 \n", "Evaluated = 72 % Time Left = 0.184068 s of 0.657408 s Par no = 51991 Density = 2.662520e+03 \n", "Evaluated = 73 % Time Left = 0.177240 s of 0.656470 s Par no = 26385 Density = 4.305937e+01 \n", "Evaluated = 74 % Time Left = 0.170738 s of 0.656709 s Par no = 37310 Density = 3.993851e+03 \n", "Evaluated = 75 % Time Left = 0.163909 s of 0.655661 s Par no = 97131 Density = 1.369731e+03 \n", "Evaluated = 76 % Time Left = 0.157313 s of 0.655499 s Par no = 3103 Density = 2.436843e+03 \n", "Evaluated = 77 % Time Left = 0.150563 s of 0.654652 s Par no = 25522 Density = 1.605624e+03 \n", "Evaluated = 78 % Time Left = 0.143859 s of 0.653935 s Par no = 910 Density = 5.151483e+01 \n", "Evaluated = 79 % Time Left = 0.137143 s of 0.653092 s Par no = 35125 Density = 2.473162e+01 \n", "Evaluated = 80 % Time Left = 0.130479 s of 0.652427 s Par no = 37920 Density = 9.703310e+02 \n", "Evaluated = 81 % Time Left = 0.123994 s of 0.652635 s Par no = 24227 Density = 3.733124e+02 \n", "Evaluated = 82 % Time Left = 0.117328 s of 0.651857 s Par no = 45909 Density = 2.877809e+03 \n", "Evaluated = 83 % Time Left = 0.110767 s of 0.651607 s Par no = 51688 Density = 4.196748e+03 \n", "Evaluated = 84 % Time Left = 0.104438 s of 0.652777 s Par no = 95852 Density = 5.509703e+02 \n", "Evaluated = 85 % Time Left = 0.097936 s of 0.652951 s Par no = 39450 Density = 1.983139e+02 \n", "Evaluated = 86 % Time Left = 0.091488 s of 0.653534 s Par no = 8808 Density = 1.392768e+02 \n", "Evaluated = 87 % Time Left = 0.084904 s of 0.653160 s Par no = 57392 Density = 3.091543e+03 \n", "Evaluated = 88 % Time Left = 0.078389 s of 0.653293 s Par no = 83393 Density = 3.849348e+03 \n", "Evaluated = 89 % Time Left = 0.071971 s of 0.654338 s Par no = 68997 Density = 4.455536e+03 \n", "Evaluated = 90 % Time Left = 0.065447 s of 0.654531 s Par no = 79817 Density = 4.377665e+01 \n", "Evaluated = 91 % Time Left = 0.058870 s of 0.654181 s Par no = 83472 Density = 2.215261e+03 \n", "Evaluated = 92 % Time Left = 0.052421 s of 0.655346 s Par no = 24893 Density = 4.393429e+03 \n", "Evaluated = 93 % Time Left = 0.045906 s of 0.655887 s Par no = 64593 Density = 1.252936e+03 \n", "Evaluated = 94 % Time Left = 0.039386 s of 0.656543 s Par no = 11536 Density = 2.153480e+03 \n", "Evaluated = 95 % Time Left = 0.032888 s of 0.657894 s Par no = 37158 Density = 3.735359e+01 \n", "Evaluated = 96 % Time Left = 0.026302 s of 0.657719 s Par no = 43390 Density = 3.406644e+01 \n", "Evaluated = 97 % Time Left = 0.019728 s of 0.657828 s Par no = 10529 Density = 3.171351e+03 \n", "Evaluated = 98 % Time Left = 0.013176 s of 0.659142 s Par no = 7315 Density = 1.658222e+03 \n", "Evaluated = 99 % Time Left = 0.006606 s of 0.661280 s Par no = 32071 Density = 5.807275e+02 \n", "Evaluated = 100 % Time Left = -0.000007 s of 0.662902 s Par no = 35331 Density = 1.974778e+01 \n", "\n", "Total Smoothing Time = 0.662931 s\n", "1 3\n", "\n", "writing snapshot file ....\n", "to_enbid_d3n64.est\n", "done with snapshot.\n", "\n", "Total Time = 0.931126 s \n", "\n" ] } ], "source": [ "name='demo'\n", "rho = EnBiD_ananke.enbid(points, name=name)" ] }, { "cell_type": "markdown", "id": "5b85b27e", "metadata": {}, "source": [ "Show result" ] }, { "cell_type": "code", "execution_count": 4, "id": "cf4db07a", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 654.357 71.6832 4765.37 ... 3993.97 4092.39 415.127 ]\n" ] } ], "source": [ "print(rho)" ] } ], "metadata": { "kernelspec": { "display_name": "base", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.12" }, "vscode": { "interpreter": { "hash": "93aa337d3be5595eae1049f286a1f242db5cd4d40fd613efbe29fb4c796799bf" } } }, "nbformat": 4, "nbformat_minor": 5 }