{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Estimate moments from Taylor expansions" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In this example we are going to estimate mean and standard deviation of an output variable of interest thanks to the Taylor variance decomposition method of order 1 or 2." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from __future__ import print_function\n", "import openturns as ot" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Create a composite random vector\n", "ot.RandomGenerator.SetSeed(0)\n", "input_names = ['x1', 'x2', 'x3', 'x4']\n", "myFunc = ot.SymbolicFunction(input_names,\n", " ['cos(x2*x2+x4)/(x1*x1+1+x3^4)'])\n", "R = ot.CorrelationMatrix(4)\n", "for i in range(4):\n", " R[i, i - 1] = 0.25\n", "distribution = ot.Normal([0.2]*4, [0.1, 0.2, 0.3, 0.4], R)\n", "distribution.setDescription(input_names)\n", "# We create a distribution-based RandomVector\n", "X = ot.RandomVector(distribution)\n", "# We create a composite RandomVector Y from X and myFunc\n", "Y = ot.CompositeRandomVector(myFunc, X)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# We create a Taylor expansion method to approximate moments\n", "taylor = ot.TaylorExpansionMoments(Y)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[0.932544]\n" ] } ], "source": [ "# get 1st order mean\n", "print(taylor.getMeanFirstOrder())" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[0.820295]\n" ] } ], "source": [ "# get 2nd order mean\n", "print(taylor.getMeanSecondOrder())" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[ 0.0124546 ]]\n" ] } ], "source": [ "# get covariance\n", "print(taylor.getCovariance())" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
[x1 : 0.181718, x2 : 0.0430356, x3 : 0.0248297, x4 : 0.750417]
" ], "text/plain": [ "class=PointWithDescription name=Unnamed dimension=4 description=[x1,x2,x3,x4] values=[0.181718,0.0430356,0.0248297,0.750417]" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# draw importance factors\n", "taylor.getImportanceFactors()" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3XlYVOXfBvB7hn2RAWUVkEUR9yXFSlLMLV8wTVErU1zSIkGztFxTs1wqzV2UVCzKTC1b3MrdTM2VXHIBBfcVFREXlnneP4j5OQ77nOEMM/fnuriKM2ee8z3MONycZzkKIYQAEREREZkNpdwFEBEREVHFYgAkIiIiMjMMgERERERmhgGQiIiIyMwwABIRERGZGQZAIiIiIjPDAEhERERkZhgAiYiIiMwMAyARERGRmWEAJCLZTZ8+HQEBAbCwsEDz5s3lLscoHT9+HAqFAmvWrJG7FJOzbt06KBQKHDx4UO5SiCoMA6CZWL58eaX/gFu4cCGWL18udxmSUSgUhX55enoa5HjHjx/HpEmTcOHCBYO0X14bNmzAmDFjEBYWhoSEBHzyySdyl1So3NzcIl+zp792794td7lEJVq9ejUaN24MW1tb+Pv7Y+rUqcjLy5O7LKoglnIXQFRaCxcuhKurK/r37y93KZLp0KEDoqKitLbZ2dkZ5FjHjx/Hxx9/jPbt26NGjRoGOUZ5bNu2DZaWlliyZAksLY33I8nCwgKJiYla2xISErBjxw58/fXXWtuDg4MrsjTSU3h4OB4+fAgbGxu5S6kwa9asQa9evdCpUyfExsbi8OHDGDduHO7cuYMvvvhC7vKoAhjvpy3Rfx48eAB7e3u5yzCI2rVro0+fPnKXoRd9X58bN27AwcGhxPCnVquRnZ0NW1vbch9LHwqFQue12r17N3bu3FkpX8OsrCw4ODjIXYZRUCqVsr2v5KBWq/HBBx/g+eefx/r166FU5ncG2tjYYPbs2YiNjYWfn5/MVZKhsQvYjPXv3x+Ojo64cOECOnfuDEdHR3h7e2PBggUAgGPHjqFt27ZwcHCAn58fVqxYofX8gm7lXbt24e2330a1atXg5OSEqKgo3LlzR+d4CxcuRP369WFjY4Pq1asjJiYGd+/e1dqnTZs2aNCgAQ4dOoTWrVvD3t4eY8eOhb+/P06cOIGdO3dqutnatGkDALh9+zZGjhyJhg0bwtHREU5OTvi///s//PPPP1pt79ixAwqFAqtWrcKUKVPg4+MDW1tbtGvXDikpKTr1/v333wgPD4eLiwscHBzQqFEjzJkzR2ufU6dOoUePHqhatSpsbW3RvHlz/Prrr2V+LYqSlJSEqKgoBAQEwMbGBp6enhg0aBBu376ts+/FixcxcOBAeHl5wdbWFoGBgYiJiUFubi6WLFmC119/HQDQqlWrQrsq582bh3r16mlen6FDhyIjI0PrGC+88AKaNGmC/fv3o1WrVrC3t8eECRMAAPv370eHDh1QrVo12NnZISAgAIMHDy7y3Aq6VBMTE5GRkaGp6dtvv9U8Nnz4cHzzzTeaurZs2QIAyMzMxHvvvQcfHx/Y2NigTp06mDVrFoQQOu0PHz4cP/zwA+rWrQs7Ozu0bNkSJ06cAJD/nqxZsyZsbW3Rtm1bSbvHs7KyMHbsWDRt2hROTk5wdHTEiy++iD179mjV6OHhgTfeeEPn+ffu3YOdnR1GjBhR7HE2btyI559/Hvb29nBxcUGPHj1w9uxZrX1GjhwJhUKBc+fOoUePHlCpVOjUqVOx7d66dQtDhgyBt7c3bGxsULt2bcyePVtrnxEjRsDCwgL79u3T2t67d2/Y29vjzJkzAP43xu7XX3/FiBEj4O7uDkdHR0RGRuLatWtaz92yZQu6d++ueW39/PwwevRoPH78WGu/Hj16wNPTE2lpaQgPD4eDgwM8PDwwfvx4rfcBkP9Z1aRJEzg4OEClUqFx48ZYvHix5vGixgAmJiaiSZMmsLW1hbu7OwYOHIjr168brA59NWvWDKGhoTrb1Wo1fH19ERkZCQA4ePAg0tLSEB0drQl/ADSfF2vXrpWsJjJigsxCQkKCACAOHDig2davXz9ha2sr6tWrJ6Kjo8WCBQtEy5YtBQCRkJAgqlevLj744AMxb948Ub9+fWFhYSHOnTun02bDhg1Fq1atxNy5c0VMTIxQKpWidevWQq1Wa/adOHGiACDat28v5s2bJ2JjY4WFhYUICQkR2dnZmv3CwsKEp6encHNzE0OHDhWLFy8WP//8s1i7dq3w8fERderUEYmJiSIxMVH88ccfQgghDhw4IGrWrClGjx4tFi9eLCZPniy8vb2FSqUSly9f1rS9fft2AUA0bdpUNGvWTMyaNUtMmjRJ2NvbixYtWmj9vP744w9hbW0t/Pz8xMSJE0VcXJwYNmyYaN++vWaf48ePC5VKJerVqyc+++wzMX/+fNG6dWuhUCjETz/9VOJrAkC8+eab4ubNm1pfjx490uwzffp0ERYWJiZPnizi4+PFsGHDhK2trXj++ee12rp48aLw9PQU9vb24r333hOLFy8W48ePF8HBwSIzM1OkpKSImJgYAUB89NFHmp/hjRs3hBBCjBs3TgAQHTt2FPPmzRNDhgwRSqVSPPfccyInJ0dznNDQUFG9enXh6uoq3n33XbF48WLx66+/iqtXrwqVSiXq1KkjZsyYIeLj48XYsWNFgwYNijx/tVotEhMTRcuWLYW9vb2mpnPnzomcnBwBQNStW1d4eHiIyZMni/nz54t//vlH5OXlaX7OgwcPFvPnzxcRERECgBg5cqSm/YI2GjduLPz8/MT06dPF1KlThZOTkwgICBBz5swRDRo0EF9++aUYN26csLKyEh06dCjxdXvS22+/LSwsLAp9LDU1Vfj4+IgPPvhALFq0SEybNk0EBgYKW1tbcerUKc1+7777rnBwcBD379/Xev6yZcsEAHH48GEhhBDHjh0TAMTq1as1+/zyyy9CqVSKBg0aiBkzZogJEyYIlUol3N3dxZUrVzT7jRgxQgAQ9erVEz179hRxcXEiPj6+yPPKyMgQwcHBwt3dXUyYMEEsWrRIvPbaawKAGDdunGa/hw8fijp16ojatWuLBw8eCCGEWLNmjQAgZs2apdnvt99+03xWNGvWTMyePVuMGDFCWFlZiUaNGml9BgwcOFB07dpVTJ8+XSxatEhERUUJpVIp+vfvr1VjZGSkUKlUolatWuLtt98WcXFxonPnzgKA+OabbzT7rV27VgAQ4eHhIi4uTsyfP19ER0eLqKgonfqe/HycN2+eACBatmwp5syZI0aOHClsbGxEcHCw1mslZR36mjVrlgCg9TkthBDbtm0TADSfS4sWLRIAxIkTJ3TacHZ2lrQmMl4MgGaiqAAIQEydOlWz7c6dO8LOzk4oFAqxcuVKzfZTp04JAGLixIk6bTZr1kzrA/zzzz8XAMQvv/wihBDixo0bwtraWnTs2FHk5eVp9ps/f74AIJYtW6bZFhYWJgCIRYsW6ZxD/fr1RVhYmM72R48eabUrRP4vXxsbGzF58mTNtoIAWLduXfH48WPN9jlz5ggA4tixY0IIIXJzc0VAQIDw8/MTd+7c0Wr3yVDbrl070bBhQ63AplarRcuWLUVQUJBOnU8DUOhXQkKCZp+CX6pPSkxMFADEnj17NNt69+4tLCwsNGGhsHq///57AUD8+eefWvtcvXpVWFpaivDwcK2f4+zZs3V+iYWGhgoAYsmSJVptrF69WgAQR44cKfG8n/bGG28IlUqlta0gvFlYWGiFJSH+FzCmT5+udZ7dunUTSqVSpKamarVha2srLly4oNl3wYIFAoCoXr26yMzM1Gz/4IMPBABx8eLFUtdeXADMzs7WCs9C5P9bcHZ2FsOGDdNsO3TokAAgEhMTtfZt27atqF+/vub7wgJgUFCQqFGjhsjIyNBs27t3rwAghgwZotlWEAAHDRpUqvMaNWqUUKlU4vz581rbY2NjhbW1teYPByGE+Pvvv4WFhYUYPny4uHHjhnBzc9P5A7AgYNWqVUtkZWVptheE3KVLl2q2FfaeHzdunLCwsBA3b97UbIuMjBQAxJdffqnZplarRXBwsNbnxJtvvim8vLy06nna0wEwKytLODk5iZCQEK3PipUrVwoA4vPPPzdIHfq6fv26sLS0FJ988onW9oEDB4qqVatqzqXgD/Lbt2/rtFG3bl2tP3TJdLELmDBo0CDN/zs7OyM4OBgODg7o1auXZntwcDCcnZ1x7tw5nee/9dZbsLKy0nz/zjvvwNLSEhs2bACQ36WTnZ2N4cOHa3U3DB48GE5OTli/fr1WezY2NhgwYECp67exsdG0m5eXh/T0dDg6OiI4OBiHDx/W2X/AgAGwtrbWfN+qVSsA0JzbkSNHkJqaiuHDh8PZ2VnruQqFAkB+t/O2bdvQq1cvZGZm4tatW7h16xbS09Px0ksvITk5GZcvXy6x9q5du2Lz5s1aXy+99JLm8ScnhDx69Ai3bt3Cc889BwCac8vNzcUvv/yCV155BU2bNi203uJs3rwZubm5Oq/P22+/DUdHR53Xx97eXmfiSsHPad26dcjNzS3xmKXVtm1bnQkVGzZsgJWVFWJjYzXbFAoF3n//fajVamzatElr/44dO8LX11fz/bPPPgsA6NmzJxwdHXW2F/YeLw8rKyvNuEa1Wo309HQoFAo0adJE6335zDPPoH79+vjuu+80265cuYIdO3YUO7YwOTkZycnJGDRoEJycnDTbn3vuOYSGhuq8bkD+v83SWL16Ndq3bw97e3vNe/vWrVvo0KEDsrOztbqxW7RogdGjR2POnDmIiIjAgwcPkJCQUOh7b8CAAVrjRfv27QtnZ2fNZwWg/Z7PysrCrVu3EBoairy8PJ1hHQqFAm+99ZbW96GhoVqvobOzM27fvo2dO3eW6twB4K+//sK9e/cwdOhQrc+KXr16wdfXV+dna6g6ysrd3R0vvfSS1nvp8ePH+PHHH9GzZ0/NuTx8+BAACp30Ymtrq3mcTBsDoJmztbWFm5ub1jaVSgUfHx+dD3CVSlXo2L6goCCt7x0dHeHl5YW0tDQAwPnz5wHozoy0trZGYGCg5vEC3t7eWh+6JVGr1Zg1axaCgoJgY2MDV1dXuLm54ejRozpj2ADozIB1cXEBAM25FYyfatCgQZHHTElJgRACH330Edzc3LS+Jk6cCCB/ckNJfHx80L59e60vLy8vzeO3bt3C0KFD4e7uDjs7O7i5uWl+3gXndv36dWRlZRVbb3GKen0KloZ4+vXx8fHRCvxAflB75ZVX8NFHH8HV1RWvvPIKli9frjNuq6wCAgIKrdfHx0dnAkPdunW1zqfA06+3SqUCAK1Q+OT2wt7j5RUfH68Zv1jwvtyxY4fO+zIqKgpbtmzRvGe+//57CCEKHRtYoKjXDcj/WVy8eFFnDFphP8+n5eXlITU1FT/++KPOe7tr164AdN/bEyZMQK1atXDgwAFMmzYNgYGBhbb99GeFpaUlAgICNJ8VQP6/v969e8PFxQWOjo5wc3NDeHg4AOj83KpVq6bzPnBxcdF6Dd999134+PjgxRdfhJ+fHwYPHoytW7cW+zMo6merUCgQHBys8x4zVB1A/s/62rVrmq8HDx4Uu39UVBROnTql+SNj3bp1yMjIQN++fTX7FITswv59Pnr0yGArEZBx4SxgM2dhYVGm7U//QjGEsn74TJ06FR999BEGDhyITz75BFWrVoVSqcTw4cOhVqt19pfi3AraHTlypNYVuyfVqlWr1O0VpUePHjhw4AA+/PBDNG7cGA4ODsjJyUFERESh51YRCnt9lEol1q5diz179mDdunX4/fffMWDAAMyaNQt79uwp92xTKX4RyfUeX7RoEd555x306tUL48aNg6urKywsLDBx4kSdINOnTx+MGTMGK1euxLBhw/Dtt98iLCxMJ6TqqzQ/T5E/NAhdunTB0KFDC92nIGwXOH36tCYUHTt2rNz1PX78GG3btkVOTg7GjRuH2rVrw97eHmfPnkV0dLTOe740r6Gvry+OHz+OjRs3YtOmTdi4cSOWLFmCIUOGaCa86cuQddSrVw/p6ema77/44guMHDmyyP27dOkCZ2dnfPvtt3jmmWfw7bffIiAgQGtySMEfmVevXtX8AVzg6tWrCAkJKf6EySQwAJLekpOT8eKLL2q+v3//Pq5evar5q71gOYHTp09rXRnIzs5Gamoq2rdvX6rjFNWduWbNGrz44otYunSp1va7d+/C1dW1TOcCADVr1gSQv25eUbUVnIeVlVWp6y+rW7duYefOnZgyZQrGjh2r2X7y5Emt/Tw8PODg4IDjx48X215RP78nX58nr5Y9fvwYaWlpiIiIKHXNLVu2RMuWLTF16lR888036NevH1avXi3p2o1+fn7YtWuXzjImp06d0jxuDNasWYOGDRvihx9+0Npe2Kze6tWro127dvjuu+/Qvn17JCUl6byfn/bk6/a0U6dOwdfXt1RDAJ5maWmJGjVq4OHDh6V6b+fk5KBfv37w8vJCjx49MHPmTPTo0QMdO3bU2Tc5OVnr+9zcXKSmpqJdu3YAgAMHDuDChQv48ccf0b17d81+mZmZZT6PJ9na2qJbt27o1q0b8vLy0L9/fyxcuBAfffRRoQuvP/mzbdGihWa7EAJnzpwp1ZVUKeoAgJ9++gnZ2dma72vXrl3iMXr27ImVK1di7Nix2LBhA0aNGqW1T5MmTQDkzwauV6+eZvuZM2dw9+5dzeNk2tgFTHqLj49HTk6O5vu4uDjk5ubi//7v/wAA7du3h7W1NebOnav1F/HSpUuRkZFR6oDh4OCgs2wMkP/X99NXbVavXl2qMXiFeeaZZxAQEIDZs2frHK/gOO7u7mjTpg0WL16Mq1ev6rRx8+bNch37SQVXFZ4+t6eX4rC0tETXrl3x888/48iRI0W2VxCWnj6nDh06wNLSUuf1iY+Px/3790v1+ty5c0enzoJfIvp2Az8tPDwcOTk5WLhwodb2WbNmQalUat53civsqtC2bdtw9OjRQvfv27cv9u/fj0mTJsHW1hY9evQotv2goCAEBQVh6dKluH//vmb7/v37sXv37jIF96f16tULW7ZsKfSOJunp6VpX4qZMmYKkpCQsW7YM06ZNQ9OmTTFo0KBCh18kJCRodWEmJibi7t27mtessPd8Xl4e5s2bV+5zefLqWcExGjZsCKDo92ZoaCicnJywYMECrc+2NWvW4MKFC+X62ZanDgBo3bq11hCR0izi3rdvX1y9ehVDhgxBdna2zljSkJAQ1KhRA4sWLdL6WcfFxcHCwgLdunUry6lRJcUrgKS37OxstGvXDr169cLp06excOFCvPDCC+jSpQsAwM3NDWPGjMHHH3+MTp06oUuXLpr9QkJCSr2IbrNmzRAXF4dPP/0UtWrVgru7O9q2bYvOnTtj8uTJGDBgAFq2bIljx47hu+++K3IcUkmUSiXi4uLw8ssvo0mTJhgwYAC8vLxw6tQpnDhxAr///jsAYMGCBXjhhRfQsGFDDB48GIGBgbh+/Tr27t2LS5cu6QxYLysXFxe0bNkS06ZNw6NHj1C9enVs2rRJZ/wRkH8v3a1bt6JVq1Z4++23ERwcjCtXrmDVqlXYv38/HB0d0bRpUyiVSkybNg3p6emwsbFB+/bt4enpiVGjRmHKlCkIDw9H586dcerUKcTFxeG5557TrB9YnKVLl2LJkiV45ZVXEBgYiHv37uGrr76Cs7NzievNlVW3bt3QunVrjBo1CmfPnkWjRo2wadMm/Pbbbxg5cqTRXAHs3Lkzhg0bhl69eqFDhw44c+YMlixZgjp16hS6f/fu3TFkyBCsXr0aPXv21JrYUZQZM2agW7duaNmyJQYMGICMjAzMmTMH7u7uGDduXLlrHz9+PDZs2IB27drhzTffROPGjXHv3j38888/+PHHH3Hnzh3Y2tri8OHDmDJlCmJiYtC2bVsAwNdff43mzZvjvffew7Jly7TatbOzQ1hYGPr27YuLFy9i7ty5aNCggWZ8WpMmTeDj44PY2FikpKTAzs4OP/zwA7Kyssp9Lq+//jrUajXCwsJQvXp1nD17FvPmzcPzzz9f5HvF3t4eU6ZMwdChQ9G2bVv06tVLU29wcDCGDBlSIXWU1wsvvICAgACsXr0aISEhOlcNlUolvvjiC7z66qsIDw9HZGQkDh06hMWLF2P48OHw9/eXtB4yUhU865hkUtQyMA4ODjr7hoWFaS0/UcDPz09ERETotLlz507x1ltvCRcXF+Ho6CjeeOMNkZ6ervP8+fPnizp16ggrKyvh4eEh3nnnHZ1lVoo6thBCXLt2TURERIgqVaoIAJolFh49eiRGjBghvLy8hJ2dnQgNDRV79+4VYWFhWsswFCwD8+QyGkLkLxmDp5ZfEUKI3bt3iw4dOogqVaoIBwcH0ahRIzFv3jytfc6ePSuioqKEp6ensLKyEt7e3qJz585izZo1hZ7DkwCImJiYYve5cOGCeOWVV4RKpRLOzs7i1VdfFZcuXRIAdJZ6SEtLE3369BFubm7CxsZGBAYGitjYWK2lSBYtWiQCAgKEhYWFzpIwc+bM0bw+np6eIiYmRty9e1frGKGhoaJx48Y6dR48eFC89tprwtfXV9jY2AgPDw/RpUsXnWVpClPcMjDvvvtuoc+5d++eePfdd4WXl5ewsrISQUFBYubMmVpLbBTVRnJyss46dUIIsXnzZgFArF27tsSaCxS3DExeXp6YOHGi8PX1Fba2tiIkJERs3rxZREZGFvkej4qKEgDEr7/+qvNYYcvACCHE+vXrxbPPPitsbW2FSqUS3bt3F8nJyVr7FCwD8/Dhw1Kf2507d8TIkSNFQECAsLKyEu7u7qJVq1Zi9uzZIi8vTzx+/Fg0aNBAZ2kXIYSYMmWKACDWr18vhPjfMiu//PKLeP/994Wrq6uwt7cX3bp101qrUwghkpKSRFhYmHBwcBDu7u4iJiZG7Nu3T+fcIyMjhYeHh07dI0aM0Ppc++6770S7du2Em5ubsLa2Fv7+/iI2NlZrSZnC1gEUQoivv/5aNGzYUFhbWwtXV1fRv39/ce3aNa19pKxDShMmTBAAxNy5c4vcZ+XKlZrzq1Gjhvjkk09Ebm6uQeoh46MQogJG9ZNJWr58OQYMGIADBw6gefPmcpdDVOkNHjwYP//8M65cuaIz07oyW7duHV5++WVs3LhR8ivCVLgpU6Zg0qRJuHLlis5KD0QAxwASERmFzMxM/PDDD3jttddMKvxRxcvLy0NCQgI6derE8EdF4hhAIiIZXblyBdu3b8eKFSvw8OHDIpdeISpJRkYGNmzYgE2bNuHs2bOIj4+XuyQyYgyAREQyOnz4MPr06QMvLy8sXry4xGU+iIpy8eJF9O7dG1WrVsWnn36qmZhDVBiOASQiIiIyMxwDSERERGRmGACJiIiIzAwDIBEREZGZYQAkIiIiMjMMgERERERmhgGQiIiIyMwwABIRERGZGQZAIiIiIjPDAEhERERkZhgAiYiIiMwMAyARERGRmWEAJCIiIjIzDIBEREREZoYBkIiIiMjMMAASERERmRkGQCIiIiIzwwBIREREZGYYAImIiIjMDAMgERERkZlhACQiIiIyMwyARERERGaGAZCIiIjIzDAAEhEREZkZBkAiIiIiM8MASERERGRmGACJiIiIzAwDIBEREZGZYQAkIiIiMjMMgERERERmhgGQiIiIyMwwABIRERGZGQZAIiIiIjPDAEhERERkZhgAiYiIiMwMAyARERGRmWEAJCIiIjIzDIBEREREZoYBkIiIiMjMMAASERERmRkGQCIiIiIzwwBIREREZGYYAImIiIjMDAMgERERkZlhACQiIiIyMwyARERERGaGAZCIqILt3r0boaGhqFatGuzs7FCnTh3MmjVL73Z37NiBrl27wsvLCw4ODmjSpAm+++67Ep+nUCh0vlauXKl5/MiRI2jatCkcHR3x8ssv4/bt25rHcnNz0axZM+zfv1/v+omo4jAAEhFVMAcHB8TGxmLXrl04efIkxo8fj/HjxyM+Pl6vdvfs2YNGjRrhxx9/xNGjRzFgwABERUVh3bp1JT43ISEBV69e1Xy98sormscGDRqEtm3b4vDhw8jIyMDUqVM1j82cOROhoaFo0aKFXrUTUcVSCCGE3EUQEZmKmzdvomHDhhg2bBjGjh0LID+YtWnTBhs3bkS7du0KfV737t3h4OCAxMRESeuJiIiAh4cHli1bVuQ+CoUCa9eu1Qp9T7K3t8fhw4dRp04dxMXFYd26dVi/fj3OnTuHTp064dChQ6hSpYqkdRORYfEKIBGRhNzc3LBs2TJMmjQJBw8eRGZmJvr27YvY2Ngiw9+RI0ewZ88ehIWFFdu2QqHA8uXLy1RPRkYGqlatWuJ+MTExcHV1RYsWLbBs2TI8eW2gcePG2Lx5M3Jzc7F161Y0atQIABAdHY3PP/+c4Y+oEmIAJCKSWHh4OAYPHow33ngD0dHRcHBwwLRp03T28/HxgY2NDZo3b46YmBgMGjSo2HaDg4OhUqlKXceqVatw4MABDBgwoNj9Jk+ejFWrVmHz5s2IjIzEkCFDMG/ePM3jS5YswZo1a1CzZk1YW1tjzJgxSExMhL29PUJCQvDSSy+hVq1aGD9+fKlrIyJ5sQuYiMgAHj58iAYNGuDixYs4dOgQGjZsqLNPamoq7t+/j3379mH06NGYP38+Xn/9dUmOv337dnTu3BlxcXGIiooq03MnTJiAhIQEXLx4sdDH09PTERISgl27dmHYsGFo3LgxRo4ciZCQEHz22Wd4+eWXpTgFIjIgXgEkIjKAs2fP4sqVK1Cr1UhLSyt0n4CAADRs2BCDBw/Ge++9h0mTJkly7J07d+Lll1/GrFmzyhz+AODZZ5/FpUuX8Pjx40Iff//99zF8+HD4+Phgx44d6NmzJxwcHBAREYEdO3boWT0RVQQGQCIiiWVnZ6NPnz549dVX8cknn2DQoEG4ceNGsc9Rq9VFBq6y2LFjByIiIvDZZ5/hrbfeKnbfn376CR06dICbmxucnJzw/PPP4/fff0dSUhJcXFxgY2Oj85ytW7fHDhopAAAgAElEQVTi5MmTiI2NBQDk5eUhJycHAJCTk4O8vDyt/VNSUlClShU4OzsXW0t6ejo6deqE6tWrw8bGBr6+voiNjcW9e/c0+3A5GiLpMAASEUls3LhxyMjIwNy5czFq1CjUrl0bAwcO1Dy+YMEC/Pbbb0hOTkZycjKWLl2KGTNmoE+fPsW2W6dOHaxdu7bIx7dv346IiAgMGzYMkZGRuHbtGq5du6YVlNauXYs6deoAAHbt2gVvb2+89dZb+P7779G4cWOEh4fj008/xdChQ3Xaf/ToEWJjYxEfHw+lMv/XR2hoKBYsWIB//vkHP/74I0JDQzX75+Tk4PXXX0erVq1K/JkplUp07doVv/76K86cOYPly5djy5YtiI6O1uzD5WiIJCSIiEgy27dvF5aWluLPP//UbEtNTRVOTk5i4cKFQggh5s6dK+rXry/s7e2Fk5OTaNq0qVi4cKHIy8srtm0AIiEhocjH+/XrJwDofIWFhQkhhPjrr7+EhYWFePKjf+PGjaJJkybC0dFRODg4CBsbG9G5c+dCaxk9erQYMWKE1rbk5GQREhIinJycxDvvvKP1vA8//FD06dNHJCQkCJVKVey5FWbOnDnCx8dH872dnZ04efKkEEKIhQsXivDwcCGEEGfPnhVBQUHi3r17ZT4GkbliACQiMlHr168XVlZW4sCBA+LevXsiMDBQvPfee0Xun5eXJ3x9fcW8efOKbbekICqEEFu3bhUBAQEiIyOjXAHw8uXLIiwsTLzxxhuabc8995yYO3euyMnJEZGRkWL06NFCCCE6dOgg1q5dW6b2icwdu4CJiExUaZejKTBjxgzcv38fvXr1KrbdkpajSU9PR//+/bF8+XI4OTmVqebXX38d9vb28Pb2hpOTE5YsWaJ5jMvREEmHy8AQEZmw0ixHAwArVqzA4MGD8csvv6B9+/Z6HbN79+6oXbs2pk+fDgBYvnw5hg8fjrt375b43GvXruHu3bs4c+YMxowZg7CwMCxcuLDQfbkcDVH5MQASEZmw48ePIyQkBDk5OVi7dm2hoWjlypUYOHAgVq9ejYiICL2P6ezsjPv372u+F0JArVbDwsIC8fHxWhNiirN79260atUKV65cgZeXl87j/fr1Q7NmzTBs2DBUrVoVu3fvRr169fDBBx9ArVZj5syZep8LkamylLsAIiIyjCeXowkODsagQYNw7NgxuLu7a/b5/vvvMXDgQKxcuVKS8AcAe/fu1VoO5pdffsFnn32GPXv2wNvbu9TtqNVqACh0eZyC5WgSEhIA6C5HQ0TF4xhAIiITVdJyNCtWrEBUVBRmzpyJZ599VrNsTEZGRrHtlrQcTd26ddGgQQPNl7e3N5RKJRo0aAAXFxcA2svRAMCGDRuQkJCA48ePIy0tDevXr0d0dDRCQ0Ph7++v1X5Zl6MhIl0MgEREJmjHjh2YPXs2EhMT4eTkBKVSicTERPz555+Ii4sDAMTHxyM3NxcxMTHw8vLSfL377rvFtn369OkSQ2JJMjIycPr0ac33dnZ2+Oqrr/DCCy+gbt26eO+999ClSxesW7dO57kff/wxIiIi0KRJE822uXPnIikpCa1bt8bLL7+MyMhIveojMnUcA0hERERkZngFkIiIiMjMMAASERERmRnOAiYis6JWA3ceALezgNsF/80CHmQDagHkqf/3X4e6Z/Gw+nlYwgJWsIAllLCEEnawhgvs4QJ7VIUDXGAPJ9jJfWpERKXGAEhEJuNaBpB8I//rzHXg0h0gPSv/qyDoZTzMD3ilERp1DX9V/6FU+1pACVc4whNO/32p4AknVIczAuGKILgjEG6wgoUeZ0hEJA0GQCKqVO5kAaevawe9gv/PfCTtsRRl2DcPalzHPVzHPfxTxD6WUML/vzD4ISzRANawQwPYoj4s4ChFyUREpcIASERGKzcP+OcSsPdc/te+c8C5WxV3fKGQdpGEXKiRghtIwQ0MxlGcx77/HlHAGjVgj6ZwQEs4IhT2aAYlbCQ9PhFRAQZAIjIa1+/9F/bOAvtSgYPn88fmmRolFPDE0Se2CGTjPLJxHnfxMwBAARvYozkcEfrfV0tYwlWegonI5DAAEpFsbmcBv58ANp0A/kwBUivw6l6pSHwFsIA/VLDGg2L3EXiMLPyFLPyF6/9ts0EwHBGKKmgLFTrBEtUMUh8RmT4GQCKqUIcvAOuPARuPA/vT8mfbGi/DBMCgcnbtPsZpPMZppGMZACUc0AIqhEOFCNjjGWmLJCKTxgBIRAaVmwfsSgbWHgF++Qe4eEfuikrPULdJqonHErSiRhb2IQv7cAUTYAUfOOMVOKMbqqA1FPx4J6Ji8BOCiCQnBLD1FPDt38BvR/O7eisjRVmmAZeBP6Tv687BJdzEfNzEfFigKpzRBVXRF1XwIhRlms9MROaAAZCIJHM+HUjYAyzfm///lZ0w0DVAX6QZpN0CebiNdCxHOpbDGv6ohn6ohv6wgb9Bj0tElQcDIBHp5VEO8NMRYNlfwLbT+Vf/TIdhTsZbawawYWUjDVfxMa5iMqqgDaqhP1zQA0rYV1gNRGR8GACJqFwOpgHL9gDfHwDuFj+htfIyQM+pFxxhb4Au4JIJZGI7MrEdFxALF/SCKwbAEaEy1EJEcmMAJKJSy84FEvcBc7cDRy/JXU0FMMAyMEFGcOVNjUykYynSsRR2aAx3DEdV9IYS1nKXRkQVhAGQiEp0/xGw+E9g1hbg8l25q6k4hhgDWAt5krepj4f4B+cxAJcxBu4YAjcM4fqCRGaAAZCIinTrPjB3G7BgR+WdyWts/GGc6+Dk4hquYAIcP1yHKo+eBUaOBGrUkLssIjIQBkAi0nHhNjDjD2DpX6Z5K7ZSM8AYwBq4KH2jErG/Xw9VvtgPYD+waBHQuzcwahRQt67cpRGRxJRyF0BExuPkVaBfAlBrPDBvu5mHPwCGmAXsgxOStykVz3in/32TkwN8/TVQvz7w6qtASop8hRGR5BgAiQhX7gIDlgMNPga+2QfkGNcwNdlIPQZQBVs4G3gNwPKyzqkB59H7dR8QAli1CqhXDxg6FLhlbDdsJqLyYAAkMmMPsoGPfwNqT8hfvFltUmv46U/qO4HURhVpG5SQx1p/KHKKuTFzTg4wfz5QsyYwdSrw8GHFFUdEkmMAJDJDQgDf7AWCJwCT1gFZUtya1gQJiZeBqSVpa9KxUFeF65ADpdv53j1g3DggKAhYtgxQFxMaichoMQASmZldZ4CQqUC/5cAl45yQarICkSl3CYVy/7MhlOllvKJ3+TLw5ptAkybAhg2GKYyIDIYBkMhMnL0JRC4CwmYChy7IXU1lIe0VwBq4LGl7UlAIO7hFHy9/A8eOARERQIcOwNmz0hVGRAbFAEhk4h7lAGPXAvUm5d+zl8pA4jGAPjgtbYMSqPZvc1idSte/oS1bgIYNgS++API4i4jI2DEAEpmwXWeAxp8A0zbl38aNykbKWcA2sIQrTknWniSEEh7DzkvX3sOHwIcfAi1aAElJ0rVLRJJjACQyQZmPgHe+A9p8CZy5Lnc1BABBUMECxpXCnS+3gO02A4wHOHwYCAkBxowBHj2Svn0i0hsDIJGJ+eNfoP4kYNGu/Nm+pAcJZwHXgoVkbUnFc6IBJ6Xk5gLTpwONGwO7dhnuOERULgyARCYi63H+Vb+X5gAXObtXElLm50A8kLA1/TlmNILDsgq4K8mZM0CbNkB0dP4SMkRkFBgAiUzA7pT8sX6LeKFFWhJeAfTDDcnakoLHXJuKO5gQwOLFQKNGwN9/V9xxiahIDIBElVieGhj3MxA2I3+ZF5KadAHQF8ZzL13bx4FQTSzlws9SOn8eaNUK+PLLij82EWlhACSqpG7cAzrOBqZu5C3cDEWqH6sSCnjiqESt6c/j++pSXtwsm5wcYMQIoGtX4A7HKhDJhQGQqBLadw54ZgqwzfiWlTMpUt0L2B8qWBvJGEDLPHdUHbZf7jKAX38FmjZllzCRTBgAiSqZ+duB1jOAy3flrsT0SbUOYBAqcLxdCdy31IUyM1vuMvKxS5hINgyARJXEg2zgjaXA0JVADm+0UEGkCYA18ViSdvSlFI5we8vIbgfDLmEiWTAAElUCydeBZ6cBK4yg586sSNYFfEuahvTkergZLC8Y6VIsBV3Cx47JXQmRWWAAJDJya48AzacCx6/IXYn5ERLNlPBFmiTt6EVYwiM2We4qinf+PPDCC8DmzXJXQmTyGACJjNj0TUD3RcA93k2rUvM2ghnAVdNawHpfJfgr4t49IDwcWLpU7kqITBoDIJEREgJ49wdgzFq5KzF3+l8B9IIj7I2gC9hjlPw1lFpuLjBoEDBuHO9nSGQgDIBERiY7F3h9CTB3m9yVkBRjAINgr38jeqqS/gzsV5+Ru4yymzoV6NMHeGwck2iITAkDIJERyXwEhM8DfjgodyUESLMMTC3IP2Xb8wu5K9DDihVAx47A7dtyV0JkUhgAiYzE9Xv5t3TbekruSkhK/pB3aRO7B8Fw+uywrDXobdcuoGVL4Nw5uSshMhkMgERGIOUG0PIz4MhFuSshLQq13k3UgLwvqsfyarIeXzKnTwPPPw8cPy53JUQmgQGQSGaHzgOhnwPnKtEYfXMhxfQDH5yQoJXysc71RtX3TehWazduAG3bMgQSSYABkEhGu1OANjOBG5lyV0KF0nMSiAq2cJZxDUD332pC8Vj+MYiSunmTIZBIAgyARDI5kJY/4eM+JzgaMf2uAdZGFYnqKDsLtQqu7xyS7fgGxRBIpDcGQCIZHL0EvDQnf9YvGS99u4BrSVJF+bjtawKL61kyVmBgDIFEemEAJKpgp68BHWYDdx7IXQmVSM9bwQVCnr59hbCBe/S/shy7QjEEEpUbAyBRBUq9BbSbxTF/lYd+AdAXlyWqo2yqngmB1bGbshy7wjEEEpULAyBRBbl0Jz/8Xb4rdyVUWvouBO0LGe6+IRTwfE+e4CmbghD4rxlc9SSSCAMgUQW4cQ9oPyv/CiBVInrMAraBJVxxUrpaSkl1PQS2G1Mr/Liyu3kTCA8Hrl2TuxKiSoEBkMjAbmcB7WcDp6/LXQmVmR4BMAgqWCBXulpKyeNTM55ZdP480LkzkGXCk1+IJMIASGRAj3KAiHnAMTPrkTMV+nQB14KFhJWUjsP9+qiy4GiFH9eoHDoEvP46oNb/Li5EpowBkMhAhAD6Lwf2mWFvnOkofwAMRMVP8/aIc6zwYxql334Dhg+Xuwoio8YASGQgk34DfjgodxWkFz26gP1wQ7o6SsEm2x/OY/ZX6DGN2rx5wOzZcldBZLQYAIkMYMV+YPJ6uasgfenTBeyLFAkrKZnHGl8o8qS4e7EJGTEC+PlnuasgMkoMgEQS+zsVGPi13FWQnJRQwBMVNxbPUu2KakN49U+HWg288Qawnz8boqcxABJJ6GoG0H0R8LjiJ3+SISjKN5HAHypYV+AYQLcd9aHM4E2lC/XgAfDyy8CFC7KVcPXqVfTu3Ru1a9eGUqnEcAnHJw4bNgzNmjWDjY0NmjRpUug+v//+O5577jlUqVIFbm5uiIyMRFpaWpFtpqWl4c0330RAQADs7OxQs2ZNTJw4EdnZ2Vr7tG7dGg4ODmjdurVOe507d8aPP/4oxSmSgTAAEkkkOxeIXARc4ULPJqO8HapBsJG0juIohT3c3zbzmb8luXEDePVVICdHlsM/fvwYbm5uGD9+PBo3bix5+wMHDsSrr75a6GOpqano2rUr2rZti6SkJPz++++4desWunfvXmR7p06dglqtxuLFi3HixAnMmjULixYtwtixYzX7jBgxAt7e3khKSoKXlxdGjhypeeyHH36AUqlEZGSkdCdJkmMAJJJI7PfA3nNyV0GSKuckkJqouKtx1Y41h2XKnQo7XqW1bx8wZozkzd68eROenp6YOnWqZtuePXtgbW2NrVu3AgD8/f0xZ84cREVFQaVSSXr8uXPnIiYmBoGBgYU+fujQIeTl5eHTTz9FzZo18cwzz2DkyJFISkpCThGBuFOnTkhISEDHjh0RGBiILl26YOTIkfjpp580+5w8eRL9+vVDUFAQ+vfvj5Mn8xc9v3v3LsaPH48FCxZIep4kPQZAIgnE7wK+2i13FSS98l0D9EcF3fJFWMBjGNcZKrUvv8xfIkZCbm5uWLZsGSZNmoSDBw8iMzMTffv2RWxsLNq1a6dX2/7+/pg0aZJebTRr1gxKpRIJCQnIy8tDRkYGEhMT0b59e1hZWZW6nYyMDFStWlXzfePGjbFlyxao1Wr88ccfaNSoEQDggw8+QExMDHx9ffWqmwyPAZBITyevAsNXyV0FGUJ5u4B9kSZlGUVyuRQCm50XK+RYJkEIoH9/yccDhoeHY/DgwXjjjTcQHR0NBwcHTJs2Te92a9asCVdXV73aCAgIwB9//IGxY8fCxsYGzs7OuHTpElatKv2HVkpKCubNm4e3335bs23GjBk4deoU/P39kZycjBkzZmDXrl1ISkpCVFQUevXqhcDAQERHR2uNHSTjwQBIpIfsXKD3UuChPEOLyNAU5YuA3hU0A9hjXEaFHMek3L4NvPYakCvtTK0ZM2YgNzcXq1evxnfffQcbG/3HgW7duhWxsbF6tXHt2jUMHjwY/fr1w4EDB7Bz505YW1ujR48eEKLk9/fly5fRqVMn9OzZE4MHD9Zs9/b2xrp163DhwgWsW7cOrq6uGDJkCBYtWoRPP/0UVapUwenTp5GcnIzFixfrdQ5kGAyARHoY/wuQxAswJqzsAdALjrCvgC5gxztN4JB40uDHMUl790o+HvDs2bO4cuUK1Gp1sTNsK9qCBQugUqnw+eefo2nTpmjdujW+/fZbbN26FX///Xexz71y5QpefPFFtGzZEvHx8cXuO3XqVHTs2BHNmjXDjh07EBkZCSsrK3Tv3h07duyQ8IxIKpZyF0BUWW0/DczcLHcVZEjlWQg6CPYGqESX52x+fOtl5kygTRsgIkLvprKzs9GnTx+8+uqrCA4OxqBBg3Ds2DG4u7vrX6eeHjx4AKVS+1qPhUX+farVxdwv+fLly3jxxRfRrFkzJCQk6LTxpJMnT2LFihVISkoCAOTl5WkmmOTk5CAvL0/f0yAD4BVAonK4kwVEJQBq3njBtJVjFnAtGP6Xne2jWnD6hPcZ1IsQQL9+wEX9L+GPGzcOGRkZmDt3LkaNGoXatWtj4MCBWvskJSUhKSkJ9+/fx82bN5GUlIR///232HbbtWuH+fPnF7tPSkoKkpKScO3aNTx8+FBznIJxdxEREThw4AAmT56M5ORkHD58GAMGDICfnx+aNm0KANi/fz/q1KmDy5cvA8gPf23atEGNGjUwY8YM3Lx5E9euXcO1a9d0ji+EwFtvvYVZs2bBwcEBABAaGoqvvvoKJ0+exDfffIPQ0NDS/SCpQvFPSKJyiP4OuMSVN0xfOcYA+sPwbwyPbz2gEBV7qzmTlJ6eHwK3bgUU5VvzZ8eOHZg9eza2b98OJycnAEBiYiIaN26MuLg4vPPOOwCgCVtA/tIsK1asgJ+fX7HdxWfPnsWtW8UPJxg0aBB27typ+b7gOKmpqfD390fbtm2xYsUKfP755/j8889hb2+P559/Hps2bYKdnR2A/KuEp0+f1ly127x5M1JSUpCSkgIfHx+t4z09bjA+Ph4eHh7o3LmzZtukSZPQu3dvPPvss+jUqRNiYmKKPQeSh0KUZhQoEWl8sxfot1zuKqgi1Bv9Nf4N2FOm5yzHZTSA4W4EbZXniQaqdCizOPNIMnFxQHS03FUQVSh2AROVQeotIHal3FVQxSn738c+OGGAOv7H/fdghj+pffihrLeKI5IDAyBRKQmRP+4v85HclVBFEWXsFVTBFs4GXANQKarALfqIwdo3W5mZwBNLnBCZAwZAolJK2APs5rArM1O2K4C1UcVAdeRzO/QMLC7eM+gxzNYffwAJCXJXQVRhGACJSuFOFjB6rdxVUMUrWwCsaaAqAEAhrOD+zmkDHoEwciRw86bcVRBVCAZAolIY+zNwM1PuKqiiiTLOAg6E4d4kLudawPqg7jIcJKHbt4H335e7CqIKwQBIVIJD54H4P+WugiqDGrhisLY9P7xhsLbpCd9+m78sDJGJYwAkKoYQwJAVXPDZbJVxEogvDNNF63SzOex+SjZI21SI6GjgEWd7kWljACQqxpLdwP40uasguQgUfausp9nAEq4wzL15PabnGqRdKkJKCvDll3JXQWRQDIBERbidBYzhxA+zVpYLv0FQwQLSBzX7B3Xg9GWS5O1SCT77DCjhLhxElRkDIFERxqwF0rPkroJkVYYu4FqwMEgJHktcDNIuleDePeCTT+SugshgGACJCnEwLb/7l8xd6a8BBuKB5Ee3zvGBy8i/JW+XSmnRIuDcObmrIDIIBkCiQoz6iRM/qGxdwH6Qfpaux6+BUOSUfhwiSSw7Gxg7Vu4qiAyCAZDoKbvOANu43i4BQBnWAfSFtLeJsVC7oNo7ByVtk8ph1SrgwAG5qyCSHAMg0VMm/CZ3BWQ0ShkAlVDAE0clPbTbnkawuCl9tzKVkRDAhx/KXQWR5BgAiZ6w7RSw84zcVZCxEKXsBPaHCtYSjgFUCFu4v31CsvZITzt2AOvXy10FkaQYAImeMJFX/0hL6QJgEGwkPWq1UyGw+pdLkBiVUaOAvDy5qyCSDAMg0X82/wvslnYYF1VyopTLwNTEYwkPqoTH8IvStUfSOHEC+P57uasgkgwDINF/ePWPnlb6LmDprtY5Xw2B7R9pkrVHEvriC7krIJIMAyARgE3Hgb1c7oueUtp1oH2RJtkxPT7m6uNG6+hR4Pff5a6CSBIMgETg1T8qnCjlLGBviWYAO9xrCMf445K0RQYyY4bcFRBJggGQzN7G48D+NLmrIGNUmi5gLzjCXqIuYM8FdpK0Qwa0ZQtw5IjcVRDpjQGQzN6XW+SugIxWKfqAg2AvyaFssv2h+ogLDlcKvApIJoABkMzav1eALSflroKMV8m3YasFaZYG8fjBB4o83n+wUli1CrhwQe4qiPTCAEhmbd52uSsgY1aaOOaPO3ofxzLPDdVi9uvdDlWQ3Fxg9my5qyDSCwMgma27D4Bv9sldBRm1UnQB14D+a/a5b6sHZWa23u1QBfrqK+DuXbmrICo3BkAyW0v/Ah7wdy4VQ5SiC9gH+t2yTSkc4Bb9j15tkAzu3wcWL5a7CqJyYwAksyQEsHiX3FWQsSupC1gFWzjruQag6z/NYHmOV5IqpcWL8z9MiCohBkAyS9tPA8k35K6CjF4JXcBBqKJf+8IC7rFn9WuD5JOaCvz5p9xVEJULAyCZJV79o1IpYSHoWno273KhBWz+uqxnKySr5cvlroCoXBgAyezczAR+5pArKoWSFoIORKZe7XuO1X8GMclszRrgwQO5qyAqMwZAMjtf7wWyc+WugiqH4gNgDVwpd8tVbjeF/YpT5X4+GYnMTOCnn+SugqjMGADJ7HzPmy1QKYkSxgD64nS52/aYWYo1ZqhyYDcwVUIMgGRWzt0EDnMBfyql4rqAbWAJV5TvNjJ2D2tDNfVwecsiY7N9O3BR//UgiSoSAyCZlVWH5K6AKpeiA2AQVLBA+cYSeHztVt6CyBip1UBiotxVEJUJAyCZldUMgFQWxfTS1oJFuZq0yvVC1ff/LmdBZLS+/lruCojKhAGQzAa7f6msiusCDkT5Zn56bKwNxUPOQjI5Z84Ae/fKXQVRqTEAktng1T8qs2LWAfRD2VcStxBOcH2bb0ST9fPPcldAVGoMgGQ2GACprIq7AuiLlDK357q/KSyu3tenJDJm69fLXQFRqTEAkllIvQUcYvcvlVFR8U8JBTxxtExtKYQ13N/hun8m7cQJ4Px5uasgKhUGQDILqw7KXQFVSkV0AftDBesyjgGsmtIC1keuS1EVGTNeBaRKggGQzAK7f6l8Cg+AQbApYzMKeIy8KkE9ZPQYAKmSYAAkk3fjHrt/qXyK6gKuicdlasfpZjPY/XpW/4LI+G3bxnsDU6XAAEgmb8cZuSugSqvILuBbZWrGc0qOFNVQZfDoUX4IJDJyDIBk8raX/3atZOZEEQHQF2mlbsP+fj1UmfuPRBVRpcBuYKoEGADJ5DEAUnkVtQyMdxlmAHvGO0lVDlUWDIBUCTAAkkm7mgGc5sRLkpAXHGFfyi5g65wacB6938AVkdG5eBE4dkzuKoiKxQBIJm0Hr/6RHgrrAg6Cfamf77HWH4octZQlUWWxc6fcFRAViwGQTBq7f0k/ugGwFvJK9UwLdVW4DjkgdUFUWfz9t9wVEBWLAZBMGmcAk9T8cadU+7n/2RDK9IcGroaM1r59cldAVCwGQDJZl+8AyTfkroIqs8K6gGvgYonPUwg7uEUfN0RJVFmkpADp6XJXQVQkBkAyWez+JX0VNgvYGydKfF61f5vD6hR/+Zs9XgUkI8YASCbrzxS5K6BKT6H9rRNs4VLSGoBCCY9h5w1WElUiHAdIRowBkExWUsk9dUTFEtCewRuMKiU+x/lyC9hu470HCbwCSEaNAZBMkhDAiatyV0GV3dMdwDVL8RzPiZmGKIUqo/378z+MiIwQAyCZpHO3gKzHcldBld5TXcCBKD7cOWY0gsOykscIkpnIyABOnZK7CqJCMQCSSTp2We4KyBQ83QVcA1eK3d9jro0hy6HKiN3AZKQYAMkkMQCSFJ7uvPNF0VPLbR8HQjWRCz/TU5KS5K6AqFAMgGSSGABJEk+sA2gDS7jiZJG7enxfHYUsG0jmLjlZ7gqICsUASCaJAZAk8USiC4IKFsgtdDfLPHdUHba/oqqiyoQBkIwUAyCZnMc5vAMISePJMYC1YFHkfu5b6kKZmV0RJVFlk5YG5Bb+hwORnBgAyeScvAbkqUvej6gkT/boBuJBofsohSPc3jpSMQVR5ZObC6Smyl0FkQ4GQDI57P4lydAVgqgAABijSURBVDyxDIwfCr+s7Hq4GSwv3KuggqhSYjcwGSEGQDI5p67JXQGZiie7gH1RyL0FhSU8YvnLnUrAAEhGiAGQTM6F23JXQKaioAtYCQU8cVTn8appLWC9r/i1AYkYAMkYMQCSybl4R+4KyGT81wXsDxWsCxkD6DHqVgUXRJUSAyAZIQZAMjmXGABJIkKR3wUcBN07fFRJfwb2q89UdElUGaUUMnyASGYMgGRShGAAJOmI/zqBa0L3xtKeX1R0NVRpnT/PpWDI6DAAkkm5dR94zM9Zkpg/tLt67R4Ew+mzwzJVQ5VOXh5wm4OTybgwAJJJucbVOEhC4r87gfgiTWu7x/JqMlRDlRoDIBkZBkAyKdcZAElCBV3A3k/MALbO9UbV9/+WqySqrNLT5a6ASAsDIJmUG5lyV0CmxguOsH+iC9j9t5pQPM6TsSKqlHgFkIwMAyCZFF4BJCkJhUAQ7DXfW6hVcH3nkIwVUaXFAEhGhgGQTAqvAJKUBARq4X9X+9z2NYHF9SwZK6JKi13AZGQYAMmk3NVdq5eo/BQC/rjz3//awD36X5kLokqLVwDJyDAAkknhEjAkJQGBGrgIAKh6JgRWx27KXBFVWrwCSEaGAZBMSjYDIElIAPDGv4BQwPO9y3KXQ5UZrwCSkWEAJJPCK4AkJZXCBi5Ihep6CGw3pspdDlVmDIBkZBgAyaQwAJKUguAIAPD49JHMlVCld/++3BUQaWEAJJPCAEhSqgkBh8wGqLLgaMk7ExUnj2tHknFhACSTwjGAJKUAxT14xDnIXQaZAgZAMjIMgGRSeAWQpFQz6z6cx+6XuwwyBQyAZGQYAMmkMACSlJquSYciT8hdBpkCBkAyMpZyF0AkJXYBk5QWpH0G5dhpcpdBJqB6FTUmyV0E0RMYAMmk8AogSWnp1WC5SyAT0cAaDIBkVNgFTCbFgu9oIjJC/GwiY8O3JJkURxu5KyAi0mXJ37ZkZPiWJJPCAEhExohXAMnY8C1JJoUBkIiMEQMgGRu+JcmkONrKXQERka4q/GwiI8MASCaFVwCJyBhV4w1lyMgwAJJJYQAkImPEAEjGhgGQTAoDIBEZo2qOcldApI0BkEwKAyARGSNeASRjwwBIJoUBkIiMEa8AkrFhACSTwpl2RGSMqtrLXQGRNgZAMileKrkrICLSxSuAZGwYAMmk+FWTuwIiIl0cA0jGhgGQTEqNqoBCIXcVRETaeAWQjA0DIJkUWyvAvYrcVRAR/Y+lElDZyV0FkTYGQDI5flXlroCI6H/8XeWugEgXAyCZHH+OAyQiI1LbXe4KiHQxAJLJ4UQQIjImtT3kroBIFwMgmRx2ARORMWEAJGPEAEgmh1cAiciYsAuYjBEDIJkcBkAiMia8AkjGiAGQTA4ngRCRsbC3Bnxc5K6CSBcDIJmcKrb5C0ITEcmtljsXpyfjxABIJqmJr9wVEBFx/B8ZLwZAMkmNfeSugIiI4//IeDEAkkliACQiY1DXU+4KiArHAEgmiQGQiIzBc4FyV0BUOAZAMkk13XjzdSKSl1uV/EkgRMaIAZBMkkIBNPeTuwoiMmfP8+ofGTEGQDJZLfzlroCIzBkDIBkzBkAyWS0C5K6AiMxZSwZAMmIMgGSyQtgFTEQysVQCIf5yV0FUNAZAMlneLrwjCBHJo4kvYGctdxVERWMAJJPWoa7cFRCROeL4PzJ2DIBk0l6qL3cFRGSOWtaUuwKi4jEAkklrXwdQ8kbsRFTBeAWQjB0DIJk0FweuB0hEFcu/GuBXTe4qiIrHAEgmj93ARFSRwhvIXQFRyRgAyeR1rCd3BURkTiIayl0BUckYAMnkPRcAONnKXQURmQM7K+DFYLmrICoZAyCZPEsLoG0duasgInPQJpjr/1HlwABIZoHdwERUETqz+5cqCQZAMgsvMQASkYEpFMArTeSuguj/27vz6KiqBI/jv8oeEhKjkARCBjKGxGDYTdhkR4OgIgiNhkWWwCAEG5RhcJQWwRbngLjRKooYFYZWaW1HsAVlsW1AWhgiKKIim+yiLGHYE+aPJ5EQgUpSya2q+/2c805I1ePlVxyS+uXdd+9zDwUQVvjXmlLTRNMpAPizzHpS7atMpwDcQwGENfq1MJ0AgD/rydk/+BAKIKxxdwZ3BQFQeXo1NZ0AcB8FENaofRXLMwCoHA0TpPpxplMA7qMAwir9Mk0nAOCPBrUynQAoGwogrHJnMyks2HQKAP4kOFAa0NJ0CqBsKICwSlQ463QB8KzbG0s1q5tOAZQNBRDWYTYwAE8a0tp0AqDsKICwTrd0Kaaa6RQA/EHCVVLW9aZTAGVHAYR1QoKkPs1NpwDgDwa1kgJ5J4UP4r8trJRzo+kEAHydyyUNaWM6BVA+FEBYKaOe1OZa0ykA+LIOKc5tJgFfRAGEtcZ2MZ0AgC8bytk/+DAKIKx1RxOp3jWmUwDwRTHVnHVFAV9FAYS1AgOk+zqZTgHAF41sz6Ly8G0UQFhtaBupepjpFAB8SbUQaQyXkMDHUQBhtahwruMBUDY5N0o1Ik2nACqGAgjr3ddJCnCZTgHAFwQHSuNuMp0CqDgKIKyXVMOZEAIAV9IvU0q82nQKoOIogIBYEgbAlQW4pP/oajoF4BkUQEDSjclS22TTKQB4s55NpeviTacAPIMCCPzisR6mEwDwZg9y9g9+hAII/KJdinRTmukUALzRTWlS87qmUwCeQwEELvDHO0wnAOCNHu5mOgHgWRRA4AIZ9aQejU2nAOBNejR2RggAf0IBBC7yxzuc28QBQHCgNO1O0ykAz+NtDrjI9bWlQa1MpwDgDXI7SPXjTKcAPM917ty5c6ZDAN5mz2Gp/kTp+GnTSQCYck2EtOUx6apqppMAnscZQOA31L5Kup/FoQGrTbqN8gf/RQEELmF8lhQXZToFABPSakkj2plOAVQeCiBwCdXDpKf6mE4BwITpd0pBgaZTAJWHAghcxt2ZUtfrTacAUJWyGkjdGppOAVQuCiBwBS9kS9VCTKcAUBUCA6QnOfMPC1AAgSuoV0N69DbTKQBUhTGdnaWgAH/HMjCAGwqLpIzHpfU/mE4CoLIkx0obJkrhnPGHBTgDCLghMEB6qT93CAH8lcslvdyf8gd78HYGuOmGetLojqZTAKgMw9tKHVJNpwCqDkPAQBn83ympwSRp58+mkwDwlMQY6ctHpKhw00mAqsMZQKAMIkKl57NNpwDgKS6XNOceyh/sQwEEyqh7Q2e4CIDvG9Ve6pJmOgVQ9RgCBsrhxGmpxRPSxt2mkwAor5Q4Kf9hJn7ATpwBBMohPER6cxgLRAO+KjBAen0w5Q/2ogAC5ZRWS5p5l+kUAMpjyu1SiyTTKQBzKIBABQxuI/VvYToFgLLo0Via0NV0CsAsCiBQQS9kO9cSAfB+9WOl1wY7s38Bm1EAgQqKDHOuBwwNMp0EwOVUC5HeGSFFs+QLQAEEPKFJovRkb9MpAFzOywOk9ATTKQDvQAEEPGRUR6l3M9MpAPyW0R2l7EzTKQDvwTqAgAcdPy21myat22k6CYDz2lwrLX9ACg40nQTwHpwBBDyoWoi0MFf6l6tNJwEgSfFR0tv/RvkDLkYBBDwsPlpalCtFhZlOAtgtJMiZoFUr2nQSwPtQAIFKkJ7gnHUI4jsMMMLlkl4bJLVLMZ0E8E68PQGV5OYG0vPZplMAdnr6d9JdGaZTAN6LAghUomFtpX+/2XQKwC4PdpXu62Q6BeDdmAUMVLJz56Q+L0l/+V/TSQD/N7SNNHug6RSA9+MMIFDJXC7pjcHceB6obLc1kmb1N50C8A0UQKAKhIdIH4yWGtUxnQTwT22udWb8BvKuBriFbxWgilwdIX08RmpQy3QSwL+k15beH+X8ogXAPRRAoArVrC4tHSulxJlOAviHutdIH94nxUSYTgL4FgogUMXio6VlY6Vra5pOAvi2+rHSp+OkhBjTSQDfwyxgwJDdh6SOM6TvDphOAviehgnSkt87v1ABKDsKIGDQ3iNSpxnS5n2mkwC+I6OeM+x7NcO+QLlRAAHDDhyVOj0lfbXHdBLA+7VNlhaNlqpzr22gQrgGEDAsNkpa8YDUNNF0EsC73dxA+vD3lD/AEyiAqDQrV65UUFCQmjRpUuFjTZo0SS6Xq9QWEfHrGFBeXl6p58PCrvxOsWLFCjVr1kyhoaFKTk5WXl5eiefnzZunxMRExcTE6P777y/x3Pbt25WSkqKjR49W6PXViJT+Pk7qll6hwwB+644mzlIv1VjqBfAICiAqxeHDhzVw4EB17tzZI8cbN26c9u7dW2Jr0KCB+vTpU2K/qKioEvvs2LHjssfdtm2bunfvro4dOyo/P19jxoxRTk6OFi9eLEk6ePCgcnJyNH36dC1ZskRz587VwoULi//+yJEj9cQTTygqKqrCrzEyTPqfUdK97St8KMCv9MuU3h4uhQSZTgL4Dwog3Pbjjz8qPj5ejz/+ePFjq1atUkhIiJYuXVpi3xEjRig7O1utWrXyyNeOjIxUfHx88bZ//35t2rRJQ4cOLbGfy+UqsV9c3OUX3HvxxReVlJSkJ598UmlpacrNzVXv3r311FNPSZK2bt2q6Oho9e3bVxkZGerYsaO+/vprSdL8+fMVHBysXr16eeQ1Ss5dDJ7Plqb3dm4hB9hubGfp9cFSUKDpJIB/oQDCbTVr1tScOXM0adIkrV27VgUFBRowYIByc3NLnOl79dVXtXXrVj3yyCNuH9vlcpUaer2c2bNnKyUlRW3bti3x+LFjx1S3bl0lJiaqR48e+uqrry57nNWrV6tLly4lHsvKytLq1aslSfXr19fx48e1fv16/fzzz/r888/VqFEjHTp0SBMnTtTMmTPdzlwWD9wkLRguhQdXyuEBrxccKM3qJ834nRTAOxXgcXxboUy6deumYcOGqV+/fhoxYoQiIiI0derU4ue/++47TZgwQXPnzlVQkPvjNampqYqOdm9Br5MnT2revHmlzv6lpqZqzpw5eu+99zR37lwVFRWpdevW2rVr1yWPtW/fvlJnCePi4nT06FGdOHFCMTExeu211zRw4EBlZmZq4MCBysrK0rhx45Sbm6tt27apadOmSk9P14IFC9x+ve7o1Uxa/oAUW92jhwW8Xkw1Z5mX4e1MJwH8F1dUoMymT5+u9PR0vf3221q3bp1CQ0MlSYWFhcrOztajjz6qlJSUMh1z8+bNbu/77rvvqqCgQPfcc0+Jx1u1alViyLl169ZKS0vTrFmzNGXKlDLluVDPnj3Vs2fP4s8/+eQTbdiwQc8995ySk5M1f/58xcfHKzMzU+3atVNsbGy5v9bFWiRJn02Quj3HWoGwQ0qcM9mD2yUClYszgCiz77//Xnv27FFRUZG2b99e/HhBQYHWrl2r3NxcBQUFKSgoSJMnT9YXX3yhoKAgLVu2zCNff/bs2br11luveH1fcHCwmjZtqi1btlxyn/PXE15o//79ioqKUnh4eKn9T506pZEjR2rWrFnasmWLzp49q/bt2ys1NVUpKSlas2ZN+V7UZSTVkFaNlzqlevzQgFfp3lD654OUP6AqUABRJqdPn1b//v3Vt29fTZkyRTk5OTpwwLmXWVRUlDZu3Kj8/PzibcSIEUpNTVV+fr5atGhR4a+/bds2LV++vNTw728pLCzUxo0bVatWrUvu06pVq1ITWD766KNLTl557LHH1LVrVzVr1kyFhYU6e/Zs8XNnzpxRYWGhm6+kbGIipCVjpIndpQAmh8DPuFzO/+33R0nRpX/vAlAJGAJGmTz00EM6cuSInn32WUVGRuqDDz7QkCFDtHDhQgUEBCg9veRCdrGxsQoLCyv1+MWuu+46TZ06tcRQ62+ZM2eOatWqpVtuuaXUc5MnT1bLli2VnJysw4cPa9q0adqxY4dycnKK93nwwQe1e/duvf7665Kc2cozZ87U+PHjNWTIEC1btkxvvfWWFi1aVOr4mzZt0ptvvqn169cXZw4ICNArr7yi+Ph4bd68WRkZGZfNXxGBAdLk26X29aV+c6T9FVt6EPAK1cOcWb53VHy5UABlQAGE21asWKGnn35ay5cvL1737o033lDjxo31wgsv6N577y33sb/55hsdOXLksvsUFRUpLy9PgwYNUmBg6TUhDh06pGHDhmnfvn2KiYlR8+bNtWrVKjVo0KB4n71792rnzp3FnyclJWnRokUaO3asnnnmGdWpU0ezZ89WVlZWiWOfO3dOw4cP14wZM4oXnw4PD1deXp5GjRqlU6dOaebMmUpISCj3v4G7OqdJX0yU+r0iLXX/0knA69xQV5o3lCFfwATuBQz4qKIi6fG/SZMWSoVFptMA7gsMkCZkSY/c5iz3AqDqUQABH/f3b6XsV6Tdh00nAa4sqYb0xmCpTbLpJIDdKICAHzh4TBr4qvS3L00nAS5tYEvpubukKCZ6AMZRAAE/ce6c9KcV0n/+VSo4aToN8KurI6QX+0l9mptOAuA8CiDgZ3Ydkkb+t/T+BtNJAKlLmpR3j5QQYzoJgAtRAAE/tWCdNPrP0j6Wi4EB0eHSlNul3I7OOn8AvAsFEPBjh49L4/8izV7pDBEDVWFAS2nanVJclOkkAC6FAghY4NPvpOFzuZ8wKlejOtKf7pZuZIYv4PUogIAlTp1x1g18YrF0+uyV9wfcFR3u3KVmVAdnjT8A3o8CCFjm2/3ShHekd/NNJ4Gvc7mcpV3+qxfDvYCvoQAClvrHFmncAmnNNtNJ4IuaJDrDva2vNZ0EQHlQAAHLvbXWWTvw+x9NJ4EvSI2T/nCrdNcNUgDDvYDPogAC0NlCac5KacoHzjqCwMVS4qQ/dJfuzqD4Af6AAgig2Kkz0gufSFM/lA4UmE4Db1A/VprYXcrOZIIH4E8ogABKOXZSevkf0rPLpO0/mU4DE5JjpYndpH4tKH6AP6IAArikwiLpnfXSjI+kz5gsYoXkWOnhblJ/ih/g1yiAANzy2VZpxsdOISwsMp0GntYlTRrdUbq1Idf4ATagAAIokx0/Sc8slV5ZKR09aToNKiIi1FnHL7eD1KC26TQAqhIFEEC5HD3hlMCXPuUWc76mcR1peFvn+r7ocNNpAJhAAQRQYet2SHPXSH/+XNp31HQa/JaIUGftvuFtpcwk02kAmEYBBOAxhUXSx187ZfCv+dKxU6YT2S08WOrWUOrdzLm2LzLMdCIA3oICCKBSHD/tlMC5a6SPNklnmThSJSJCpW7pUp/mzseIUNOJAHgjCiCASnfgqPRuvrT4K2npZiaPeFpkqHRrI6lPM+mWdCk8xHQiAN6OAgigSp0tlFZvdcrg4k3Sup0SP4XKLjVO6pAqdb3e2cKCTScC4EsogACMOnjMGSJevElasknae8R0Iu90vvB1SHG2+GjTiQD4MgogAK+yYZe06ntpzTbpn9udJWaKLPwpReEDUJkogAC82tET0todzpb/g5S/S/p2v//cjSQyVLq+tpR+fkuQGiVIsVGmkwHwZxRAAD7nxGlp427pyz3StoPS9p+cbcdP0u7D3nnGMCxYSom7oOjVlhomSHWvkVwu0+kA2IYCCMCvnCmUfvj510J4vhzuPSIVnJQKTv3y8ZetosvTRIRKsdV/3WpWl+KjpISrpMQYqc4vW83qnnl9AOAJFEAAVjt1pnQpPFMohQRJIYHOx9Cgkp9f+OfAANOvAADKjgIIAABgGX53BQAAsAwFEAAAwDIUQAAAAMtQAAEAACxDAQQAALAMBRAAAMAyFEAAAADLUAABAAAsQwEEAACwDAUQAADAMhRAAAAAy1AAAQAALEMBBAAAsAwFEAAAwDIUQAAAAMtQAAEAACxDAQQAALAMBRAAAMAyFEAAAADLUAABAAAsQwEEAACwDAUQAADAMhRAAAAAy1AAAQAALEMBBAAAsAwFEAAAwDIUQAAAAMtQAAEAACxDAQQAALAMBRAAAMAyFEAAAADLUAABAAAsQwEEAACwDAUQAADAMhRAAAAAy1AAAQAALEMBBAAAsAwFEAAAwDIUQAAAAMtQAAEAACxDAQQAALAMBRAAAMAyFEAAAADLUAABAAAsQwEEAACwDAUQAADAMhRAAAAAy1AAAQAALEMBBAAAsAwFEAAAwDIUQAAAAMtQAAEAACxDAQQAALAMBRAAAMAyFEAAAADLUAABAAAsQwEEAACwDAUQAADAMhRAAAAAy1AAAQAALEMBBAAAsAwFEAAAwDIUQAAAAMtQAAEAACxDAQQAALAMBRAAAMAyFEAAAADLUAABAAAsQwEEAACwDAUQAADAMhRAAAAAy1AAAQAALEMBBAAAsAwFEAAAwDIUQAAAAMtQAAEAACxDAQQAALAMBRAAAMAyFEAAAADLUAABAAAsQwEEAACwDAUQAADAMhRAAAAAy1AAAQAALEMBBAAAsAwFEAAAwDIUQAAAAMtQAAEAACzz/5QgD9ApLyjkAAAAAElFTkSuQmCC\n", "text/plain": [ "class=Graph name=Unnamed implementation=class=GraphImplementation name=Unnamed title=Importance Factors from Taylor expansions - y0 xTitle= yTitle= axes=OFF grid=OFF legendposition= legendFontSize=1 drawables=[class=Drawable name=Unnamed implementation=class=Pie name=Unnamed labels=[x1 : 18.2%,x2 : 4.3%,x3 : 2.5%,x4 : 75.0%] radius=1 center=class=Point name=Unnamed dimension=2 values=[0,0] color palette=[#ff0000,#ccff00,#00ff66,#0066ff] derived from class=DrawableImplementation name=Unnamed legend= data=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=4 dimension=1 data=[[0.181718],[0.0430356],[0.0248297],[0.750417]] color=blue fillStyle=solid lineStyle=solid pointStyle=plus lineWidth=1]" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# draw importance factors\n", "taylor.drawImportanceFactors()" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/html": [ "[0.932544]
" ], "text/plain": [ "class=Point name=Unnamed dimension=1 values=[0.932544]" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Get the value of the output at the mean point\n", "taylor.getValueAtMean()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/html": [ "[[ -0.35812 ]
\n",
" [ -0.0912837 ]
\n",
" [ -0.0286496 ]
\n",
" [ -0.228209 ]]
sheet #0
\n",
"[[ -1.51555 0.0350552 0.0220043 0.0876379 ]
\n",
" [ 0.0350552 -0.605625 0.00280441 -0.373018 ]
\n",
" [ 0.0220043 0.00280441 -0.427983 0.00701103 ]
\n",
" [ 0.0876379 -0.373018 0.00701103 -0.932544 ]]