|
11 | 11 | }, |
12 | 12 | { |
13 | 13 | "cell_type": "code", |
14 | | - "execution_count": null, |
| 14 | + "execution_count": 1, |
15 | 15 | "metadata": {}, |
16 | 16 | "outputs": [], |
17 | 17 | "source": [ |
|
43 | 43 | }, |
44 | 44 | { |
45 | 45 | "cell_type": "code", |
46 | | - "execution_count": null, |
| 46 | + "execution_count": 2, |
47 | 47 | "metadata": {}, |
48 | | - "outputs": [], |
| 48 | + "outputs": [ |
| 49 | + { |
| 50 | + "name": "stdout", |
| 51 | + "output_type": "stream", |
| 52 | + "text": [ |
| 53 | + "A coefficients: [1.04581094e-03-2.45891920e-02j 5.61754632e-06-3.39444972e-04j\n", |
| 54 | + " 3.54498410e-08-2.25608771e-06j 1.33314410e-10-8.61596363e-09j\n", |
| 55 | + " 3.25668489e-13-2.12414217e-11j]\n", |
| 56 | + "B coefficients: [2.22095481e-05-8.65285036e-04j 1.48987829e-07-6.11710129e-06j\n", |
| 57 | + " 5.86860069e-10-2.42490110e-08j 1.47867923e-12-6.12742871e-11j\n", |
| 58 | + " 2.58354155e-15-1.07228248e-13j]\n" |
| 59 | + ] |
| 60 | + } |
| 61 | + ], |
49 | 62 | "source": [ |
50 | 63 | "from deeptrack.backend import mie\n", |
51 | 64 | "\n", |
|
70 | 83 | }, |
71 | 84 | { |
72 | 85 | "cell_type": "code", |
73 | | - "execution_count": null, |
| 86 | + "execution_count": 3, |
74 | 87 | "metadata": {}, |
75 | | - "outputs": [], |
| 88 | + "outputs": [ |
| 89 | + { |
| 90 | + "name": "stdout", |
| 91 | + "output_type": "stream", |
| 92 | + "text": [ |
| 93 | + "A coefficients: [4.51815568e-03-4.14667025e-02j 9.00848393e-05-8.90449307e-04j\n", |
| 94 | + " 1.30191942e-06-9.70936049e-06j 1.04767670e-08-6.35738175e-08j\n", |
| 95 | + " 5.30930040e-11-2.77931029e-10j]\n", |
| 96 | + "B coefficients: [2.08384699e-04-1.97056962e-03j 2.98780924e-06-2.15999592e-05j\n", |
| 97 | + " 2.41258487e-08-1.40444143e-07j 1.22573034e-10-6.07956361e-10j\n", |
| 98 | + " 4.27642496e-13-1.88504769e-12j]\n" |
| 99 | + ] |
| 100 | + } |
| 101 | + ], |
76 | 102 | "source": [ |
77 | 103 | "particle_radii = [0.5, 0.6, 0.7]\n", |
78 | | - "relative_refract_index = 1.5 + 0.01j\n", |
| 104 | + "relative_refract_index = [1.5 + 0.01j, 1.3 + 0.02j, 1.1 + 0.03j]\n", |
79 | 105 | "max_order = 5\n", |
80 | 106 | "\n", |
81 | | - "A, B = mie.stratified_coefficients(relative_refract_index, particle_radii, max_order)\n", |
| 107 | + "A, B = mie.stratified_coefficients(\n", |
| 108 | + " relative_refract_index, particle_radii, max_order,\n", |
| 109 | + ")\n", |
82 | 110 | "\n", |
83 | 111 | "print(\"A coefficients:\", A)\n", |
84 | 112 | "print(\"B coefficients:\", B)" |
|
93 | 121 | }, |
94 | 122 | { |
95 | 123 | "cell_type": "code", |
96 | | - "execution_count": null, |
| 124 | + "execution_count": 4, |
97 | 125 | "metadata": {}, |
98 | 126 | "outputs": [ |
99 | 127 | { |
|
110 | 138 | "(-0.5, 99.5, 4.5, -0.5)" |
111 | 139 | ] |
112 | 140 | }, |
113 | | - "execution_count": 21, |
| 141 | + "execution_count": 4, |
114 | 142 | "metadata": {}, |
115 | 143 | "output_type": "execute_result" |
116 | 144 | }, |
117 | 145 | { |
118 | 146 | "data": { |
119 | | - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxoAAAA7CAYAAAD8WPxcAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAC59JREFUeJzt3U9sHNUdwPHf7Mzuep0YxyQEQUriEBRDcuihFKknpAQQSIAqcSlSDjkEgiAVqqCip4J868U5ISNyCAJx4wKpXFk+ABGRkapKRUL9oypYHJqQKhQ7WW/Wuzszve289xvv213zNttsvh8p0rz9vZl5ntkZz8v7vXGQpmkqAAAAAOBRYdgNAAAAADB66GgAAAAA8I6OBgAAAADv6GgAAAAA8I6OBgAAAADv6GgAAAAA8I6OBgAAAADv6GgAAAAA8I6OBgAAAADv6GgAAAAM2fHjx2V6enrYzQC8oqOBW0oQBD39++yzz4bdVAAARETkvffes35HjY2NycGDB+XUqVNy5cqVYTcPGJho2A0A+vHBBx9Y5ffff1+WlpZynz/00EM3s1kAAHQ1Ozsr+/fvl3q9Ll988YXMz8/LwsKCfP3113LmzBlJkmTYTQS8oqOBW8qxY8es8pdffilLS0u5zwEA+H/z1FNPycMPPywiIidOnJCdO3fK3NycfPzxx/L8888PuXWAf6ROYeScPXtWjhw5Irt375ZyuSyHDh2S+fn5XL0gCOStt97KfT49PS3Hjx8ffEMBALe1I0eOiIjIysoKczQwkhjRwMiZn5+Xw4cPy7PPPitRFMm5c+fk5ZdfliRJ5JVXXhl28wAAEBGRixcviojIzp07h9wSYDDoaGDkfP7551KpVNrlU6dOyZNPPilzc3N0NAAAQ7O2tiZXr16Ver0uFy5ckNnZWalUKvL000/L8vLysJsHeEdHAyPH7GSsra1Js9mURx99VBYXF2VtbU0mJyeH2DoAwO3qscces8r79u2TDz/8UPbs2TOkFgGDRUcDI+fChQvy5ptvyvLystRqNStGRwMAMCxvv/22HDx4UKIokrvvvltmZmakUGC6LEYXHQ2MlIsXL8rRo0flwQcflLm5ObnvvvukVCrJwsKCnD59uqdXB8ZxfBNaCgC43TzyyCPtt04BtwM6Ghgp586dk42NDfnkk09k79697c8//fTTXN2pqSlZXV21Pms0GnL58uVBNxMAAGDkMV6HkRKGoYiIpGna/mxtbU3Onj2bq3vgwAE5f/689dm7777LiAYAAIAHjGhgpDzxxBNSKpXkmWeekZMnT0q1WpUzZ87I7t27cyMVJ06ckJdeekmee+45efzxx+Wrr76SxcVF2bVr15BaDwAAMDoY0cBImZmZkY8++kiCIJDXX39d3nnnHXnxxRfl1VdfzdV94YUX5I033pDz58/La6+9JisrK7K0tCTbtm0bQssBAABGS5CaOSYAAAAA4AEjGgAAAAC8o6MBAAAAwDs6GgAAAAC8o6MBAAAAwDs6GgAAAAC8o6MBAAAAwLue/2DfA3+Y29oegq2tpg3kHbye2jbobaYDaaenI3rL/LwD2OaAtuvt5/d1jq1t+t/krbLNdMjXzLcnf+tn/yNmZvb0sJsAALetf/7+N844IxoAAAAAvKOjAQAAAMA7OhoAAAAAvOt5jka674ZdNhLJ00TVTYJNl0VExCzrWNqhno6pVOnA2qYdMxPeg1xMHDF7/1Z6tqpr1XTtQ6d4J46Y2W5XanhuPdc2zYLj+Pa7zx5tOcV9IBN0tmiL+fVbnnfhWC+3Tdc+9OkOOsd63qc6oVZM/RdGr/vLnWpjO/mfV+3f3Kf+ea2Y2ou5j4K+EDpvU8y6rpjaZmCUg1zMWB7EHJsRNPUP+6brvNbM71pf14ujcs/XS4/1ural97o+9nkzjlPP++un7qCOaaf1bvb+uhnmcerSTvd20h7r9R7z0W7n3fjHXHeu+/xN/+712JY+n2cY0QAAAADgHR0NAAAAAN71nDp1558qVrllFOMxexwlHjOWS/Z24nI2NJMU7VhSMobMiiotwiwXVX6SUS5E9nqFQhYLI3u9MEw2rSciEqlyaKY7qOGl0KhbUDHXCJOu20niGOvKZWOZKW1qPTMWq9Q0XTdOsj6o3n9irJskhc6x2I6Z+0jizil1qY6Z5VwqXrYYqPXMchBLx5iISp3TqXHGuoXcPhzrJZsv6/VyMdd6rnamaee6rtQ4Z4qdYz2HraZ3pI6Uq7SgvrOhKpspUKHajnGnS3TMKOt7Uho57lfGPSkt6ftVdvADdd+JStnJL481rdi2cqO9PFmuC7qr/mrNKsfGfSdW9yDznqTvXWYKcKrvXa57kBWzQ660Xuv+pFNuk84xZwqu817iuHf2sc1O289tJ5fi7Ij1cT8y63Y9Np2247zHOYK3UjajM+2mxxS3Labudk3rddY1ngX0f4U71svV7RRzpfXqmPHcl/vd5Nqmc/+dU37TsI+0XnObufWM51UVC4zn3kC123wOLoTqGdgoh6HjprAJRjQAAAAAeEdHAwAAAIB3dDQAAAAAeNfzHI1jv1uwysur97eX//XDXVbs2n+3Z4VVO7G5eC3r25RX7SS7YjUrhzfsvLKonpWjDXu9QrNgLKt8tDjddLkrnQ9ulHX+dxIFmy6LiCRFM2avF5eMWNER0/NcjLJrnotez47peS7quJnzXkr2BIcoysrFoh0rRS1j2Y6VQyM3PWxZsaIRi1Tibcmoq2PmPJeCK5m4iyTtPCelZcQacdQx1oztL8aGUXdDxRqt0Fi2t9loZrFWy14vaWTltKn+n6CprouGcV00dEw6x4xpA2FD7FjDjKn5UMZ60YaKtcx6OpZuuiyi57Ko61eXHcz5HHpuh3WNFtW5L2fllp6LVjHO/XYrJM3tRuwO+zy1dmTLkxP2a8Pvn/y+vfyLHd8Iurvnl3+3yoWJiWx527hdeTybXJiMj1mhZDy7mSYV+5qMy6GxrL4HRjkudv5uxbk5Ptlyqn9vGLvP/76xy1aOt84VN9fVMXP+UR+58K4c/kFMYXC+SXOrr3d3zhEJOsdcc0L6ma/ifLW+quua29JjzHUs8vfV3tbTnPNsBsH12tZ+YuazXW4ehlFZxcy5frk5gblrzVhPXevWdeiKqee51PEKdfOZLVLPb+Vy9ot6vGTPEZwob7SXp8o1KzZVysp3larSD0Y0AAAAAHhHRwMAAACAdz2nTv3x8JRVbh3NUqc2flq2Kx/Ihmom9l6zQofv+q69fGjishV7oHylvXxv8QcrdkeQDekU1VhjbIyF1VQu0WqSDZ1/H9v5DZeaO9rL321MWrF/13dY5Su1bDj+h3V7OH59PRuCj6v2IQ2r2ZhaWFOpYutZOVq3QhLVsqGw4ved08jCunoF2YaR8tSwY0HLeK1ZrGL9pKGYQ42RfbzTYvZdSCK7HxuXs3KtrGOdU1Ra5uuSdcx6zbLdzriSGsvqrwdX7OHEsGK+ctTOFxo3hhq3G0OLIiKTpewVpHeO2yfxjmjDWLZTZMaNnKTxgr2/spGDVNTv5XVoqvFbs1yL7Wu0apSvtexXV682s/JqQ8XqWfla3d5mrZaV46r9vSjUsrZEVfvcR8Z1ULpuhaRYzc5hqWqfw+L1WNXN8rPCqn2eghvZMQ42VD5Ywxg+1q+2DI1jWrJ/pmQiOxatKfs43diVjXOv32Ofl+pPsp+/WrGHrqf3ZKlTv576VtDd4qW/WuW/GOf3zzf2W7G/1e5tL39T3WXF/nM9u8dfV9/R5np2DgN1Hw9vZHXDuooZX8OCfart9EU7k1QiI3U4/2puVd5iqstW/+qwM1XLGestDcUZU/H+9t95H/ZrTB2vH+0jjczaRudQVz2njrnSwRxpXc7XHnt69XpB/xpLeovp5xLX/guOV8Y7U+rMYB/Xi/3qdRVzvUI9lybZOU3eSr1Uj9nxmPl8o/+UgbFv9Xp1M019l3pmuX/71fbyofFLVuznlZX28s/KqqFdMKIBAAAAwDs6GgAAAAC8o6MBAAAAwLsgTXVSMgAAAAD8OIxoAAAAAPCOjgYAAAAA7+hoAAAAAPCOjgYAAAAA7+hoAAAAAPCOjgYAAAAA7+hoAAAAAPCOjgYAAAAA7+hoAAAAAPDuf7bZzAng6bn3AAAAAElFTkSuQmCC", |
| 147 | + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxoAAAA7CAYAAAD8WPxcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAALn0lEQVR4nO3dT2wc1R3A8d/szO56nRjHJARBSuIQFENy6KEUqSekBBBIgCpxKVIOOQSCIBWqoKKngnzrxTkhI3IIAnHjAqlcWT4AEZGRqkpFQv2jKlgcmpAqFDtZb9a7OzO97bz3G+/bXfM222y+HynSvP29mXme2RnPy/u9cZCmaSoAAAAA4FFh2A0AAAAAMHroaAAAAADwjo4GAAAAAO/oaAAAAADwjo4GAAAAAO/oaAAAAADwjo4GAAAAAO/oaAAAAADwjo4GAAAAAO/oaAAAAAzZ8ePHZXp6etjNALyio4FbShAEPf377LPPht1UAABEROS9996zfkeNjY3JwYMH5dSpU3LlypVhNw8YmGjYDQD68cEHH1jl999/X5aWlnKfP/TQQzezWQAAdDU7Oyv79++Xer0uX3zxhczPz8vCwoJ8/fXXcubMGUmSZNhNBLyio4FbyrFjx6zyl19+KUtLS7nPAQD4f/PUU0/Jww8/LCIiJ06ckJ07d8rc3Jx8/PHH8vzzzw+5dYB/pE5h5Jw9e1aOHDkiu3fvlnK5LIcOHZL5+flcvSAI5K233sp9Pj09LcePHx98QwEAt7UjR46IiMjKygpzNDCSGNHAyJmfn5fDhw/Ls88+K1EUyblz5+Tll1+WJEnklVdeGXbzAAAQEZGLFy+KiMjOnTuH3BJgMOhoYOR8/vnnUqlU2uVTp07Jk08+KXNzc3Q0AABDs7a2JlevXpV6vS4XLlyQ2dlZqVQq8vTTT8vy8vKwmwd4R0cDI8fsZKytrUmz2ZRHH31UFhcXZW1tTSYnJ4fYOgDA7eqxxx6zyvv27ZMPP/xQ9uzZM6QWAYNFRwMj58KFC/Lmm2/K8vKy1Go1K0ZHAwAwLG+//bYcPHhQoiiSu+++W2ZmZqRQYLosRhcdDYyUixcvytGjR+XBBx+Uubk5ue+++6RUKsnCwoKcPn26p1cHxnF8E1oKALjdPPLII+23TgG3AzoaGCnnzp2TjY0N+eSTT2Tv3r3tzz/99NNc3ampKVldXbU+azQacvny5UE3EwAAYOQxXoeREoahiIikadr+bG1tTc6ePZure+DAATl//rz12bvvvsuIBgAAgAeMaGCkPPHEE1IqleSZZ56RkydPSrValTNnzsju3btzIxUnTpyQl156SZ577jl5/PHH5auvvpLFxUXZtWvXkFoPAAAwOhjRwEiZmZmRjz76SIIgkNdff13eeecdefHFF+XVV1/N1X3hhRfkjTfekPPnz8trr70mKysrsrS0JNu2bRtCywEAAEZLkJo5JgAAAADgASMaAAAAALyjowEAAADAOzoaAAAAALyjowEAAADAOzoaAAAAALyjowEAAADAu57/YN8Df5jb2h6Cra2mDeQdvJ7aNuhtpgNpp6cjesv8vAPY5oC26+3n93WOrW363+Stss10yNfMtyd/62f/I2Zm9vSwmwAAt61//v43zjgjGgAAAAC8o6MBAAAAwDs6GgAAAAC863mORrrvhl02EsnTRNVNgk2XRUTELOtY2qGejqlU6cDaph0zE96DXEwcMXv/Vnq2qmvVdO1Dp3gnjpjZbldqeG491zbNguP49rvPHm05xX0gE3S2aIv59Vued+FYL7dN1z706Q46x3repzqhVkz9F0av+8udamM7+Z9X7d/cp/55rZjai7mPgr4QOm9TzLqumNpmYJSDXMxYHsQcmxE09Q/7puu81szvWl/Xi6Nyz9dLj/W6tqX3uj72eTOOU8/766fuoI5pp/Vu9v66GeZx6tJO93bSHuv1HvPRbufd+Mdcd677/E3/7vXYlj6fZxjRAAAAAOAdHQ0AAAAA3vWcOnXnnypWuWUU4zF7HCUeM5ZL9nbicjY0kxTtWFIyhsyKKi3CLBdVfpJRLkT2eoVCFgsje70wTDatJyISqXJopjuo4aXQqFtQMdcIk67bSeIY68plY5kpbWo9Mxar1DRdN06yPqjef2KsmySFzrHYjpn7SOLOKXWpjpnlXCpethio9cxyEEvHmIhKndOpcca6hdw+HOslmy/r9XIx13qudqZp57qu1Dhnip1jPYetpnekjpSrtKC+s6EqmylQodqOcadLdMwo63tSGjnuV8Y9KS3p+1V28AN134lK2ckvjzWt2LZyo708Wa4Luqv+as0qx8Z9J1b3IPOepO9dZgpwqu9drnuQFbNDrrRe6/6kU26TzjFnCq7zXuK4d/axzU7bz20nl+LsiPVxPzLrdj02nbbjvMc5grdSNqMz7abHFLctpu52Tet11jWeBfR/hTvWy9XtFHOl9eqY8dyX+93k2qZz/51TftOwj7Rec5u59YznVRULjOfeQLXbfA4uhOoZ2CiHoeOmsAlGNAAAAAB4R0cDAAAAgHd0NAAAAAB41/McjWO/W7DKy6v3t5f/9cNdVuzaf7dnhVU7sbl4LevblFftJLtiNSuHN+y8sqielaMNe71Cs2Asq3y0ON10uSudD26Udf53EgWbLouIJEUzZq8Xl4xY0RHT81yMsmuei17Pjul5Luq4mfNeSvYEhyjKysWiHStFLWPZjpVDIzc9bFmxohGLVOJtyairY+Y8l4IrmbiLJO08J6VlxBpx1DHWjO0vxoZRd0PFGq3QWLa32WhmsVbLXi9pZOW0qf6foKmui4ZxXTR0TDrHjGkDYUPsWMOMqflQxnrRhoq1zHo6lm66LKLnsqjrV5cdzPkcem6HdY0W1bkvZ+WWnotWMc79diskze1G7A77PLV2ZMuTE/Zrw++f/L69/Isd3wi6u+eXf7fKhYmJbHnbuF15PJtcmIyPWaFkPLuZJhX7mozLobGsvgdGOS52/m7FuTk+2XKqf28Yu8//vrHLVo63zhU319Uxc/5RH7nwrhz+QUxhcL5Jc6uvd3fOEQk6x1xzQvqZr+J8tb6q65rb0mPMdSzy99Xe1tOc82wGwfXa1n5i5rNdbh6GUVnFzLl+uTmBuWvNWE9d69Z16Iqp57nU8Qp185ktUs9v5XL2i3q8ZM8RnChvtJenyjUrNlXKyneVqtIPRjQAAAAAeEdHAwAAAIB3PadO/fHwlFVuHc1SpzZ+WrYrH8iGaib2XrNCh+/6rr18aOKyFXugfKW9fG/xByt2R5AN6RTVWGNsjIXVVC7RapINnX8f2/kNl5o72svfbUxasX/Xd1jlK7VsOP6HdXs4fn09G4KPq/YhDavZmFpYU6li61k5WrdCEtWyobDi953TyMK6egXZhpHy1LBjQct4rVmsYv2koZhDjZF9vNNi9l1IIrsfG5ezcq2sY51TVFrm65J1zHrNst3OuJIay+qvB1fs4cSwYr5y1M4XGjeGGrcbQ4siIpOl7BWkd47bJ/GOaMNYtlNkxo2cpPGCvb+ykYNU1O/ldWiq8VuzXIvta7RqlK+17FdXrzaz8mpDxepZ+Vrd3matlpXjqv29KNSytkRV+9xHxnVQum6FpFjNzmGpap/D4vVY1c3ys8KqfZ6CG9kxDjZUPljDGD7Wr7YMjWNasn+mZCI7Fq0p+zjd2JWNc6/fY5+X6k+yn79asYeup/dkqVO/nvpW0N3ipb9a5b8Y5/fPN/Zbsb/V7m0vf1PdZcX+cz27x19X39HmenYOA3UfD29kdcO6ihlfw4J9qu30RTuTVCIjdTj/am5V3mKqy1b/6rAzVcsZ6y0NxRlT8f7233kf9mtMHa8f7SONzNpG51BXPaeOudLBHGldztcee3r1ekH/Gkt6i+nnEtf+C45XxjtT6sxgH9eL/ep1FXO9Qj2XJtk5Td5KvVSP2fGY+Xyj/5SBsW/1enUzTX2Xema5f/vV9vKh8UtW7OeVlfbyz8qqoV0wogEAAADAOzoaAAAAALyjowEAAADAuyBNdVIyAAAAAPw4jGgAAAAA8I6OBgAAAADv6GgAAAAA8I6OBgAAAADv6GgAAAAA8I6OBgAAAADv6GgAAAAA8I6OBgAAAADv6GgAAAAA8O5/ttnMCeDpufcAAAAASUVORK5CYII=", |
120 | 148 | "text/plain": [ |
121 | 149 | "<Figure size 1000x1000 with 2 Axes>" |
122 | 150 | ] |
|
151 | 179 | ], |
152 | 180 | "metadata": { |
153 | 181 | "kernelspec": { |
154 | | - "display_name": ".venv", |
| 182 | + "display_name": "py_env_book", |
155 | 183 | "language": "python", |
156 | 184 | "name": "python3" |
157 | 185 | }, |
|
165 | 193 | "name": "python", |
166 | 194 | "nbconvert_exporter": "python", |
167 | 195 | "pygments_lexer": "ipython3", |
168 | | - "version": "3.10.12" |
| 196 | + "version": "3.10.15" |
169 | 197 | } |
170 | 198 | }, |
171 | 199 | "nbformat": 4, |
|
0 commit comments