ERpy.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    }
   ],
   "source": [
    "%pylab inline\n",
    "import sys\n",
    "sys.path.append(\"/Users/becqg/Documents/code/Robinson\")\n",
    "#sys.path.append(\"C:\\Users\\Guillaume\\Documents\\projets\\code\\Robinson\")\n",
    "import ER\n",
    "import ERpy\n",
    "import imp"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<h2> getTraceMode, getMultiplexedMode, TMTOMM, MMTOTM </h2>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[  1.  11.  21.   2.  12.  22.   3.  13.  23.   4.  14.  24.   5.  15.  25.]\n",
      "[  1.   2.   3.   4.   5.  11.  12.  13.  14.  15.  21.  22.  23.  24.  25.]\n",
      "[  1.  11.  21.   2.  12.  22.   3.  13.  23.   4.  14.  24.   5.  15.  25.]\n",
      "[  1.  11.  21.   2.  12.  22.   3.  13.  23.   4.  14.  24.   5.  15.  25.]\n",
      "[  1.   2.   3.   4.   5.  11.  12.  13.  14.  15.  21.  22.  23.  24.  25.]\n",
      "[  1.  11.  21.   2.  12.  22.   3.  13.  23.   4.  14.  24.   5.  15.  25.]\n"
     ]
    }
   ],
   "source": [
    "imp.reload(ERpy)\n",
    "X = array([1., 11, 21, 2, 12, 22, 3, 13, 23, 4, 14, 24, 5, 15, 25])\n",
    "N = 3\n",
    "LX = 5\n",
    "Y = ERpy.getTraceMode(N, LX, X)\n",
    "Z = ERpy.getMultiplexedMode(N, LX, Y)\n",
    "print(numpy.round(X, 2))\n",
    "print(numpy.round(Y, 2))\n",
    "print(numpy.round(Z, 2))\n",
    "Y = ERpy.MMTOTM(N, LX, X)\n",
    "Z = ERpy.TMTOMM(N, LX, Y)\n",
    "print(numpy.round(X, 2))\n",
    "print(numpy.round(Y, 2))\n",
    "print(numpy.round(Z, 2))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[  1.   2.   3.   4.   5.]\n",
      " [ 11.  12.  13.  14.  15.]\n",
      " [ 21.  22.  23.  24.  25.]]\n",
      "[  1.   2.   3.   4.   5.  11.  12.  13.  14.  15.  21.  22.  23.  24.  25.]\n",
      "(15,)\n",
      "[  1.  11.  21.   2.  12.  22.   3.  13.  23.   4.  14.  24.   5.  15.  25.]\n"
     ]
    }
   ],
   "source": [
    "imp.reload(ERpy)\n",
    "X = array([[1., 2, 3, 4, 5], [11, 12, 13, 14, 15], [21, 22, 23, 24, 25]])\n",
    "print(X)\n",
    "Y = X.flatten()\n",
    "print(Y)\n",
    "print(Y.shape)\n",
    "Z = ERpy.getMultiplexedMode(3, 5, Y)\n",
    "print(Z)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<h2> NDTOMM, NDTOTM, MMTOND, TMTOND </h2>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[  1.   2.   3.   4.   5.]\n",
      " [ 11.  12.  13.  14.  15.]\n",
      " [ 21.  22.  23.  24.  25.]]\n",
      "[  1.   2.   3.   4.   5.  11.  12.  13.  14.  15.  21.  22.  23.  24.  25.]\n",
      "[  1.  11.  21.   2.  12.  22.   3.  13.  23.   4.  14.  24.   5.  15.  25.]\n",
      "[  1.  11.  21.   2.  12.  22.   3.  13.  23.   4.  14.  24.   5.  15.  25.]\n",
      "[[  1.   2.   3.   4.   5.]\n",
      " [ 11.  12.  13.  14.  15.]\n",
      " [ 21.  22.  23.  24.  25.]]\n",
      "[  1.   2.   3.   4.   5.  11.  12.  13.  14.  15.  21.  22.  23.  24.  25.]\n",
      "[  1.   2.   3.   4.   5.  11.  12.  13.  14.  15.  21.  22.  23.  24.  25.]\n",
      "[  1.   2.   3.   4.   5.  11.  12.  13.  14.  15.  21.  22.  23.  24.  25.]\n",
      "[[  1.   2.   3.   4.   5.]\n",
      " [ 11.  12.  13.  14.  15.]\n",
      " [ 21.  22.  23.  24.  25.]]\n",
      "[  1.   2.   3.   4.   5.  11.  12.  13.  14.  15.  21.  22.  23.  24.  25.]\n"
     ]
    }
   ],
   "source": [
    "imp.reload(ERpy)\n",
    "X = array([[1., 2, 3, 4, 5], [11, 12, 13, 14, 15], [21, 22, 23, 24, 25]])\n",
    "print(X)\n",
    "print(numpy.ravel(X, 'K'))\n",
    "(N, LX, X) = ERpy.NDTOMM(X)\n",
    "print(X)\n",
    "print(numpy.ravel(X, 'K'))\n",
    "X = ERpy.MMTOND(N, LX, X)\n",
    "print(X)\n",
    "print(numpy.ravel(X, 'K'))\n",
    "(N, LX, X) = ERpy.NDTOTM(X)\n",
    "print(X)\n",
    "print(numpy.ravel(X, 'K'))\n",
    "X = ERpy.TMTOND(N, LX, X)\n",
    "print(X)\n",
    "print(numpy.ravel(X, 'K'))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "ename": "ValueError",
     "evalue": "cannot slice a 0-d array",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[1;31mValueError\u001b[0m                                Traceback (most recent call last)",
      "\u001b[1;32m<ipython-input-3-dfd4097f1dea>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m     10\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mNZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mZ\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mERpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mNDTOMM\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mZ0\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;36m2\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m:\u001b[0m\u001b[1;36m90\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     11\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 12\u001b[1;33m \u001b[1;33m(\u001b[0m\u001b[0mF\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mY\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mE\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mERpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mWIENER\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mNX\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLX\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mX\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mNZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m10\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m0\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m0\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m     13\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     14\u001b[0m \u001b[1;32mprint\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mY\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32mC:\\Users\\Guillaume\\Documents\\projets\\code\\Robinson\\ERpy.pyc\u001b[0m in \u001b[0;36mWIENER\u001b[1;34m(N, LX, X, M, LZ, Z, LR, LW, FLOOR)\u001b[0m\n\u001b[0;32m     74\u001b[0m     \u001b[0mLS\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mN\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mN\u001b[0m \u001b[1;33m*\u001b[0m \u001b[1;33m(\u001b[0m\u001b[1;36m5\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mLR\u001b[0m \u001b[1;33m+\u001b[0m \u001b[1;36m6\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m+\u001b[0m \u001b[0mM\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mN\u001b[0m \u001b[1;33m*\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mLR\u001b[0m \u001b[1;33m+\u001b[0m \u001b[1;36m2\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m+\u001b[0m \u001b[1;36m2\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mM\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mM\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     75\u001b[0m     \u001b[0mS\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0marray\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mLS\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 76\u001b[1;33m     \u001b[1;33m(\u001b[0m\u001b[0mLF\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mF\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mE\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLY\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mY\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mS\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mER\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mWIENER_1\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mN\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLX\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mX\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mM\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLR\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLW\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mFLOOR\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLF\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mF\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mE\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLY\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mY\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mS\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m     77\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     78\u001b[0m     \u001b[1;32mreturn\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mF\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mY\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mE\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32mC:\\Users\\Guillaume\\Documents\\projets\\code\\Robinson\\ER.pyc\u001b[0m in \u001b[0;36mWIENER_1\u001b[1;34m(N, LX, X, M, LZ, Z, LR, LW, FLOOR, LF, F, E, LY, Y, S)\u001b[0m\n\u001b[0;32m    273\u001b[0m     \u001b[0mIH\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mIGAM\u001b[0m \u001b[1;33m+\u001b[0m \u001b[0mMN\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    274\u001b[0m     \u001b[0mIFGT\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mIH\u001b[0m \u001b[1;33m+\u001b[0m \u001b[0mM\u001b[0m \u001b[1;33m*\u001b[0m \u001b[0mM\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 275\u001b[1;33m     \u001b[0mS\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mIH\u001b[0m \u001b[1;33m-\u001b[0m \u001b[1;36m1\u001b[0m\u001b[1;33m:\u001b[0m \u001b[1;33m]\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mHEAT\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mM\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m1\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mM\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m1\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mZ\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m1\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mS\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mIH\u001b[0m \u001b[1;33m-\u001b[0m \u001b[1;36m1\u001b[0m\u001b[1;33m:\u001b[0m \u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    276\u001b[0m     \u001b[1;31m# print(Z, S[IH: IH + M * M])\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    277\u001b[0m     \u001b[1;32mif\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mLW\u001b[0m \u001b[1;33m<=\u001b[0m \u001b[1;36m1\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;31mValueError\u001b[0m: cannot slice a 0-d array"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEACAYAAABS29YJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xtwk3XaPvCrh6THgGKlrbahWEqSItKiTdIqNcCOy6q1\nOMyOyyzI2O77c8u6KlL+EJkVdhYYWUYExmWrTnfeWQZ3Vmd8d9VVsOtGQdqkcpKtLVDaQotgqQik\n9JSW5/dHSZq0OTXn5Lk+MxlTmsOTR+fi9s79/T4xgiAIICKiqBYb6gMgIqLAY9gTEYkAw56ISAQY\n9kREIsCwJyISAYY9EZEIuAz7gYEBaDQaFBQUQKvVYseOHQAAk8mE8vJyyOVyLF26FL29vdbn7Nq1\nC3l5ecjPz8ehQ4cCe/REROSRGHdz9n19fUhOTsbg4CDuv/9+fPDBB/jggw/Q2dmJ7du3Y+3atcjJ\nyUF1dTW6u7tRWlqKAwcOoL29HWvWrMHRo0eD9VmIiMgJt22c5ORkAEBvby+Gh4eRkJAAo9GIyspK\nJCQkoKKiAgaDAQBgMBiwZMkSyOVyPPzwwxAEASaTKbCfgIiI3HIb9jdv3sS8efOQnp6O5557DnK5\nHI2NjVAqlQAApVIJo9EIYDTsVSqV9bkKhcL6OyIiCp14dw+IjY3FiRMn0NHRgUcffRQPPvggJrPD\nQkxMjE8HSEREvnMb9hY5OTl49NFHYTAYUFRUhObmZhQWFqK5uRlFRUUAAI1Gg7q6OutzWlparL+z\nNWvWLJw9e9YPh09EJB65ublobW316rku2zg9PT24evUqAOCHH37AgQMHUF5eDo1Gg9raWvT396O2\nthZarRYAoFarsX//fpw/fx56vR6xsbGQyWQTXvfs2bMQBIE3QcCrr74a8mMIlxvPBc+FaM7FwACE\n+noIr78O4ec/h5CdDSEtDcLjj0PYvBnC559DMJkmPM+XItllZX/x4kWsWrUKIyMjyMjIQHV1NTIz\nM1FVVYUVK1ZAoVBg/vz5eO211wAA6enpqKqqwqJFiyCVSlFTU+P1gRERRY0LF4D6+rHbN98AeXlA\ncTHw+OPA5s3ArFlAANveLsN+7ty5DkcnZTIZ/vGPfzh8zgsvvIAXXnjBP0dHRBRpBgeBY8fsw31g\nANBqR8N982agqAhITQ3qYXncs6fA0Ol0oT6EsMFzMYbnYkzYn4vxVfuJE8Ds2aPBXlYGbNkC5OYG\ntGr3hNtFVQF505gYhOBtiYh8465qLy4OaNXuS3Yy7ImInOnqGg30hoaxqt3Sa7fcAtxrt8WwJyLy\nlW3Vbgn3vj77YA9Br90Ww56IaLKc9dptWzJBrNo9wbAnInLFUa+9vz+sqnZPMOyJiGxFYNXuCYY9\nEYnX4CBw/Lh9uNv22rVaQK0O+6rdEwx7IhIPV3PtllsYzLUHAsOeiKJTiOfaww3Dnoiig6O5dpFU\n7Z5g2BNR5GHVPmkMeyIKf1E6IRNMDHsiCi9RMtcebhj2RBRarNqDgmFPRMHDXnvIMOyJKHAcVe0h\n3PlRzBj2ROQflqrdMvpo6bWzag8LDHsi8g7n2iMKw56I3GOvPeIx7IloIkvVbqncbSdkSkpYtUcg\nhj2R2I2v2hsawu4qS+Q7hj2R2Lja+dHSluGETNRh2BNFM65GpVsY9kTRRMT7tZNrDHuiSDX+KkuH\nD7NqJ6cY9kSRgnvIkA8Y9kThyLbXblm0xAkZ8oEv2Rnr6pednZ1YuHAh5syZA51Oh3379gEANm7c\niKysLBQWFqKwsBCffPKJ9Tm7du1CXl4e8vPzcejQIa8OiigiXbgAvP8+sHbt6Bz7tGnAr38NnD4N\nPPYY8O9/A5cvAx9+CKxfDyxcyKCnoHFZ2V+6dAmXLl1CQUEBenp6oFarceLECbz++uuQyWR46aWX\n7B7f3d2N0tJSHDhwAO3t7VizZg2OHj068U1Z2VOkczTXzj1kKMB8yc54V7/MyMhARkYGACAtLQ1z\n5sxBY2MjADh8Q4PBgCVLlkAul0Mul0MQBJhMJshkMq8OjihsuJqQefxxYPNm9toprLkMe1utra1o\namqCRqPBwYMHsXv3brz33nt48sknsXr1ashkMhiNRqhUKutzFAoFjEYjFi9eHJCDJwoId3Ptmzez\naqeI47Jnb2EymfDUU09hx44dSElJQVVVFdrb27F//36cPXsWNTU1ABxX+zGsdCjcOeq1V1UBZ84A\nZWXAf/7DXjtFPLeVvdlsxrJly7By5UqUl5cDAKZPnw4AmDp1Kn7zm99g9erVqK6uhkajQV1dnfW5\nLS0tKCoqcvi6GzdutN7X6XTQ6XQ+fAwiD7nb+ZFVO4URvV4PvV7vl9dy+QWtIAhYtWoV0tLS8Prr\nr1v//OLFi8jMzMTw8DBeeeUVTJkyBa+88gq+//57PPzwwzhw4ADa2trw0ksv8QtaCi1eZYmiSMC+\noP3qq6+wd+9e3HfffSgsLAQAbNmyBe+++y6OHz8OqVSK0tJSVFVVAQDS09NRVVWFRYsWQSqVWts7\nREHhbq6dVTuJGBdVUeTialQSGa6gpejHnR+JGPYUhVi1E03AsKfIxqssEXmEYU+Rhfu1E3mFYU/h\ny91cO6t2Io8x7Cl8dHWNjT2yaifyK4Y9hYalareE++HDrNqJAohhT8HBCRmikGLYk/9xrp0o7DDs\nyXes2onCHsOeJsdV1W4Jd7WaVTtRmGHYk2vu5tq1WlbtRBGAYU9jxk/IWKp2TsgQRTyGvZhxrp1I\nNBj2YsHVqESixrCPVl1d9puD2U7IlJSwaicSGYZ9NODOj0TkBsM+ErmakLG0ZTghQ0Q2GPbhjqtR\nicgPGPbhxlHVnpdnH+6s2olokhj2oTQ4CBw/Phbshw9PrNofeACQyUJ9pEQU4Rj2wcQ9ZIgoRBj2\ngcLVqEQURhj2/uJsrp2rUYkoDDDsvcHVqEQUYRj2nuBcOxFFOIb9eJxrJ6IoxLDnfu1EJALiCnt3\nvXatlldZIqKo5Et2xrr6ZWdnJxYuXIg5c+ZAp9Nh3759AACTyYTy8nLI5XIsXboUvb291ufs2rUL\neXl5yM/Px6FDh7w6KDsXLgDvvw+sXTu60+O0acCvfw2cPg08/jjw+edAdzfw4YfA+vXAokUMeiKi\ncVxW9pcuXcKlS5dQUFCAnp4eqNVqnDhxAnv27EFnZye2b9+OtWvXIicnB9XV1eju7kZpaSkOHDiA\n9vZ2rFmzBkePHp34ps7+dnI0186dH4mIAPhW2ce7+mVGRgYyMjIAAGlpaZgzZw4aGxthNBqxYcMG\nJCQkoKKiAlu3bgUAGAwGLFmyBHK5HHK5HIIgwGQyQeZsqwBXq1Efewz4wx/Yayci8gOXYW+rtbUV\nTU1NUKvVeOaZZ6BUKgEASqUSRqMRwGjYq1Qq63MUCgWMRiMWL1488QXlcvvVqJs3s2onIgoQj8Le\nZDLhqaeewo4dO5Camjqp/42IcVKVb1y6FLj9diAmBrqSEuh0Oo9fk4goGgmCgNM/nEZDVwPqu+rR\neaITReYiv7y227A3m81YtmwZVq5cifLycgBAUVERmpubUVhYiObmZhQVjR6MRqNBXV2d9bktLS3W\n3423cdcufxw/EVHEuj54HcYLRmu4N3Q1IFWaiuKsYhRnFeN/nvkf3H/X/dbHb9q0yev3chn2giCg\nsrIS9957L1588UXrn2s0GtTW1mLbtm2ora2FVqsFAKjVaqxbtw7nz59HW1sbYmNjnffriYhExFK1\nW0L9cOdhnP3xLAozClGcVYzKwkq8XfY27pLdFZD3dzmNc+jQIZSWluK+++6ztmO2bt2KBx98ECtW\nrMCxY8cwf/587N27F6m3eu07d+7E7t27IZVKUVNTgwULFkx803DYG4eIKIBMgyYYLxhR31VvDXiZ\nVIbi7GJo79aiOLsYBRkFkMZJPX5NcS2qIiIKM7ZVe31nPRouNKD1Squ1ai/OLoY2S+tz1c6wJyIK\nIldVe3HWaLBPtmr3BMOeiChAxlft9V31dr12f1XtnmDYExH5iaOq3XZCxpteu78w7ImIvDB+rr2+\nqx5nr5xFQUZB0Kt2TzDsiYg84KzXrs3Shrxq9wTDnohonPFz7eFetXuCYU9EoudoNaqlai/JLkFx\nVjHmZcwL26rdEwx7IhIVT6r24qxiZMoyQ32ofsWwJ6KoFqq59nDDsCeiqOHJXHs0Vu2eYNgTUcQK\n57n2cMOwJ6KIYDvXfrjzcEhXo0Yihj0RhaVIn2sPNwx7Igo5RxMygdj5UcwY9kQUdIHYr51cY9gT\nUUAJgoAzV85Yp2Msc+2FmYXWYBfrhEwwMeyJyK9sq/aGroaxCRmRzbWHG4Y9EXnN0VWWxLAaNRIx\n7InIY9cHr6PxQqPT1ajRsIdMtGLYE5FD4XSVJfIdw56IAHA1arRj2BOJUDTu106uMeyJRICrUYlh\nTxRlPJmQYdUuPgx7ogjHXjt5gmFPFEGcrUa1VO3arNEVqazaaTyGPVEY41w7+QvDnihM8CpLFEgB\nDfuKigp8/PHHmD59Ok6ePAkA2LhxI9555x3ceeedAIAtW7bgZz/7GQBg165d2L17NyQSCd566y08\n9NBDfj1gonByffA6jBeM1tFH2157SXYJ95Ahvwpo2B88eBCpqal4+umnrWG/adMmyGQyvPTSS3aP\n7e7uRmlpKQ4cOID29nasWbMGR48e9esBE4XK+Ln2w52HuRqVgsqX7Ix394AFCxago6Njwp87ekOD\nwYAlS5ZALpdDLpdDEASYTCbIZDKvDo4olNzNtVcUVrBqp4jhNuyd2b17N9577z08+eSTWL16NWQy\nGYxGI1QqlfUxCoUCRqMRixcv9svBEgWKo7l226ssVRZW4u2yt1m1U8TyKuyrqqrwu9/9DtevX8e6\ndetQU1OD6upqh9V+TEyMw9fYuHGj9b5Op4NOp/PmUIi84uoqS8VZxaicX8mqPYrdNN/E0HdDGOwa\nxEDnAAY7BzHYNTj2z65BzK6ZjbSytJAep16vh16v98treTSN09HRgbKyMmvP3taJEyewevVqfPXV\nV/jwww9RV1eHnTt3AgAKCgpw8ODBCW0c9uwpmNzNtbPXHl0sQT7QOWAN7vFhbu4xQ5ouRUJWAhKy\nE0b/aXs/OwHSDCli42ND/XHsBLRn78jFixeRmZmJ4eFh7Nu3D48++igAQK1WY926dTh//jza2toQ\nGxvLfj0FnbvVqOy1R64JFbmTIJdMlyAxO9Ea3ok5iZj60NSwDvJAcxv2y5cvxxdffIGenh5kZ2dj\n06ZN0Ov1OH78OKRSKUpLS1FVVQUASE9PR1VVFRYtWgSpVIqampqAfwASN0uv3TL6OL5qZ689cnjS\nWplQkWfbBPmtYBdjkHuCi6ooojiaa+dq1PDnsrVy636ktlaCiStoKSrxKkuRwWlrxaYyN/8wsbVi\nqcwtoS72IPcEw56igrOq3TLXXpJdwqo9yCZU5OOqcbvWSnYCEu4eF+JsrfgVw54iDq+yFHo+t1YY\n5EHHsKew526unXvI+Je3UytsrYQ3hj2FFU+ussSdH73naZA7ba3wy86IxbCnkBpftRu6DEiRpvAq\nS15w2iN3NH44ri9uuZ+YnQhJuoRBHoUY9hQ04+faufOj5yb9ZWdWAscPyQ7DngLG3c6PrNpH+dQj\nZ2uFPMSwJ79wNCFju/OjWKv2SffInSwKYpCTrxj25BV3e8iIoWp32VoZP37orLWSlQBpJoOcAo9h\nT26Jca7dq4qcrRUKYwx7msDVXLv2bm3EV+2TCnIHC4HYWqFIxLAXuWiba/d290O7qvxutlYo+jDs\nRcbdatRw3vnRNsidVeUTWisOFgUxyEmMGPZRzNHOj20/toVlr93nlZ1srRC5xLCPIuE6IeMyyF1t\nmmXzT67sJPINwz5ChcuEjNv9yLlpFlFYYNhHiFCsRvV6iT7HD4nCDsM+DDmbkJmXMQ8lWSV+qdq5\nHzmRuDDsw4C/59o9ba04rcgZ5ERRh2EfZL7OtbusyLkfORE5wbAPsMnMtbsdP7x18WVOrRDRZDHs\n/cjRXLtlv/aSjBKUSEsw9+ZcpF5J9fzCEraVOVd2EpGXGPY+uD54HcYLRhg6DPjvyf/ifMt53H3j\nbhSiEHmDeUg3pSPpchKGLgy5niNna4WIAoxh74Ztj7yjuQNtzW243HoZ/Z39SLqchMzeTKTeSIWQ\nJiAlOwWyGTKu7CSisCPqsHc1R953vg83zt+AcEVA75ReXEy9iGu3X0NCVgLSctMwM38mFPcqkCJP\nYWuFiMJe1Ib9ZKZWpFlSmKeb8b3se5xNPIsTMSdwMu4k0vPSMTd/LrQ52rDZQ4aIyBsRGfYjQyOe\nb2ProD8+kj6Cb+O+Rf1gPeov8dqoRBT9Ahr2FRUV+PjjjzF9+nScPHkSAGAymbBixQocO3YM8+fP\nx969e5GamgoA2LVrF3bv3g2JRIK33noLDz30kMMD1kv0zld2jtvGVhAEnLlyxjodM34PGW3W6KIl\nVu1EFM0CGvYHDx5Eamoqnn76aWvYb9u2DZ2dndi+fTvWrl2LnJwcVFdXo7u7G6WlpThw4ADa29ux\nZs0aHD161OEBj5hHnPbI3c21a7O0rNqJSHR8Cft4dw9YsGABOjo67P7MaDRiw4YNSEhIQEVFBbZu\n3QoAMBgMWLJkCeRyOeRyOQRBgMlkgkwmm/C6lqB3V7VXFlbinbJ3IuYqS0RE4cht2DvS2NgIpVIJ\nAFAqlTAajQBGw16lUlkfp1AoYDQasXjx4gmv8Ycv/4D6rnoYugyj+7XfqtorCyvD9ipLRESRyquw\nn8z/RsTExDj883+9/S9kT8nGqimrUPbTMuh0Om8OhYgo6pjNwMmTQF2dHn19er+8pldhX1RUhObm\nZhQWFqK5uRlFRUUAAI1Gg7q6OuvjWlparL8b7/BfD3vz1kREUae7G6ivH7sdOQLMmAGUl+uwZYvO\n+rhNmzZ5/R5ehb1Go0FtbS22bduG2tpaaLVaAIBarca6detw/vx5tLW1ITY21mG/nohIrCxVu224\n//ADoNEAxcXA+vWj92+7zb/v6zbsly9fji+++AI//PADsrOz8fvf/x5VVVVYsWIFFAoF5s+fj9de\new0AkJ6ejqqqKixatAhSqRQ1NTX+PVoiogjT3Q00NIwF+9dfj1btxcXAwoXAyy8DKhUQG+AF/GG9\ngpaIKJIMDwPffOO8ai8u9q1qj8gVtAx7Iop046v2I0cAuXws2IuLAaXSf1U7w56IKMACXbV7gmFP\nRORnthMyDQ1jVbtWOxbuwei122LYExH5YHzV3tAwVrVbwj3QVbsnGPZERJPgbK49UL12f2HYExE5\n4Wiu/cqV8KvaPcGwJyK6ZXzVfvTo2ISMVguUlIRn1e4Jhj0RiVI4TMgEE8OeiEQh2HPt4YZhT0RR\nx1HVHqm9dn9h2BNRxIvUCZlgYtgTUUQRW6/dXxj2RBTWXPXaLS2ZYK9GjUQMeyIKG6567ZZwZ9Xu\nHYY9EYWM2CdkgolhT0RB4ewqS7abg6nVrNoDhWFPRAHh6ipLrNqDj2FPRD7jXPsos/lH9PW1IClp\nFqTSO0N9OHYY9kQ0ac72axdD1S4IIxgYOI++vpYJt5s3+5CcrMSsWW9g6tQHQ32odhj2ROSSs6pd\nrY7uufbh4V7095+6FeSnrIHe398KiSQNyckKJCerkJystN6XSjMRExMT6kN3iGFPRHYsvfbDhx3v\n/BhNc+2CIGBw8IJNqI/dzOYfkJSUdyvMLYGuRFLSbMTHp4b60CeNYU8kYmJZjToyMoD+/tYJgd7f\nfwqxsSk2gT5WpScmyhETEwV/o93CsCcSEXdz7Vpt5FbtgiDAbL5sE+Zj1frg4AUkJuYgJUV1qzpX\nICVFhaQkBSSSCP+bzEMMe6IoFa1V+82bZgwMtE0I9L6+FgC41UdX2PXUExNnIjZWEuIjDy2GPVGU\niLbVqKNjjBN76QMDHUhIyLL20G1DXSJJC+kXpP0jI/hxeBhT4+OREhcXsuNwhGFPFIHcTciUlETG\natTRMcZz40L91K0xxhvWPnpSksLm/izExSUG7Jgsgf3j8DB+NJvH7nvw801BwDSJBP+rVOKRadMC\ndozeYNgTRQDLXLulco+0qn10jPG0gy9ILWOMSptKffQmld7ldZU+cCuwr9gE8hUPw1sAcHt8PKZJ\nJLg9Pn70Znvfxc9JsbEcvfQXhj1FO2e9dq02vFejWsYYbSddJo4x2s6mK5GcPBtxcSkOX8/XCttZ\nKE+L4MD2RcjCPicnB1OmTEFcXBwkEgmMRiNMJhNWrFiBY8eOYf78+di7dy9SU+3nWRn2FG3GX2XJ\ndq49HKv20THGMxP66fZjjKOhHp+YhyFpLkxxd+Hq8E2XAX3F9mebCtsSwtOipMIOlZCF/cyZM3Hk\nyBFMs+lrbdu2DZ2dndi+fTvWrl2LnJwcVFdX++2AiULN2c6P4TYhYzvGeK23GVdvNONGXwuG+k9B\nMF/EsCQbfZJ7cDXuHlyOzcF3MTNwTsjCpZvJdoENwOOAZmAHli/ZGe/rm49/Y6PRiA0bNiAhIQEV\nFRXYunWrr29BFFKudn5cuBB4+eXgzLU7bYmY+zHQ3wZh4DTihlqRZD6LKcPtuGOkAwIEdCIb5zED\nPbE5uB7/E/RL/h+ElBmYKkmyC+n74uOhG9fjnhYfj6Qwm0gh7/hU2d9zzz2QyWSYOXMmKioq8MQT\nT2DGjBk4deoUEhMT0dfXB5VKhXPnztm/KSt7ClOBnmv3pIft7EvIFMEEVdwF5MVdQA46cZdwDnfe\nPAfZSBcG4zNhluYiJiEPkiQFUpJn47aUfExLysQ0iYQVdpQIWWX/1VdfITMzE83NzSgrK4NarWaI\nU0RxNde+cCGwfv3EXnv/yAi+G5z8l45XPJgSmZkggTrhR9w+0oYUcxsShloRM3gGZuE0hJt9duOL\nycmPIjlZFfAxRooOPoV9ZmYmAEClUuGJJ57Ahx9+iKKiIjQ3N6OwsBDNzc0oKipy+NyNGzda7+t0\nOuh0Ol8Ohcgt26r9oHEE9f8dxpWhYajUZuQWDKPolWEsyRnGQPxoQH8xPIz/Gx7Gj8cdB7azPnVO\nYiLmu+lhj+3G+N+xL0ivnkJ//xlIJHeMTbvc/gCSk3/p8xgjRSa9Xg+9Xu+X1/K6jdPX14eRkRHI\nZDJcvnwZOp0On376Kd599110dnZi27ZtqK6uxsyZM/kFLQXEgG1LxMEkSJdpGGcvDaPzqhmXB4Zh\nwjDipg5DSDUjJhaYGhePOxMkmCYJzJeO7ndjnDVue93I3Y2RgiMk0zjt7e148sknAQB33HEHfvnL\nX6KiooKjlzQp/prDvi0uHnH98Rjqice1CxJcbotHf3c8ZqVLcO+MeDygiEfx3HjMmOb/KRHXuzEm\nTwj00X1eZkTVbowUHFxURSE1MDIy+qWiBz1rZysdJzuHbf4xHscNsWiojwnKXLu73RiTku6xqc7F\ntxsjBQfDnnzmy5QIENg5bEdz7YG6NurYbowTWy8AJuyXzt0YKZgY9gTA8wp7skvTbX+2nSKx3Pf3\nHHYwrrLkaDfG/v5T6O9vv7Ubo22gK8JiN0Yihn0Ucfelo6sK+yYwoQ0yLcxXOgZyrt2Ti0rbjzIG\nfjdGIl8w7MOML186CoIQ1UvTA7Ffu/uLStu3XZKTFRxjpIjEsA8AbwPbkzlsR19CBqolEkqOqnZv\ne+22uzGOH2V0vBujgmOMFHUY9k442g/bNqBd9bcFwOMWSCRW2IHgj/3a3e/GqBg3yqi4NcYYPX9J\nEjkT1WEfqP2w3Y35iTWwPWVbtVvCvafHs177xDHGsfbL6BjjTOv4om1fnWOMJHYRGfb/e/GiV3PY\nDr94ZGAHnKteu7MJGfuLStv30wHLGKP9haU5xkjkXESG/Ypvv2VLJEy5mmu3VO2210Y1m686nHgZ\nvaj03XbTLpZQ5xgj0eRFZNiH+xe0YjL+KktHjozt1265KRQjGBpyflFpR20XjjES+RfDnjzm/ipL\nvZg79zTi4sbv83Lm1hjjxH1eOMZIFBwMe3LK8VWWBCxa9B2Ki1swe3YLpkxpQX//+DFGpV0vnWOM\nRKHHsCcAE+faGxsHkJh4BjpdCwoKTkEub0FqaguGhhztxqi6NcYo5xgjUZhi2IvUaK9dwJEjl3Hu\nXAt6e08hP78FSmULpk9vgURyAcnJ99xaYDS2E2NysgISye2hPnwimiSGvQgMDZlx4kQbmptbcOlS\nC4aGWpCWdgozZrQgPh6Ii1MiLU2B225TYWzPdI4xEkUThn0UsezG+P33LWhvb8HVq6cQG9uCqVPb\nce1aFgYHlUhJUSA7W4nc3NH7Esmd/IKUSAQY9hFm/G6MN2604PLlFgwOju7GeOGCEufOKRATo0Ra\nmhJ5eQo88EAepk3jGCORmDHsw9ToboynHSw4asXISBquXFHg7Fkljh5VwmxWIitLgXnz7oZWG+Pz\nfu1EFH0Y9iHk+KLSp26NMfYgMTEPQ0MKfPedCk1NChw8qEBzswLz5qU6XI1KROQMwz4IXF9Uemw3\nxps3FWhrU+Drr1X4z3/k+PrrOLvVqFqt71dZIiJxYtj7ifuLSs+0u0ydVKpAe7sSDQ232e0ho1b7\nfpUlIqLxGPaTZL8bo/2+6QAcbtyVmDgTPT0Sn/drJyLyFsPeCfe7MU7c58WyG6Oza6NqtZO/yhIR\nkT+IOuwdX1R6bDdG+4tKW/Z5mbgb4/idH48eZdVOROFFFGHveIzx1K3dGO+Y1G6M7nd+ZNVOROEn\nasJeEAQMDX3nsPUycTdGlccXlXa88yMnZIgoskRk2JtMJ12MMdoGunJSuzE667WzaieiSBeRYW8w\nKK39dNsrHE32otKuro3KXjsRRZOwC/svv/wSzz77LIaHh/H888/jt7/9rf2bennAjqr28XPtajVw\nO3fvJaIoFHZhX1hYiJ07d2LGjBn46U9/ikOHDiEtLW3sTT08YMuETDTPtev1euh0ulAfRljguRjD\nczGG52KML2Hv95i8du0aAKC0tBQzZszAI488AoPB4PZ5w8Oj445vvgmsWAHk5gIKBbBnD5CQALz8\nMtDZCTQas8C4AAAFYElEQVQ1Ae+8A1RWAvn5kR30wOh/yDSK52IMz8UYngv/iPf3CzY2NkKpVFp/\nzs/PR0NDAx577DG7x9nOtTc02FftCxcC69dHftVORBQu/B72nlIoRqditNrRqp0TMkREAST42dWr\nV4WCggLrz88995zw0Ucf2T0mNzdXAMAbb7zxxtskbrm5uV5ns98r+6lTpwIYnciRy+X47LPP8Oqr\nr9o9prW11d9vS0RELgSkjfPGG2/g2WefhdlsxvPPP283iUNERMEXkkVVREQUXEGddfnyyy+hUqmQ\nl5eH3bt3B/OtQ66zsxMLFy7EnDlzoNPpsG/fPgCAyWRCeXk55HI5li5dit7e3hAfafCMjIygsLAQ\nZWVlAMR7Lm7cuIFVq1Zh9uzZyM/Ph8FgEO25ePvtt1FSUoL7778fL774IgDx/HdRUVGB9PR0zJ07\n1/pnrj77rl27kJeXh/z8fBw6dMjt6wc17F944QXU1NSgrq4Ob775Jnp6eoL59iElkUiwY8cONDU1\n4f3338eGDRtgMpmwZ88eyOVynDlzBllZWfjzn/8c6kMNmp07dyI/P9+6M6lYz8Wrr74KuVyOb775\nBt988w2USqUoz8WVK1ewZcsWfPbZZ2hsbMTp06exf/9+0ZyLZ555Bp9++qndnzn77N3d3fjTn/6E\nf//739izZw+ef/55t68ftLD3drFVtMjIyEBBQQEAIC0tDXPmzEFjYyOMRiMqKyuRkJCAiooK0ZyT\nrq4u/Otf/8KvfvUr64pAsZ6Luro6rF+/HomJiYiPj8fUqVNFeS6SkpIgCAKuXbuG/v5+9PX14bbb\nbhPNuViwYAFuH7fXi7PPbjAYsGTJEsjlcjz88MMQBAEmk8nl6wct7J0tthKj1tZWNDU1Qa1W250X\npVIJo9EY4qMLjjVr1uCPf/wjYm1WzYnxXHR1dWFgYABVVVXQaDR47bXX0N/fL8pzkZSUhD179iAn\nJwcZGRl48MEHodFoRHkuLJx9doPBAJVKZX2cQqFwe164PjXITCYTnnrqKezYsQOpqalheeH1QPvo\no48wffp0FBYW2n1+MZ6LgYEBnD59GsuWLYNer0dTUxP+/ve/i/JcXL58GVVVVfj222/R0dGB+vp6\nfPTRR6I8FxaT+eyOLtRkK2hhX1RUhJaWFuvPTU1N0Gq1wXr7sGA2m7Fs2TKsXLkS5eXlAEbPS3Nz\nMwCgubkZRUVFoTzEoDh8+DD++c9/YubMmVi+fDk+//xzrFy5UpTnYtasWVAoFCgrK0NSUhKWL1+O\nTz/9VJTnwmg0QqvVYtasWbjjjjvw85//HAcPHhTlubBw9tk1Gg2+/fZb6+NaWlrcnpeghb3tYquO\njg589tln0Gg0wXr7kBMEAZWVlbj33nutUwbA6L+02tpa9Pf3o7a2VhR/AW7ZsgWdnZ1ob2/H3/72\nNyxatAh//etfRXkuACAvLw8GgwE3b97Exx9/jJ/85CeiPBcLFizA119/jStXrmBwcBCffPIJHnnk\nEVGeCwtnn12tVmP//v04f/489Ho9YmNjIZPJXL+Y12tvvaDX6wWlUink5uYKO3fuDOZbh9zBgweF\nmJgYYd68eUJBQYFQUFAgfPLJJ8L169eFJ554QsjOzhbKy8sFk8kU6kMNKr1eL5SVlQmCIIj2XJw6\ndUrQaDTCvHnzhLVr1wq9vb2iPRd/+ctfhNLSUuGBBx4QNmzYIIyMjIjmXPziF78QMjMzBalUKmRl\nZQm1tbUuP/sbb7wh5ObmCiqVSvjyyy/dvj4XVRERiQC/oCUiEgGGPRGRCDDsiYhEgGFPRCQCDHsi\nIhFg2BMRiQDDnohIBBj2REQi8P8BISEiomnU6eoAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x78217f0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "imp.reload(ERpy)\n",
    "X0 = arange(0, 300).reshape(3, -1)\n",
    "Z0 = empty((3, 90))\n",
    "plot(X0.T)\n",
    "Z0[0, :] = 0.3 * X0[1, 10:]\n",
    "Z0[1, :] = 0.5 * X0[2, 10:]\n",
    "Z0[2, :] = 0.7 * X0[0, 10:]\n",
    "plot(Z0.T)\n",
    "(NX, LX, X) = ERpy.NDTOMM(X0[:, :90])\n",
    "(NZ, LZ, Z) = ERpy.NDTOMM(Z0[:2, :90])\n",
    "\n",
    "(F, Y, E) = ERpy.WIENER(NX, LX, X, NZ, LZ, Z, 10, 0, 0)\n",
    "\n",
    "print(Y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Help on function MACRO_partial in module ERpy:\n",
      "\n",
      "MACRO_partial(N, LG, G, I, J, K)\n",
      "    MACRO_partial multichannel cross partial correlation\n",
      "    \n",
      "    G = MACRO(X, Y, LG)\n",
      "    \n",
      "    Example\n",
      "    \n",
      "    LG = 4\n",
      "    N = 3\n",
      "    G = array([\n",
      "        1., 0.9, 0.8, 0.7, 0.2, 0.3, 0.2, 0.1, 0.1, 0.1, 0.1, 0.1, \n",
      "        0.1, 0.1, 0.2, 0.1, 1, 0.8, 0.6, 0.4, 0.1, 0., 0., 0., \n",
      "        0., 0., 0., 0., 0.2, 0.3, .4, .5, 1., 0.7, 0.4, 0.1])\n",
      "    GP = MACRO_partial(N, LG, G, 0, 1, 2)\n",
      "    print(GP)\n",
      "\n"
     ]
    }
   ],
   "source": [
    "help(ERpy.MACRO_partial)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[1.0  0.2  0.1] \n",
      " [0.9  0.3  0.1] \n",
      " [0.8  0.2  0.1] \n",
      " [0.7  0.1  0.1]]\n",
      "\n",
      "[[0.1  1.0  0.1] \n",
      " [0.1  0.8  0.0] \n",
      " [0.2  0.6  0.0] \n",
      " [0.1  0.4  0.0]]\n",
      "\n",
      "[[0.0  0.2  1.0] \n",
      " [0.0  0.3  0.7] \n",
      " [0.0  0.4  0.4] \n",
      " [0.0  0.5  0.1]]\n",
      "\n",
      "None\n",
      "(0, 0, 0, array([ nan,  0.9,  0.8,  0.7]))\n",
      "(0, 0, 1, array([ 1.1       ,  0.98888889,  0.87777778,  0.76666667]))\n",
      "(0, 0, 2, array([ 1.6       ,  1.23333333,  0.87777778,  0.7       ]))\n",
      "(0, 1, 0, array([       -inf, -0.28284271, -0.37416574, -0.08164966]))\n",
      "(0, 1, 1, array([ 0.1      ,  0.1      ,  0.0942809,  0.1      ]))\n",
      "(0, 1, 2, array([ 0.15811388,  0.12909944,  0.10540926,  0.1       ]))\n",
      "(0, 2, 0, array([        inf,  1.03333333,  0.04714045, -0.13471506]))\n",
      "(0, 2, 1, array([ 0.65555556,         inf,  0.04714045, -0.1       ]))\n",
      "(0, 2, 2, array([ 0.9486833 ,  0.51639778,  0.10540926,  0.        ]))\n",
      "(1, 0, 0, array([        nan,  0.07071068,  0.16035675,  0.12247449]))\n",
      "(1, 0, 1, array([ 0.76666667,  0.21111111,  0.32998316,  0.21111111]))\n",
      "(1, 0, 2, array([ 1.10679718,  0.25819889,  0.31622777,  0.2       ]))\n",
      "(1, 1, 0, array([-1.3       ,  0.075     ,  0.15714286,  0.075     ]))\n",
      "(1, 1, 1, array([ 0.1,  0.1,  0.2,  0.1]))\n",
      "(1, 1, 2, array([ 0.1,  0.1,  0.2,  0.1]))\n",
      "(1, 2, 0, array([-0.2857738 , -0.02357023, -0.03779645, -0.02357023]))\n",
      "(1, 2, 1, array([-0.01111111,        -inf, -0.4       , -0.1       ]))\n",
      "(1, 2, 2, array([ 0.,  0.,  0.,  0.]))\n",
      "(2, 0, 0, array([        nan,  0.03333333,  0.04714045,  0.05773503]))\n",
      "(2, 0, 1, array([ 0.21111111,         nan,  0.04714045,  0.06666667]))\n",
      "(2, 0, 2, array([ 0.31622777,  0.12909944,  0.10540926,  0.1       ]))\n",
      "(2, 1, 0, array([-0.08164966,  1.15494108,  0.97010881,  0.68353656]))\n",
      "(2, 1, 1, array([ 0.1,  inf,  1.6,  0.9]))\n",
      "(2, 1, 2, array([ 0.1,  1. ,  0.8,  0.6]))\n",
      "(2, 2, 0, array([-0.05      , -0.01111111, -0.01111111, -0.01111111]))\n",
      "(2, 2, 1, array([-0.01111111,        -inf, -3.2       , -0.9       ]))\n",
      "(2, 2, 2, array([ 0.,  0.,  0.,  0.]))\n"
     ]
    }
   ],
   "source": [
    "imp.reload(ERpy)\n",
    "LG = 4\n",
    "N = 3\n",
    "G = array([\n",
    "    1., 0.9, 0.8, 0.7, 0.2, 0.3, 0.2, 0.1, 0.1, 0.1, 0.1, 0.1, \n",
    "    0.1, 0.1, 0.2, 0.1, 1, 0.8, 0.6, 0.4, 0.1, 0., 0., 0., \n",
    "    0., 0., 0., 0., 0.2, 0.3, .4, .5, 1., 0.7, 0.4, 0.1])\n",
    "print(ER.RITE(1, LG, N, N, G))\n",
    "for i in range(3): \n",
    "    for j in range(3): \n",
    "        for k in range(3): \n",
    "            GP = ERpy.MACRO_partial(N, LG, G, i, j, k)\n",
    "            print(i, j, k, GP)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Help on function RITE in module ER:\n",
      "\n",
      "RITE(NB, M, N, L, A)\n",
      "    RITE \n",
      "    \n",
      "    p. 55\n",
      "    \n",
      "    This don't follow the code given in p. 55 but should produce the same result. \n",
      "    \n",
      "    RITE(NB, M, N, L, A)\n",
      "    \n",
      "    NB = number of matrices printed per line = u \n",
      "     M = number of rows in each matrix a_k\n",
      "     N = number of columns in each matrix a_k\n",
      "     L = number of matrices = q\n",
      "     A = (a_1, a_2, \\dots, a_q) stored in multiplexed mode, i.e., stored as A(I,J,K) where I denotes the row, J denotes te column, and K denotes the matrix index.\n",
      "\n"
     ]
    }
   ],
   "source": [
    "help(ER.RITE)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}