diff --git a/.gitignore b/.gitignore
index adb55f2..955ed56 100644
--- a/.gitignore
+++ b/.gitignore
@@ -183,4 +183,4 @@ _deps
# Blender files
*.blend
-*.blend1
\ No newline at end of file
+*.blend1!/utk_experiments/utk/cmake-build-debug/
diff --git a/.idea/.name b/.idea/.name
index e5bf71f..deb6568 100644
--- a/.idea/.name
+++ b/.idea/.name
@@ -1 +1 @@
-psa
\ No newline at end of file
+utk_tests
\ No newline at end of file
diff --git a/.idea/codeStyles/Project.xml b/.idea/codeStyles/Project.xml
deleted file mode 100644
index eaf5af8..0000000
--- a/.idea/codeStyles/Project.xml
+++ /dev/null
@@ -1,19 +0,0 @@
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
\ No newline at end of file
diff --git a/.idea/misc.xml b/.idea/misc.xml
index 9f709d9..ac6945a 100644
--- a/.idea/misc.xml
+++ b/.idea/misc.xml
@@ -3,7 +3,12 @@
-
+
+
+
+
+
+
\ No newline at end of file
diff --git a/.idea/vcs.xml b/.idea/vcs.xml
index a000c97..8f3769a 100644
--- a/.idea/vcs.xml
+++ b/.idea/vcs.xml
@@ -4,8 +4,5 @@
-
-
-
\ No newline at end of file
diff --git a/raytracer/media/scenes/grid.blend1 b/raytracer/media/scenes/grid.blend1
new file mode 100644
index 0000000..131623d
Binary files /dev/null and b/raytracer/media/scenes/grid.blend1 differ
diff --git a/raytracer/media/scenes/sportscar.blend1 b/raytracer/media/scenes/sportscar.blend1
new file mode 100644
index 0000000..280edbe
Binary files /dev/null and b/raytracer/media/scenes/sportscar.blend1 differ
diff --git a/utk_experiments/main.cpp b/utk_experiments/main.cpp
index 299bdc3..328849b 100644
--- a/utk_experiments/main.cpp
+++ b/utk_experiments/main.cpp
@@ -1,4 +1,5 @@
#include "utk/metrics/PCF.hpp"
+#include "utk/metrics/RDF_Heck.hpp"
#include
#include
#include
@@ -57,6 +58,28 @@ void writeRadspecToFile(const std::string &filename,
}
}
+template
+void writeRDFtoFile(const std::string &filename,
+ std::pair, std::vector> radspec) {
+ std::ofstream file;
+ file.open(filename);
+
+ auto xs = radspec.first;
+ auto ys = radspec.second;
+
+ if (xs.size() != ys.size()) {
+ std::cerr << "Dimensions of radial spactrum are unequal: xDim: "
+ << xs.size() << " yDim: " << ys.size() << std::endl;
+ std::terminate();
+ }
+
+ for (int i = 0; i < xs.size(); ++i) {
+ file << std::setprecision(std::numeric_limits::digits10 + 2)
+ << std::fixed;
+ file << xs[i] << ", " << ys[i] << std::endl;
+ }
+}
+
int main() {
// load Points from file and generate metrics
@@ -75,10 +98,10 @@ int main() {
for (float i = 0.015625f; i <= 0.5f; i += 0.015625f) {
// PCF
- std::cout << "Calculating PCF with " << NBINS << " bins for max radius " << i << "…" << std::endl;
- auto pcf = utk::PCF{true, 0.01, i, NBINS, 0.001}.compute(pointView);
+ std::cout << "Calculating RDF with " << NBINS << " bins for max radius " << i << "…" << std::endl;
+ auto rdf = utk::RDF{true, NBINS, i}.compute(pointView);
- writePCFToFile("pcffiles/pcf_" + std::to_string(i) + ".txt", pcf);
+ writeRDFtoFile("rdffiles/rdf_" + std::to_string(i) + ".txt", rdf);
}
// FFT based spectrum
diff --git a/utk_experiments/plotting/plotter.ipynb b/utk_experiments/plotting/plotter.ipynb
index 5a9d541..52ac494 100644
--- a/utk_experiments/plotting/plotter.ipynb
+++ b/utk_experiments/plotting/plotter.ipynb
@@ -2,7 +2,7 @@
"cells": [
{
"cell_type": "code",
- "execution_count": 1,
+ "execution_count": 6,
"id": "d0850083-6fec-4a51-8128-219c643397ce",
"metadata": {},
"outputs": [],
@@ -86,7 +86,7 @@
},
{
"cell_type": "code",
- "execution_count": 5,
+ "execution_count": 3,
"id": "4bc08371-8322-4f1a-b78a-8a0f94dd8a7f",
"metadata": {},
"outputs": [],
@@ -261,6 +261,30 @@
"source": [
"plt.close()"
]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "id": "f57e77eb-695c-4c98-9a80-281011682727",
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjEAAAGgCAYAAABbvTaPAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAABULUlEQVR4nO3deVxU9f4/8NcwMwyLzMg2LIqIiiiipmiAlbuopVR21bIv6b2mtmmk3tJu/a62aMstbzfLa1ZSptmmbRouN9cQFxQ33EUFZdiEGdZhO78/yAPDDMsAAzPwej4e83jM+Zz3OfM5p8l58zmfRSIIggAiIiIiG2PX1hUgIiIiagomMURERGSTmMQQERGRTWISQ0RERDaJSQwRERHZJCYxREREZJOYxBAREZFNYhJDRERENolJDBEREdkkJjFERERkk8xKYtasWYMBAwZAqVRCqVQiIiICv/32m7h/1qxZkEgkBq/w8HCDc+j1esyfPx8eHh5wdnZGVFQU0tLSDGJyc3MRHR0NlUoFlUqF6Oho5OXlNf0qiYiIqN2RmLN20i+//AKpVIpevXoBAL744gu8++67OHHiBPr164dZs2YhIyMD69evF4+xt7eHm5ubuP3000/jl19+QWxsLNzd3bFo0SLcvn0biYmJkEqlAICJEyciLS0Nn3zyCQBg7ty56N69O3755ZdGX1hlZSVu3boFFxcXSCSSRh9HREREbUcQBOTn58PX1xd2dg20tQjN5OrqKnz66aeCIAjCzJkzhQcffLDO2Ly8PEEulwubN28Wy27evCnY2dkJcXFxgiAIQnJysgBASEhIEGMOHTokABDOnz/f6HqlpqYKAPjiiy+++OKLLxt8paamNvhbL0MTVVRU4LvvvkNhYSEiIiLE8r1790KtVqNz584YMWIE3nzzTajVagBAYmIiysrKEBkZKcb7+voiJCQE8fHxGD9+PA4dOgSVSoWwsDAxJjw8HCqVCvHx8QgKCjJZH71eD71eL24LfzYwpaamQqlUNvUyiYiIqBXpdDr4+fnBxcWlwVizk5jTp08jIiICJSUl6NSpE7Zu3Yrg4GAAVY+Bpk6dCn9/f6SkpODVV1/F6NGjkZiYCIVCAY1GA3t7e7i6uhqc08vLCxqNBgCg0WjEpKcmtVotxpiycuVKLF++3Kj8Tv8dIiIish2N6QpidhITFBSEpKQk5OXl4YcffsDMmTOxb98+BAcHY/r06WJcSEgIhgwZAn9/f2zbtg1Tpkyp85yCIBhU1lTFa8fUtnTpUixcuFDcvpPJERERUftkdhJjb28vduwdMmQIjh49ig8++ABr1641ivXx8YG/vz8uXboEAPD29kZpaSlyc3MNWmMyMzMxbNgwMSYjI8PoXFlZWfDy8qqzXgqFAgqFwtzLISIiIhvV7HliBEEw6ItSU05ODlJTU+Hj4wMACA0NhVwux65du8SY9PR0nDlzRkxiIiIioNVqceTIETHm8OHD0Gq1YgwRERGRWS0xL7/8MiZOnAg/Pz/k5+dj8+bN2Lt3L+Li4lBQUIBly5bhkUcegY+PD65du4aXX34ZHh4eePjhhwEAKpUKs2fPxqJFi+Du7g43NzcsXrwY/fv3x9ixYwEAffv2xYQJEzBnzhyxdWfu3LmYNGlSnZ16iYiIqOMxK4nJyMhAdHQ00tPToVKpMGDAAMTFxWHcuHEoLi7G6dOn8eWXXyIvLw8+Pj4YNWoUvvnmG4MexqtWrYJMJsO0adNQXFyMMWPGIDY2VpwjBgA2btyIBQsWiKOYoqKisHr16ha6ZCIiImoPzJrszpbodDqoVCpotVqOTiIiIrIR5vx+c+0kIiIisklMYoiIiMgmMYkhIiIim8QkhoiIiGwSkxgiIiKySUxiiIiIyCYxiSEiIiJRcWkF1u67gqtZBW1dlQYxiSEiIiLRqt0XsfK38xj7/r62rkqDmMQQERGR6HDKbQBAZY2pcK9lF2Lafw/h9/PGCzS3JSYxREREJJKYKHv912QcuXYbf4s9Bo22BD8l3UR5RWWr1602s9ZOIiIiovZLEAQkpeYZlReVVojvJ36wH7lFZThxIw/Lovq1Yu2MsSWGiIiIAAD7LmaZLPdwUYjvc4vKAACx8ddw5qa2VepVF7bEEBEREQDg9/OZBtvT/nsIbs72iDurMRn/w/E0hHRRtUbVTGISQ0RERAAAO4lhj5gj127XG19S1rb9Yvg4iYiIiAAAMjtT3XrrVlJW0XCQBTGJISIi6sBe+yUZD3/8B0rKKiA1M4kpLm3bJIaPk4iIiDqYy5kF+PTAVYwMUuPzP1IAANtPp5udxBS1cUsMkxgiIqIOpEBfLs7Gu/loqliu0ZWY/Tipn6+yRetmLj5OIiIi6kA02mKT5e/EXWiwI29tL03o0xJVajImMURERB1IUT39WBKumpfEtDUmMURERB1I3p+T1TXFvb08xPcT+nm3RHWahUkMERFRB5JbVNrkY12d7cX3a/5vcEtUp1nYsZeIiKgDKC6tgEJm16yWmHt6usPNSY4+PkpIJOZ1ArYEJjFERETtXLq2GGPf24deXi6I6OFutP+BAT7YdirdoGxxZG/Y2UlwMjUPO85mAAAEAMsfDGmNKjcKkxgiIiIbJwhCvS0jH++5gsLSCpxMzcNJE6tU9/LshAf6+2Db6epEZt6InpBLq3qddF+yDQDgpVQYHduW2CeGiIjIhuUU6BGx8ne88WuyQfn20+n4x9bTKKuoxOXMgnrP4eZsj/88Nkjc/us93cUEBgA+fWIInh8TiFFB6patfDOxJYaIiMiGff5HCjS6Enx6MAWvTAoWy5/ZeBwAEODhjMtZDScxNWfrLa8QDPaPDfbC2GCvFqx1y2ASQ0REZMMK9YbzvhSXViA5XStuv7HtXIPncK8x6ggAyivbdnXqxmISQ0REZKN0JWXYcjzNoGzJllP4KemWWedx62SYxJTVaomxVuwTQ0REZKOe23QCupJygzJzExgAcHOq1RJTYRstMUxiiIiIrNi5dB2+PZYKQTBuHdl/Mctgu8SMVaVfeaCv+N7V6HGSbbTE8HESERGRFZv4wQEAQGdHOSIbmOo/K1/f6PPWTHjujEQa1K0zTtzIw9Qhfk2oaetjSwwREZGVOXtLi1nrj+DsreoOuufS8w1iUm8XGR03efXBOs8ps5Ogq6ujuO0glxrFfD0nHP9bNAIjens2pdqtji0xREREVuaxTxKgKynH8eu5YplMajiZ3bObjhsdV3tJgT7eLjivqUp+yisFRA30xcd7r6CXuhMeu7sbfj+fifE1Wncc5FL09OzUkpdiUUxiiIiILOzw1Rx4KR3Q3cO5UfF3OuvW7rRbc2beU2lao+NqWhsdir7eSszdcAznNfkY21eN58cGoodnJwwP9ICzQoZNc8KbcDXWQyKY6inUDuh0OqhUKmi1WiiVyrauDhERdVCXMvIxbtV+AMC1tx5o1DF3pvmvzV5qh3enDsDnB1NwsoEk5s5npWuLseX4TTx2dze41erAa43M+f1mSwwREZEFJafrWuxcpRWVeH5zklnH+Kgc8eyoXi1WB2vCjr1EREQWVN/CjKboyxs/TLqjYxJDRERkQfWlMOc1OhTqq/u9nEvXof+ynWadf+OTYUZlrz/Yz6xz2Cqzkpg1a9ZgwIABUCqVUCqViIiIwG+//SbuFwQBy5Ytg6+vLxwdHTFy5EicPXvW4Bx6vR7z58+Hh4cHnJ2dERUVhbQ0wymTc3NzER0dDZVKBZVKhejoaOTl5TX9KomIiNqIXY2WmIpKAYIg4NUfz2DGugRM+PcBPPjRH+L+5zefQGl542fLldlJcE8vD3H7vkAPHHtlLKIjurdI3a2dWUlM165d8dZbb+HYsWM4duwYRo8ejQcffFBMVN555x28//77WL16NY4ePQpvb2+MGzcO+fnVY9tjYmKwdetWbN68GQcPHkRBQQEmTZqEiorq5rMZM2YgKSkJcXFxiIuLQ1JSEqKjo1vokomIiFpPzadJpeWVOHtLhw0J1xF/JQcAcDmzaoXpA5eycDGj/tWma7OzM2zn6eujhEcnRfMqbEPM6tg7efJkg+0333wTa9asQUJCAoKDg/Hvf/8b//jHPzBlyhQAwBdffAEvLy9s2rQJ8+bNg1arxWeffYYNGzZg7NixAICvvvoKfn5+2L17N8aPH49z584hLi4OCQkJCAuraiJbt24dIiIicOHCBQQFBbXEdRMREbWKmmlGYWk5dCVlJuOW/HDa7HM/MrgLAOC1B/vh56Rb7bYDb12a3CemoqICmzdvRmFhISIiIpCSkgKNRoPIyEgxRqFQYMSIEYiPjwcAJCYmoqyszCDG19cXISEhYsyhQ4egUqnEBAYAwsPDoVKpxBhT9Ho9dDqdwYuIiMiaDHljN3YnZ5rcpzfjMRIAfDRjMJZMqFr/6ImI7vj+6WFQOcqbXUdbYnYSc/r0aXTq1AkKhQJPPfUUtm7diuDgYGg0GgCAl5eXQbyXl5e4T6PRwN7eHq6urvXGqNVqo89Vq9VijCkrV64U+9CoVCr4+dnGug9ERNS+ldVaTPHzP1KMYsJX/A/ZBY1f96i7uxMeGOADlVPHSlpqMzuJCQoKQlJSEhISEvD0009j5syZSE5OFvfXHkpWc3bButSOMRXf0HmWLl0KrVYrvlJTUxt7SURERBYhCAIWfH2iwTiNrsSs8/qoHBsO6gDMTmLs7e3Rq1cvDBkyBCtXrsTAgQPxwQcfwNu7au2F2q0lmZmZYuuMt7c3SktLkZubW29MRkaG0edmZWUZtfLUpFAoxFFTd15ERERt6U6n3ZZWVmHeo6f2qtnzxAiCAL1ej4CAAHh7e2PXrl3ivtLSUuzbtw/Dhg0DAISGhkIulxvEpKen48yZM2JMREQEtFotjhw5IsYcPnwYWq1WjCEiIrIF2mLTnXibq7CUE+IBZo5OevnllzFx4kT4+fkhPz8fmzdvxt69exEXFweJRIKYmBisWLECgYGBCAwMxIoVK+Dk5IQZM2YAAFQqFWbPno1FixbB3d0dbm5uWLx4Mfr37y+OVurbty8mTJiAOXPmYO3atQCAuXPnYtKkSRyZREREVq+yUkBuUSncOynqHInUVLPvDcBnB1OwZGKfFj2vrTIricnIyEB0dDTS09OhUqkwYMAAxMXFYdy4cQCAF198EcXFxXjmmWeQm5uLsLAw7Ny5Ey4uLuI5Vq1aBZlMhmnTpqG4uBhjxoxBbGwspFKpGLNx40YsWLBAHMUUFRWF1atXt8T1EhERNdvFjHy89MMpPD8mECODDAejvPLTGWw6fAOPDvWDWulg9rm7dHbEzbxicfuxu/2QfEuHv90bgMkDfDHnvh7wVpl/3vaIq1gTERGZKXLVPnFiuporUwuCgICl25t17uOvjsPx67l48stjAIC3H+mP6UO7NeuctsSc32+unURERGSmW3mmRxOZO8qotokh3nBztkdf3+of72AfVbPO2Z6Z9TiJiIiIgPJK06ODbtV4DGSuJ+8NQMy43gAAZ/vqLhaBXp2afM72jkkMERGRmSoqjXtilJZX4pE1h5p0vvsCPfDKpGBxu7OTPdY9MQTOCikc5NJ6juzYmMQQERGZqbxGEhN3RgN9eQWkdvVP7Fqf8B7uRmXjguueG42qMIkhIiKqR4auBK/9moyZEd1xd4AbAKDmkJinvkps9mfYNTCzPZnGJIaIiKgeS344hT0XsrDtVDoGdlVZZHhzsC9H0TYFkxgiIqJ6XMspEt+fTNPiZJq2Rc//r6kDMTzQo0XP2VFwiDUREVE9WupJz6I/Rx7VFODhjL+Edm1woWQyjUkMERFRHW4XluJqVmGjYp8d1bPOfZueDMPE/t5G5UxdmodJDBERUR0e/Ohgo2MHd3PF/r+Pgr+7k9E+Fwc5FDLjodLREf7Nql9Hxz4xREREtVRUCkjJLkTq7cZPXueldEA3dye8+kCwuGRATQp5dbvBJ9GhcO+kwCC/zi1R3Q6LSQwREVEtb/12DusOpNQbYycBVI5y5BZVrVStVioAGCYrd3ipFAaT1jnIpQj1d23BGndMfJxEREQdXn5JmcF2QwkMAMikdnh6ZHU/GHfnqiTGo5NCLPvn5GB8Oy8CahcHONR4nFTRPtdebnVsiSEiog5tQ8J1vPrjGbz7lwGYOsSv8QcKQC919bpGd2bs7ePtgo8fHwxtcRkeHeonjjySS6u78VZUMIlpCUxiiIioQ3v1xzMAgL9/fwof/n4Z0eGN62xbKQgYFaTG4sjeBpPVSSQS3N/fxyi+5jBqtsS0DCYxRETUIVVWCrCrtd7RjdtFeHP7uUYdL6AqMXludKDZn92ls6PZx5Ax9okhIqIOZ+2+Kxj42k6cS9eZddx/Hhskvhea0Jry3VMReH/aQIR0UZl9LBljSwwREXU4K387DwCY+MEBs47r7CgX3zflgdDQ7m4Y2t2tCUeSKWyJISIiaiR/dyfc/WcScn+Icb8Xal1siSEionajUF8OZ4XhT1tJWQWkdhLIpVV/t5eWVzbp3BP6ecPf3Rn/jQ7FkZTbGNbLvdn1peZhEkNERO3Ct0dT8eIPp/CvqQPxl9CuAABtcRlG/2sv/N2dsOWZewAAr/+a3Kjz3RfogWlD/JChK0HC1RyxP4ybsz0mhBivg0Stj0kMERG1Cy/+cAoAsPi7k2IS8/PJW8gpLEVOYSmuZReiq6sjNiRcb9T5XhjXG4O7Vc2q++R9PSxTaWoW9okhIqJ2K/mWVnw/8l97cS2ncStSA4C/m/FCjmRdmMQQEVG7k6ErwROfH8EvJ9MNyse+v198P2943a0ruxcOh3uN5QPIOjGJISIim7U7OQNPfH4EGboSg/J3d1zA/otZKNCXmzyur48Skf286jxvL7VLi9aTLIN9YoiIyGY9+eUxAMDyX86KZfZSO+iKy+o6BADg0ckeQd5Kk/vm1tNCQ9aFSQwREdm87ac14nu5VILOTvJ6ogF3Z3t0UsgQ3sMNCVdvi+V7Fo9EgIezxepJLYuPk4iIqF2xl9lB6dBAEvNnf5cNs8Ow6ckwsdxRLrVo3ahlMYkhIiKrVqAvx4ZD15CpK8GVrAJsOnwDpeWVSLyeazLeTiJBUVlFved072QPAJBL7aCssZSAQsafRVvCx0lERGRV0rXF8FY6QCKpWmH6//14BltO3MS3x9Jw9pYWlQLwdtx5aOvo95JTWIpNh2/U+xkezqZHHtkzibEp/K9FRERW46ekm4hY+Tte+fGMWLblxE0AwOmbVQkMgDoTmMa60xJTG1tibAv/axERkdV4d8cFAMDGBlpSmqJmZ19VjUdIapfqVhmZlD+LtoSPk4iIyGrY/fkICQAEQcD+S9ktdm7PTgrkFVW14HjUmMhOrXTAhtl3Gy0cSdaP/8WIiMhq1MhhsO9iFmatP9pi547o6Y6X7++LrAI9utcaRn1foGeLfQ61HiYxRERkNWq2xMRfyWmRczrKpXh2VE/89Z4Atra0M/yvSUREbUYQBCz89iTsJBK8N20gajTEIO6Mps7jGuPBu3zx9iMDIJfaQWonafgAsjnswURERG3mdmEptp64iR+Op2H/xSxcza5eZfrG7aJGneOFsb2NyqaGdsUHjw6Cg1zKBKYdY0sMERG1uopKAYu/O4nfzlSvMv2PH0+bfR5XJzlC/V3F7adG9MT1nEL8fXxQi9STrBuTGCIialFFpeV4+7fzmBDig4ie7gb7yisqcV6Tjw/+dwm7kjMM9qXeLjbrc54Z2RPPje5lcFxYgBuWTOzT9MqTTTHrcdLKlSsxdOhQuLi4QK1W46GHHsKFCxcMYmbNmgWJRGLwCg8PN4jR6/WYP38+PDw84OzsjKioKKSlpRnE5ObmIjo6GiqVCiqVCtHR0cjLy2vaVRIRUav5aM9lfHHoOh5bl2C0b+Vv5zHpw4NGCUxTuDnbw8leBjfn6onrZFI+OupIzEpi9u3bh2effRYJCQnYtWsXysvLERkZicLCQoO4CRMmID09XXxt377dYH9MTAy2bt2KzZs34+DBgygoKMCkSZNQUVG91sWMGTOQlJSEuLg4xMXFISkpCdHR0c24VCIiag2XMgrq3PfZwZQW+5w7I5lca0xiV1Ra/5pJ1L6Y9TgpLi7OYHv9+vVQq9VITEzE8OHDxXKFQgFvb2+T59Bqtfjss8+wYcMGjB07FgDw1Vdfwc/PD7t378b48eNx7tw5xMXFISEhAWFhVauLrlu3DhEREbhw4QKCgvisk4jIWlXcWRugldScZVdo3Y+mNtas0UlarRYA4ObmZlC+d+9eqNVq9O7dG3PmzEFmZqa4LzExEWVlZYiMjBTLfH19ERISgvj4eADAoUOHoFKpxAQGAMLDw6FSqcSY2vR6PXQ6ncGLiIhaX7kFkpj7Aj3Q10eJT58YIpbV/JTnxwTinl7uGN1H3eKfTdaryR17BUHAwoULce+99yIkJEQsnzhxIqZOnQp/f3+kpKTg1VdfxejRo5GYmAiFQgGNRgN7e3u4uroanM/LywsaTdWcABqNBmq18RdRrVaLMbWtXLkSy5cvb+rlEBFRCymvrGz2OdZGh2LehkQAwNi+anw6c2i98S+MMx5mTe1fk5OY5557DqdOncLBgwcNyqdPny6+DwkJwZAhQ+Dv749t27ZhypQpdZ5PEARx2XUABu/riqlp6dKlWLhwobit0+ng5+fX6OshIqKWUVZR3UaSfEuHIG8Xca4WuVRisN+U8f28MKZGi4qni4NlKko2r0lJzPz58/Hzzz9j//796Nq1a72xPj4+8Pf3x6VLlwAA3t7eKC0tRW5urkFrTGZmJoYNGybGZGQY91zPysqCl5eXyc9RKBRQKBQm9xERUfOVV1TinR0XMKynO0YG1f3YpryiuiXm/v8cwCODu+LsLS18VA4NJjAAsObxUNjZSbDuiSH49lhqnXO+COwA0+GZ1SdGEAQ899xz2LJlC37//XcEBAQ0eExOTg5SU1Ph4+MDAAgNDYVcLseuXbvEmPT0dJw5c0ZMYiIiIqDVanHkyBEx5vDhw9BqtWIMERG1rh+Op+GT/VcbXJSxdsfeH46n4bwmH3suZDXqc+z+bLUZF+yFdU8MMRhCDQB9vF0AAOP7mR5AQh2HWS0xzz77LDZt2oSffvoJLi4uYv8UlUoFR0dHFBQUYNmyZXjkkUfg4+ODa9eu4eWXX4aHhwcefvhhMXb27NlYtGgR3N3d4ebmhsWLF6N///7iaKW+fftiwoQJmDNnDtauXQsAmDt3LiZNmsSRSUREbeRWXkmj4prSsffZUT3x0Z4rjYr9+bl7kV9SBvdObH3v6MxKYtasWQMAGDlypEH5+vXrMWvWLEilUpw+fRpffvkl8vLy4OPjg1GjRuGbb76Bi4uLGL9q1SrIZDJMmzYNxcXFGDNmDGJjYyGVSsWYjRs3YsGCBeIopqioKKxevbqp10lERM0kb8REcvryCpy9Zf7oUKld4x8M2MvsmMAQADOTmIaePzo6OmLHjh0NnsfBwQEffvghPvzwwzpj3Nzc8NVXX5lTPSIisiC5wXwshgMtCvXl+O2MBo5yqalDAVQtE/DxXtOtLQoZ1yMm8/FbQ0REjVIziSkpMxxG/ca2ZCz+7iSe3XTc5LEenewN5nBZOrEPgn2U4vbjYd3Q1dURT97bcF9Loju4ACQRETXKnWHSAFCgL8fV7AJIIMEn+6/gx6Rb9R47PNATXsrqodJRd/lix9nqeb86O9njwIuj6pxGg8gUJjFERNQoZTWGTqdrixG1+o9GH/t4uD88Xar7sTjIpHjz4f6I/uwInh8bCMD0/GBE9WESQ0RE9Tp8NQd///4UvGu0pFzLKWr08V/+7W6E+lfNC/bBo3ehpKwCrs72cHW2x9F/jGHyQk3GJIaIiIzsSs7Auv1X8d60gXjyi2PI15fjxu3qxCW3sLTR5+rjUz069cG7uhjsYwJDzcEkhoiIjMz58hgA4P/9dAZFZRVG+29pixt1nqiBvlBz2QCyEI5OIiKiOuUUlpocNn2jgcdJo4I84efmiDcfDqk3jqg52BJDRER1ktlJ4CCXokBfblB+QZNf73GfPDEEMjsJHxeRRTGJISKiOsmkdnC0N260v5pdWOcxQ7u7GswpQ2QpTGKIiMhAfkmZ+F4uldQ7C29t51+fAJkdW1+odTCJISIiA6//miy+l9nZQdbIdY2WR/WDgxkJD1Fzsb2PiIgMfHssTXwvtZOgopGrUs8I62apKhGZxCSGiKid+/18BhZ9e9Kgc25peSXW/5GCy5kF+CrhOoJe+Q2PfnIIh67koJubkxiXlJqHCxn1d+Lt4+2CQ0tHsx8MtTo+TiIiaocK9OVwtpdCIpHgb7FVc774uzthwZiqKf6/iL+GN7efMzgm4eptJFxNgL97dRJz28SkdjtfGI5PD1zFt8fSsDyqH2YO6265CyGqB9NmIqJ2JvH6bYT8cweW/5JsUJ5doBffH7+RW+fx1xuYA6aHhzPemjIAexaPxBMR/s2rLFEzMIkhImpn3t1xAQAQG38NglDdn8VZUd343pwOuDKpHezsJAjwcOY8MNSmmMQQEbUzNfumaIurh0s71UhcHOT1//Mf7KPEY3ezoy5ZN/aJISJqZ6Q15mnJyq9+hFReKaCotBz7LmQBqL8FZaBfZ3R2kluqikQtgkkMEVE7U3Oyuawa/WCKSssxY91hJKXmNXiOzk5yFJRUj2ZylEtRXFaBgV1VLVpXoubg4yQiIhu053wmxr2/DydNJCQ1J6er2RJzu7DMZAITNdDXqKyTQmbQh+abeeGYPNAXq2cMbl7FiVoQkxgiIhv019ijuJRZgOe+Pm5QXlkpIO6sRtyumcT8cuqW0XneeWQA/vPYIKNyFwcZ5g7vgWE93fHe1IEY0LUzPnxsEPxqzCFD1Nb4OImIyMYUl1aYfA8AW0/cNNi+lVcivi8trzQ6l4O96VFKLg4yuDnbY9Oc8OZUlcii2BJDRGRDDl/NQd//FyduO9ZKQv64km2wfbGB2Xa7ujqaLHeU829csn78lhIRWaGSsgrY/zkfS03//PmswbbTn8mGIAiY//UJ/Hoq3WD/wcuGSU1NSgcZBvl1NrmPC1GTLWBLDBGRldEWl2Hg8p2Y8WlCg7FOiqqWmHRtiVEC05Cou3zFyerCe7gZ7OMkdmQLmMQQEVmZvRcyoS+vRMLV22LZhkPX8NvpdAi1FpQ+cSMPy34+i5t5xXWez8/NEQtG9wIA/F94N6yc0h89PJ0xb3hPMeaLv92N3xeNgEJW9bMw0I9Dqcn68XESEZGVqayVqVzPKcSrP1U9RgpUdzKKj42/htj4a3Web4i/G14Y1xvDenmgn68SLg5yo9l4FTIpenh2wpF/jEWBvhxqF4fmXwiRhTGJISKyMrVbW2ouHXAps8Ds843uo4ZEIkF4D/cGY1WOcqgcOVMv2QYmMUREVqZmElNRKaCo1jDq+kjtJKiorDqBo1yKWfd0x6QBPi1dRSKrwCSGiMjK1GyIKS2vhK5GS0xD7iQwAPCvqQPxABMYaseYxBARWZmafWJKyyuhq7GGUUNCuijx+oMhOH4jDxNDvC1RPSKrwSSGiMjK1GxNeXnraZgz2vmjGYPh7+6MQd1cLVAzIuvCJIaIqI39cvIWPtl/Fc+M7IkRQZ7Ql1X3gdl2uvFzv7g6yeHv7myJKhJZJSYxREStSKMtwe/nM/HwoC7ikgHzvz4BAHh643GE93DDqCC1Wef8+bl7sO5ACv52T/eWri6RVWMSQ0TUiqaujUfq7WJcyszHPyf3M9qfcPW2wSR3Denu7iSuME3U0XDGXiKiVpR6u2pm3Z1nM8w6rvZCjf8X3g1J/28cdrwwvMXqRmRrmMQQEbWB0opKs+JrTkB3X6AH3nioPzo72UMhk9ZzFFH7xsdJRESt4P2dF5CUphW3s/L1+Pfui5g2xK9Rx3d2qk5ipFximggAkxgiIosqKatASnYh/vP7ZaN9a/Zewb93X2rUeZQO1UmMzI6N6ESAmY+TVq5ciaFDh8LFxQVqtRoPPfQQLly4YBAjCAKWLVsGX19fODo6YuTIkTh79qxBjF6vx/z58+Hh4QFnZ2dERUUhLS3NICY3NxfR0dFQqVRQqVSIjo5GXl5e066SiKiNzPz8CCZ+cMDkPn154x8pOSuq/+aUsSWGCICZScy+ffvw7LPPIiEhAbt27UJ5eTkiIyNRWFgoxrzzzjt4//33sXr1ahw9ehTe3t4YN24c8vPzxZiYmBhs3boVmzdvxsGDB1FQUIBJkyahoqJ6boQZM2YgKSkJcXFxiIuLQ1JSEqKjo1vgkomIWs/hlMaPNLrD2d64n0uXztUde2VSJjFEACARhNrrpTZeVlYW1Go19u3bh+HDh0MQBPj6+iImJgYvvfQSgKpWFy8vL7z99tuYN28etFotPD09sWHDBkyfPh0AcOvWLfj5+WH79u0YP348zp07h+DgYCQkJCAsLAwAkJCQgIiICJw/fx5BQUEN1k2n00GlUkGr1UKpVDb1EomImqX7km1mH+PRSYGX7++D2PhrOPVnP5qPHx+Mk6l5WB9/DdsX3Ide6k4tXVUiq2DO73ezHqxqtVX/c7m5uQEAUlJSoNFoEBkZKcYoFAqMGDEC8fHxAIDExESUlZUZxPj6+iIkJESMOXToEFQqlZjAAEB4eDhUKpUYU5ter4dOpzN4ERHZIrlUgimDu+Ln5+7FKw/0xV9Cu2J8P28svb8vTi+LZAJD9Kcmd+wVBAELFy7Evffei5CQEACARqMBAHh5eRnEenl54fr162KMvb09XF1djWLuHK/RaKBWG89YqVarxZjaVq5cieXLlzf1coiIWlxTG7rl0uq/L5+8r4fBPg6pJqrW5JaY5557DqdOncLXX39ttE9Sa7UyQRCMymqrHWMqvr7zLF26FFqtVnylpqY25jKIiCxm/R/XmnQc+7wQNU6Tkpj58+fj559/xp49e9C1a1ex3Nu7atn32q0lmZmZYuuMt7c3SktLkZubW29MRobxbJZZWVlGrTx3KBQKKJVKgxcRUWsSBAGfH0zB0WtVnXlf+zW5Ucedf30C/lpj3SOOPiJqHLOSGEEQ8Nxzz2HLli34/fffERAQYLA/ICAA3t7e2LVrl1hWWlqKffv2YdiwYQCA0NBQyOVyg5j09HScOXNGjImIiIBWq8WRI0fEmMOHD0Or1YoxRETWZt/FLLz2azKm/vcQfkhMg6O84Uc/9/f3hoNcarCOEleiJmocs/rEPPvss9i0aRN++uknuLi4iC0uKpUKjo6OkEgkiImJwYoVKxAYGIjAwECsWLECTk5OmDFjhhg7e/ZsLFq0CO7u7nBzc8PixYvRv39/jB07FgDQt29fTJgwAXPmzMHatWsBAHPnzsWkSZMaNTKJiKgtpN4uEt8v+u5kvbF/vac7Lmbk482H+otlkwf64vDVHCyPMl4YkoiMmZXErFmzBgAwcuRIg/L169dj1qxZAIAXX3wRxcXFeOaZZ5Cbm4uwsDDs3LkTLi4uYvyqVasgk8kwbdo0FBcXY8yYMYiNjYVUWv1Xy8aNG7FgwQJxFFNUVBRWr17dlGskImoV5kxe9+yoXvDopDAo+/CxQSirqDTo2EtEdWvWPDHWjPPEEFFr2J2cgUNXc1BcVoFNh280+riT/4w0WNSRiKqY8/vNtZOIiOohCAIqBdOLLm45noaF39b/2Kgu9mxtIWo2/l9ERFSPmeuPYtz7+1Bq4lHRl4euN/m89jL+80vUXGyJISKqQ0WlgP0XswAA8VeysWL7Obg522POfT3w0Z7LSErNMzpmw+y7Ef3ZEaPy2ky17BCReZjEEBHV4cClLPH9rPVHxfcJV+te1LF/F5VR2fvTBsJeZofnNp0AAIQFuLVgLYk6LiYxRES1aIvK8MeVbDyz8bjZx5rqrDtlcFdk6krE7X9NHdis+hFRFSYxREQ1bEi4jld/PNPk4+taGqWzk7343qERk+ARUcOYxBAR1dCcBOaO+aN74cPfLxuU2cvs8O/pd6GotAKeLoo6jiQic7B7PBF1aPsvZmH62kO4mlXQ6GNmRviL753spVA6GP49uCgyCP95bBAAYOWU6hl5HxrUBTPCujWzxkR0B1tiiKhD0pdX4OvDN7Dsl6pFGhd/dxJbnrmnUcfWfGSkkNmhvMJ4ztCogb4Y00cNZwX/mSWyFP7fRUQd0t+/O4WfT94St4/fyMO/dlxo1LFyaXUS4+IgR06B3mQcExgiy+LjJCLqEH49dQufHrgqbtdMYO5YveeyUZkpT97XAyun9IfSQYZV0+9CZbtcvIXI+vHPBCLqEO7M0dLV1QmVzVgyblhPd3gpHfDY3d3w6FA/SCSSZp2PiJqOSQwRtXslZRXi+6e+SjTr2HnDe2Dt/qsm993pG6N0lCMr3/QjJSKyHD5OIqJ2T1dc1uRjF48Pwjdzw+uNWffEEPT0dMbns4Y0+XOIyHxsiSGidkcQBIMRRNomJDFj+6oxrKcH5FI7hPVwr3Fu49i7/Drjf4tGNqWqRNQMbIkhonblUkY+Br++C+tqPAKqK4mpY3JdAMCnM4fib/cGGJULYP8XImvBJIaI2pUlW04jt6gMb24/h51nNQCAvCLTScyArp2x6ckwcbvm0Om6sA8vkfVgEkNE7cb1nEIkXs8Vt+duSMR7Oy/gyS+PmYx3dZIjvIc7Hg/rhrcf6Q9/d+c6zx0dXjVL76LIoJatNBE1GfvEEFG7ERt/zais9hpGNbk62cPOToI3H65aGmBQN1e88uMZxIwNNIp97cF++PuEICgdjFepJqK2wSSGiNoNcx/1dHYyTEh6e7ng23kRJmMlEgkTGCIrwySGiGzekh9OITNfD9/ODmYd5+pkb6EaEVFrYBJDRDatpKwCm4+mAgCGdnc169jaLTFEZFuYxBCRTaqsFPDx3stQK6tbX45eyzWI8ehkD2eFDNdzikyeozNbYohsGpMYIrJJu85l4F87L9YbsyyqH+4L9MTHey/DTiLBF/HX8HhYN6w7kAIA6OFR92gkIrJ+TGKIyCY1Zq0iCSRQOcqxdGJfAMDiyCDYSYC7/FzhpVQgpIvK0tUkIgtiEkNENuPMTS0SruYgt6gUPT07NRhvV2vuOumfBQ8M8LFE9YiolTGJISKb8MflbDz+6WFx+7G7/Ro8JthXackqEVEb44y9RGQTfj2VbrC982xGvfEjenvWOwMvEdk+tsQQkVX7+eQt7DyrMUpicgpLxfcTQ7yxavpdUMjsELB0OwAg1N+84dZEZHuYxBCR1SqvqMSCr0/UGzN/dC+D9YwWjAnEr6duYWZEdwvXjojaGpMYIrJaRWUVDca4OBj+M7ZwXG8sHNfbUlUiIivCPjFEZHEXNPmY/OFB/H6+/n4sgiDg2Y3H8dL3pwAAuTUeGdWlk4Kz7hJ1VExiiMjintmYiNM3tfhb7LF649Jyi7HtdDq+OZaKAn05Rr+3r8Fz9/ZqeKg1EbVPTGKIyOIyGzExHVC1DtIdcWc0qKhseFlqTlhH1HGxTwwRWVx5RcPJCABoi8vE94u/O1ln3FtT+iMpNQ89PJ3hIJc2u35EZJuYxBCRxTXUohJ3Jh3r/7iGwym3G3W+qUP88Ojd3VqiakRkw5jEEJHFlVdW1rlvV3IGnvrqeKPP9fCgLuLyAUTUsbFPDBFZXH0NMV8fuWGyvJfadIfdSqFxj6aIqP1jEkNEraqkrAKCICD2jxT0enk7fj+faTJu98IRcJAb/xM1tLubpatIRDbC7CRm//79mDx5Mnx9fSGRSPDjjz8a7J81axYkEonBKzw83CBGr9dj/vz58PDwgLOzM6KiopCWlmYQk5ubi+joaKhUKqhUKkRHRyMvL8/sCySi1nc6TYv4K9km9w1/Zw/6vBqHZb8ko7yOJprXH+wHAOjt5SKW7V08Eu/8ZQAeY18YIvqT2UlMYWEhBg4ciNWrV9cZM2HCBKSnp4uv7du3G+yPiYnB1q1bsXnzZhw8eBAFBQWYNGkSKiqqh1fOmDEDSUlJiIuLQ1xcHJKSkhAdHW1udYmoDUxefRAz1h3GnC+P4aekmwb7MvP10Jcb9pGJDPbCL8/dC08XBZ4fE4j/C/cHAHz42CCM6aPGd09FoLuHM6YN8WN/GCISSQSh6Q+YJRIJtm7dioceekgsmzVrFvLy8oxaaO7QarXw9PTEhg0bMH36dADArVu34Ofnh+3bt2P8+PE4d+4cgoODkZCQgLCwMABAQkICIiIicP78eQQFBZk8d006nQ4qlQparRZKpbKpl0hEZiopq0CfV+PMOmbBmEAuFUBEAMz7/bZIn5i9e/dCrVajd+/emDNnDjIzq595JyYmoqysDJGRkWKZr68vQkJCEB8fDwA4dOgQVCqVmMAAQHh4OFQqlRhTm16vh06nM3gRUeurOddLY0X0cLdATYiovWvxIdYTJ07E1KlT4e/vj5SUFLz66qsYPXo0EhMToVAooNFoYG9vD1dXV4PjvLy8oNFoAAAajQZqtdro3Gq1WoypbeXKlVi+fHlLXw4R1UNXUoZ/7biAAn05+ndR4a/3BODfuy+ZdY6+PkoM6e7acCARUS0tnsTceUQEACEhIRgyZAj8/f2xbds2TJkypc7jBEGARFL9rLvm+7pialq6dCkWLlwobut0Ovj5+TXlEoiokVZuPy8Okd5y/CZ8VI51DpmuaffC4cgtKsPQ7m6orBRgx34uRNQEFp/szsfHB/7+/rh0qeqvM29vb5SWliI3N9egNSYzMxPDhg0TYzIyjFe7zcrKgpeXl8nPUSgUUCgUFrgCIqrL/otZBtu7z9W/SjUAzBvRA73U1aOOmMAQUVNZfJ6YnJwcpKamwsfHBwAQGhoKuVyOXbt2iTHp6ek4c+aMmMRERERAq9XiyJEjYszhw4eh1WrFGCJqfRm6Ejy76TiO38gFABToyw32f5+YZuowAw/d1cUidSOijsfslpiCggJcvnxZ3E5JSUFSUhLc3Nzg5uaGZcuW4ZFHHoGPjw+uXbuGl19+GR4eHnj44YcBACqVCrNnz8aiRYvg7u4ONzc3LF68GP3798fYsWMBAH379sWECRMwZ84crF27FgAwd+5cTJo0qVEjk4jIMj7ZfxXbTqVj26l0/H18kNmdeH1UDujpaXomXiIic5mdxBw7dgyjRo0St+/0Q5k5cybWrFmD06dP48svv0ReXh58fHwwatQofPPNN3BxqW4+XrVqFWQyGaZNm4bi4mKMGTMGsbGxkEqrV6PduHEjFixYII5iioqKqnduGiKyvOKy6rmc3t1xoc64u7u7oaS8Aj09O2Hriap5YlZNH4hJA3whl3KicCJqGc2aJ8aacZ4Youa5mlUA386OcJBX/3Hx0ven8M2x1AaP/eHpCIT6Vy0PcF6jQ9KNPEwf6ldnx3wiojvM+f3mKtZEZCT+SjZmrDuMYT3dER3uj8//SMHDg7riZFqeyXg7CeDv7oyU7EIAgI/KUdzXx1uJPt78Q4KIWh6TGCIyEvvHNQBA/JUcxF/JAQAcvZZrMnb6ED/8MyoYmw7fwBvbzgEA1C4cKUhElseH00RkRFfS+A67j4R2hZO9DGqlg1gmY78XImoFbIkhIiP5JeUNBwH4ek447g6o6vtyf4g3job7i9tERJbGJIaIAADxl7NxM68YU4f44XZhaaOOCfap7usik9rh9YdCLFU9IiIjTGKICAAw49PDAIC8ojKka0vqjHtrSn8s2XIaANDJgf+EEFHb4b9ARB3ct0dTEeDpLG6/uf1cnbGbngzDsF4euJVXjM5O9pByyQAiakNMYog6sMTruXjxh1ONjh/WywMAsDCSM2cTUdtjEkPUARXqy6ErKcO1P+d1acic+wLwl1CuCk9E1oVJDFEHNPnDg7iaXYhF43qb3D8zwh9fHLoOAFg/ayhG9VG3ZvWIiBqFSQxRB1BcWoFd5zIworcnVI5yXP2zBea9XReNYu8OcMNjYd3w8OCuOJeuw8ggz9auLhFRozCJIeoAVv52Dl8euo77Aj2w7okh9cZ+Oy9CfH+XX2cL14yIqOk4rSZRB/DdsTQAwIFL2Ui9XdTGtSEiahlMYog6AKVjdaPruFX727AmREQth0kMUQegcpS3dRWIiFockxiidq60vBKdHe2NypU1Ztt1czbeT0Rk7dixl6gduplXjL9/dxIp2YV1LiEwvp83DqfchouDDO/+ZSCmf3IIz47q1co1JSJqOiYxRO3Qx3suI/5KTr0xPp0dsXvhCACAvcwOJ/9fJOy4jAAR2RA+TiJqZzYevo6Nh2+Y3OdkLxXf39/fG/YyO9jLqv4ZYAJDRLaGLTFENubOEGk/NyeT+/+x9Uydx34SPQSZ+SUY2t2tzuOJiGwFkxgiG6Ivr8B97+wBAFx4YwIUMimKSsux/o9rcJRL8WPSzXqP7+erxL3OHq1RVSIii2MSQ2RDtMVl4vvcwjKc1+Tgf+cysSHheoPHDvTrDFeOQiKidoR9YohsSJG+Qny/7XQ6Zq0/WmcCs2xysMH2IC4hQETtDJMYIhtSoC8X3284dK3eWGmtjro1O/USEbUHTGKIrNy17EK8t/MCCvXlBknMtZy610BydZLDW+VoUOYoZxJDRO0L+8QQWTFBEDD2/X0orxQgkUgwoIuqwWPu7++NF8b2Rk/PTnh2VE98tOcKAODhwV0sXV0iolbFJIbISgmCgLTcYpRXCgCAhKs50BaVNnjcy/f3RVfXquHTfx/fBy+M7Y1CfQVUTlw/iYjaFyYxRFamolLAOzvOY+vxm7irRmfcIym3cSTldoPH114HSSa1g8qJT46JqP1hEkPUhkrKKnApowD+Hk747EAKxvRVY86Xx5Ch0wMAdiZnmH1O9n0hoo6CSQxRKzpzU4s95zMxZ3gPOMiliNmchLizGnH/B/+71KjzPHlvAD49mAIA6KXuhEkDfPDv3Zfgo3KARMLlA4ioY2ASQ9SKJn14EEDVOkUPD+pikMA0VvyS0fBROeD742nIKyrDq5OCMaK3Jx66qws6s98LEXUgTGKIWkFeUSm2nU4Xt9/dcQHv7rjQpHP5dq4aOr3zheG4lFGAYT3dAQDdPZybX1EiIhvCJIaoFSzYnIT9F7OafHwfbxec1+QblKldHKB2cWhu1YiIbBaHLBC1guYkMADwwaOD0KWzI14Y27uFakREZPuYxBBZ2NWsArPiFTI7/OexQQCAAV1V+OW5exHk7YI/lozG82MDLVFFIiKbxMdJRBaSri3GkZTbSLja8Nwu7/5lAP7+/SkAwLYF96GXuhMierjD3dkednYcbUREZAqTGCILefqr40hKzWtU7CODu0LpKEduYSl6qTsBADxdFBasHRGR7WMSQ2QhjU1ggKoh1+P7eVuuMkRE7RD7xBBZQOWf6x01xoCuDS/qSERExpjEEFnA+7suNiqur48SXz0ZZuHaEBG1T2YnMfv378fkyZPh6+sLiUSCH3/80WC/IAhYtmwZfH194ejoiJEjR+Ls2bMGMXq9HvPnz4eHhwecnZ0RFRWFtLQ0g5jc3FxER0dDpVJBpVIhOjoaeXl5Zl8gUWvQlZRh74VMVFQKeOu381i957JRjK/KAVMGdxG3V88YhN+evw9KB86yS0TUFGYnMYWFhRg4cCBWr15tcv8777yD999/H6tXr8bRo0fh7e2NcePGIT+/eqKumJgYbN26FZs3b8bBgwdRUFCASZMmoaKiQoyZMWMGkpKSEBcXh7i4OCQlJSE6OroJl0hkeYu+PYlZ649i6ZZT+O++KyZjXprYB+9NHShu68sqW6t6RETtkkQQhMY/vK99sESCrVu34qGHHgJQ1Qrj6+uLmJgYvPTSSwCqWl28vLzw9ttvY968edBqtfD09MSGDRswffp0AMCtW7fg5+eH7du3Y/z48Th37hyCg4ORkJCAsLCqpvaEhARERETg/PnzCAoKMqqLXq+HXq8Xt3U6Hfz8/KDVaqFUKpt6iUQmVVQK+PxgCsJ6uGFA187ovmSbUczKKf2xdMtpcXttdCjG9/PGR3su43/nMvDl7DB0UrBvPRFRTTqdDiqVqlG/3y3aJyYlJQUajQaRkZFimUKhwIgRIxAfHw8ASExMRFlZmUGMr68vQkJCxJhDhw5BpVKJCQwAhIeHQ6VSiTG1rVy5Unz0pFKp4Ofn15KXRmRg64mbeHP7OUSt/gMlZRUmYx67u5vB9p3Ovs+O6oUtz9zDBIaIqJlaNInRaKpW5PXy8jIo9/LyEvdpNBrY29vD1dW13hi1Wm10frVaLcbUtnTpUmi1WvGVmpra7OshqsuZm1rx/b1v/96oY0or+PiIiKglWWR0kkRiOMOoIAhGZbXVjjEVX995FAoFlEqlwYuopaTeLsLId/fg84MpAAC7Gt/D7ILSRp1jkJ9rw0FERNRoLZrEeHtXTdZVu7UkMzNTbJ3x9vZGaWkpcnNz643JyMgwOn9WVpZRKw9Ra1i16yKu5RThtV+TAQACjLuS2cvsMMTfMFHp0tkRAPDmwyHo5u5k+YoSEXUgLZrEBAQEwNvbG7t27RLLSktLsW/fPgwbNgwAEBoaCrlcbhCTnp6OM2fOiDERERHQarU4cuSIGHP48GFotVoxhqg1FZVW93v58cRN3MorNop5/cF+eOuRAejj7YLVM6oWcNz+/H34df69eDzMv9XqSkTUUZjds7CgoACXL1fPgZGSkoKkpCS4ubmhW7duiImJwYoVKxAYGIjAwECsWLECTk5OmDFjBgBApVJh9uzZWLRoEdzd3eHm5obFixejf//+GDt2LACgb9++mDBhAubMmYO1a9cCAObOnYtJkyaZHJlEZGkyafXjo5hvkoz2z743ANOHVnXkjYsZLparHOVQdeGMvERElmB2EnPs2DGMGjVK3F64cCEAYObMmYiNjcWLL76I4uJiPPPMM8jNzUVYWBh27twJFxcX8ZhVq1ZBJpNh2rRpKC4uxpgxYxAbGwupVCrGbNy4EQsWLBBHMUVFRdU5Nw2RpdnV0RfroxmDMbqPGo72UpP7iYjIcpo1T4w1M2ecOVF9CvXleHrjcey/mGW0748lo8V+L0RE1Hzm/H5zogqiGm7kFCHubDoeD/OHs0KGTw9cxRvbztUZ7610aMXaERFRTUxiiGqI+ugg8orKkK4twasPBNebwIzpo4bUrv6pA4iIyHKYxBDVkFdUBgD443I2UnOL6oz7JDoU44I53J+IqC1ZZLI7Ilt3PacII97dW+d+9072DU7gSERElsUkhsgEfXn9SwQ4c90jIqI2x3+JqcO5nJmP1NvFGNWnan2utNwiVFQKJhMXH5UD0rUlRuVStsIQEbU5JjHU4Yx9fz8A4Nt5EZj5+REU17EK9cODumDV9LugLSrDk18exeg+Xoi/ko3sglIEeDi3ZpWJiMgEJjHUodScFun1X5PrTGAAYGDXqpl2VU5yfPdU1XIXT43oAcD0AqVERNS62CeG2i1tcRme+PwIfkhME8t0JeXi+2vZhfUePzTAzahMIpEwgSEishJsiaF2a93+q9h/MQv7L2bhkdCuSL6lw+o9l8T9+fryeo4Ggn040zMRkTVjEkPtRklZBW4XlsL3z2UAcgr1Bvvf2XEeey8YLx1wxz293LFscj9k6PTo6urIFhciIivHJIbajUc/SUBSah6CfZTw7ewAT5fqJQHScovqTWAA4K0pA+Dn5oRAL5d644iIyDowiaF2Iyk1DwCQnK5DcrrOYN9DH8VDLpWgrML0eqcjenvCz83J0lUkIqIWxCSG2oXKyvoXY88u0Ne739NF0ZLVISKiVsAkhtoFbXFZk4577O5uuJZdiBcnBLVwjYiIyNI4xJpslkZbgtmxRxF/JRtZ9bS0dKqxREA3NyfELxktbk8b0hVfzw2Hukb/GSIisg1siSGboispgyAAKkc5Fn6bhPgrOfjf+UyEmZjT5Y4DL47CoNd3AQByi0rho3KAnQSoFICe6k6tVXUiImphTGLIJgiCgLIKAQ/85wBSbxdjzn0BiL+SI+4/nHK7zmNdne2xYEwg/vO/S3j9wRBIJBIkvjIO+vJKKB3krVF9IiKyAIlQcx72dkSn00GlUkGr1UKp5KRltkyjLcGDHx2Et9IBJ9O0DcbPvjcAFzPykXA1BxufDMfdAW6orBRwNbsQPT2dOf8LEZEVM+f3my0xZNUqKwVM+vAAsgtKkaGrf4QRAMQvGS1OdleTnZ0EvfjoiIioXWHHXrJq+y5mIbugtM792xbci39NHShuqzlUmoiow2BLDFmlL+Kv4Vy6Dg5yab1x/XxV8FU5opNChmAfJWRS5uVERB0FkxiyOmUVlfjnz2fr3P/ag/2w82wGhvf2AFDVcTd+6Wg4yOpPeIiIqH1hEkNW5Vy6Du/tvFhvTICHM756MsygjKOMiIg6HiYx1ObO3NTiSlYBHryrCyZ+cKDBeB8VJ6YjIiImMWQFJn14EACwMzmjzpjdC0fgXLoOGboS9FJzlWkiImISQ1Zk26l0g+1AdSd4qxzw+ayhkEvtOESaiIgMMImhVnc1qwC3C0uhcpTjx6SbdcZ991QEOjvZt2LNiIjIljCJoVYlCAJGv7evUbEqR3bWJSKiunFSDWpVqbeLGx3L5QGIiKg+TGKoVZ2+Wf/aR/9bNAIAcF+gR2tUh4iIbBgfJ1Gr0uhK6t3f07MTEl8Zy0dJRETUICYx1KryiupeB2lwt84AAPdOXP+IiIgaxiSGLE5bVIZfT9/C/otZ2HHWcC6YNx4KwaND/XDgUjYGd3NtoxoSEZEtYhJDFvf6tmR8n5hmct//hfsDAEb1UbdmlYiIqB1gx16yuLoSGCIiouZgEkMWdSuv8UOqiYiIzMEkhiwq+rPDbV0FIiJqp9gnhiyipKwCe85n4kpWYZ0x/XyVrVgjIiJqb1q8JWbZsmWQSCQGL29vb3G/IAhYtmwZfH194ejoiJEjR+Ls2bMG59Dr9Zg/fz48PDzg7OyMqKgopKWxX4UtEAQBp9Ly0OfVODy98bhY7qNyQC91Jwzxd8WBF0fhjYdCEPvXu9uwpkREZOss0hLTr18/7N69W9yWSqXi+3feeQfvv/8+YmNj0bt3b7zxxhsYN24cLly4ABcXFwBATEwMfvnlF2zevBnu7u5YtGgRJk2ahMTERINzkXXQlZRhz/lMuDrZ47/7riD+So5RzKGlYwy274xKIiIiaiqLJDEymcyg9eUOQRDw73//G//4xz8wZcoUAMAXX3wBLy8vbNq0CfPmzYNWq8Vnn32GDRs2YOzYsQCAr776Cn5+fti9ezfGjx9v8jP1ej30er24rdPpLHBlVNuRlNuYtvZQvTHrnhjSSrUhIqKOxCIdey9dugRfX18EBATg0UcfxdWrVwEAKSkp0Gg0iIyMFGMVCgVGjBiB+Ph4AEBiYiLKysoMYnx9fRESEiLGmLJy5UqoVCrx5efnZ4lLoxpKyysbTGCurrgf44K9WqlGRETUkbR4EhMWFoYvv/wSO3bswLp166DRaDBs2DDk5ORAo9EAALy8DH/UvLy8xH0ajQb29vZwdXWtM8aUpUuXQqvViq/U1NQWvjKqLbeeJQQA4INH74KdHVeiJiIiy2jxx0kTJ04U3/fv3x8RERHo2bMnvvjiC4SHhwMAJBLDHzZBEIzKamsoRqFQQKHgmjut5VZeMcav2m9yXzc3J/y64F4oHbiIIxERWY7F54lxdnZG//79cenSJbGfTO0WlczMTLF1xtvbG6WlpcjNza0zhtpW4vXbGPbW78jXlxvtmzWsO3YtHM4EhoiILM7iSYxer8e5c+fg4+ODgIAAeHt7Y9euXeL+0tJS7Nu3D8OGDQMAhIaGQi6XG8Skp6fjzJkzYgxZ1oFLWfhoz2VUVgpiWUWlgF9O3kJeUSle+OZknccui+oHhYwjyIiIyPJa/HHS4sWLMXnyZHTr1g2ZmZl44403oNPpMHPmTEgkEsTExGDFihUIDAxEYGAgVqxYAScnJ8yYMQMAoFKpMHv2bCxatAju7u5wc3PD4sWL0b9/f3G0EllW9GdHAABKBxnkUjus2XcFd3d3w3cNrIE0oZ/xiDQiIiJLafEkJi0tDY899hiys7Ph6emJ8PBwJCQkwN+/al6QF198EcXFxXjmmWeQm5uLsLAw7Ny5U5wjBgBWrVoFmUyGadOmobi4GGPGjEFsbCzniGlhe85nwqezA/p4m545d9vpdCRcvQ0AuJ5TZLCvh4czrmZXz8b747P3YEAXleUqS0REVItEEASh4TDbo9PpoFKpoNVqoVRyevvazqXrMPGDAwCA356/DxpdCUYFqaEvr0DQK3ENHr/xyTA8/mnVukj/uL8v5gzvYdH6EhFRx2DO7zfXTuqgLmbki+/vJDMenRQIbsR6RsN6uiOihztCuihx5qYOkwf6WqyeREREdWES00GZGq6eXaDH/otZJuN7e3VCdkEp/rdwBFyd7QEAm+dGoFBfDi+lg0XrSkREZAqTmA6gvKISS7acRnd3JzjZy7Bq90XklxgPj67P9gX3QSKRQFpj8rpOChk6KfgVIiKitsFfoA7g9/OZ+L6BkUW1Deiqwqk0rbgtk1p8ND4REZFZ+MvUARy8nG32MS4OMrj/+dhoYFeOOiIiIuvDlph2SBAE/HzyFraeuIlRQWp8eeh6vfELxgTiSEoO5o8OxE9JN/HtsTQsjgyCXGqHT/Zfxd/HB7VSzYmIiBqPQ6zbEUEQ8HbcBXx+MAWlFZUNxrs722PD7DCDEUklZRXILSqFj8rRklUlIiIyiUOsO6iP917Bf/ddqXP/I4O7YvvpdNzTyx1KBznemzbQaJSSg1zKBIaIiGwCk5h2oLyiEunaEry740KdMQdeHIUunR3x5sMhcJBz5mMiIrJ9TGLagRe/P4UtJ26a3CeRADtihsPPzQkA4GDHBIaIiNoHJjHtQF0JzOa54eju7gxvFSejIyKi9odJTDsW5OUizq5LRETU3nCeGBv36YGrde7r7CRvxZoQERG1LiYxNmzzkRt4Y9s5o/IgLxccf3WcyfWRiIiI2gs+TrJRd9ZDusNH5YB0bQmWTOyDp0b0bMOaERERtQ4mMTYoQ1eCv8UeNSjbs3gkh04TEVGHwiTGxuw8q8HcDYnittpFgXVPDGECQ0REHQ6TGBuSri02SGAA4Lfn74N7J0Ub1YiIiKjtsGOvDfnvXuMlBZjAEBFRR8WWGBux5XgavvhzNeqIHu7ILtDj5Qf6tnGtiIiI2g6TGBtQoC/Hwm9PittrnwiF0oFzwBARUcfGJMaKCYKAzw6mGM0FwwSGiIiIfWKs2tlbOpOT2RERERGTGKt2IjXPqOytKf1bvyJERERWiI+TrFBlpYDfzmjw6o9nDMrXPD4YE/v7tFGtiIiIrAuTGCu0es9lvL/rorg9xN8VDw/uggkh3m1YKyIiIuvCJMYK1UxgAOCjxwfDS+nQRrUhIiKyTuwTY2Vu5BQZbL83dSATGCIiIhOYxFiRxOu5GP7uHoOyAE/nNqoNERGRdWMSYyW2n07HI2viDcoCPJzRz1fZRjUiIiKybuwTYwWOpNzGMxuPG5Qtj+qHx8O6QSZlnklERGQKkxgrMGNdgvh+oF9nOMmleHhwFyYwRERE9WAS08bStcUorxTE7TcfCkFIF1Ub1oiIiMg28E/9NrTjrAYRK383KOvS2bGNakNERGRbmMS0ka8SrmPehkSDskcGd0VnJy7uSERE1Bh8nNQG4i9n45UaSwooHWTYvXAE1JwPhoiIqNGYxLSyY9duY1bsUXF7ZoQ//jm5H+zsJG1YKyIiItvDJKaVlJZX4stD1/DGtnNiWdRAXywcF8QEhoiIqAmYxLSCzw+m4LVfkw3K9v19JPzdORsvERFRU1l9x96PP/4YAQEBcHBwQGhoKA4cONDWVWqU3MJS/HE5G09+cdQggXk8rBt2LxzOBIaIiKiZrLol5ptvvkFMTAw+/vhj3HPPPVi7di0mTpyI5ORkdOvWra2rZ9LVrAK89dt57EzOMNq3NjoU4/t5t0GtiIiI2h+JIAhCw2FtIywsDIMHD8aaNWvEsr59++Khhx7CypUrDWL1ej30er24rdPp4OfnB61WC6Wy5dYfupyZj42Hb6C0vLLqVVGJ24WlyMrXo7C0HKm3i42OGdvXCy9NCEKgl0uL1YOIiKg90ul0UKlUjfr9ttqWmNLSUiQmJmLJkiUG5ZGRkYiPjzeKX7lyJZYvX27xet3MK8H6P641GOetdMCT9wVg6hA/qBw59wsREVFLs9okJjs7GxUVFfDy8jIo9/LygkajMYpfunQpFi5cKG7faYlpaf5uTnh6ZE/YS+2gkNvBXmoHZ4UMMjsJvJQOCPJ2QWcnORQyaYt/NhEREVWz2iTmDonEcPixIAhGZQCgUCigUCgsXp/uHs54aUIfi38OERER1c9qRyd5eHhAKpUatbpkZmYatc4QERFRx2O1SYy9vT1CQ0Oxa9cug/Jdu3Zh2LBhbVQrIiIishZW/Thp4cKFiI6OxpAhQxAREYFPPvkEN27cwFNPPdXWVSMiIqI2ZtVJzPTp05GTk4PXXnsN6enpCAkJwfbt2+Hv79/WVSMiIqI2ZtXzxDSHOePMiYiIyDqY8/tttX1iiIiIiOrDJIaIiIhsEpMYIiIisklMYoiIiMgmMYkhIiIim8QkhoiIiGwSkxgiIiKySUxiiIiIyCZZ9Yy9zXFnDj+dTtfGNSEiIqLGuvO73Zi5eNttEpOfnw8A8PPza+OaEBERkbny8/OhUqnqjWm3yw5UVlbi1q1bcHFxgUQiadFz63Q6+Pn5ITU1lUsaWBDvc+vgfW49vNetg/e5dVjqPguCgPz8fPj6+sLOrv5eL+22JcbOzg5du3a16GcolUr+D9IKeJ9bB+9z6+G9bh28z63DEve5oRaYO9ixl4iIiGwSkxgiIiKySUximkChUOCf//wnFApFW1elXeN9bh28z62H97p18D63Dmu4z+22Yy8RERG1b2yJISIiIpvEJIaIiIhsEpMYIiIisklMYoiIiMgmMYkhIiIim8QkBsDHH3+MgIAAODg4IDQ0FAcOHKg3ft++fQgNDYWDgwN69OiB//73v0YxP/zwA4KDg6FQKBAcHIytW7daqvo2o6Xvc2xsLCQSidGrpKTEkpdhE8y51+np6ZgxYwaCgoJgZ2eHmJgYk3H8Thtr6fvM77Rp5tznLVu2YNy4cfD09IRSqURERAR27NhhFMfvs2ktfa8t/p0WOrjNmzcLcrlcWLdunZCcnCw8//zzgrOzs3D9+nWT8VevXhWcnJyE559/XkhOThbWrVsnyOVy4fvvvxdj4uPjBalUKqxYsUI4d+6csGLFCkEmkwkJCQmtdVlWxxL3ef369YJSqRTS09MNXh2dufc6JSVFWLBggfDFF18Id911l/D8888bxfA7bcwS95nfaWPm3ufnn39eePvtt4UjR44IFy9eFJYuXSrI5XLh+PHjYgy/z6ZZ4l5b+jvd4ZOYu+++W3jqqacMyvr06SMsWbLEZPyLL74o9OnTx6Bs3rx5Qnh4uLg9bdo0YcKECQYx48ePFx599NEWqrXtscR9Xr9+vaBSqVq8rrbO3Htd04gRI0z+uPI7bcwS95nfaWPNuc93BAcHC8uXLxe3+X02zRL32tLf6Q79OKm0tBSJiYmIjIw0KI+MjER8fLzJYw4dOmQUP378eBw7dgxlZWX1xtR1zvbOUvcZAAoKCuDv74+uXbti0qRJOHHiRMtfgA1pyr1uDH6nDVnqPgP8TtfUEve5srIS+fn5cHNzE8v4fTZmqXsNWPY73aGTmOzsbFRUVMDLy8ug3MvLCxqNxuQxGo3GZHx5eTmys7PrjanrnO2dpe5znz59EBsbi59//hlff/01HBwccM899+DSpUuWuRAb0JR73Rj8Thuy1H3md9pQS9zn9957D4WFhZg2bZpYxu+zMUvda0t/p2UtchYbJ5FIDLYFQTAqayi+drm55+wIWvo+h4eHIzw8XNx/zz33YPDgwfjwww/xn//8p6WqbZMs8f3jd9pYS98TfqdNa+p9/vrrr7Fs2TL89NNPUKvVLXLO9q6l77Wlv9MdOonx8PCAVCo1yjIzMzONstE7vL29TcbLZDK4u7vXG1PXOds7S93n2uzs7DB06NAO+1cr0LR73Rj8Thuy1H2uraN/p5tzn7/55hvMnj0b3333HcaOHWuwj99nY5a617W19He6Qz9Osre3R2hoKHbt2mVQvmvXLgwbNszkMREREUbxO3fuxJAhQyCXy+uNqeuc7Z2l7nNtgiAgKSkJPj4+LVNxG9SUe90Y/E4bstR9rq2jf6ebep+//vprzJo1C5s2bcIDDzxgtJ/fZ2OWute1tfh32mJdhm3EnSFln332mZCcnCzExMQIzs7OwrVr1wRBEIQlS5YI0dHRYvydob8vvPCCkJycLHz22WdGQ3//+OMPQSqVCm+99ZZw7tw54a233urww/cscZ+XLVsmxMXFCVeuXBFOnDgh/PWvfxVkMplw+PDhVr8+a2LuvRYEQThx4oRw4sQJITQ0VJgxY4Zw4sQJ4ezZs+J+fqeNWeI+8zttzNz7vGnTJkEmkwkfffSRwZDevLw8MYbfZ9Msca8t/Z3u8EmMIAjCRx99JPj7+wv29vbC4MGDhX379on7Zs6cKYwYMcIgfu/evcKgQYMEe3t7oXv37sKaNWuMzvndd98JQUFBglwuF/r06SP88MMPlr4Mq9fS9zkmJkbo1q2bYG9vL3h6egqRkZFCfHx8a1yK1TP3XgMwevn7+xvE8DttrKXvM7/Tpplzn0eMGGHyPs+cOdPgnPw+m9bS99rS32mJIPzZW5KIiIjIhnToPjFERERku5jEEBERkU1iEkNEREQ2iUkMERER2SQmMURERGSTmMQQERGRTWISQ0RERDaJSQwRERHZJCYxREREZJOYxBAREZFNYhJDRERENun/A1AIg0vB5RtiAAAAAElFTkSuQmCC",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "(xs, ys) = load_radspec(\"../cmake-build-debug/rdffiles/rdf_0.250000.txt\")\n",
+ "\n",
+ "plt.plot(xs,ys)\n",
+ "plt.show()"
+ ]
}
],
"metadata": {
diff --git a/utk_experiments/utk b/utk_experiments/utk
index a12c408..3a248e3 160000
--- a/utk_experiments/utk
+++ b/utk_experiments/utk
@@ -1 +1 @@
-Subproject commit a12c408ddac03425444c857a725572bfe6fa4837
+Subproject commit 3a248e347abe6c377a0e87a9434fcec2782a6545