Skip to content

Commit 788ed54

Browse files
nikiquarknikiquark
committed
fix plasma particles probleb
1 parent 1a8f1a9 commit 788ed54

File tree

4 files changed

+273
-9
lines changed

4 files changed

+273
-9
lines changed

‎beam-compare.ipynb

+17-7
Large diffs are not rendered by default.

‎lcode2dPy/config/default_config.py

-1
Original file line numberDiff line numberDiff line change
@@ -7,7 +7,6 @@
77

88
# Parameters of simulation window:
99
default_config.set('window-width', 5.0) # Has different meanings in 2d and 3d!
10-
default_config.set('r-step', 5.0)
1110
default_config.set('window-width-step-size', 0.05)
1211
default_config.set('window-width-steps', 101) # Only needed in 3d
1312
# TODO: get rid of this parameter in 3d

‎lcode2dPy/config/template.py

+3-1
Original file line numberDiff line numberDiff line change
@@ -2,7 +2,7 @@
22
lcode_template = '''
33
# Simulation area:
44
geometry = c
5-
window-width = {window-width}; r-step = {r-step}
5+
window-width = {window-width}; r-step = {window-width-step-size}
66
window-length = {window-length}; xi-step = {xi-step}
77
time-limit = {time-limit}; time-step = {time-step}
88
continuation = n # Plasma continuation (no/beam/longplasma, n/y/Y)
@@ -22,6 +22,8 @@
2222
plasma-profile = {plasma-profile}
2323
plasma-zshape = {plasma-zshape}
2424
plasma-width = {plasma-width}
25+
plasma-width-2 = {plasma-width-2}
26+
plasma-density-2 = {plasma-density-2}
2527
plasma-temperature = {plasma-temperature}
2628
ion-model = {ion-model}
2729
ion-mass = 157000

‎plasma_checker.ipynb

+253
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,253 @@
1+
{
2+
"cells": [
3+
{
4+
"cell_type": "code",
5+
"execution_count": 1,
6+
"metadata": {
7+
"Collapsed": "false"
8+
},
9+
"outputs": [],
10+
"source": [
11+
"from lcode2dPy.simulation.interface import Simulation\n",
12+
"from lcode2dPy.diagnostics.targets import BeamDiagnostics, PlasmaDiagnostics\n",
13+
"from lcode2dPy.config.default_config import default_config\n",
14+
"from lcode2dPy.beam.beam_generator import make_beam, Gauss, rGauss\n",
15+
"import numpy as np\n",
16+
"import matplotlib.pyplot as plt\n",
17+
"import subprocess"
18+
]
19+
},
20+
{
21+
"cell_type": "code",
22+
"execution_count": 2,
23+
"metadata": {
24+
"Collapsed": "false"
25+
},
26+
"outputs": [],
27+
"source": [
28+
"# Параметры\n",
29+
"\n",
30+
"time_step = 200\n",
31+
"time_limit = 200.5\n",
32+
"window_length = 0.2\n",
33+
"window_width = 2\n",
34+
"r_step = 0.05\n",
35+
"xi_step = 0.05"
36+
]
37+
},
38+
{
39+
"cell_type": "code",
40+
"execution_count": 3,
41+
"metadata": {
42+
"Collapsed": "false"
43+
},
44+
"outputs": [],
45+
"source": [
46+
"!rm ../c_code/*"
47+
]
48+
},
49+
{
50+
"cell_type": "code",
51+
"execution_count": 4,
52+
"metadata": {
53+
"Collapsed": "false"
54+
},
55+
"outputs": [],
56+
"source": [
57+
"# Config\n",
58+
"config = default_config\n",
59+
"config.set('time-step', time_step)\n",
60+
"config.set('time-limit', time_limit)\n",
61+
"config.set('window-length', window_length)\n",
62+
"config.set('window-width', window_width)\n",
63+
"config.set('window-width-step-size', r_step)\n",
64+
"config.set('xi-step', xi_step)\n",
65+
"\n",
66+
"# Beam\n",
67+
"gamma = 426\n",
68+
"angspread = 1e-5\n",
69+
"m_proton = 958/0.51\n",
70+
"\n",
71+
"beam_pars = dict(xi_distr=Gauss(sigma=100, vmin=-window_length, vmax=0),\n",
72+
" r_distr=rGauss(vmin=0, vmax=window_width),\n",
73+
" ang_distr=Gauss(sigma=angspread, vmin=None, vmax=None),\n",
74+
" pz_distr=Gauss(gamma*m_proton, gamma*m_proton*1e-4, vmin=None, vmax=None),\n",
75+
" Ipeak_kA=2*40/1000,\n",
76+
" q_m=1/m_proton,\n",
77+
" saveto=\"../c_code\")\n",
78+
"\n",
79+
"diagnostics = [\n",
80+
" BeamDiagnostics(period=time_limit//time_step * time_step),\n",
81+
" PlasmaDiagnostics(period=time_limit//time_step * time_step)\n",
82+
"]\n",
83+
"\n",
84+
"sim = Simulation(beam_pars=beam_pars, diagnostics=diagnostics, config=config)"
85+
]
86+
},
87+
{
88+
"cell_type": "code",
89+
"execution_count": 5,
90+
"metadata": {
91+
"Collapsed": "false"
92+
},
93+
"outputs": [
94+
{
95+
"name": "stdout",
96+
"output_type": "stream",
97+
"text": [
98+
"Number of particles: 10093\n",
99+
"Number of particles in the middle layer: 1261\n"
100+
]
101+
}
102+
],
103+
"source": [
104+
"sim.step(int(time_limit//time_step))"
105+
]
106+
},
107+
{
108+
"cell_type": "code",
109+
"execution_count": 6,
110+
"metadata": {
111+
"Collapsed": "false"
112+
},
113+
"outputs": [
114+
{
115+
"data": {
116+
"text/plain": [
117+
"0"
118+
]
119+
},
120+
"execution_count": 6,
121+
"metadata": {},
122+
"output_type": "execute_result"
123+
}
124+
],
125+
"source": [
126+
"# Готовим C диагностики\n",
127+
"!rm ../c_code/*\n",
128+
"!cp beamfile.bin ../c_code/beamfile.bin\n",
129+
"config.c_config('../c_code/lcode.cfg')\n",
130+
"!echo 0 > ../c_code/beamfile.bit \n",
131+
"\n",
132+
"p = subprocess.Popen([\"/home/global/lcode/lcode\", \"lcode.cfg\"], cwd=\"../c_code\")\n",
133+
"p.wait()"
134+
]
135+
},
136+
{
137+
"cell_type": "code",
138+
"execution_count": 7,
139+
"metadata": {
140+
"Collapsed": "false"
141+
},
142+
"outputs": [],
143+
"source": [
144+
"c_plasma=np.loadtxt('../c_code/00200.pls')\n",
145+
"xis = np.arange(-xi_step,-window_length,-xi_step)\n",
146+
"xis = np.unique(c_plasma.T[0])[::-1]\n",
147+
"c_rs = np.array([c_plasma[c_plasma.T[0]==xi].T[2] for xi in xis])\n",
148+
"plasma = diagnostics[1].data\n",
149+
"py_pls = plasma[time_step]\n",
150+
"py_rs = np.array([layer.r for layer in py_pls])"
151+
]
152+
},
153+
{
154+
"cell_type": "code",
155+
"execution_count": 8,
156+
"metadata": {
157+
"Collapsed": "false"
158+
},
159+
"outputs": [
160+
{
161+
"data": {
162+
"text/plain": [
163+
"(0, 0.1)"
164+
]
165+
},
166+
"execution_count": 8,
167+
"metadata": {},
168+
"output_type": "execute_result"
169+
},
170+
{
171+
"data": {
172+
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAWjUlEQVR4nO3dfZBddX3H8fcnu8lmk0CSDRsISSCRBGMi+LQEqSDWKAS1xgdSE61GixN0iJ360AKj9QHRSnGkrWJrFBxEO4A4tqlEUwV8qLWYjSgxYGBBa1Yc2ZAIBAhhk2//uCfT7fVu9uzdc3NOfvN5zezk3vPwu5+b7P3szTnn/lYRgZmZpWtc2QHMzKy1XPRmZolz0ZuZJc5Fb2aWOBe9mVniXPRmZonLVfSSlkvaLqlP0iUN1r9Y0k8kDUo6v27dGkn3ZV9rigpuZmb5aKTr6CW1AfcCLwf6gc3A6oi4e8g284CjgfcBGyLi5mx5F9AL9AABbAFeEBG7i34iZmbWWJ539EuBvoh4ICL2ATcAK4ZuEBG/ioi7gAN1+54LfDsidmXl/m1geQG5zcwsp/Yc28wGdgy53w+cnnP8RvvOrt9I0lpgLcDkyZNfsGjRopzDm5kZwJYtW3ZGRHejdXmKXg2W5Z03Ide+EbEeWA/Q09MTvb29OYc3MzMASf8z3Lo8h276gblD7s8BHsz52GPZ18zMCpCn6DcDCyXNlzQBWAVsyDn+JuAcSdMlTQfOyZaZmdlhMmLRR8QgsI5aQd8D3BQR2yRdJunVAJJOk9QPrAQ+J2lbtu8u4KPUflhsBi7LlpmZ2WEy4uWVh1uzx+jveuwJVv70/hYkGp3Iffqitaryz1qRGIXwc7FWe2X3VD79rBOb2lfSlojoabQuz8nYI8L08e2sPG562TGAxmegy6CKJKlGioIk9GQSeirJePaUzpaMm0zRz504gcsXzik7hplZ5XiuGzOzxLnozcwS56I3M0uci97MLHEuejOzxLnozcwS56I3M0tcMtfRHzgwyP79j5cdoyApfW4xnedStU+Rj01KzyUd48ZNoL39qMLHTabo9+y5m829ry07hplZ02bOfCWnPPsfCx83maLv6JjFwoUfKDtGYaoyfUExEnouSui5pPTvkohJnc3NczOShIq+mxPmvq3sGGZmleOTsWZmiXPRm5klzkVvZpY4F72ZWeJc9GZmiXPRm5klzkVvZpa4ZK6j37dvHzt37iw7htVJa9oAs9bq7Oykq6ur8HGTKfqBgQE+//nPlx3DzKxpS5YsYeXKlYWPm0zRd3V1sWrVqrJjFEZJfdQ+Hf53sVaaMmVKS8ZNpug7OztZtGhR2THMzCrHJ2PNzBLnojczS5yL3swscS56M7PEuejNzBLnojczS5yL3swscclcR3/giafZe/8jZccoSELTBiT0VMxarW1qBx0nHl34uMkU/eDDe9n1lXvKjmFm1rTOU49x0R9K+7GTOPYvn192DGvEswaY5aKJrankZIp+3IQ2xh03uewYZmaVk+tkrKTlkrZL6pN0SYP1HZJuzNbfIWletny8pOskbZV0j6RLi41vZmYjGbHoJbUBVwPnAYuB1ZIW1212AbA7IhYAVwFXZMtXAh0RcQrwAuDCgz8EzMzs8Mjzjn4p0BcRD0TEPuAGYEXdNiuA67LbNwPLVJvPNYDJktqBTmAf8Gghyc3MLJc8RT8b2DHkfn+2rOE2ETEIPALMoFb6jwO/BX4NfDIidtU/gKS1knol9Q4MDIz6SZiZ2fDyFH2jaybqr44ebpulwH7geGA+8F5Jz/iDDSPWR0RPRPR0d3fniGRmZnnlKfp+YO6Q+3OAB4fbJjtMMxXYBbwR+FZEPB0RDwE/BHrGGtrMzPLLU/SbgYWS5kuaAKwCNtRtswFYk90+H7gtar8V+tfAS1UzGXgh8ItiopuZWR4jXkcfEYOS1gGbgDbg2ojYJukyoDciNgDXANdL6qP2Tv7gL2+9Gvgi8HNqh3e+GBF3teB58Pjvd3P3D25vxdCHX6Qzb0Ak9FzMWq3r+DksOO2FhY+b6wNTEbER2Fi37INDbu+ldill/X57Gi1vhcd2DvD9L197OB7KzKwlnnnGWeUV/ZFg5jNO4l3XfbXsGIVRSvMGJPRUUEpPxqpm3LjWTCicTNGPG9fGhImdZccwM6scz0dvZpY4F72ZWeJc9GZmiXPRm5klzkVvZpY4F72ZWeJc9GZmiUvmOvpHH36SO//j12XH+MN5PUtSkRhWVZ6aopJmnng0i888vvBxkyn6px4fpG/LQ2XHAPzhSTtC+Bu1cgJc9IfSfcJRXPDJs8qOYWZWOT5Gb2aWOBe9mVniXPRmZolz0ZuZJc5Fb2aWOBe9mVniXPRmZolz0ZuZJS6ZD0w99ctf8ruP/23ZMYrhj6dbq/l7rJImnXYax7zjwsLHTabo2b+f/Y88UnaK4vjT6dZiSf0C+kQceGpvS8ZNpug7Fixg/k03lh3DzKxyfIzezCxxLnozs8S56M3MEueiNzNLnIvezCxxLnozs8S56M3MEueiNzNLXDIfmLp3972897vvLTuGmVnTzpx9JhcvvbjwcZMp+s62ThZ1LSo7hjXgj9qb5TPnqDktGTeZop979FyuPPvKsmOYmVVOrmP0kpZL2i6pT9IlDdZ3SLoxW3+HpHlD1p0q6UeStknaKmlicfHNzGwkIxa9pDbgauA8YDGwWtLius0uAHZHxALgKuCKbN924MvAOyJiCfAS4OnC0puZ2YjyvKNfCvRFxAMRsQ+4AVhRt80K4Lrs9s3AMkkCzgHuioifAUTEwxGxv5joZmaWR56inw3sGHK/P1vWcJuIGAQeAWYAJwMhaZOkn0j660YPIGmtpF5JvQMDA6N9DmZmdgh5ir7RJRP1v55muG3agTOBN2V/vlbSsj/YMGJ9RPRERE93d3eOSGZmlleeou8H5g65Pwd4cLhtsuPyU4Fd2fLvRcTOiHgC2Ag8f6yhzcwsvzxFvxlYKGm+pAnAKmBD3TYbgDXZ7fOB2yIigE3AqZImZT8AzgbuLia6mZnlMeJ19BExKGkdtdJuA66NiG2SLgN6I2IDcA1wvaQ+au/kV2X77pb0KWo/LALYGBG3tOi5mJlZA4qK/Tb4np6e6O3tHf2Ov/kJXHtu8YFsbCr2/WVWaUteA6//QlO7StoSET2N1iXzyVimzIQzLio7hTXkKRDMcjl2SUuGTafop86Bl3247BRmZpXjaYrNzBLnojczS5yL3swscS56M7PEuejNzBLnojczS5yL3swscclcR79v8AC7n9hXdgyr4w/GmuU3cfw4pk2aUPi4yRT9Pb99lBVX/7DsGGZmTXvVqbP4zBuLn+A3maKfPb2Tj7/2lLJjWAPyDAhmuZzQNakl4yZT9MdM6eCNp59Qdgwzs8rxyVgzs8S56M3MEueiNzNLnIvezCxxLnozs8S56M3MEueiNzNLXDLX0e8Z3M/de54sO4bV8QwIZvl1jW9n4eSJhY+bTNHf98RTvPrOvrJjmJk1bcXMaXxuybzCx02m6E+a1MGNzzmp7BjWgGdAMMune0JrKjmZoj+6vY2zu44qO4aZWeX4ZKyZWeJc9GZmiXPRm5klzkVvZpY4F72ZWeJc9GZmiXPRm5klLpnr6J9+ejc7H/5e2THMjgzhySmqaOLE2UyfvrTwcZMp+ief3MHdd7+37BhmZk2bOfOVLvpDmTz5mZzxwlvLjmF2BPHkFFXT1j65JeMmU/RtbR1MmjSv7BhmZpWT62SspOWStkvqk3RJg/Udkm7M1t8haV7d+hMk7ZH0vmJim5lZXiMWvaQ24GrgPGAxsFrS4rrNLgB2R8QC4Crgirr1VwHfHHtcMzMbrTzv6JcCfRHxQETsA24AVtRtswK4Lrt9M7BMkgAkvQZ4ANhWTGQzMxuNPEU/G9gx5H5/tqzhNhExCDwCzJA0GbgY+MihHkDSWkm9knoHBgbyZjczsxzyFH2jU/P1F+EOt81HgKsiYs+hHiAi1kdET0T0dHd354hkZmZ55bnqph+YO+T+HODBYbbpl9QOTAV2AacD50v6O2AacEDS3oj4zJiTm5lZLnmKfjOwUNJ84DfAKuCNddtsANYAPwLOB26LiADOOriBpA8De1zyZmaH14hFHxGDktYBm4A24NqI2CbpMqA3IjYA1wDXS+qj9k5+VStDN/Loo49y5513Hu6HNTMrTHd3N4sX11/UOHa5PjAVERuBjXXLPjjk9l5g5QhjfLiJfLk99thj3H777a18CDOzllqyZElLil5RscmNenp6ore3d9T7RQRVey5mZqM1blxzkwpL2hIRPY3WJTMFgiSyS/fNzGwIz0dvZpY4F72ZWeJc9GZmiXPRm5klzkVvZpY4F72ZWeJc9GZmiUvmOvrBXXt59NZflx3DzKxpE+YexZQXzip83GSK/sDeQZ7q+33ZMczMmqbxrTnIkkzRTzh+CrMuXVp2DDOzyvExejOzxLnozcwS56I3M0uci97MLHEuejOzxLnozcwS56I3M0uci97MLHHJfGDq4d/s4NYvfLbsGGZmTZu75FTOOH914eMmU/QE/uXgZnZEa1WHJVP0M+bM5Q0f/kTZMczMKsfH6M3MEueiNzNLnIvezCxxLnozs8S56M3MEueiNzNLnIvezCxxLnozs8Ql84Gpnf2PsfGzW8uOYWbWtHmnHsOLV51c+LjJFP34jnZmnzyt7BhmZk3rmjWpJeMmU/RTuztZ9tbFZccwM6ucXMfoJS2XtF1Sn6RLGqzvkHRjtv4OSfOy5S+XtEXS1uzPlxYb38zMRjJi0UtqA64GzgMWA6sl1b91vgDYHRELgKuAK7LlO4E/iYhTgDXA9UUFNzOzfPK8o18K9EXEAxGxD7gBWFG3zQrguuz2zcAySYqIOyPiwWz5NmCipI4igpuZWT55in42sGPI/f5sWcNtImIQeASYUbfN64E7I+Kp+geQtFZSr6TegYGBvNnNzCyHPEWvBsvqZ8c/5DaSllA7nHNhoweIiPUR0RMRPd3d3TkimZlZXnmKvh+YO+T+HODB4baR1A5MBXZl9+cAXwfeEhH3jzWwmZmNTp6i3wwslDRf0gRgFbChbpsN1E62ApwP3BYRIWkacAtwaUT8sKjQZmaW34hFnx1zXwdsAu4BboqIbZIuk/TqbLNrgBmS+oD3AAcvwVwHLAD+RtJPs6+ZhT8LMzMblqr2C7V7enqit7d31Ps9uXUrv1r5py1IZGZ2eBz9ivOY/alPNbWvpC0R0dNoXTKfjG2fOZNjLrqo7BhmZk3rWLiwJeMmU/Tjjz2W7netKzuGmVnleJpiM7PEuejNzBLnojczS5yL3swscS56M7PEuejNzBLnojczS1wy19Hv27+Ph554qOwYZmZN62zvZEZn/QzvY5dM0d+7+15W37K67BhmZk1bPm85V559ZeHjJlP0s6fM5vIXXV52DDOzph0/5fiWjJtM0U+fOJ0VC+p/w6GZmflkrJlZ4lz0ZmaJc9GbmSXORW9mljgXvZlZ4lz0ZmaJc9GbmSUumevo2fsoPHhn2SnMzJo3uRuOXVz4sOkU/cP3wZdeXXYKM7PmLXkdrPxi4cOmU/THnAxv3Vh2CjOz5k0+piXDplP0HUfBvBeVncLMrHJ8MtbMLHEuejOzxLnozcwS56I3M0uci97MLHEuejOzxLnozcwSl8x19A/veYrv3PO7smOYmTVtbtck/uik4j80lUzR9+9+kou/trXsGGZmTXvVqbNc9IeyaNZR/NclLy07hplZ0yaOb2vJuMkUfUd7G8dP6yw7hplZ5eQ6GStpuaTtkvokXdJgfYekG7P1d0iaN2Tdpdny7ZLOLS66mZnlMWLRS2oDrgbOAxYDqyXVT5h8AbA7IhYAVwFXZPsuBlYBS4DlwGez8czM7DDJ845+KdAXEQ9ExD7gBmBF3TYrgOuy2zcDyyQpW35DRDwVEb8E+rLxzMzsMMlzjH42sGPI/X7g9OG2iYhBSY8AM7Ll/1237+z6B5C0Flib3d0jaXuu9I0dA+wcw/6tVvV8UP2MVc8H1c9Y9XzgjKN14nAr8hS9GiyLnNvk2ZeIWA+sz5FlRJJ6I6KniLFaoer5oPoZq54Pqp+x6vnAGYuU59BNPzB3yP05wIPDbSOpHZgK7Mq5r5mZtVCeot8MLJQ0X9IEaidXN9RtswFYk90+H7gtIiJbviq7Kmc+sBD4cTHRzcwsjxEP3WTH3NcBm4A24NqI2CbpMqA3IjYA1wDXS+qj9k5+VbbvNkk3AXcDg8BFEbG/Rc/loEIOAbVQ1fNB9TNWPR9UP2PV84EzFka1N95mZpYqz15pZpY4F72ZWeKOuKKX1CXp25Luy/6c3mCb50r6kaRtku6S9IYh6+Zn0zTcl03bMKGMjNl235L0e0nfqFu+TNJPJP1U0n9KWlDBjJL0MUn3SrpH0l9UKd+Q9Z+WtKfIbEVllPSVbGqQn0u6VtL4iuWr0mtlTbbNfZLWDFm+WtLW7HX+LUmFTv1YQL4JktZnr5NfSHp9kfnyOuKKHrgEuDUiFgK3ZvfrPQG8JSIOTr3w95KmZeuuAK7K9t9NbfqGMjICXAm8ucHyfwLeFBHPBf4F+EAFM76V2qWziyLiWdQ+MV2lfEjqAaY1WleQsWb8CrAIOAXoBN5esXyVeK1I6gI+RO2DmkuBD0manl3K/Q/AH0fEqcBdwLqq5MtWvx94KCJOpjaFzPcKzpdPRBxRX8B2YFZ2exawPcc+P6N2aaeofYqtPVt+BrCpzIzAS4BvNNj/9Oz2pcDHK5jxx8CCKvw7D5OvDbg923dPFTPWrX838LGq5KvSawVYDXxuyP3PZcvGAwPUPhEq4J+BtVXJl93eAUxuxfffaL6OxHf0x0bEbwGyP2ceamNJS4EJwP3UpmX4fUQMZqsbTslwuDM28HZgo6R+au+0PlFwPhh7xpOAN0jqlfRNSQsrlm8dsOHgGC0y1owAZIds3gx8q8BsMLZ8VXqtNJqGZXZEPA28E9hK7YOYi6ld6l2JfEOOInw0OxT7VUnHFpwvl0rORy/pO8BxDVa9f5TjzAKuB9ZExAFJuaZkyDl2IRmH8W7gFRFxh6S/Aj5FE/+tb3HGDmBvRPRIeh1wLXBWFfJJOh5YSe1d6pi0+O/woM8C34+IH4x2xxbmq9JrpWGW7AfkO4HnAQ8An6b2P+DLq5CPWr/OAX4YEe+R9B7gkwxzqLGVKln0EfGy4dZJ+p2kWRHx26zIHxpmu6OBW4APRMTBidV2AtMktWfvVJqekqGIjMPs2w08JyLuyBbdSJPv9FqVMdMPfC27/XXgixXK9zxgAdCX/WyfJKkvatNoVyXjwTE+BHQDF4523xbnq9JrpZ///0N7DvBd4LnZ+PdnY93E8Ochysj3MLXzhV/Pln+V1pznGNGReOhm6HQLa4B/q98guzrg68CXIuKrB5dH7aDZ7dSmaRh2/8OR8RB2A1MlnZzdfzlwT4HZDhpLRoB/BQ7+7sazgXsLynVQ0/ki4paIOC4i5kXEPOCJZkq+lRkBJL0dOJfa8dwDBWeDsf0dVum1sgk4JzsBOx04J1v2G2Bx9uYIWvNaaTpf9nf47/zfD4Fl1GYJOPzKPkkw2i9qxw5vBe7L/uzKlvcAX8hu/xnwNPDTIV/PzdY9g9qJxD5qP2E7ysiY3f8BtZNJT1J7V3Butvy11I47/ozaO4NnVDDjNGr/Y9oK/Ija/0Iqk69urFadjB3r3+EgtXNHB79HP1ixfFV6rfx5lqMPeNuQ5e+gVu53USvVGRXLdyLw/SzfrcAJrfheHOnLUyCYmSXuSDx0Y2Zmo+CiNzNLnIvezCxxLnozs8S56M3MEueiNzNLnIvezCxx/wsiAPmSNcJIWQAAAABJRU5ErkJggg==",
173+
"text/plain": [
174+
"<Figure size 432x288 with 1 Axes>"
175+
]
176+
},
177+
"metadata": {
178+
"needs_background": "light"
179+
},
180+
"output_type": "display_data"
181+
}
182+
],
183+
"source": [
184+
"plt.plot(xis, c_rs);\n",
185+
"plt.ylim(0, 0.1)"
186+
]
187+
},
188+
{
189+
"cell_type": "code",
190+
"execution_count": 9,
191+
"metadata": {
192+
"Collapsed": "false"
193+
},
194+
"outputs": [],
195+
"source": [
196+
"plt.plot(xis, py_rs);\n",
197+
"plt.ylim(0, 0.1)"
198+
]
199+
},
200+
{
201+
"cell_type": "code",
202+
"execution_count": null,
203+
"metadata": {},
204+
"outputs": [
205+
{
206+
"ename": "Error",
207+
"evalue": "Kernel is dead",
208+
"output_type": "error",
209+
"traceback": [
210+
"Error: Kernel is dead",
211+
"at g._sendKernelShellControl (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:52:1006305)",
212+
"at g.sendShellMessage (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:52:1006074)",
213+
"at g.requestExecute (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:52:1008616)",
214+
"at d.requestExecute (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:37:328037)",
215+
"at S.requestExecute (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:32:19306)",
216+
"at w.executeCodeCell (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:52:300924)",
217+
"at w.execute (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:52:300551)",
218+
"at w.start (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:52:296215)",
219+
"at async t.CellExecutionQueue.executeQueuedCells (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:52:310950)",
220+
"at async t.CellExecutionQueue.start (/home/okhotnikov/.vscode-server/extensions/ms-toolsai.jupyter-2021.9.1101343141/out/client/extension.js:52:310490)"
221+
]
222+
}
223+
],
224+
"source": [
225+
"diff = c_rs - py_rs\n",
226+
"plt.plot(xis, diff);"
227+
]
228+
}
229+
],
230+
"metadata": {
231+
"interpreter": {
232+
"hash": "40d3a090f54c6569ab1632332b64b2c03c39dcf918b08424e98f38b5ae0af88f"
233+
},
234+
"kernelspec": {
235+
"display_name": "Python 3.7.4 64-bit ('base': conda)",
236+
"name": "python3"
237+
},
238+
"language_info": {
239+
"codemirror_mode": {
240+
"name": "ipython",
241+
"version": 3
242+
},
243+
"file_extension": ".py",
244+
"mimetype": "text/x-python",
245+
"name": "python",
246+
"nbconvert_exporter": "python",
247+
"pygments_lexer": "ipython3",
248+
"version": "3.7.4"
249+
}
250+
},
251+
"nbformat": 4,
252+
"nbformat_minor": 4
253+
}

0 commit comments

Comments
 (0)